The KM3NeT research infrastructure is currently under construction at two locations in the Mediterranean Sea. The KM3NeT/ORCA water-Cherenkov neutrino detector off the French coast will instrument several megatons of seawater with photosensors. Its main objective is the determination of the neutrino mass ordering. This work aims at demonstrating the general applicability of deep convolutional neural networks to neutrino telescopes, using simulated datasets for the KM3NeT/ORCA detector as an example. To this end, the networks are employed to achieve reconstruction and classification tasks that constitute an alternative to the analysis pipeline presented for KM3NeT/ORCA in the KM3NeT Letter of Intent. They are used to infer event reconstruction estimates for the energy, the direction, and the interaction point of incident neutrinos. The spatial distribution of Cherenkov light generated by charged particles induced in neutrino interactions is classified as shower- or track-like, and the main background processes associated with the detection of atmospheric neutrinos are recognized. Performance comparisons to machine-learning classification and maximum-likelihood reconstruction algorithms previously developed for KM3NeT/ORCA are provided. It is shown that this application of deep convolutional neural networks to simulated datasets for a large-volume neutrino telescope yields competitive reconstruction results and performance improvements with respect to classical approaches.

The gSeaGen code is a GENIE-based application developed to efficiently generate high statistics samples of events, induced by neutrino interactions, detectable in a neutrino telescope. The gSeaGen code is able to generate events induced by all neutrino flavours, considering topological differences between track-type and shower-like events. Neutrino interactions are simulated taking into account the density and the composition of the media surrounding the detector. The main features of gSeaGen are presented together with some examples of its application within the KM3NeT project.

Jets in blazars are an excellent forum for studying acceleration at relativistic shocks using the highly-variable emission seen across the electromagnetic spectrum. Our recent work on combining multi-wavelength leptonic emission models with simulated thermal+non-thermal distributions from shock acceleration theory has resulted in new insights into plasma conditions in blazars. This has demonstrated the ability to infer the cyclotron frequency, the plasma density and thus also the Alfven speed, thereby determining the rapidity of particle energization. An important inference was that MHD turbulence levels decline with remoteness from jet shocks. This paper outlines new results from our recent extension of this program to a two-zone, time-evolving construction, modeling together both extended, enhanced emission states from larger radiative regions, and prompt flare events from compact acceleration zones. These are applied to flares in the FSRQ blazar 3C 279 monitored by Fermi-LAT in gamma-rays in late 2013. With impulsive injection episodes from the shock zone, as the acceleration first proceeds and then abates, the radiative simulations obtain a pronounced spectral hardening in the optical and gamma-ray bands as the flare grows, followed by a softening during the decay phase. For 3C 279, while model radio and X-ray synchrotron flares are temporally correlated, there is a lag in both bands relative to GeV gamma rays and optical emission on timescales of a number of hours. This delay is governed by the short cooling time associated with the bright external Compton signal.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Gamma-ray bursts (GRBs) are brief flashes of gamma rays, considered to be the most energetic explosive phenomena in the Universe. The emission from GRBs comprises a short (typically tens of seconds) and bright prompt emission, followed by a much longer afterglow phase. During the afterglow phase, the shocked outflow -- produced by the interaction between the ejected matter and the circumburst medium -- slows down, and a gradual decrease in brightness is observed. GRBs typically emit most of their energy via gamma-rays with energies in the kiloelectronvolt-to-megaelectronvolt range, but a few photons with energies of tens of gigaelectronvolts have been detected by space-based instruments. However, the origins of such high-energy (above one gigaelectronvolt) photons and the presence of very-high-energy (more than 100 gigaelectronvolts) emission have remained elussive. Here we report observations of very-high-energy emission in the bright GRB 180720B deep in the GRB afterglow -ten hours after the end of the prompt emission phase, when the X-ray flux had already decayed by four orders of magnitude. Two possible explanations exist for the observed radiation: inverse Compton emission and synchrotron emission of ultrarelativistic electrons. Our observations show that the energy fluxes in the X-ray and gamma-ray range and their photon indices remain comparable to each other throughout the afterglow. This discovery places distinct constraints on the GRB environment for both emission mechanisms, with the inverse Compton explanation alleviating the particle energy requirements for the emission observed at late times. The late timing of this detection has consequences for the future observations of GRBs at the highest energies.

The FSRQ CTA 102 (z=1.032) has been tremendously active over the last few years. During its peak activity lasting several months in late 2016 and early 2017, the gamma-ray and optical fluxes rose by up to a factor 100 above the quiescence level. We have interpreted the peak activity as the ablation of a gas cloud by the relativistic jet, which can nicely account for the months-long lightcurve in 2016 and 2017. The peak activity was in the middle of a 2-year-long high-state, which was characterized by increased fluxes and increased rms variability compared to the previous low-states, and which was flanked by two bright flares. In this presentation, we put the cloud-ablation scenario into the broader context of the 2-year-long high-state.

Mildly-relativistic shocks that are embedded in colliding magnetohydrodynamic flows are prime sites for relativistic particle acceleration and production of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars and gamma-ray bursts. The principal energization mechanisms at these shocks are diffusive shock acceleration and shock drift acceleration. In recent work, we had self-consistently coupled shock acceleration and radiation transfer simulations in blazar jets in a basic one-zone scenario. These one-zone models revealed that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of the shock layer hydromagnetic turbulence. In this paper, we expand our previous work by including full time dependence and treating two zones, one being the site of acceleration, and the second being a larger emission zone. This construction is applied to multiwavelength flares of the flat spectrum radio quasar 3C 279, fitting snap-shot SEDs and generating light curves that are consistent with observed variability timescales. We also present a generic study for the typical flaring behavior of the BL Lac object Mrk 501. The model predicts correlated variability across all wavebands, but cross-band time lags depending on the type of blazar (FSRQ vs. BL Lac), as well as distinctive spectral hysteresis patterns in all wavelength bands, from mm radio waves to gamma-rays. These evolutionary signatures serve to provide diagnostics on the competition between acceleration and radiative cooling.

We present a procedure to generally constrain the environments of neutrino-producing sites in photomeson production models of jetted Active Galactic Nuclei (AGN) where any origin of the dominant target photon field can be accommodated. For this purpose we reconstruct the minimum target photon spectrum required to produce the (observed) neutrino spectrum, and derive the distributions of all corresponding secondary particles. These initiate electromagnetic cascades with an efficiency that is linked to the neutrino production rate. The derived photon spectra represent the minimum radiation emerging from the source that is strictly associated with the photo-hadronically produced neutrinos. Using the 2014-15 neutrino spectrum observed by IceCube from TXS 0506+056, we conduct a comprehensive study of these cascade spectra and compare them to the simultaneous multi-wavelength emission. For this set of observations, photopion production from a co-spatially produced (co-moving) photon target can be ruled out as well as a setup where synchrotron or Compton-synchrotron supported cascades on a stationary (AGN rest frame) target photon field operate in this source. However, a scenario where Compton-driven cascades develop in the stationary soft-X-ray photon target which photo-hadronically produced the observed neutrinos appears feasible with required proton kinetic jet powers near the Eddington limit. The source is then found to produce neutrinos inefficiently, and emits GeV photons significantly below the observed Fermi-LAT-flux. Hence, the neutrinos and the bulk of the gamma rays observed in 2014/2015 from TXS 0506+056 cannot have been initiated by the same process.

In this work, we present a spectral and temporal analysis of {\it Swift}/XRT and {\it NuSTAR} observations of GRS 1716--249 during its recent 2016--2017 outburst. This low mass X-ray binary underwent an extraordinary outburst after a long quiescence of 23 years, since its last major outburst in 1993. The source was observed over two different epochs during 2017 April, 07 and 10. The best fit joint spectral fitting in the energy range 0.5 $-$ 79.0 keV indicates that the spectrum is best described by relatively cold, weak disk blackbody emission, dominant thermal Comptonization emission, and a relativistically broadened fluorescent iron K$\alpha$ emission line. We observed a clear indication of a Compton hump around 30 keV. We also detected an excess feature of $\sim1.3$ keV. Assuming a lamp-post geometry of the corona, we constrained the inner disk radius for both observations to 11.92$^{+8.62}_{-11.92}$ R$_{ISCO}$ (i.e., an upper limit) and 10.39$^{+9.51}_{-3.02}$ R$_{ISCO}$ (where R$_{ISCO}\equiv$ radius of the innermost stable circular orbit) for the first epoch (E1) and second epoch (E2), respectively. A significant ($\sim5\sigma$) type$-$C quasi-periodic oscillation (QPO) at $1.20\pm0.04$ Hz is detected for the first time for GRS 1716--249, which drifts to $1.55\pm0.04$ Hz ($\sim6\sigma$) at the end of the second observation. The derived spectral and temporal properties show a positive correlation between the QPO frequency and the photon index.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

PKS 1830-211 is a known macrolensed quasar located at a redshift of z=2.5. Its high-energy gamma-ray emission has been detected with the Fermi-LAT instrument and evidence for lensing was obtained by several authors from its high-energy data. Observations of PKS 1830-211 were taken with the H.E.S.S. array of Imaging Atmospheric Cherenkov Telescopes in August 2014, following a flare alert by the Fermi- LAT collaboration. The H.E.S.S observations were aimed at detecting a gamma-ray flare delayed by 20-27 days from the alert flare, as expected from observations at other wavelengths. More than twelve hours of good quality data were taken with an analysis threshold of $\sim67$ GeV. The significance of a potential signal is computed as a function of the date as well as the average significance over the whole period. Data are compared to simultaneous observations by Fermi-LAT. No photon excess or significant signal is detected. An upper limit on PKS 1830-211 flux above 67 GeV is computed and compared to the extrapolation of the Fermi-LAT flare spectrum.

Mildly relativistic, oblique shocks are frequently invoked as possible sites of relativistic particle acceleration and production of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars, via diffusive shock acceleration (DSA). In recent work, we had self-consistently coupled DSA and radiation transfer simulations in blazar jets. These one-zone models determined that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of the hydromagnetic turbulence responsible for pitch-angle scattering. In this paper, we expand our previous work by including full time dependence and treating two emission zones, one being the site of acceleration. This modeling is applied to a multiwavelength flare of the flat spectrum radio quasar 3C~279, fitting snap-shot SEDs and light curves. We predict spectral hysteresis patterns in various energy bands as well as cross-band time lags with optical and GeV gamma-rays as well as radio and X-rays tracing each other closely with zero time lag, but radio and X-rays lagging behind the optical and gamma-ray variability by several hours.

This review provides an overview of recent advances in multi-wavelength and multi-messenger observations of blazars, the current status of theoretical models for blazar emission, and prospects for future facilities. The discussion of observational results will focus on advances made possible through the Fermi Gamma-Ray Space Telescope and ground-based gamma-ray observatories (H.E.S.S., MAGIC, VERITAS) as well as the recent first evidence for a blazar being a source of IceCube neutrinos. The main focus of this review will be the discussion of our current theoretical understanding of blazar multi-wavelength and multi-messenger emission, in the spectral, time, and polarization domains. Future progress will be expected in particular through the development of the first X-ray polarimeter, IXPE, and the installation of the Cherenkov Telescope Array (CTA), both expected to become operational in the early to mid 2020s.

We test different physically motivated models for the spectral shape of the $\gamma$-ray emission in a sample of 128 blazars with known redshifts detected by the Fermi Large Area Telescope (LAT) at energies above 50 GeV. The first nine years of LAT data in the energy range from 300 MeV to 2 TeV are analyzed in order to extend the spectral energy coverage of the 2FHL blazars in our sample. We compare these spectral data to four leptonic models for the production of $\gamma$-rays through Compton scattering by a population of electrons with different spectral shapes. In the first three models we consider Compton scattering in the Thomson regime with different acceleration mechanisms for the electrons. In the fourth model we consider Compton scattering by a pure power law distribution of electrons with spectral curvature due to scattering in the Klein-Nishina regime. The majority of blazar $\gamma$-ray spectra are preferentially fit with either a power law with exponential cut-off in the Thomson regime or a power law electron distribution with Compton scattering in the Klein-Nishina regime, while a log-parabola with a low-energy power-law and broken power-law spectral shape in the Thomson regime appears systematically disfavoured, which is likely a consequence of the restriction to pure Thomson scattering which we imposed on those models. This finding may be an indication that the $\gamma$-ray emission from FSRQs in the 2FHL catalog is dominated by Compton scattering of radiation from the dusty torus, while in the case of BL Lac objects, it is dominated by synchrotron self-Compton radiation.

In this paper, we consider the impact of Lorentz Invariance Violation (LIV) on the $\gamma-\gamma$ opacity of the Universe to VHE-gamma rays, compared to the effect of local under-densities (voids) of the Extragalactic Background Light, and on the Compton scattering process. Both subluminal and superluminal modifications of the photon dispersion relation are considered. In the subluminal case, LIV effects may result in a significant reduction of the $\gamma-\gamma$ opacity for photons with energies $\gtrsim 10$ TeV. However, the effect is not expected to be sufficient to explain the apparent spectral hardening of several observed VHE $\gamma$-ray sources in the energy range from 100 GeV - a few TeV, even when including effects of plausible inhomogeneities in the cosmic structure. Superluminal modifications of the photon dispersion relation lead to a further enhancement of the EBL $\gamma\gamma$ opacity. We consider, for the first time, the influence of LIV on the Compton scattering process. We find that this effect becomes relevant only for photons at ultra-high energies, $E \gtrsim 1$ PeV. In the case of a superluminal modification of the photon dispersion relation, both the kinematic recoil effect and the Klein-Nishina suppression of the cross section are reduced. However, we argue that the effect is unlikely to be of astrophysical significance.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

We report on recent results from a target-of-opportunity program to obtain spectropolarimetry observations with the Southern African Large Telescope (SALT) on flaring gamma-ray blazars. SALT spectropolarimetry and contemporaneous multi-wavelength spectral energy distribution (SED) data are being modelled self-consistently with a leptonic single-zone model. Such modeling provides an accurate estimate of the degree of order of the magnetic field in the emission region and the thermal contributions (from the host galaxy and the accretion disk) to the SED, thus putting strong constraints on the physical parameters of the gamma-ray emitting region. For the specific case of the $\gamma$-ray blazar 4C+01.02, we demonstrate that the combined SED and spectropolarimetry modeling constrains the mass of the central black hole in this blazar to $M_{\rm BH} \sim 10^9 \, M_{\odot}$.

We compare Particle-in-Cell simulation results of relativistic electron-ion shear flows with different bulk Lorentz factors, and discuss their implications for spine-sheath models of blazar versus gamma-ray burst (GRB) jets. Specifically, we find that most properties of the shear boundary layer scale with the bulk Lorentz factor: the lower the Lorentz factor, the thinner the boundary layer, and the weaker the self-generated fields. Similarly, the energized electron spectrum peaks at an energy near the ion drift energy, which increases with bulk Lorentz factor, and the beaming of the accelerated electrons gets narrower with increasing Lorentz factor. This predicts a strong correlation between emitted photon energy, angular beaming and temporal variability with the bulk Lorentz factor. Observationally, we expect systematic differences between the high-energy emissions of blazars and GRB jets.

The unexpectedly hard very-high-energy (VHE; $E > 100$ GeV) $\gamma$-ray spectra of a few distant blazars have been interpreted as evidence for a reduction of the $\gamma\gamma$ opacity of the Universe due to the interaction of VHE $\gamma$-rays with the extragalactic background light (EBL) compared to the expectation from our current knowledge of the density and cosmological evolution of the EBL. One of the suggested solutions to this problem consisted of the inhomogeneity of the EBL. In this paper, we study the effects of such inhomogeneities on the energy density of the EBL (which then also becomes anisotropic) and the resulting $\gamma\gamma$ opacity. Specifically, we investigate the effects of cosmic voids along the line of sight to a distant blazar. We find that the effect of such voids on the $\gamma\gamma$ opacity, for any realistic void size, is only of the order of $\lesssim 1$ % and much smaller than expected from a simple linear scaling of the $\gamma\gamma$ opacity with the line-of-sight galaxy under-density due to a cosmic void.

Collimated outflows (jets) appear to be a ubiquitous phenomenon associated with the accretion of material onto a compact object. Despite this ubiquity, many fundamental physics aspects of jets are still poorly understood and constrained. These include the mechanism of launching and accelerating jets, the connection between these processes and the nature of the accretion flow, and the role of magnetic fields; the physics responsible for the collimation of jets over tens of thousands to even millions of gravitational radii of the central accreting object; the matter content of jets; the location of the region(s) accelerating particles to TeV (possibly even PeV and EeV) energies (as evidenced by gamma-ray emission observed from many jet sources) and the physical processes responsible for this particle acceleration; the radiative processes giving rise to the observed multi-wavelength emission; and the topology of magnetic fields and their role in the jet collimation and particle acceleration processes. This chapter reviews the main knowns and unknowns in our current understanding of relativistic jets, in the context of the main model ingredients for Galactic and extragalactic jet sources. It discusses aspects specific to active Galactic nuclei (especially blazars) and microquasars, and then presents a comparative discussion of similarities and differences between them.

The knowledge of the structure of the magnetic field inside a blazar jet, as deduced from polarization observations at radio to optical wavelengths, is closely related to the formation and propagation of relativistic jets that result from accretion onto supermassive black holes. However, a largely unexplored aspect of the theoretical understanding of radiation transfer physics in blazar jets has been the magnetic field geometry as revealed by the polarized emission and the connection between the variability in polarization and flux across the spectrum. Here, we explore the effects of various magnetic geometries that can exist inside a blazar jet: parallel, oblique, toroidal, and tangled. We investigate the effects of changing the orientation of the magnetic field, according to the above-mentioned geometries, on the resulting high-energy spectral energy distributions (SEDs) and spectral variability patterns (SVPs) of a typical blazar. We use the MUlti-ZOne Radiation Feedback (MUZORF) model of Joshi et al. (2014) to carry out this study and to relate the geometry of the field to the observed SEDs at X-ray and gamma-ray energies. One of the goals of the study is to understand the relationship between synchrotron and inverse Compton peaks in blazar SEDs and the reason for the appearance of gamma-ray "orphan flares" observed in some blazars. This can be associated with the directionality of the magnetic field, which creates a difference in the radiation field as seen by an observer versus that seen by the electrons in the emission region.

Blazars radiate from radio through gamma-ray frequencies thereby being ideal targets for multifrequency studies. Such studies allow constraining the properties of the emitting jet. 3C 279 is among the most notable blazars and therefore subject to extensive multifrequency campaigns. We report the results of a campaign ranging from near-IR to gamma-ray energies of an outburst of 3C 279 in June 2015. The overall spectral energy distribution from near-IR to gamma rays can be well represented by either a leptonic or a lepto-hadronic radiation transfer model. Even though the data are equally well represented by the two models, their inferred parameters challenge the physical conditions in the jet. In fact, the leptonic model requires parameters with a magnetic field far below equipartition with the relativistic particle energy density. On the contrary, equipartition may be achieved with the lepto-hadronic model, which however implies an extreme total jet power close to Eddington luminosity.

Relativistic shocks are one of the most plausible sites of the emission of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars, via diffusive shock acceleration (DSA) of relativistic particles. This paper summarizes recent results on a self-consistent coupling of diffusive shock acceleration and radiation transfer in blazar jets. We demonstrate that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of hydromagnetic turbulence responsible for pitch-angle scattering by requiring a strongly energy-dependent pitch-angle mean free path. The prominent soft X-ray excess ("Big Blue Bump") in the SED of the BL Lac object AO 0235+164 can be modelled as the signature of bulk Compton scattering of external radiation fields by the thermal electron population, which places additional constraints on the level of hydromagnetic turbulence. It has further been demonstrated that internal shocks propagating in a jet pervaded by a helical magnetic field naturally produce polarization-angle swings by 180$^o$, in tandem with multi-wavelength flaring activity, without requiring any helical motion paths or other asymmetric jet structures. The specific application of this model to 3C279 presents the first consistent, simultaneous modeling of snap-shot SEDs, multi-wavelength light curves and time-dependent polarization signatures of a blazar during a polarization-angle (PA) rotation. This model has recently been generalized to a lepto-hadronic model, in which the high-energy emission is dominated by proton synchrotron radiation. It is shown that in this case, the high-energy (X-ray and $\gamma$-ray) polarization signatures are expected to be significantly more stable (not showing PA rotations) than the low-energy (electron-synchrotron) signatures.

The magnification effect due to gravitational lensing enhances the chances of detecting moderate-redshift ($z \sim 1$) sources in very-high-energy (VHE; $E > 100$ GeV) $\gamma$-rays by ground-based Atmospheric Cherenkov Telescope facilities. It has been shown in previous work that this prospect is not hampered by potential $\gamma-\gamma$ absorption effects by the intervening (lensing) galaxy, nor by any individual star within the intervening galaxy. In this paper, we expand this study to simulate the light-bending effect of a realistic ensemble of stars. We first demonstrate that, for realistic parameters of the galaxy's star field, it is extremely unlikely (probability $\lesssim 10^{-6}$) that the direct line of sight between the $\gamma$-ray source and the observer passes by any star in the field close enough to be subject to significant $\gamma\gamma$ absorption. Our simulations then focus on the rare cases where $\gamma\gamma$ absorption by (at least) one individual star might be non-negligible. We show that gravitational light bending will have the effect of avoiding the $\gamma-\gamma$ absorption spheres around massive stars in the intervening galaxy. This confirms previous results by Barnacka et al. and re-inforces arguments in favour of VHE $\gamma$-ray observations of lensed moderate-redshift blazars to extend the redshift range of objects detected in VHE $\gamma$-rays, and to probe the location of the $\gamma$-ray emission region in those blazars.

We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3{\sigma} level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.

The BL Lacertae object 1ES 1440+122 was observed in the energy range from 85 GeV to 30 TeV by the VERITAS array of imaging atmospheric Cherenkov telescopes. The observations, taken between 2008 May and 2010 June and totalling 53 hours, resulted in the discovery of $\gamma$-ray emission from the blazar, which has a redshift $z$=0.163. 1ES 1440+122 is detected at a statistical significance of 5.5 standard deviations above the background with an integral flux of (2.8$\pm0.7_{\mathrm{stat}}\pm0.8_{\mathrm{sys}}$) $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ (1.2\% of the Crab Nebula's flux) above 200 GeV. The measured spectrum is described well by a power law from 0.2 TeV to 1.3 TeV with a photon index of 3.1 $\pm$ 0.4$_{\mathrm{stat}}$ $\pm$ 0.2$_{\mathrm{sys}}$. Quasi-simultaneous multi-wavelength data from the Fermi Large Area Telescope (0.3--300 GeV) and the Swift X-ray Telescope (0.2--10 keV) are additionally used to model the properties of the emission region. A synchrotron self-Compton model produces a good representation of the multi-wavelength data. Adding an external-Compton or a hadronic component also adequately describes the data.

In this study, we use a one-zone leptonic and a lepto-hadronic model to investigate the multi-wavelength emission and the prominent flare of the flat spectrum radio quasar 3C 454.3 in Nov 2010. We perform a parameter study with both models to obtain broadband fits to the spectral energy distribution of 3C 454.3. Starting with the baseline parameters obtained from the fits, we then investigate different flaring scenarios for both models to explain an extreme outburst and spectral hardening of 3C 454.3 that occurred in Nov 2010. We find that the one zone lepto-hadronic model can successfully explain both the broadband multi-wavelength spectral energy distribution and light curves in the optical R, Swift XRT and Fermi gamma-ray bandpasses for 3C 454.3 during quiescence and the peak of the Nov. 2010 flare. We also find that the one-zone leptonic model produces poor fits to the broadband spectra in the X-ray and high energy gamma-ray band passes for the Nov. 2010 flare.

The expected level of gamma-gamma absorption in the Broad Line Region (BLR) radiation field of gamma-ray loud Flat Spectrum Radio Quasars (FSRQs)is evaluated as a function of the location of the gamma-ray emission region. This is done self-consistently with parameters inferred from the shape of the spectral energy distribution (SED) in a single-zone leptonic EC-BLR model scenario. We take into account all geometrical effects both in the calculation of the gamma-gamma opacity and the normalization of the BLR radiation energy density. As specific examples, we study the FSRQs 3C279 and PKS 1510-089, keeping the BLR radiation energy density at the location of the emission region fixed at the values inferred from the SED. We confirm previous findings that the optical depth due to $\gamma\gamma$ absorption in the BLR radiation field exceeds unity for both 3C279 and PKS 1510-089 for locations of the gamma-ray emission region inside the inner boundary of the BLR. It decreases monotonically, with distance from the central engine and drops below unity for locations within the BLR. For locations outside the BLR, the BLR radiation energy density required for the production of GeV gamma-rays rapidly increases beyond observational constraints, thus making the EC-BLR mechanism implausible. Therefore, in order to avoid significant gamma-gamma absorption by the BLR radiation field, the gamma-ray emission region must therefore be located near the outer boundary of the BLR.

A large number of new "black widow" and "redback" energetic millisecond pulsars with irradiated stellar companions have been discovered through radio and optical searches of unidentified \textit{Fermi} sources. Synchrotron emission, from particles accelerated up to several TeV in the intrabinary shock, exhibits modulation at the binary orbital period. Our simulated double-peaked X-ray light curves modulated at the orbital period, produced by relativistic Doppler-boosting along the intrabinary shock, are found to qualitatively match those observed in many sources. In this model, redbacks and transitional pulsar systems where the double-peaked X-ray light curve is observed at inferior conjunction have intrinsically different shock geometry than other millisecond pulsar binaries where the light curve is centered at superior conjunction. We discuss, and advocate, how current and future optical observations may aid in constraining the emission geometry, intrabinary shock and the unknown physics of pulsar winds.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

LS 5039 is a gamma-ray binary system observed in a broad energy range, from radio to TeV energies. The binary system exhibits both flux and spectral modulation as a function of its orbital period. The X-ray and very-high-energy (VHE, E > 100 GeV) gamma-ray fluxes display a maximum/minimum at inferior/superior conjunction, with spectra becoming respectively harder/softer, a behaviour that is completely reversed in the high-energy domain (HE, 0.1 < E < 100 GeV). The HE spectrum cuts off at a few GeV, with a new hard component emerging at E > 10 GeV that is compatible with the low-energy tail of the TeV emission. The low 10 - 100 GeV flux, however, makes the HE and VHE components difficult to reconcile with a scenario including emission from only a single particle population. We report on new observations of LS 5039 conducted with the High Energy Stereoscopic System (H.E.S.S.) telescopes from 2006 to 2015. This new data set enables for an unprecedentedly-deep phase-folded coverage of the source at TeV energies, as well as an extension of the VHE spectral range down to ~120 GeV, which makes LS 5039 the first gamma-ray binary system in which a spectral overlap between satellite and ground-based gamma-ray observatories is obtained.

Blazars are astrophysical sources whose emission is dominated by non-thermal processes, typically interpreted as synchrotron and inverse Compton emission. Although the general picture is rather robust and consistent with observations, many aspects are still unexplored. Polarimetric monitoring can offer a wealth of information about the physical processes in blazars. Models with largely different physical ingredients can often provide almost indistinguishable predictions for the total flux, but usually are characterized by markedly different polarization properties. We explore, with a pilot study, the possibility to derive structural information about the emitting regions of blazars by means of a joint analysis of rapid variability of the total and polarized flux at optical wavelengths. Short timescale (from tens of seconds to a couple of minutes) optical linear polarimetry and photometry for two blazars, BL Lacertae and PKS 1424+240, was carried out with the PAOLO polarimeter at the 3.6m Telescopio Nazionale Galileo. Several hours of almost continuous observations were obtained for both sources. Our intense monitoring allowed us to draw strongly different scenarios for BL Lacertae and PKS 1424+240, with the former characterized by intense variability on time-scales from hours to a few minutes and the latter practically constant in total flux. Essentially the same behavior is observed for the polarized flux and the position angle. The variability time-scales turned out to be as short as a few minutes, although involving only a few percent variation of the flux. The polarization variability time-scale is generally consistent with the total flux variability. Total and polarized flux appear to be essentially uncorrelated. However, even during our relatively short monitoring, different regimes can be singled out. (abridged)

We study hadronic models of broad-band emission of jets in radio-loud active galactic nuclei, and their implications for the accretion in those sources. We show that the models that account for broad-band spectra of blazars emitting in the GeV range in the sample of Boettcher et al. have highly super-Eddington jet powers. Furthermore, the ratio of the jet power to the radiative luminosity of the accretion disc is $\sim 3000$ on average and can be as high as $\sim 10^5$. We then show that the measurements of the radio core shift for the sample imply low magnetic fluxes threading the black hole, which rules out the Blandford-Znajek mechanism to produce powerful jets. These results require that the accretion rate necessary to power the modelled jets is extremely high, and the average radiative accretion efficiency is $\sim 4 \times 10^{-5}$. Thus, if the hadronic model is correct, the currently prevailing picture of accretion in AGNs needs to be significantly revised. Also, the obtained accretion mode cannot be dominant during the lifetimes of the sources, as the modelled very high accretion rates would result in too rapid growth of the central supermassive black holes. Finally, the extreme jet powers in the hadronic model are in conflict with the estimates of the jet power by other methods.

The polarization signatures of the blazar emissions are known to be highly variable. In addition to small fluctuations of the polarization angle around a mean value, sometimes large (> 180^o) polarization angle swings are observed. We suggest that such p henomena can be interpreted as arising from light-travel-time effects within an underlying axisymmetric emission region. We present the first simultaneous fitting of the multi-wavelength spectrum, variability and time-dependent polarization features of a correlated optical and gamma-ray flaring event of the prominent blazar 3C279, which was accompanied by a drastic change of its polarization signatures. This unprecedented combination of spectral, variability, and polarization information in a coherent physical model allows us to place stringent constraints on the particle acceleration and magnetic-field topology in the relativistic jet of a blazar, strongly favoring a scenario in which magnetic energy dissipation is the primary driver of the flare event.

We introduce a new time-dependent lepto-hadronic model for blazar emission that takes into account the radiation emitted by secondary particles, such as pions and muons, from photo hadronic interactions. Starting from a baseline parameter set guided by a fit to the spectral energy distribution of the blazar 3C 279, we perform a parameter study to investigate the effects of perturbations of the input parameters to mimic different flaring events to study the resulting lightcurves in the optical, X-ray, high energy (HE: E > 100 MeV) and very-high-energy (VHE: E > 100 GeV) gamma-rays as well as the neutrino emission associated with charged-pion and muon decay. We find that flaring events from an increase in the efficiency of Fermi II acceleration will produce a positive correlation between all bandpasses and a marked plateau in the HE gamma-ray lightcurve. We also predict a distinctive dip in the HE lightcurve for perturbations caused by a change in the proton injection spectral index. These plateaus / dips could be a tell tale signature of hadronic models for perturbations that lead to more efficient acceleration of high energy protons in parameter regimes where pion and muon synchrotron emission is non-negligible.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

We study the acceleration, transport, and emission of particles in relativistic jets. Localized stochastic particle acceleration, spatial diffusion, and synchrotron as well as synchrotron self-Compton emission are considered in a leptonic model. To account for inhomogeneity, we use a 2D axi-symmetric cylindrical geometry for both relativistic electrons and magnetic field. In this first phase of our work, we focus on steady-state spectra that develop from a time-dependent model. We demonstrate that small isolated acceleration region in a much larger emission volume are sufficient to accelerate particles to high energy. Diffusive escape from these small regions provides a natural explanation for the spectral form of the jet emission. The location of the acceleration regions within the jet is found to affect the cooling break of the spectrum in this diffusive model. Diffusion-caused energy-dependent inhomogeneity in the jets predicts that the SSC spectrum is harder than the synchrotron spectrum. There can also be a spectral hardening towards the high-energy section of the synchrotron spectrum, if particle escape is relatively slow. These two spectral hardening effects indicate that the jet inhomogeneity might be a natural explanation for the unexpected hard {\gamma}-ray spectra observed in some blazars.

This paper presents an analysis of the dependence of the opacity for high-energy gamma-rays to gamma-gamma absorption by low-energy photons, on the polarization of the gamma-ray and target photons. This process has so far only been considered using the polarization-averaged gamma-gamma absorption cross section. It is demonstrated that in the case of polarized gamma-ray emission, subject to source-intrinsic gamma-gamma absorption by polarized target photons, this may lead to a slight over-estimation of the gamma-gamma opacity by up to ~ 10 % in the case of a perfectly ordered magnetic field. Thus, for realistic astrophysical scenarios with partially ordered magnetic fields, the use of the polarization-averaged gamma-gamma cross section is justified for practical purposes, such as estimates of minimum Doppler factors inferred for gamma-ray bursts and blazars, based on gamma-gamma transparency arguments, and this paper quantifies the small error incurred by the unpolarized-radiation approximation. Furthermore, it is shown that polarization-dependent gamma-gamma absorption of initially polarized gamma-rays can lead to a slight increase in the polarization beyond the spectral break caused by gamma-gamma absorption, to an amount distinctly different from the change in polarization expected if the same spectral break was produced by a break in the underlying electron distribution. This may serve as a diagnostic of whether gamma-gamma absorption is relevant in sources such as gamma-ray bursts and blazars where the gamma-ray emission may be intrinsically highly polarized.

In this paper, we discuss the light-curve features of various flaring scenarios in a time-dependent leptonic model for low-frequency-peaked blazars. The quasar 3C273 is used as an illustrative example. Our code takes into account Fermi-II acceleration and all relevant electron cooling terms, including the external radiation fields generally found to be important in the modeling of the SEDs of FSRQs, as well as synchrotron self absorption and gamma-gamma pair-production. General parameters are constrained through a fit to the average spectral energy distribution (SED) of the blazar by numerically solving the time-dependent Fokker-Planck equation for the electron evolution in a steady-state situation. We then apply perturbations to several input parameters (magnetic field, particle injection luminosity, acceleration time scale) to simulate flaring events and compute time-dependent SEDs and light curves in representative energy bands (radio, optical, X-rays, gamma-rays). Time lags between different bands are evaluated using a discrete cross correlation analysis. We find that Fermi-II acceleration has a significant effect on the distributions and that flaring events caused by increased acceleration efficiency of the Fermi II process will produce a correlation between the radio, optical and gamma-ray bandpasses, but an anti-correlation between these three bandpasses and the X-ray band, with the X-rays lagging behind the variations in other bands by up to several hours.

We perform time-dependent, spatially-resolved simulations of blazar emission to evaluate several flaring scenarios related to magnetic-field amplification and enhanced particle acceleration. The code explicitly accounts for light-travel-time effects and is applied to flares observed in the flat spectrum radio quasar (FSRQ) PKS 0208-512, which show optical/{\gamma}-ray correlation at some times, but orphan optical flares at other times. Changes in both the magnetic field and the particle acceleration efficiency are explored as causes of flares. Generally, external Compton emission appears to describe the available data better than a synchrotron self-Compton scenario, and in particular orphan optical flares are difficult to produce in the SSC framework. X-ray soft-excesses, {\gamma}-ray spectral hardening, and the detections at very high energies of certain FSRQs during flares find natural explanations in the EC scenario with particle acceleration change. Likewise, optical flares with/without {\gamma}-ray counterparts can be explained by different allocations of energy between the magnetization and particle acceleration, which may be related to the orientation of the magnetic field relative to the jet flow. We also calculate the degree of linear polarization and polarization angle as a function of time for a jet with helical magnetic field. Tightening of the magnetic helix immediately downstream of the jet perturbations, where flares occur, can be sufficient to explain the increases in the degree of polarization and a rotation by >= 180 degree of the observed polarization angle, if light-travel-time effects are properly considered.

We extend our approach of modeling spectral energy distribution (SED) and lightcurves of blazars to include external Compton (EC) emission due to inverse Compton scattering of an external anisotropic target radiation field. We describe the time-dependent impact of such seed photon fields on the evolution of multifrequency emission and spectral variability of blazars using a multi-zone time-dependent leptonic jet model, with radiation feedback, in the internal shock model scenario. We calculate accurate EC-scattered high-energy spectra produced by relativistic electrons throughout the Thomson and Klein-Nishina regimes. We explore the effects of varying the contribution of (1) a thermal Shakura-Sunyaev accretion disk, (2) a spherically symmetric shell of broad-line clouds, the broad line region (BLR), and (3) a hot infrared emitting dusty torus (DT), on the resultant seed photon fields. We let the system evolve to beyond the BLR and within the DT and study the manifestation of the varying target photon fields on the simulated SED and lightcurves of a typical blazar. The calculations of broadband spectra include effects of gamma-gamma absorption as gamma-rays propagate through the photon pool present inside the jet due to synchrotron and inverse Compton processes, but neglect gamma-gamma absorption by the BLR and DT photon fields outside the jet. Thus, our account of gamma-gamma absorption is a lower limit to this effect. Here, we focus on studying the impact of parameters relevant for EC processes on high-energy (HE) emission of blazars.

We present a detailed analysis of time- and energy-dependent synchrotron polarization signatures in a shock-in-jet model for gamma-ray blazars. Our calculations employ a full 3D radiation transfer code, assuming a helical magnetic field throughout the jet. The code considers synchrotron emission from an ordered magnetic field, and takes into account all light-travel-time and other relevant geometric effects, while the relevant synchrotron self-Compton and external Compton effects are taken care of with the 2D MCFP code. We consider several possible mechanisms through which a relativistic shock propagating through the jet may affect the jet plasma to produce a synchrotron and high-energy flare. Most plausibly, the shock is expected to lead to a compression of the magnetic field, increasing the toroidal field component and thereby changing the direction of the magnetic field in the region affected by the shock. We find that such a scenario leads to correlated synchrotron + SSC flaring, associated with substantial variability in the synchrotron polarization percentage and position angle. Most importantly, this scenario naturally explains large PA rotations by > 180 deg., as observed in connection with gamma-ray flares in several blazars, without the need for bent or helical jet trajectories or other non-axisymmetric jet features.

Diffusive shock acceleration (DSA) at relativistic shocks is likely to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud AGN. An important recent development for blazar science is the ability of Fermi-LAT data to pin down the power-law index of the high energy portion of emission in these sources, and therefore also the index of the underlying non-thermal particle population. This diagnostic potential was not possible prior to Fermi launch, when gamma-ray information was dominated by the highly-absorbed TeV band. This paper highlights how multiwavelength spectra including X-ray band and Fermi data can be used to probe diffusive acceleration in relativistic, oblique, MHD shocks in blazar jets. The spectral index of the non-thermal particle distributions resulting from Monte Carlo simulations of DSA, and the fraction of thermal particles accelerated to non-thermal energies, depend sensitively on the particles' mean free path scale, and also on the magnetic field obliquity to the shock normal. We investigate self-consistently the radiative synchrotron/Compton signatures of the resulting thermal and non-thermal particle distributions. Important constraints on the frequency of particle scattering and the level of field turbulence are identified for the blazar AO 0235+164. The possible interpretation that turbulence levels decline with remoteness from jet shocks, and a significant role for non-gyroresonant diffusion, are discussed.

We present Particle-in-Cell simulation results of relativistic shear flows for hybrid positron-electron-ion plasmas and compare to those for pure e+e- and pure e-ion plasmas. Among the three types of relativistic shear flows, we find that only hybrid shear flow is able to energize the electrons to form a high-energy spectral peak plus a hard power-law tail. Such electron spectra are needed to model the observational properties of gamma-ray bursts.

We report on the VERITAS observations of the high-frequency peaked BL Lac object 1ES 1959+650 in the period 2007-2011. This source is detected at TeV energies by VERITAS at 16.4 standard deviation (sigma) significance in 7.6 hours of observation in a low flux state. A multiwavelength spectral energy distribution (SED) is constructed from contemporaneous data from VERITAS, Fermi-LAT, RXTE PCA, and Swift UVOT. Swift XRT data is not included in the SED due to a lack of simultaneous observations with VERITAS. In contrast to the orphan gamma-ray flare exhibited by this source in 2002, the X-ray flux of the source is found to vary by an order of magnitude, while other energy regimes exhibit less variable emission. A quasi-equilibrium synchrotron self-Compton model with an additional external radiation field is used to describe three SEDs corresponding to the lowest, highest, and average X-ray states. The variation in the X-ray spectrum is modeled by changing the electron injection spectral index, with minor adjustments of the kinetic luminosity in electrons. This scenario produces small-scale flux variability of order >~2 in the HE (E>1 MeV) and VHE (E>100 GeV) gamma-ray regimes, which is corroborated by the Fermi-LAT, VERITAS, and Whipple 10m telescope light curves.

We present a theoretical analysis of the expected X-ray and gamma-ray polarization signatures resulting from synchrotron self-Compton emission in leptonic models, compared to the polarization signatures from proton synchrotron and cascade synchrotron emission in hadronic models for blazars. Source parameters resulting from detailed spectral-energy-distribution modeling are used to calculate photon-energy-dependent upper limits on the degree of polarization, assuming a perfectly organized, mono-directional magnetic field. In low-synchrotron-peaked blazars, hadronic models exhibit substantially higher maximum degrees of X-ray and gamma-ray polarization than leptonic models, which may be within reach for existing X-ray and gamma-ray polarimeters. In high-synchrotron-peaked blazars (with electron-synchrotron-dominated X-ray emission), leptonic and hadronic models predict the same degree of X-ray polarization, but substantially higher maximum gamma-ray polarization in hadronic models than leptonic ones. These predictions are particularly relevant in view of the new generation of balloon-borne X-ray polarimeters (and possibly GEMS, if revived), and the ability of Fermi-LAT to measure gamma-ray polarization at < 200 MeV. We suggest observational strategies combining optical, X-ray, gamma-ray polarimetry to determine the degree of ordering of the magnetic field and to distinguish between leptonic and hadronic high-energy emission.

We describe new implementations of leptonic and hadronic models for the broadband emission from relativistic jets in AGN in a temporary steady state. For the leptonic model, a temporary equilibrium between particle injection/acceleration, radiative cooling, and escape from a spherical emission region is evaluated, and the self-consistent radiative output is calculated. For the hadronic model, a temporary equilibrium between particle injection/acceleration, radiative and adiabatic cooling, and escape is evaluated for both primary electrons and protons. A new, semi-analytical method to evaluate the radiative output from cascades initiated by internal gamma-gamma pair production is presented. We use our codes to fit snap-shot spectral energy distributions of a representative set of Fermi-LAT detected blazars. We find that the leptonic model provides acceptable fits to the SEDs of almost all blazars with parameters close to equipartition between the magnetic field and the relativistic electron population. However, the hard gamma-ray spectrum of AO 0235+164, in contrast to the very steep IR-optical-UV continuum, poses a severe problem for the leptonic model. If charge neutrality in leptonic models is provided by cold protons, the kinetic energy carried by the jet should be dominated by protons. We find satisfactory representations of the snapshot SEDs of most blazars in our sample with the hadronic model presented here. However, in the case of two quasars the characteristic break at a few GeV energies can not be well modelled. All of our hadronic model fits require powers in relativistic protons in the range L_p ~ 1e47 - 1e49 erg/s.

The synergy between the Fermi-LAT and ground-based Cherenkov telescope arrays gives us the opportunity for the first time to characterize the high-energy emission from blazars over 5 decades in energy, from 100 MeV to 10 TeV. In this study, we perform a Fermi-LAT spectral analysis for TeV-detected blazars and combine it with archival TeV data. We examine the observational properties in the gamma-ray band of our sample of TeV-detected blazars and compare the results with X-ray and GeV-selected populations. The spectral energy distributions (SEDs) that result from combining Fermi-LAT and ground-based spectra are studied in detail. Simple parameterizations such as a power-law function do not always reproduce the high-energy SEDs, where spectral features that could indicate intrinsic absorption are observed.

Densely time sampled multi-frequency flux measurements of the extreme BL Lac object S5 0716+714 over the past three years allow us to study its broad-band variability, and the detailed underlying physics, with emphasis on the location and size of the emitting regions and the evolution with time. We study the characteristics of some prominent mm-/gamma-ray flares in the context of the shock-in-jet model and investigate the location of the high energy emission region. The rapid rise and decay of the radio flares is in agreement with the formation of a shock and its evolution, if a geometrical variation is included in addition to intrinsic variations of the source. We find evidence for a correlation between flux variations at gamma-ray and radio frequencies. A two month time-delay between gamma-ray and radio flares indicates a non-cospatial origin of gamma-rays and radio flux variations in S5 0716+714.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

Using 2.5-dimensional Particle-in-Cell simulations, we study the kinetic physics of relativistic shear flow boundary in collisionless electron-positron (e+e-) plasmas. We find efficient magnetic field generation and particle energization at the shear boundary, driven by streaming instabilities across the shear interface and sustained by the shear flow. Nonthermal, anisotropic high-energy particles are accelerated across field lines to produce a power-law tail, truncated at energies below the shear Lorentz factor. These results have important implications for the dissipation and radiation of jets in blazars, gamma-ray bursts and other relativistic outflows.

We report on the development of a numerical code to calculate the angle-dependent synchrotron + synchrotron self-Compton radiation from relativistic jet sources with partially ordered magnetic fields and anisotropic particle distributions. Using a multi-zone radiation transfer approach, we can simulate magnetic-field configurations ranging from perfectly ordered (unidirectional) to randomly oriented (tangled). We demonstrate that synchrotron self-Compton model fits to the spectral energy distributions (SEDs) of extragalactic jet sources may be possible with a wide range of magnetic-field values, depending on their orientation with respect to the jet axis and the observer. This is illustrated with the example of a spectral fit to the SED of Mrk~421 from multi-wavelength observations in 2006, where acceptable fits are possible with magnetic-field values varying within a range of an order of magnitude for different degrees of B-field alignment and orientation.

In this paper, we review the prospects for studies of active galactic nuclei (AGN) using the envisioned future Cherenkov Telescope Array (CTA). This review focuses on jetted AGN, which constitute the vast majority of AGN detected at gamma-ray energies. Future progress will be driven by the planned lower energy threshold for very high energy (VHE) gamma-ray detections to ~10 GeV and improved flux sensitivity compared to current-generation Cherenkov Telescope facilities. We argue that CTA will enable substantial progress on gamma-ray population studies by deepening existing surveys both through increased flux sensitivity and by improving the chances of detecting a larger number of low-frequency peaked blazars because of the lower energy threshold. More detailed studies of the VHE gamma-ray spectral shape and variability might furthermore yield insight into unsolved questions concerning jet formation and composition, the acceleration of particles within relativistic jets, and the microphysics of the radiation mechanisms leading to the observable high-energy emission. The broad energy range covered by CTA includes energies where gamma-rays are unaffected from absorption while propagating in the extragalactic background light (EBL), and extends to an energy regime where VHE spectra are strongly distorted. This will help to reduce systematic effects in the spectra from different instruments, leading to a more reliable EBL determination, and hence will make it possible to constrain blazar models up to the highest energies with less ambiguity.

Recent detections of very-high-energy (VHE, E > 100 GeV) gamma-ray blazars which do not belong to the high frequency peaked BL Lac (HBL) class, suggest that gamma-gamma absorption and pair cascading might occur in those objects. In the presence of even weak magnetic fields, these Compton-supported pair cascades will be deflected and contribute to the Fermi gamma-ray flux of radio galaxies. We demonstrate that, in this scenario, the magnetic field can not be determined from a fit of the cascade emission to the gamma-ray spectrum alone, and the degeneracy can only be lifted if the synchrotron emission from the cascades is observed as well. We illustrate this fact with the example of NGC 1275. We point out that the cascade synchrotron emission may produce spectral features reminiscent of the big blue bump observed in the spectral energy distributions of several blazars, and illustrate this idea for 3C 279.

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.

[abridged] We present results of modeling the SED and multiwavelength variability of the bright FSRQ PKS1510-089 with our time-dependent multizone Monte Carlo/Fokker-Planck code (Chen et al. 2001). As primary source of seed photons for inverse Compton scattering, we consider radiation from the broad line region (BLR), from the molecular torus, and the local synchrotron radiation (SSC). Different scenarios are assessed by comparing simulated light curves and SEDs with one of the best flares by PKS1510-089, in March 2009. The time-dependence of our code and its correct handling of light travel time effects allow us to fully take into account the effect of the finite size of the active region, and in turn to fully exploit the information carried by time resolved observed SEDs, increasingly available since the launch of Fermi. We confirm that the spectrum adopted for the external radiation has an important impact on the modeling of the SED, in particular for the lower energy end of the Compton component, observed in the X-ray band, which in turn is one of the most critical bands to assess the differences between EC and SSC emission. In the context of the scenario presented here, where the flaring is caused by the increase of the number of relativistic electrons ascribed to the effect of the interaction of a portion of the jet (blob) with a shock, we can not firmly discriminate the three main scenarios for gamma-ray emission. However, results show clearly the differences produced by a more realistic treatment of the emitting source in the shape of SEDs and their time variability over relevant, observable time-scales, and demonstrate the crucial importance of time-dependent multi-zone models to advance our understanding of the physics of these sources, by taking full advantage of the wealth of information offered by the high quality data of current multiwavelength campaigns.

In this review, recent progress in theoretical models for the broadband (radio through gamma-ray) emission from blazars are summarized. The salient features of both leptonic and hadronic models are reviewed. I present sample modeling results of spectral energy distributions (SEDs) of different types of Fermi-detected blazars along the traditional blazar sequence, using both types of models. In many cases, the SEDs of high-frequency peaked blazars (HBLs) have been found to be well represented by simple synchrotron + synchrotron self-Compton (SSC) models. However, a few HBLs recently discovered as very-high-energy (VHE) gamma-ray emitters by VERITAS are actually better represented by either external-Compton or hadronic models. Often, spectral modeling with time-independent single-zone models alone is not sufficient to constrain models, as both leptonic and lepto-hadronic models are able to provide acceptable fits to the overall SED. This degeneracy can be lifted by considering further constraints from spectral variability. Recent developments of time-dependent and inhomogeneous blazar models will be discussed, including detailed numerical simulations as well as a semi-analytical approach to the time-dependent radiation signatures of shock-in-jet models.

We report on the discovery of high-energy (HE; E > 0.1 GeV) and very high-energy (VHE; E > 100 GeV) gamma-ray emission from the high-frequency-peaked BL Lac object RBS 0413. VERITAS, a ground-based gamma-ray observatory, detected VHE gamma rays from RBS 0413 with a statistical significance of 5.5 standard deviations (sigma) and a gamma-ray flux of (1.5 \pm 0.6stat \pm 0.7syst) \times 10^(-8) photons m^(-2) s^(-1) (\sim 1% of the Crab Nebula flux) above 250 GeV. The observed spectrum can be described by a power law with a photon index of 3.18 \pm 0.68stat \pm 0.30syst. Contemporaneous observations with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope detected HE gamma rays from RBS 0413 with a statistical significance of more than 9 sigma, a power-law photon index of 1.57 \pm 0.12stat +0.11sys -0.12sys and a gamma-ray flux between 300 MeV and 300 GeV of (1.64 \pm 0.43stat +0.31sys -0.22sys) \times 10^(-5) photons m^(-2) s^(-1). We present the results from Fermi-LAT and VERITAS, including a spectral energy distribution modeling of the gamma-ray, quasi-simultaneous X-ray (Swift-XRT), ultraviolet (Swift-UVOT) and R-band optical (MDM) data. We find that, if conditions close to equipartition are required, both the combined synchrotron self-Compton/external-Compton and the lepto-hadronic models are preferred over a pure synchrotron self-Compton model.

The discovery of very-high-energy (VHE, E > 100 GeV) gamma-ray emission from intermediate- and low-frequency peaked blazars suggests that gamma-gamma absorption and pair cascading might occur in those objects. In previous papers, we investigated the Compton emission from VHE gamma-ray induced pair cascades, deflected by moderate magnetic fields, in a largely model-independent way, and demonstrated that this emission can explain the Fermi fluxes and spectra of the radio galaxies Cen A and NGC 1275. In this paper, we describe a generalization of our Monte-Carlo cascade code to include the angle-dependent synchrotron output from the cascades, allowing for the application to situations with non-negligible magnetic fields, leading to potentially observable synchrotron signatures, but still not dominating the radiative energy loss of cascade particles. We confirm that the synchrotron radiation from the cascades in NGC 1275 and Cen A are negligible for the parameters used in our previous works. We demonstrate that the magnetic field can not be determined from a fit of the cascade emission to the gamma-ray spectrum alone, and the degeneracy can only be lifted if the synchrotron emission from the cascades is observed as well. We illustrate this fact with the example of NGC 1275. We point out that the cascade synchrotron emission may produce spectral features in the same energy range in which the big blue bump is observed in the spectral energy distributions of several blazars, and may make a non-negligible contribution to this feature. We illustrate this idea with the example of 3C 279.

We present the results of our optical monitoring of the BL Lac object S5 0716+714 on seven nights in 2006 December. The monitoring was carried out simultaneously at three optical wavelengths with a novel photometric system. The object did not show large-amplitude internight variations during this period. Intranight variations were observed on four nights and probably on one more. Strong bluer-when-brighter chromatism was detected on both intranight and internight timescales. The intranight variation amplitude decreases in the wavelength sequence of B', R', and V'. Cross correlation analyses revealed that the variability at the $B'$ and $V'$ bands lead that at the $R'$ band by about 30 minutes on one night.

Origin of extragalactic gamma-ray background (EGRB) is still a matter of debate. EGRB can either have truly diffuse or unresolved discrete point sources origin. Majority of the Fermi and EGRET detected identified sources are blazar. So, they are expected to be a significant contributors to the EGRB. In order to estimate their contribution to the EGRB one needs to construct their luminosity functions. Here, we examine and construct the evolution and luminosity function of blazars using first LAT AGN catalog. We consider both pure luminosity and pure density evolution models for flat spectrum radio quasars (FSRQs). We also describe the methodology to estimate the contribution from misaligned blazars to the EGRB.

We present multiwavelength spectral analyses of two Fermi-LAT blazars, OJ 287 and 3C 279, that are part of the Boston University multiwaveband polarization program. The data have been compiled from observations with Fermi, RXTE, the VLBA, and various ground-based optical and radio telescopes. We simulate the dynamic spectral energy distributions (SEDs) within the framework of a multi-slice, time-dependent leptonic jet model for blazars, with radiation feedback, in the internal shock scenario. We use the physical jet parameters obtained from the VLBA monitoring to guide our modeling efforts. We discuss the role of intrinsic parameters and the interplay between synchrotron and inverse Compton radiation processes responsible for producing the resultant SEDs.

The synergy between the Fermi-LAT and ground-based Cherenkov telescope arrays gives us the opportunity for the first time to characterize the high energy emission (100 MeV - 10 TeV) from more than 30 blazars. In this study we performed a Fermi-LAT spectral analysis for all TeV-detected blazars and combined it with archival TeV spectra. Our results for low synchrotron-peaked BL Lacs (LBL) show hints of absorption features in the GeV band that could be interpreted as internal opacity at the source. We note that simple or broken power laws cannot describe all the observed GeV-TeV spectra and more complex spectral shapes seem required.

Aims. We use a sample of 83 core-dominated active galactic nuclei (AGN) selected from the MOJAVE (Monitoring of Jets in AGN with VLBA Experiments) radio-flux-limited sample and detected with the Fermi Large Area Telescope (LAT) to study the relations between non-simultaneous radio, optical, and gamma-ray measurements. Methods. We perform a multi-band statistical analysis to investigate the relations between the emissions in different bands and reproduce these relations by modeling of the spectral energy distributions of blazars. Results. There is a significant correlation between the gamma-ray luminosity and the optical nuclear and radio (15 GHz) luminosities of blazars. We report a well defined positive correlation between the gamma-ray luminosity and the radio-optical loudness for quasars and BL Lacertae type objects (BL Lacs). A strong positive correlation is found between the radio luminosity and the gamma-ray-optical loudness for quasars, while a negative correlation between the optical luminosity and the gamma-ray-radio loudness is present for BL Lacs. Modeling of these correlations with a simple leptonic jet model for blazars indicates that variations of the accretion disk luminosity (and hence the jet power) is able to reproduce the trends observed in most of the correlations. To reproduce all observed correlations, variations of several parameters, such as the accretion power, jet viewing angle, Lorentz factor, and magnetic field of the jet, are required.

(abridged) We present a new time-dependent multi-zone radiative transfer code and its application to study the SSC emission of Mrk 421. The code couples Fokker-Planck and Monte Carlo methods, in a 2D geometry. For the first time all the light travel time effects (LCTE) are fully considered, along with a proper treatment of Compton cooling, which depends on them. We study a set of simple scenarios where the variability is produced by injection of relativistic electrons as a `shock front' crosses the emission region. We consider emission from two components, with the second one either being pre-existing and co-spatial and participating in the evolution of the active region, or spatially separated and independent, only diluting the observed variability. Temporal and spectral results of the simulation are compared to the multiwavelength observations of Mrk 421 in March 2001. We find parameters that can adequately fit the observed SEDs and multiwavelength light curves and correlations. There remain however a few open issues, most notably: i) systematic soft intra-band X-ray lags. ii) The quadratic correlation between the TeV and X-ray flux during the flare decay has not been reproduced. These features are among those affected by the spatial extent and geometry of the source. The difficulty of producing hard X-ray lags is exacerbated by a bias towards soft lags caused by the combination of energy dependent radiative cooling time-scales and LCTE. About the second emission component, our results strongly favor the scenario where it is co-spatial and it participates in the flare evolution, suggesting that different phases of activity may occur in the same region. The cases presented in this paper represent only an initial study, and despite their limited scope they make a strong case for the need of true time-dependent and multi-zone modeling.

We present multiwavelength observations of the ultraluminous blazar-type radio loud quasar PKS 0528+134 in quiescence during the period July to December 2009. Significant flux variability on a time scale of several hours was found in the optical regime, accompanied by a weak trend of spectral softening with increasing flux. We suggest that this might be the signature of a contribution from the accretion disk at the blue end of the optical spectrum. The optical flux is weakly polarized with rapid variations of the degree and direction of polarization, while the polarization of the 43 GHz radio core remains steady. Optical spectropolarimetry suggests a trend of increasing degree of polarization with increasing wavelength, providing additional evidence for an accretion disc contribution towards the blue end of the optical spectrum. We constructed four SEDs indicating that even in the quiescent state, the bolometric luminosity of PKS 0528+134 is dominated by its gamma-ray emission. A leptonic single-zone jet model produced acceptable fits to the SEDs with contributions to the high-energy emission from synchrotron self-Compton radiation and Comptonization of direct accretion disk emission. Fit parameters close to equipartition were obtained. The moderate variability on long time scales implies the existence of on-going particle acceleration, while the observed optical polarization variability seems to point towards a turbulent acceleration process. Turbulent particle acceleration at stationary features along the jet therefore appears to be a viable possibility for the quiescent state of PKS 0528+134.

The growing number of extragalactic high-energy (HE, E > 100 MeV) and very-high-energy (VHE, E > 100 GeV) gamma-ray sources that do not belong to the blazar class suggests that VHE gamma-ray production may be a common property of most radio-loud Active Galactic Nuclei (AGN). In a previous paper, we have investigated the signatures of Compton-supported pair cascades initiated by VHE gamma-ray absorption in monochromatic radiation fields, dominated by Ly-alpha line emission from the Broad Line Region. In this paper, we investigate the interaction of nuclear VHE gamma-rays with the thermal infrared radiation field from a circumnuclear dust torus. Our code follows the spatial development of the cascade in full 3-dimensional geometry. We provide a model fit to the broadband SED of the dust-rich, gamma-ray loud radio galaxy Cen A and show that typical blazar-like jet parameters may be used to model the broadband SED, if one allows for an additional cascade contribution to the Fermi gamma-ray emission.

We describe the time-dependent radiation transfer in blazar jets, within the internal shock model. We assume that the central engine, which consists of a black hole and an accretion disk, spews out relativistic shells of plasma with different velocity, mass, and energy. We consider a single inelastic collision between a faster (inner) and a slower (outer) moving shell. We study the dynamics of the collision and evaluate the subsequent emission of radiation via the synchrotron and synchrotron self Compton (SSC) processes after the interaction between the two shells has begun. The collision results in the formation of a forward shock (FS) and a reverse shock (RS) that convert the ordered bulk kinetic energy of the shells into magnetic field energy and accelerate the particles, which then radiate. We assume a cylindrical geometry for the emission region of the jet. We treat the self-consistent radiative transfer by taking into account the inhomogeneity in the photon density throughout the region. In this paper, we focus on understanding the effects of varying relevant input parameters on the simulated spectral energy distribution (SED) and spectral variability patterns.

We present results of five years of optical (UBVRI) observations of the very-high-energy gamma-ray blazar 1ES 1011+496 at the MDM Observatory. We calibrated UBVRI magnitudes of five comparison stars in the field of the object. Most of our observations were done during moderately faint states of 1ES 1011+496 with R > 15.0. The light curves exhibit moderate, closely correlated variability in all optical wavebands on time scales of a few days. A cross-correlation analysis between optical bands does not show significant evidence for time lags. We find a positive correlation (Pearson's r = 0.57; probability of non-correlation P(>r) ~ 4e-8) between the R-band magnitude and the B - R color index, indicating a bluer-when-brighter trend. Snap-shot optical spectral energy distributions (SEDs) exhibit a peak within the optical regime, typically between the V and B bands. We find a strong (r = 0.78; probability of non-correlation P (>r) ~ 1e-15) positive correlation between the peak flux and the peak frequency, best fit by a relation $\nu F_{\nu}^{\rm pk} \propto \nu_{\rm pk}^k$ with k = 2.05 +/- 0.17. Such a correlation is consistent with the optical (synchrotron) variability of 1ES 1011+496 being primarily driven by changes in the magnetic field.

In this review, recent progress in theoretical models for blazar emission will be summarized. The salient features of both leptonic and lepto-hadronic approaches to modeling blazar spectral energy distributions will be reviewed. I will present sample modeling results of spectral energy distributions (SEDs) of different types of blazars along the blazar sequence, including Fermi high-energy gamma-ray data, using both types of models. Special emphasis will be placed on the implications of the recent very-high-energy (VHE) gamma-ray detections of non-traditional VHE gamma-ray blazars, including intermediate and low-frequency-peaked BL Lac objects and even flat-spectrum radio quasars. Due to the featureless optical spectra of BL Lac objects, the redshifts of several BL Lacs remain unknown. I will briefly discuss possible constraints on their redshift using spectral modeling of their SED including Fermi + ground-based VHE gamma-ray data. It will be shown that in some cases, spectral modeling with time-independent single-zone models alone is not sufficient to constrain models, as both leptonic and lepto-hadronic models are able to provide acceptable fits to the overall SED. Subsequently, recent developments of time-dependent and inhomogeneous blazar models will be discussed, including detailed numerical simulations as well as a semi-analytical approach to the time-dependent radiation signatures of shock-in-jet models.

Recent blazar detections by HESS, MAGIC, and VERITAS suggest that very-high-energy (VHE, E > 100 GeV) gamma-rays may be produced in most, if not all, types of blazars, including those that possess intense circumnuclear radiation fields. In this paper, we investigate the interaction of nuclear VHE gamma-rays with the circumnuclear radiation fields through gamma-gamma absorption and pair production, and the subsequent Compton-supported pair cascades. We have developed a Monte-Carlo code to follow the spatial development of the cascade in full 3-dimensional geometry, and calculate the radiative output due to the cascade as a function of viewing angle with respect to the primary VHE gamma-ray beam (presumably the jet axis of the blazar). We show that even for relatively weak magnetic fields, the cascades can be efficiently isotropized, leading to substantial off-axis cascade emission peaking in the Fermi energy range at detectable levels for nearby radio galaxies. We demonstrate that this scenario can explain the Fermi flux and spectrum of the radio galaxy NGC 1275.

We report on 18 months of multiwavelength observations of the blazar 3C 454.3 (Crazy Diamond) carried out in July 2007-January 2009. We show the results of the AGILE campaigns which took place on May-June 2008, July-August 2008, and October 2008-January 2009. During the May 2008-January 2009 period, the source average flux was highly variable, from an average gamma-ray flux F(E>100MeV) > 200E-8 ph/cm2/s in May-June 2008, to F(E>100MeV)~80E-8 ph/cm2/s in October 2008-January 2009. The average gamma-ray spectrum between 100 MeV and 1 GeV can be fit by a simple power law (Gamma_GRID ~ 2.0 to 2.2). Only 3-sigma upper limits can be derived in the 20-60 keV energy band with Super-AGILE. During July-August 2007 and May-June 2008, RXTE measured a flux of F(3-20 keV)= 8.4E-11 erg/cm2/s, and F(3-20 keV)=4.5E-11 erg/cm2/s, respectively and a constant photon index Gamma_PCA=1.65. Swift/XRT observations were carried out during all AGILE campaigns, obtaining a F(2-10 keV)=(0.9-7.5)E-11 erg/cm2/s and a photon index Gamma_XRT=1.33-2.04. BAT measured an average flux of ~5 mCrab. GASP-WEBT monitored 3C 454.3 during the whole 2007-2008 period from the radio to the optical. A correlation analysis between the optical and the gamma-ray fluxes shows a time lag of tau=-0.4 days. An analysis of 15 GHz and 43 GHz VLBI core radio flux observations shows an increasing trend of the core radio flux, anti- correlated with the higher frequency data. The modeling SEDs, and the behavior of the long-term light curves in different energy bands, allow us to compare the jet properties during different emission states, and to study the geometrical properties of the jet on a time-span longer than one year.

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.

The MAGIC collaboration recently reported the detection of the quasar 3C279 at > 100 GeV gamma-ray energies. Here we present simultaneous optical (BVRI) and X-ray (RXTE PCA) data from the day of the VHE detection and discuss the implications of the snap-shot spectral energy distribution for jet models of blazars. A one-zone synchrotron-self-Compton origin of the entire SED, including the VHE gamma-ray emission can be ruled out. The VHE emission could, in principle, be interpreted as Compton upscattering of external radiation (e.g., from the broad-line regions). However, such an interpretation would require either an unusually low magnetic field of B ~ 0.03 G or an unrealistically high Doppler factor of Gamma ~ 140. In addition, such a model fails to reproduce the observed X-ray flux. This as well as the lack of correlated variability in the optical with the VHE gamma-ray emission and the substantial gamma-gamma opacity of the BLR radiation field to VHE gamma-rays suggests a multi-zone model. In particular, an SSC model with an emission region far outside the BLR reproduces the simultaneous X-ray -- VHE gamma-ray spectrum of 3C279. Alternatively, a hadronic model is capable of reproducing the observed SED of 3C279 reasonably well. However, the hadronic model requires a rather extreme jet power of L_j ~ 10^{49} erg s^{-1}, compared to a requirement of L_j ~ 2 X 10^{47} erg s^{-1} for a multi-zone leptonic model.

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.

Prompted by a high optical state in September 2007, the Whole Earth Blazar Telescope (WEBT) consortium organized an intensive optical, near-IR (JHK) and radio observing campaign on the intermediate BL Lac object 3C 66A throughout the fall and winter of 2007 -- 2008. The source remained in a high optical state throughout the observing period and exhibited several bright flares on time scales of ~ 10 days. This included an exceptional outburst around September 15 - 20, 2007, reaching a peak brightness at R ~ 13.4. Our campaign revealed microvariability with flux changes up to |dR/dt| ~ 0.02 mag/hr. Our observations do not reveal evidence for systematic spectral variability or spectral lags. We infer a value of the magnetic field in the emission region of B ~ 19 e_B^{2/7} \tau_h^{-6/7} D_1^{13/7} G. From the lack of systematic spectral variability, we can derive an upper limit on the Doppler factor, D <= 28 \tau_h^{-1/8} e_B^{3/16}. This is in agreement with superluminal motion measurements of \beta_{app} \le 27 and argues against models with very high Lorentz factors of \Gamma > 50, required for a one-zone SSC interpretation of some high-frequency-peaked BL Lac objects detected at TeV gamma-ray energies.

A recent optical monitoring campaign on the prominent quasar 3C279 revealed a period of a remarkably clean exponential decay of BVRI fluxes with time, with a time constant of 12.8 d, over about 14 days. This is clearly too long to be associated with radiative cooling. Here we propose that this may be the signature of deceleration of the synchrotron emitting jet component. We develop a model analogous to the relativistic blast wave model for GRBs, including radiative energy losses and radiation drag, to simulate the deceleration of a relativistically moving plasmoid in the moderately dense AGN environment. Synchrotron, SSC and external Compton emission are evaluated self-consistently. We show that the observed optical light curve decay can be successfully reproduced with this model. The decelerating plasmoid model predicts a delayed X-ray flare, about 2 - 3 weeks after the onset of the quasi-exponential light curve decay in the optical. A robust prediction of this model, which can be tested with Fermi and simultaneous optical monitoring, is that the peak in the gamma-ray light curve at ~ 100 MeV is expected to be delayed by a few days with respect to the onset of the optical decay, while the VHE gamma-rays are expected to track the optical light curve closely with a delay of at most a few hours.

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 detection of very high-energy gamma-ray emission from the intermediate-frequency-peaked BL Lacertae object W Comae (z=0.102) by VERITAS. The source was observed between January and April 2008. A strong outburst of gamma-ray emission was measured in the middle of March, lasting for only four days. The energy spectrum measured during the two highest flare nights is fit by a power-law and is found to be very steep, with a differential photon spectral index of Gamma = 3.81 +- 0.35_stat +- 0.34_syst. The integral photon flux above 200GeV during those two nights corresponds to roughly 9% of the flux from the Crab Nebula. Quasi-simultaneous Swift observations at X-ray energies were triggered by the VERITAS observations. The spectral energy distribution of the flare data can be described by synchrotron-self-Compton (SSC) or external-Compton (EC) leptonic jet models, with the latter offering a more natural set of parameters to fit the data.

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}$.

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.

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.

Context: BL Lacertae objects have been the targets for numerous recent multiwavelength campaigns, continuum spectral variability studies, and theoretical spectral and variability modeling. A meaningful interpretation of the results of such studies requires a reliable knowledge of the objects' redshifts; however, the redshifts for many are still unknown or uncertain. Aims: Therefore, we hope to determine or constrain the redshifts of six BL Lac objects with unknown or poorly known redshifts. Methods: Observations were made of these objects with the MDM 2.4 m Hiltner telescope. Although no spectral features were detected, and thus no redshifts could be measured, lower redshift limits were assigned to the objects based on the expected equivalent widths of absorption features in their host galaxies. Redshifts were also estimated for some objects by assuming the host galaxies are standard candles and using host galaxy apparent magnitudes taken from the literature. Results: The commonly used redshift of $z=0.102$ for 1219+285 is almost certainly wrong, while the redshifts of the other objects studied remain undetermined.

The quasar 3C279 was the target of an extensive multiwavelength monitoring campaign from January through April 2006, including an optical-IR-radio monitoring campaign by the Whole Earth Blazar Telescope (WEBT) collaboration. In this paper we focus on the results of the WEBT campaign. The source exhibited substantial variability of optical flux and spectral shape, with a characteristic time scale of a few days. The variability patterns throughout the optical BVRI bands were very closely correlated with each other. In intriguing contrast to other (in particular, BL Lac type) blazars, we find a lag of shorter- behind longer-wavelength variability throughout the RVB ranges, with a time delay increasing with increasing frequency. Spectral hardening during flares appears delayed with respect to a rising optical flux. This, in combination with the very steep IR-optical continuum spectral index of ~ 1.5 - 2.0, may indicate a highly oblique magnetic field configuration near the base of the jet. An alternative explanation through a slow (time scale of several days) acceleration mechanism would require an unusually low magnetic field of < 0.2 G, about an order of magnitude lower than inferred from previous analyses of simultaneous SEDs of 3C279 and other FSRQs with similar properties.

We discuss our radio (Australia Telescope Compact Array and Australian Long Baseline Array) and X-ray (XMM-Newton) monitoring observations of the unusual ultraluminous supernova SN 1978K in NGC 1313 at ~25 years after the explosion. SN 1978K is a rare example of a Type IIn supernova that has remained bright enough to have long-term X-ray and radio observations. The observations probe the dense medium that was ejected by the progenitor star prior to its explosion; the star might have been a Luminous Blue Variable. The radio imaging shows that the source remains compact, but it may be marginally resolved. The radio monitoring shows deviations from a smooth decay suggesting that inhomogeneities are present in the radio emitting region. It appears that a major change occurred in the mass-loss rate of the progenitor star shortly before the supernova event. The X-ray emission between 2000 and 2006 is consistent with the radiation coming from two shocks. All the X-ray data can be fit using the same model (with no systematic evolution or short-term variability) but this has a surprising requirement: the X-ray emitting regions have a very large abundance of helium. This would be consistent with the X-ray emitting shocks being located in a helium-rich layer that was ejected by the progenitor star, or helium-rich material was ejected in the supernova explosion. The unusual properties of the supernova motivated a search for an associated gamma-ray burst (GRB). We show that SN 1978K was inside the ~4 sigma error box of GRB 771029. If this association is correct, the GRB was exceptionally underluminous. However, the quality of the gamma-ray burst locations at that time was poor, and this is likely just a chance alignment.

The nature of ultraluminous X-ray sources in nearby galaxies is one of the major open questions in modern X-ray astrophysics. One possible explanation for these objects is an inhomogeneous, radiation dominated accretion disk around a $\sim 10 M_{\odot}$ black hole -- the so-called ``photon bubble'' model. While previous studies of this model have focused primarily on its radiation-hydrodynamics aspects, in this paper, we provide an analysis of its X-ray spectral (continuum and possible edge and line) characteristics. Compton reflection between high and low density regions in the disk may provide the key to distinguishing this model from others, such as accretion onto an intermediate mass black hole. We couple a Monte Carlo/Fokker-Planck radiation transport code with the XSTAR code for reflection to simulate the photon spectra produced in a photon bubble model for ULXs. We find that reflection components tend to be very weak and in most cases not observable, and make predictions for the shape of the high-energy Comptonizing spectra. In many cases the Comptonization dominates the spectra even down to $\sim$ a few keV. In one simulation, a $\sim 9 \kev$ feature was found, which may be considered a signature of photon bubbles in ULXs; furthermore, we make predictions of high energy power-laws which may be observed by future instruments.

The BL Lac object 3C 66A was observed in an extensive multiwavelength monitoring campaign from July 2003 till April 2004. The spectral energy distribution (SED) was measured over the entire electromagnetic spectrum, with flux measurements from radio to X-ray frequencies and upper limits in the very high energy (VHE) gamma-ray regime. Here, we use a time-dependent leptonic jet model to reproduce the SED and optical spectral variability observed during our multiwavelength campaign. Our model simulations could successfully reproduce the observed SED and optical light curves and predict an intrinsic cutoff value for the VHE gamma-ray emission at ~ 4 GeV. The effect of the optical depth due to the intergalactic infrared background radiation (IIBR) on the peak of the high-energy component of 3C 66A was found to be negligible. Also, the presence of a broad line region (BLR) in the case of 3C 66A may play an important role in the emission of gamma-ray photons when the emission region is very close to the central engine, but further out, the production mechanism of hard X-ray and gamma-ray photons becomes rapidly dominated by synchrotron self-Compton emission. We further discuss the possibility of an observable X-ray spectral variability pattern. The simulated results do not predict observable hysteresis patterns in the optical or soft X-ray regimes for major flares on multi-day time scales.

A detailed parameter study of the gamma-gamma absorption effects in LS 5039 is presented. For a range of plausible locations of the VHE gamma-ray emission region and the allowable range of viewing angles, the de-absorbed, intrinsic VHE gamma-ray spectra and total VHE photon fluxes and luminosities are calculated and compared to luminosity constraints based on wind accretion onto the compact object in LS 5039. It is found that (1) it is impossible to choose the viewing angle and location of the VHE emission region in a way that the intrinsic (deabsorbed) fluxes and spectra in superior and inferior conjunction are identical; consequently, the intrinsic VHE luminosities and spectral shapes must be fundamentally different in different orbital phases, (2) if the VHE luminosity is limited by wind accretion from the companion star and the system is viewed at an inclination angle of i > 40 deg., the emission is most likely beamed by a larger Doppler factor than inferred from the dynamics of the large-scale radio outflows, (3) the still poorly constrained viewing angle between the line of sight and the jet axis is most likely substantially smaller than the maximum of ~ 64 deg. inferred from the lack of eclipses.(4) Consequently, the compact object is more likely to be a black hole rather than a neutron star. (5) There is a limited range of allowed configurations for which the expected > 1 TeV neutrino flux would actually anti-correlate with the observed VHE gamma-ray emission. If hadronic models for the gamma-ray production in LS 5039 apply, a solid detection of the expected VHE neutrino flux and its orbital modulation with km^3 scale water-Cherenkov neutrino detectors might require the accumulation of data over more than 3 years.

The Galactic high-mass X-ray binary and jet source (microquasar) LSI +61 303 has recently been detected at TeV gamma-ray energies by the MAGIC telescope. We have applied a time-dependent leptonic jet model to the broadband spectral energy distribution and suggested (though not unambiguously detected) orbital modulation of the very high energy gamma-ray emission of this source. Our model takes into account time dependent electron injection and acceleration, and the adiabatic and radiative cooling of non-thermal electrons. It includes synchrotron, synchrotron self-Compton and external inverse Compton (with seed photons from the companion star), as well as gamma-gamma absorption of gamma-rays by starlight photons. The model can successfully reproduce the available multiwavelength observational data. Our best fit to the SED indicates that a magnetic field of B_0 ~ 5 X 10^3 G at ~ 10^3 R_g is required, and electrons need to be accelerated out to TeV energies (gamma_2 = 10^6) with a nonthermal injection spectrum with a spectral index of q = 1.7, indicating the operation of acceleration mechanisms beyond the standard first-order Fermi mechanism at relativistic or non-relativistic shocks. The orbital modulation of the VHE gamma-ray emission can be explained solely by the geometrical effect of changes in the relative orientation of the stellar companion with respect to the compact object and jet as it impacts the position and depth of the gamma-gamma absorption trough. Such a scenario predicts a trend of spectral hardening during VHE gamma-ray low orbital phases.

Blazars are the most violent steady/recurrent sources of high-energy gamma-ray emission in the known Universe. They are prominent emitters of electromagnetic radiation throughout the entire electromagnetic spectrum. The observable radiation most likely originates in a relativistic jet oriented at a small angle with respect to the line of sight. This review starts out with a general overview of the phenomenology of blazars, including results from a recent multiwavelength observing campaign on 3C279. Subsequently, issues of modeling broadband spectra will be discussed. Spectral information alone is not sufficient to distinguish between competing models and to constrain essential parameters, in particular related to the primary particle acceleration and radiation mechanisms in the jet. Short-term spectral variability information may help to break such model degeneracies, which will require snap-shot spectral information on intraday time scales, which may soon be achievable for many blazars even in the gamma-ray regime with the upcoming GLAST mission and current advances in Atmospheric Cherenkov Telescope technology. In addition to pure leptonic and hadronic models of gamma-ray emission from blazars, leptonic/hadronic hybrid models are reviewed, and the recently developed hadronic synchrotron mirror model for TeV gamma-ray flares which are not accompanied by simultaneous X-ray flares (``orphan TeV flares'') is revisited.

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)

In this paper we analyse five observations of the BL Lac object AO 0235+16 performed with the Chandra and XMM-Newton satellites during the years 2000-2005. In the February 2002 observation the source is found in a bright state and presents a steep X-ray spectrum, while in all the other epochs it is faint and the spectrum is hard. The soft X-ray spectrum appears to be strongly absorbed, likely by the intervening system at z=0.524, which also absorbs the optical-UV radiation. We find that models that consider spectral curvature are superior to single power law ones in fitting the X-ray spectrum. In particular, we favour a double power law model, which agrees with the assumption of a superposition of two different components in the X-ray domain. Both in the Chandra and in one of the XMM-Newton observations, a tentative detection of the redshifted Fe Kalpha emission line may suggest its origin from the inner part of an accretion disc. Thermal emission from this accretion disc might explain the UV-soft-X-ray bump that appears in the spectral energy distributions, when the X-ray spectra are complemented with the optical-UV data from the Optical Monitor onboard XMM-Newton. More likely, the bump can be interpreted in terms of an additional synchrotron component emitted from an inner region of the jet with respect to that where the lower-energy emission comes from. An inspection of the X-ray light curves reveals that intraday variability occurs only when the source is in a bright state.

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.

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.

We present lightcurves of the afterglow of GRB050502A, including very early data at t-t_{GRB} < 60s. The lightcurve is composed of unfiltered ROTSE-IIIb optical observations from 44s to 6h post-burst, R-band MDM observations from 1.6 to 8.4h post-burst, and PAIRITEL J H K_s observations from 0.6 to 2.6h post-burst. The optical lightcurve is fit by a broken power law, where t^{alpha} steepens from alpha = -1.13 +- 0.02 to alpha = -1.44 +- 0.02 at \~5700s. This steepening is consistent with the evolution expected for the passage of the cooling frequency nu_c through the optical band. Even in our earliest observation at 44s post-burst, there is no evidence that the optical flux is brighter than a backward extrapolation of the later power law would suggest. The observed decay indices and spectral index are consistent with either an ISM or a Wind fireball model, but slightly favor the ISM interpretation. The expected spectral index in the ISM interpretation is consistent within 1 sigma with the observed spectral index beta = -0.8 +- 0.1; the Wind interpretation would imply a slightly (~2 sigma) shallower spectral index than observed. A small amount of dust extinction at the source redshift could steepen an intrinsic spectrum sufficiently to account for the observed value of beta. In this picture, the early optical decay, with the peak at or below 4.7e14 Hz at 44s, requires very small electron and magnetic energy partitions from the fireball.

A challenge to standard leptonic SSC models are so-called orphan TeV flares, i.e. enhanced very high energy (VHE) gamma-ray emission without any contemporaneous X-ray flaring activity, that have recently been observed in TeV-blazars (e.g., 1ES 1959+650). In order to explain the orphan TeV flare of 1ES 1959+650 observed in June 2002, the co-called hadronic synchrotron mirror model has been developed. Here, relativistic protons are proposed to exist in the jet, and interact with reflected electron-synchrotron radiation of the precursor SSC flare. If the reflector is located in the cloud region, time shifts of several days are possible between the precursor and the orphan flare. The external photons, blueshifted in the comoving jet frame, are able to excite the \Delta(1232)-resonance when interacting with protons of Lorentz factors \gamma'_p~10^3-10^4. The decay products of this resonance include charged pions which, upon decay, give rise to neutrino production during the orphan flare. In this paper we calculate the expected neutrino emission for the June 4, 2002, orphan TeV flare of 1ES 1959+650. We compare our results with the recent observations of AMANDA-II of a neutrino event in spatial and temporal coincidence with the orphan flare of this blazar. We find that the expected neutrino signal from the hadronic synchrotron mirror model is insufficient to explain the observed neutrino event from the direction of 1ES 1959+650.

The BL Lac object 3C 66A was the target of an extensive multiwavelength monitoring campaign from July 2003 through April 2004, involving observations throughout the electromagnetic spectrum. Radio, infrared, and optical observations were carried out by the WEBT-ENIGMA collaboration. At higher energies, 3C 66A was observed in X-rays (RXTE), and at very-high-energy (VHE) gamma-rays (STACEE, VERITAS). In addition, the source has been observed with the VLBA at 9 epochs throughout the period September 2003 - December 2004, including 3 epochs contemporaneous with the core campaign.The WEBT campaign revealed microvariability with flux changes of ~ 5 % on time scales as short as \~ 2 hr. The source was in a relatively bright state, with several bright flares on time scales of several days. The SED indicates a $\nu F_{\nu}$ peak in the optical regime. A weak trend of optical spectral hysteresis with a trend of spectral softening throughout both the rising and decaying phases, has been found. The 3 - 10 keV X-ray flux of 3C 66A during the core campaign was historically high and its spectrum very soft, indicating that the low-frequency component of the broadband SED extends beyond ~ 10 keV. No significant X-ray flux and/or spectral variability was detected. STACEE and Whipple observations provided upper flux limits at > 150 GeV and > 390 GeV, respectively. The 22 and 43 GHz data from the 3 VLBA epochs between September 2003 and January 2004 indicate a smooth jet with only moderate internal structure. Evidence for superluminal motion (8.5 +/- 5.6 h^{-1} c) was found in only one out of 6 components. The radial radio brightness profile suggests a magnetic field decay $\propto r^{-1}$ and, thus, a predominantly perpendicular magnetic field.

Very-high-energy gamma-ray flares of TeV blazars are generally accompanied by simultaneous flaring activity in X-rays. The recent observations by the Whipple collaboration of an ``orphan'' TeV flare of 1ES 1959+650 (without simultaneous X-ray flare) is very hard to reconcile with the standard leptonic SSC model which is routinely very successfully employed to explain the SED and spectral variability of TeV blazars. In this paper, an alternative scenario is suggested in which the ``orphan'' TeV flare may originate from relativistic protons, interacting with an external photon field supplied by electron-synchrotron radiation reflected off a dilute reflector. While the external photons will be virtually ``invisible'' to the co-moving ultrarelativistic electrons in the jet due to Klein-Nishina effects, their Doppler boosted energy is high enough to excite the $\Delta$ resonance from relativistic protons with Lorentz factors of ~ 10^3 - 10^4. This model is capable of explaining the ``orphan'' TeV flare of 1ES 1959+650 with plausible parameters, thus constraining the number and characteristic energy of relativistic protons in the jet of this blazar.

BL Lacertae was the target of an extensive multiwavelength monitoring campaign in the second half of 2000. In this paper, we are using leptonic and hadronic jet models to fit the observed broadband spectra and spectral variability patterns. We start out with global spectral models. Subsequently, we investigate various flaring scenarios for comparison with the observed short-term variability. For our leptonic jet model, we find that the short-term variability, in particular the optical and X-ray spectral variability, can be best represented with a flaring scenario dominated by a spectral-index change of the spectrum of ultrarelativistic electrons injected into the jet. Based on this result, a detailed model simulation, reproducing the observed optical and X-ray spectral variability and broadband SED of BL Lacertae simultaneously, is presented. Our leptonic modeling results are compared to fits using the hadronic synchrotron-proton blazar (SPB) model. That model can reproduce the observed SEDs of BL Lacertae in a scenario with muon-synchrotron dominated high-energy emission. It requires a significantly higher magnetic field than the leptonic model (~ 40 G vs. ~ 2 G in the leptonic model) and a lower Doppler factor associated with the bulk motion of the emission region (D ~ 8 vs. D ~ 18 in the leptonic model). The hadronic model predicts a significantly larger > 100 GeV flux than the leptonic models, well within the anticipated capabilities of VERITAS and MAGIC.

Our recent progress on time-dependent modeling of the multiwavelength spectra and variability of blazars with leptonic and hadronic jet models is reviewed. Special emphasis is placed on X-ray spectral variability of low-frequency peaked (LBLs) and intermediate BL Lac objects (IBLs). As an example, recent observational and modeling results of an extensive multiwavelength campaign of BL Lacertae in 2000 are presented. It is demonstrated how combined spectral and variability modeling of LBLs and IBLs can significantly constrain emission models and potential variability scenarios. In the case of BL Lacertae, the variability appears to be driven primarily by fluctuations in the spectral index of the non-thermal, ultrarelativistic electron population in the jet. Such constraints allow us to refine predictions of the intrinsic GeV gamma-ray emission and the dominant electron cooling mechanism in these objects.

Absorption and reprocessing of Gamma-ray burst radiation in the environment of cosmological GRBs can be used as a powerful probe of the elusive nature of their progenitors. Although it is widely accepted that long-duration GRBs are associated with the deaths of massive stars, at least two fundamentally different scenarios concerning the final collapse are currently being considered. Delayed reddened excesses in the optical afterglows of several GRBs indicate that a supernova, possibly of Type Ic, takes place within a few days of a GRB. This supports the collapsar model, where the core of a massive star collapses promptly to a black hole. Variable X-ray features observed in the prompt and afterglow spectra of several GRBs, suggest that a highly metal enriched and dense shell of material surrounds the sources of GRBs. In some cases, evidence for expansion of these shells with velocities of a substantial fraction of the speed of light has been claimed. These observations have been interpreted as support for the supranova model, where a massive star collapses first to a supramassive neutron star, which later collapses to a black hole following loss of rotational support. In this review paper, I will present a brief overview of the current status of the observational evidence for X-ray spectral features in GRBs, and discuss their implications for both the collapsar and the supranova model.

BL Lacertae was the target of an extensive multiwavelength monitoring campaign in the second half of 2000. Simultaneous or quasi-simultaneous observations were taken at radio (UMRAO and Metsaehovi) and optical(WEBT collaboration) frequencies, in X-rays (BeppoSAX and RXTE), and at VHE gamma-rays (HEGRA). The WEBT optical campaign achieved an unprecedented time coverage, virtually continuous over several 10 - 20 hour segments. It revealed intraday variability on time scales of ~ 1.5 hours and evidence for spectral hardening associated with increasing optical flux. During the campaign, BL Lacertae underwent a major transition from a rather quiescent state prior to September 2000, to a flaring state for the rest of the year. This was also evident in the X-ray activity of the source. BeppoSAX observations on July 26/27 revealed a rather low X-ray flux and a hard spectrum, while a BeppoSAX pointing on Oct. 31 - Nov. 2, 2000, indicated significant variability on time scales of < a few hours, and provided evidence for the synchrotron spectrum extending out to ~ 10 keV during that time. During the July 26/27 observation, there is a tantalizing, though not statistically significant, indication of a time delay of ~ 4 - 5 hr between the BeppoSAX and the R-band light curve. Also, a low-significance detection of a time delay of 15 d between the 14.5 GHz and the 22 GHz radio light curves is reported. Several independent methods to estimate the co-moving magnetic field in the source are presented, suggesting a value of ~ 2 e_B^{2/7} G, where e_B is the magnetic-field equipartition factor w.r.t. the electron energy density in the jet.

We present two BeppoSAX observations of BL Lac (2200+420) as part of a multiwavelength campaign performed in 2000. The source was in different states of activity: in June, the X-ray spectrum was faint and hard (alpha sim 0.8), with positive residuals towards low energies. In October we detected the highest [2-10] keV flux ever measured for the source. During this observation, the spectrum was soft (alpha sim 1.56) up to 10 keV, while above this energy a hard component was dominating. The BeppoSAX data are confirmed by simultaneous RXTE short observations. During the first observation the soft X-ray flux was variable on timescales of a few hours, while the hard X-ray flux was almost constant. During the second observation, instead, the soft spectrum displayed an erratic behaviour with large variations (up to factors 3-4) on timescales smaller than 2 hrs. The analysis of the multiwavelength SED of October evidenced an intriguing feature: the optical and X-ray sections of the SED are misaligned, while in the prevailing standard picture, they are both thought to be produced via synchrotron emission. We suggested four scenarios to account for this discrepancy: a higher than galactic dust-to-gas ratio towards the source, the first detection of bulk Compton emission in the X-ray band, the presence of two synchrotron emitting regions located at different distances from the nucleus, the detection of a Klein-Nishina effect on the synchrotron spectrum. We evidenced the favorable and critical points of each scenario, but, at present, we cannot discriminate between them.

Absorption and Reprocessing of Gamma-ray burst radiation in the environment of cosmological GRBs can be used as a powerful probe of the elusive nature of their progenitors. In particular, transient X-ray emission line and absorption features in the prompt and early afterglows of GRBs are sensitive to details of the location and density structure of the reprocessing and/or absorbing material. To date, there have been only rather few detections of such features, and the significance is marginal in most individual cases. However, transient X-ray emission lines in GRB afterglows have now been found by four different X-ray satellites, which may justify a more detailed theoretical investigation of their origin. In this paper, I will first present a brief review of the status of observations of transient X-ray emission line and absorption features. I will then discuss general physics constraints which those results impose on isotropy, homogeneity, and location of the reprocessing material with respect to the GRB source, and review the various currently discussed, specific models of GRBs and their environments in which the required conditions could arise.

We present two BeppoSAX observations of BL Lacertae as part of a multiwavelength radio-to-TeV campaign. During the first observation we observe a faint Compton spectrum, while during the second, we detect a synchrotron spectrum with the highest [2-10] keV flux ever measured; above 10 keV an inverse Compton component begin to dominate. The synchrotron flux is very variable with time scales of 1 hr. We describe four different SED shifting the synchrotron peak both in frequency and flux intensity and we sketch a scenario in which a blob moves along a jet and can be located in or outside the BLR. This implies different radiative mechanism at work, SSC or external Compton, producing different high energy spectra.

We are presenting a detailed parameter study of the time-dependent electron injection and kinematics and the self-consistent radiation transport in jets of intermediate and low-frequency peaked BL Lac objects. Using a time-dependent, combined synchrotron-self-Compton and external-Compton jet model, we study the influence of variations of several essential model parameters, such as the electron injection compactness, the relative contribution of synchrotron to external soft photons to the soft photon compactness, the electron-injection spectral index, and the details of the time profiles of the electron injection episodes giving rise to flaring activity. In the analysis of our results, we focus on the expected X-ray spectral variability signatures in a region of parameter space particularly well suited to reproduce the broadband spectral energy distributions of intermediate and low-frequency peaked BL Lac objects. We demonstrate that SSC- and external-Compton dominated models for the gamma-ray emission from blazars are producing significantly different signatures in the X-ray variability, in particular in the soft X-ray light curves and the spectral hysteresis at soft X-ray energies, which can be used as a powerful diagnostic to unveil the nature of the high-energy emission from BL Lac objects.

Spectral fitting of the radio through hard X-ray emission of BL Lac objects has previously been used to predict their level of high-energy (GeV - TeV) emission. In this paper, we point out that such spectral fitting can have very large uncertainties with respect to predictions of the VHE emission. This is demonstrated with the example of W Comae. We show that the best currently available contemporaneous optical - X-ray spectrum of W Comae still allows for a large range of possible parameters, resulting in drastically different > 40 GeV fluxes. We find that all acceptable leptonic-model fits to the optical - X-ray emission of W Comae predict a cut-off of the high-energy emission around ~ 100 GeV. We suggest that detailed measurements and analysis of the soft X-ray variability of W Comae may be used to break the degeneracy in the choice of possible fit parameters, and thus allow a more reliable prediction of the VHE emission from this object. Using the available soft X-ray variability measured by BeppoSAX, we predict a > 40 GeV flux from W Comae of ~ (0.4 - 1) * 10^{-10} photons/(cm^2 s) with no significant emission at E > 100 GeV for a leptonic jet model. We compare our results concerning leptonic jet models with detailed predictions of the hadronic Synchrotron-Proton Blazar model. This hadronic model predicts > 40 GeV fluxes very similar to those found for the leptonic models, but results in > 100 GeV emission which should be clearly detectable with future high-sensitivity instruments like VERITAS. Thus, we suggest this object as a promising target for VHE gamma-ray and co-ordinated broadband observations to distinguish between leptonic and hadronic jet models for blazars.

It is generally accepted that relativistic jet outflows power the nonthermal emission from active galactic nuclei (AGN). The composition of these jets -- leptonic versus hadronic -- is still under debate. We investigate the microphysical details of the conversion process of the kinetic energy in collimated relativistic pair outflows into radiation through interactions with the ambient interstellar medium. Viewed from the coordinate system comoving with the pair outflow, the interstellar protons and electrons represent a proton-electron beam propagating with relativistic speed in the pair plasma. We demonstrate that the beam excites both electrostatic and low-frequency magnetohydrodynamic Alfven-type waves via a two-stream instability in the pair background plasma, and we calculate the time evolution of the distribution functions of the beam particles and the generated plasma wave turbulence power spectra. For standard AGN jet outflow and environment parameters we show that the initial beam distributions of interstellar protons and electrons quickly relax to plateau-distributions in parallel momentum, transferring thereby one-half of the initial energy density of the beam particles to electric field fluctuations of the generated electrostatic turbulence. On considerably longer time scales, the plateaued interstellar electrons and protons will isotropise by their self-generated transverse turbulence and thus be picked-up in the outflow pair plasma. These longer time scales are also characteristic for the development of transverse hydromagnetic turbulence from the plateaued electrons and protons. This hydromagnetic turbulence upstream and downstream is crucial for diffusive shock acceleration to operate at external or internal shocks associated with pair outflows.

We have organized an extensive multiwavelength campaign on BL Lacertae in the second half of 2000. Simultaneous or quasi-simultaneous observations were taken at radio frequencies, in the optical - carried out by the WEBT collaboration -, in X-rays - using BeppoSAX and RXTE -, and at VHE gamma-rays with the CAT and HEGRA Cherenkov telescope facilities. In this paper, we are presenting first results from this campaign. The WEBT optical campaign achieved an unprecedented time coverage, virtually continuous over several 10 - 20 hour segments, and revealed intraday variability on time scales of ~ 1.5 hours. The X-ray observations of Nov. 1 - 2, 2000, revealed significant variability on similar time scales, and provided evidence for the synchrotron spectrum extending out to ~ 10 keV during that time. At higher energies, the onset of the hard power-law commonly observed in radio-loud quasars and radio-selected BL Lac objects is seen. From the energy-resolved X-ray variability patterns, we find evidence that electron cooling might not be dominated by one single radiation mechanism; rather, a combination of synchrotron and Compton cooling seems to be required.

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.

We find that energy losses due to synchrotron self-Compton (SSC) emission in blazar jets can produce distinctive signatures in the time-averaged synchrotron and SSC spectra of these objects. For a fairly broad range of particle injection distributions, SSC-loss dominated synchrotron emission exhibits a spectral dependence $F_\nu \sim \nu^{-3/2}$. The presence or absence of this dependence in the optical and ultraviolet spectra of flat spectrum radio quasars such as 3C~279 and in the soft X-ray spectra of high frequency BL Lac objects such as Mrk 501 gives a robust measure of the importance of SSC losses. Furthermore, for partially cooled particle distributions, spectral breaks of varying sizes can appear in the synchrotron and SSC spectra and will be related to the spectral indices of the emission below the break. These spectral signatures place constraints on the size scale and the non-thermal particle content of the emitting plasma as well as the observer orientation relative to the jet axis.

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.

It is generally accepted that gamma-ray bursts (GRBs) are initiated by a relativistic pair fireball, converting its internal energy into kinetic energy of a relativistically moving plasmoid and subsequently into radiation. Here, we investigate the early stages of this evolution, after the pair fireball has become optically thin to gamma-gamma pair production. We show that for a short period of time, ~ 0.1 - a few seconds after the initial explosion, the pair plasmoid evolution might be dominated by collisional processes prior to the formation of a collisionless shock. We simulate these processes during the early pair plasmoid evolution and calculate the expected radiative signatures. We show that the radiation from the collision-dominated pair plasmoid phase results in a short (~ a few ms) flash of thermal soft X-ray emission, followed by a transition phase of < 1 s during which the fireball turns Thomson thin, but its radiation remains dominated by thermal Comptonization, peaking at around E_pk ~ 100 MeV - a few GeV. While the very early thermal emission could be associated with the quasi-thermal radiation signatures found in the very early phases of several bright BATSE GRBs, the predicted subsequent flash of high-energy emission should be easily detectable with the GLAST satellite.

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.

The current status of leptonic jet models for blazars is reviewed. Differences between the quasar and BL-Lac subclasses of blazars may be understood in terms of the dominance of different radiation mechanisms in the gamma-ray regime. Spectral variability patterns of different blazar subclasses appear to be significantly different and require different intrinsic mechanisms causing gamma-ray flares. As examples, recent results of long-term multiwavelength monitoring of PKS 0528+134, 3C 279, and Mrk 501 are interpreted in the framework of leptonic jet models. Short-term variability patterns give important additional clues about the source geometry and the relevant radiation mechanisms in blazars. Challenges for future observational efforts are discussed.

Light curves of 3C 279 are presented in optical (R-band), X-rays (RXTE/PCA), and gamma rays (CGRO/EGRET) for 1999 Jan-Feb and 2000 Jan-Mar. During both of those epochs the gamma-ray levels were high, and all three observed bands demonstrated substantial variation, on time scales as short as one day. Correlation analyses provided no consistent pattern, although a rather significant optical/gamma-ray correlation was seen in 1999, with a gamma-ray lag of ~2.5 days, and there are other suggestions of correlations in the light curves. For comparison, correlation analysis is also presented for the gamma-ray and X-ray light curves during the large gamma ray flare in 1996 Feb and the two gamma-bright weeks leading up to it; the correlation at that time was strong, with a gamma-ray/X-ray offset of no more than 1 day.

We suggest a simple fitting formula to represent Comptonized X- and gamma-ray spectra from a hot ($kT_e \gtrsim 10$ keV), Thomson thick ($\tau_T \gtrsim 5$) hybrid thermal/nonthermal plasma in spherical geometry with homogeneous soft-photon injection throughout the Comptonizing region. We have used this formula to fit a large data base of Monte-Carlo generated photon spectra, and provide correlations between the physical parameters of the plasma and the fit parameters of our analytic fit function. These allow us to construct Thomson thick Comptonization spectra without performing computer intensive Monte Carlo simulations of the high-Thomson-depth hybrid-plasma Comptonization problem. Our formulae can easily be used in data analysis packages such as XSPEC, thus rendering rapid $\chi^2$ fitting of such spectra to real data feasible.

Of the blazars detected by EGRET in GeV gamma rays, 3C 279 is not only the best-observed by EGRET, but also one of the best-monitored at lower frequencies. We have assembled eleven spectra, from GHz radio through GeV gamma rays, from the time intervals of EGRET observations. Although some of the data have appeared in previous publications, most are new, including data taken during the high states in early 1999 and early 2000. All of the spectra show substantial gamma-ray contribution to the total luminosity of the object; in a high state, the gamma-ray luminosity dominates over that at all other frequencies by a factor of more than 10. There is no clear pattern of time correlation; different bands do not always rise and fall together, even in the optical, X-ray, and gamma-ray bands. The spectra are modeled using a leptonic jet, with combined synchrotron self-Compton + external Compton gamma-ray production. Spectral variability of 3C 279 is consistent with variations of the bulk Lorentz factor of the jet, accompanied by changes in the spectral shape of the electron distribution. Our modeling results are consistent with the UV spectrum of 3C 279 being dominated by accretion disk radiation during times of low gamma-ray intensity.

The most popular models for the complex phase and time lags in the rapid aperiodic variability of Galactic X-ray binaries are based Comptonization of soft seed photons in a hot corona, where small-scale flares are induced by flares of the soft seed photon input (presumably from a cold accretion disc). However, in their original version, these models have neglected the additional cooling of the coronal plasma due to the increased soft seed photon input, and assumed a static coronal temperature structure. In this paper, our Monte-Carlo/Fokker-Planck code for time-dependent radiation transfer and electron energetics is used to simulate the self-consistent coronal response to the various flaring scenarios that have been suggested to explain phase and time lags observed in some Galactic X-ray binaries. It is found that the predictions of models involving slab-coronal geometries are drastically different from those deduced under the assumption of a static corona. However, with the inclusion of coronal cooling they may even be more successful than in their original version in explaining some of the observed phase and time lag features. The predictions of the model of inward-drifting density perturbations in an ADAF-like, two-temperature flow also differ from the static-corona case previously investigated, but may be consistent with the alternating phase lags seen in GRS 1915+105 and XTE J1550-564. Models based on flares of a cool disc around a hot, inner two-temperature flow may be ruled out for most objects where significant Fourier-frequency-dependent phase and time lags have been observed.

We discuss the space-and-time-dependent Monte Carlo code we have developed to simulate the relativistic radiation output from compact astrophysical objects, coupled to a Fokker-Planck code to determine the self-consistent lepton populations. We have applied this code to model the emission from a magnetized neutron star accretion shell near the Alfven radius, reprocessing the radiation from the neutron sar surface. We explore the parameter space defined by the accretion rate, stellar surface field and the level of wave turbulence in the shell. Our results are relevant to the emission from atoll sources, soft-X-ray transient X-ray binaries containing weakly magnetized neutron stars, and to recently suggested models of accretion-powered emission from anomalous X-ray pulsars.

We investigate the potentially observable prompt or delayed X-ray spectral features from the currently popular gamma-ray burst (GRB) models. During the evolution of many GRB progenitors, a disk around the central GRB source is produced. Shock heating as the GRB ejecta collide with the disk may produce observable X-ray features. We first summarize predictions deduced from previous calculations which invoke photoionization and relativistic blast waves. We then calculate the quasi-thermal X-ray line features produced assuming the ejecta are nonrelativistic (which is more likely for the disk interactions of many GRB models). In the framework of the Hypernova/Collapsar model, delayed (a few days - several months after the GRB) bursts of line-dominated, thermal X-ray emission may be expected. The He-merger scenario predicts similar X-ray emission line bursts <~ a few days after the GRB. These X-ray signatures should be observable with Chandra and XMM-Newton out to at least z ~ 1. Weak emission line features <~ a few days after the GRB may also result from the supranova GRB scenario. In all three cases, significant X-ray absorption features, in particular during the prompt GRB phase, are expected. No significant X-ray spectral features might result from compact-object binary mergers.

The time and angle dependent line and continuum emission from a dense torus around a cosmological gamma-ray burst source is simulated, taking into account photoionization, collisional ionization, recombination, and electron heating and cooling due to various processes. The importance of the hydrodynamical interaction between the torus and the expanding blast wave is stressed. Due to the rapid deceleration of the blast wave as it interacts with the dense torus, the material in the torus will be illuminated by a drastically different photon spectrum than observable through a low-column-density line of sight, and will be heated by the hydrodynamical interaction between the blast wave and the torus. A model calculation to reproduce the Fe K-alpha line emission observed in the X-ray afterglow of GRB 970508 is presented. The results indicate that ~ 10^{-4} solar masses of iron must be concentrated in a region of less than 10^{-3} pc. The illumination of the torus material due to the hydrodynamic interaction of the blast wave with the torus is the dominant heating and ionization mechanism leading to the formation of the iron line. These results suggest that misaligned GRBs may be detectable as X-ray flashes with pronounced iron emission line features.

We present a theoretical model for the alternating phase lags associated with QPO fundamental and harmonic frequencies observed in some Galactic black-hole candidates. We assume that the accretion flow exhibits a transition from an outer cool, optically thick accretion disk to a hot, inner advection-dominated accretion flow (ADAF), and that the QPOs are related to small-scale oscillations of the accretion rate and the transition radius. We present an analytical estimate of the expected phase lags at the fundamental and first harmonic frequencies of the QPOs and perform detailed time-dependent Monte-Carlo simualtions of the radiation transport in the oscillating ADAF / cool disk system. We find that this model is well suited to reproduce alternating phase lags between the fundamental and the first harmonic. It also naturally explains the trend observed in GRS 1915+105 that, as the soft X-ray luminosity increases, the QPO frequency increases and the phase lag associated with the QPO fundamental frequency changes sign from positive to negative. The relation between the disk temperature and the QPO frequency observed in GRS 1915+105 is consistent with a secular instability modulating the disk evaporation at the transition radius.

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.

During the observation period 1997, the nearby Blazar Mkn 501 showed extremely strong emission and high variability. We examine multiwavelength aspects of this event using radio, optical, soft and hard X-ray and TeV data. We concentrate on the medium-timescale variability of the broadband spectra, averaged over weekly intervals. We confirm the previously found correlation between soft and hard X-ray emission and the emission at TeV energies, while the source shows only minor variability at radio and optical wavelengths. The non-linear correlation between hard X-ray and TeV fluxes is consistent with a simple analytic estimate based on an SSC model in which Klein-Nishina effects are important for the highest-energy electrons in the jet, and flux variations are caused by variations of the electron density and/or the spectral index of the electron injection spectrum. The time-averaged spectra are fitted with a Synchrotron Self-Compton (SSC) dominated leptonic jet model, using the full Klein-Nishina cross section and following the self-consistent evolution of relativistic particles along the jet, accounting for gamma-gamma absorption and pair production within the source as well as due to the intergalactic infrared background radiation. The contribution from external inverse-Compton scattering is tightly constrained by the low maximum EGRET flux and found to be negligible at TeV energies. We find that high levels of the X-ray and TeV fluxes can be explained by a hardening of the energy spectra of electrons injected at the base of the jet, in remarkable contrast to the trend found for gamma-ray flares of the flat-spectrum radio quasar PKS 0528+134.

The Compton Gamma-Ray Observatory (CGRO) experiments EGRET and COMPTEL observed the Virgo sky region continuously for 7 weeks between December 10, 1996 and January 28, 1997. The prominent quasar 3C~273 was found to be the brightest source in gamma-rays and was significantly detected by EGRET and COMPTEL. The EGRET experiment observed a time-variable flux at energies above 100 MeV, which reached in a 2-week flaring period (December 30, 1996 to January 14, 1997) its highest flux level observed during the CGRO-era. COMPTEL, however, does not observe obvious time variability at energies below ~30 MeV contemporaneous to EGRET. In particular, no flare was observed, indicating that this outburst is solely a high-energy (>100 MeV) phenomenon. The energy spectrum between 3 MeV and 10 GeV is well represented by a simple power-law model. Below 3 MeV a spectral turnover is indicated. Performing spectral analysis for different time periods, we found evidence for a spectral hardening during the flaring period, which is consistent with the flare occurring mainly at the higher energies and with its absence at COMPTEL energies of a few MeV. This may be interpreted as an indication that the emission in the EGRET energy range is dominated by a different radiation mechanism than the MeV emission. We argue that the most likely mechanism for the high-energy flare is inverse-Compton scattering of reprocessed accretion-disk radiation.

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.

The multiwavelength spectrum of BL Lacertae during its July 1997 outburst is analyzed in terms of different variations of the homogeneous leptonic jet model for the production of high-energy radiation from blazars. We find that a two-component gamma-ray spectrum, consisting of a synchrotron self-Compton and an external Compton component, is required in order to yield an acceptable fit to the broadband spectrum. Our analysis indicates that in BL Lac, unlike other BL Lac objects, the broad emission line region plays an important role for the high-energy emission. Several alternative blazar jet models are briefly discussed. In the appendix, we describe the formalism in which the process of Comptonization of reprocessed accretion disk photons is treated in the previously developed blazar jet simulation code which we use.

The current status of leptonic jet models for gamma-ray blazars is reviewed. Differences between the quasar and BL-Lac subclasses of blazars may be understood in terms of the dominance of different radiation mechanisms in the gamma-ray regime. Spectral variability patterns of different blazar subclasses appear to be significantly different and require different intrinsic mechanisms causing gamma-ray flares. As examples, recent results of long-term multiwavelength monitoring of PKS 0528+134 and Mrk 501 are interpreted in the framework of leptonic jet models. A simple quasi-analytic toy model for broadband spectral variability of blazars is presented.

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.

The multiwavelength spectra of PKS 0528+134 during six years of observations by EGRET have been analyzed using synchrotron self-Compton (SSC) and external radiation Compton (ERC) models. We find that a two-component model, in which the target photons are produced externally to the gamma-ray emitting region, but also including an SSC component, is required to suitably reproduce the spectral energy distributions of the source. Our analysis indicates that there is a trend in the observed properties of PKS 0528+134, as the source goes from a gamma-ray low state to a flaring state. We observe that during the higher gamma-ray states, the bulk Lorentz factor of the jet increases and the ERC component dominates the high-energy emission. Our model calculations indicate the trend that the energies of the electrons giving rise to the synchrotron peak decreases, and the power-ratio of the gamma-ray and low energy spectral components increases, as the source goes from a low to a high gamma-ray state.

The recent simultaneous detection of optical, X-ray and gamma-ray photons from GRB990123 during the burst provides the first broadband multi-wavelength characterization of the burst spectrum and evolution. Here we show that a direct correlation exists between the time-varying gamma-ray spectral shape and the prompt optical emission. This combined with the unique signatures of the time-resolved spectra of GRB990123 convincingly supports earlier predictions of the saturated Comptonization model. Contrary to other suggestions, we find that the entire continuum from optical to gamma-rays can be generated from a single source of leptons (electrons and pairs). The optical flux only appears to lag the gamma-ray flux due to the high initial Thomson depth of the plasma. Once the plasma has completely thinned out, the late time afterglow behavior of our model is the same as in standard models based on the Blandford-McKee (1976) solution.

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.

A prediction for the energy shift of the synchrotron spectrum of flat-spectrum radio quasars (FSRQs) during high-energy flares is presented. If the $\gamma$-ray emission of FSRQs is produced by Comptonization of external radiation, then the peak of the synchrotron spectrum is predicted to move to lower energies in the flare state. This is opposite to the well-known broadband spectral behavior of high-frequency peaked BL-Lac objects where the external radiation field is believed to be weak and synchrotron-self Compton scattering might be the dominant $\gamma$-ray radiation mechanism. The synchrotron peak shift, if observed in FSRQs, can thus be used as a diagnostic to determine the dominant radiation mechanism in these objects. I suggest a few FSRQs as promising candidates to test the prediction of the external-Comptonization 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.

The time-dependent Comptonized output of a cool soft X-ray source drifting inward through an inhomogeneous hot inner disk or corona is numerically simulated. We propose that this scenario can explain from first principles the observed trends in the hard time lags and power spectra of the rapid aperiodic variability of the X-ray emission of Galactic black-hole candidates.

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$.

We present a detailed study of the observable effects of photoionization and Comptonization of line and continuum radiation from a cold accretion disk with a thin, warm, photoionized transition layer in the framework of self-consistent accretion-disk corona models for Galactic black-hole candidates. We use an iterative method to treat the non-linear radiation feedback between the transition layer and the hot corona numerically using a Monte-Carlo Comptonization code in combination with the photoionization and line transfer code XSTAR. The subset of the parameter space allowed in self-consistent accretion-disk corona systems on energetic grounds is explored, checking for the emergence of emission lines and/or absorption edges, with special emphasis on the spectral range $E \lesssim 1$ keV where such features might become observable with the advent of the AXAF satellite. Comparing our model calculations to the broadband spectrum of GX~339-4, we find good agreement with the observed spectral features and discuss how the future detection of the predicted features at lower energies can be used to constrain parameters.

Models of relativistic jets filled with ultrarelativistic pair plasma are very successful in explaining the broadband radiation of $\gamma$-ray blazars. Assuming that the initial injection and cooling of ultrarelativistic pair plasma in an AGN jet has occurred, producing the observed high-energy $\gamma$-ray radiation, we investigate the further evolution of the pair plasma as it continues to move out from the central engine. The effects of thermalization and reacceleration, the emission of pair bremsstrahlung and annihilation radiation and the bulk Compton process, and the possible application to MeV blazars are discussed. A model calculation to the special case of PKS~0208-512 is presented.

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 effect of inverse-Compton scattering of flares of soft radiation in different geometries of a hot, Comptonizing region and a colder accretion disk around a solar-mass black hole. The photon-energy dependent light curves, their Fourier transforms, power spectra and Fourier-period dependent time lags of hard photons with respect to softer photons are discussed. On the basis of a comparison with existing data we find arguments against Comptonization of external soft radiation as well as Comptonization in a homogeneous medium as dominant mechanisms for the rapid aperiodic variability in Galactic black-hole candidates. Possible further observational tests for the influence of Comptonization on the rapid aperiodic variability of Galactic black-hole candidates are suggested.

We present model calculations to the observed spectral variability in the gamma-ray spectrum of the ultraluminous blazar-type quasar PKS 0528+134. We argue that the observed appearance of a spectral break between the COMPTEL and the EGRET energy range can plausibly be explained by a variation of the Doppler beaming factor in the framework of a relativistic jet model for AGNs, combining the ERC and the SSC model. This may be an instructive example for the different dependence of these two spectral components on the Doppler factor, as predicted in earlier work.

We present a model for the inverted radio spectra of active active galactic nuclei as well as the central regions of normal galaxies. The model is based on the unified scenario for active galaxies, stating that the central engines of active galaxies consists of a supermassive black hole surrounded by an accretion disk and a radio jet. The nuclei of normal (i. e. less active) galaxies are supposed to be scaled-down versions of the same phenomenon. We show that the radio emission of a jet component, becoming optically thin to the radio emission of a monoenergetic pair plasma at decreasing frequencies as it moves outward and expands, is well suited to explain the observational results. We present a model calculation for the special case of the nucleus of M~81.

We present model calculations reproducing broadband spectra of $\gamma$-ray blazars by a relativistic leptonic jet, combining the EIC and the SSC model. To this end, the evolution of the particle distribution functions inside a relativistic pair jet and of the resulting photon spectra is investigated. Inverse-Compton scattering of both external (EIC) as well as synchrotron photons (SSC) is treated using the full Klein-Nishina cross-section and the full angle-dependence of the external photon source. We present model fits to the broadband spectra of Mrk~421 and 3C279 and the X-ray and $\gamma$-ray spectrum of PKS~1622-297. We find that the most plausible way to explain both the quiescent and the flaring states of these objects consists of a model where EIC and SSC dominate the observed spectrum in different frequency bands. For both Mrk~421 and 3C279 the flaring states can be reproduced by a harder spectrum of the injected pairs.

We present the first completely analytical computation of the full differential \gamma-\gamma pair production rate from compact radiation fields, exact to 2nd order QED, and use this result to investigate the validity of previously known approximations.

The evolution of the particle distribution functions inside a relativistic jet containing an e^-e^+ pair plasma and of the resulting $\gamma$-ray and X-ray spectra is investigated. The first phase of this evolution is governed by heavy radiative energy losses. For this phase, approximative expressions for the energy-loss rates due to inverse-Compton scattering, using the full Klein-Nishina cross section, are given as one-dimensional integrals for both cooling by inverse-Compton scattering of synchrotron photons (SSC) and of accretion disk photons (EIC). We calculate instantaneous and time-integrated $\gamma$-ray spectra emitted by such a jet for various sets of parameters, deducing some general implications on the observable broadband radiation. Finally, we present model fits to the broadband spectrum of Mrk 421. We find that the most plausible way to explain both the quiescent and the flaring state of Mrk 421 consists of a model where EIC and SSC dominate the observed spectrum in different frequency bands. In our model the flaring state is primarily related to an increase of the maximum Lorentz factor of the injected pairs.

The contribution of the pair annihilation process in relativistic electron-positron jets to the gamma-ray emission of blazars is calculated. Under the same assumptions as for the calculation of the yield of inverse Compton scattered accretion disk radiation (Dermer and Schlickeiser 1993) we calculate the emerging pair annihilation radiation taking into account all spectral broadening effects due to the energy spectra of the annihilating particles and the bulk motion of the jet. It is shown that the time-integrated pair annihilation spectrum appears almost like the well-known gamma-ray spectrum from decaying $\pi^o$-mesons at rest, yielding a broad bumpy feature located between 50 and 100 MeV. We also demonstrate that for pair densities $> 10^9$ cm$^{-3}$ in the jet the annihilation radiation will dominate the inverse Compton radiation, and indeed may explain reported spectral bumps at MeV energies. The refined treatment of the inverse Compton radiation leads to spectral breaks of the inverse Compton emission in the MeV energy range with a change in spectral index $\Delta \alpha $ larger than 0.5 as detected in PKS 0528+134 and 3C273.