We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

The majority of the activity around nearby (z ~ 0) supermassive black holes is found in low-luminosity active galactic nuclei (LLAGN), the most of them being classified as low ionization nuclear emission regions. Although these sources are well studied from radio up to X-rays, they are poorly understood in gamma-rays. In this work we take advantage of the all sky-surveying capabilities of the Large Area Telescope on board Fermi Gamma ray Space Telescope to study the whole Palomar sample of LLAGN in gamma-rays. Precisely, the four radio-brightest LLAGN in the sample are identified as significant gamma-ray emitters, all of which are recognized as powerful Fanaroff-Riley I galaxies. These results suggest that the presence of powerful radio jets is of substantial importance for observing a significant gamma-ray counterpart even if these jets are misaligned with respect to the line of sight. We also find that most of the X-ray-brightest LLAGN do not have a significant gamma-ray and strong radio emission, suggesting that the X-rays come mainly from the accretion flow in these cases. A detailed analysis of the spectral energy distributions (SEDs) of NGC 315 and NGC 4261, both detected in gamma-rays, is provided where we make a detailed comparison between the predicted hadronic gamma-ray emission from a radiatively inefficient accretion flow (RIAF) and the gamma-ray emission from a leptonic jet-dominated synchrotron self-Compton (SSC) model. Both SEDs are better described by the SSC model while the RIAF fails to explain the gamma-ray observations.

Multiwavelength monitoring of the blazar 3C 279 observed a very bright, 12-hour, orphan gamma-ray flare on 20 Dec 2013 with a uniquely hard Fermi-LAT spectrum and high Compton dominance. We work with a one-zone, leptonic model with both first- and second-order Fermi acceleration, which now reproduces the unique flaring behavior. We present a simplified analytic electron energy distribution to provide intuition about how particle acceleration shapes multi-wavelength blazar jet emission spectra. The contributions of individual processes in relativistic jets is fundamental to understanding the particle energy budget in the formation and propagation of astrophysical jets. We show that first- and second-order Fermi acceleration are sufficient to explain the flare, and that magnetic reconnection is not needed. Our analysis suggests that the flare is initiated by an increase in the particle energies due to shock acceleration, which also increases the stochastic acceleration. The higher energy particle preferentially occupy the outer jet, along the sheath, which decreases the apparent magnetic field and synchrotron radiation, while increasing electron exposure to the broad line region photon fields, driving up the external Compton emission.

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.

The Hubble constant $H_{0}$ and matter density $\Omega_{m}$ of the Universe are measured using the latest $\gamma$-ray attenuation results from Fermi-LAT and Cherenkov telescopes. This methodology is based upon the fact that the extragalactic background light supplies opacity for very high energy photons via photon-photon interaction. The amount of $\gamma$-ray attenuation along the line of sight depends on the expansion rate and matter content of the Universe. This novel strategy results in a value of $H_{0}=67.4_{-6.2}^{+6.0}$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.14_{-0.07}^{+0.06}$. These estimates are independent and complementary to those based on the distance ladder, cosmic microwave background (CMB), clustering with weak lensing, and strong lensing data. We also produce a joint likelihood analysis of our results from $\gamma$ rays and these from more mature methodologies, excluding the CMB, yielding a combined value of $H_{0}=66.6\pm 1.6$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.29\pm 0.02$.

The broadband spectrum from the 2013 December 20 $\gamma$-ray flare from 3C~279 is analyzed with our previously-developed one-zone blazar jet model. We are able to reproduce two SEDs, a quiescent and flaring state, the latter of which had an unusual SED, with hard $\gamma$-ray spectrum, high Compton dominance, and short duration. Our model suggests that there is insufficient energy for a comparable X-ray flare to have occurred simultaneously, which is an important constraint given the lack of X-ray data. We show that first- and second-order Fermi acceleration are sufficient to explain the flare, and that magnetic reconnection is not needed. The model includes particle acceleration, escape, and adiabatic and radiative energy losses, including the full Compton cross-section, and emission from the synchrotron, synchrotron self-Compton, and external Compton processes. We provide a simple analytic approximation to the electron distribution solution to the transport equation that may be useful for simplified modeling in the future.

In this paper, we characterize the first $\gamma$-ray flaring episode of the FSRQ PKS 0346-27 (z=0.991), as revealed by Fermi-LAT monitoring data, and the concurrent multi-wavelength variability observed from radio through X-rays. The quasi-simultaneous multi-wavelength coverage allowed us to construct time-resolved spectral energy distributions (SEDs). PKS 0346-27 entered an elevated $\gamma$-ray activity state starting from the beginning of 2018. The high-state continued throughout the year, displaying the highest fluxes in May 2018. We find evidence of short-time scale variability down to $\sim$1.5 hours, which constrains the $\gamma$-ray emission region to be compact. The extended flaring period was characterized by a persistently harder spectrum with respect to the quiescent state, indicating changes in the broadband spectral properties of the source. This was confirmed by the multi-wavelength observations, which show a shift in the position of the two SED peaks by $\sim$2 orders of magnitude in energy and peak flux value. As a result, during the high state the non-thermal jet emission completely outshines the thermal contribution from the dust torus and accretion disk. The broadband SED of PKS 0346-27 transitions from a typical Low-Synchrotron-Peaked (LSP) to the Intermediate-Synchrotron-Peaked (ISP) class, a behavior previously observed in other flaring $\gamma$-ray sources. Our one-zone leptonic emission model of the high-state SEDs constrains the $\gamma$-ray emission region to have a lower magnetic field, larger radius, and higher maximum electron Lorentz factors with respect to the quiescent SED. Finally, we note that the bright and hard $\gamma$-ray spectrum observed during the peak of flaring activity in May 2018 implies that PKS 0346-27 could be a promising target for future ground-based Cherenkov observatories such as the CTA.

This White Paper explores advances in the study of Active Galaxies which will be enabled by new observing capabilities at MeV energies (hard X-rays to gamma-rays; 0.1-1000 MeV), with a focus on multi-wavelength synergies. This spectral window, covering four decades in energy, is one of the last frontiers for which we lack sensitive observations. Only the COMPTEL mission, which flew in the 1990s, has significantly probed this energy range, detecting a handful of AGN. In comparison, the currently active Fermi Gamma-ray Space Telescope, observing at the adjacent range of 0.1-100 GeV, is 100-1000 times more sensitive. This White Paper describes advances to be made in the study of sources as diverse as tidal disruption events, jetted AGN of all classes (blazars, compact steep-spectrum sources, radio galaxies and relics) as well as radio-quiet AGN, most of which would be detected for the first time in this energy regime. New and existing technologies will enable MeV observations at least 50-100 times more sensitive than COMPTEL, revealing new source populations and addressing several open questions, including the nature of the corona emission in non-jetted AGN, the precise level of the optical extragalactic background light, the accretion mode in low-luminosity AGN, and the structure and particle content of extragalactic jets.

Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?

The Extragalactic Background Light (EBL) can be probed via the absorption imprint it leaves in the spectra of gamma-ray sources ($\gamma\gamma \rightarrow e^-e^+$). We recently developed a dedicated technique to reconstruct the EBL, and its evolution with redshift, from $\gamma$ ray optical depth data using a large sample of blazars detected by the $Fermi$ Large Area Telescope. Here, we extend this dataset to the TeV regime using ground-based Cherenkov observations of 38 blazars and report the first homogeneous measurement of the EBL spectral intensity covering the ultraviolet to infrared wavelengths ($\sim$0.1-100$\mathrm{\mu m}$). A minimal EBL throughout the wavelength range with respect to integrated galaxy light is found, allowing little additional unresolved emission from faint or truly diffuse populations setting an upper limit of $\lesssim 4~{\rm nW\cdot m^{-2}sr^{-1}}$ at 1.4 ${\rm \mu m}$. In particular, the cosmic optical background (COB) at $z=0$ is found to be $27.8_{-2.0}^{+2.1}~{\rm nW\cdot m^{-2}sr^{-1}}$. This work lays the foundation for accurate gamma-ray measurements of the EBL across its whole spectral range using a combination of GeV and TeV data.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is active galactic nuclei and gamma-ray bursts.

In 2007 September, the Crab Nebula exhibited a bright gamma-ray flare in the GeV energy range that was detected by AGILE. The observed emission at >160 MeV indicates that the radiating electrons had energies above the classical synchrotron radiation-reaction limit, thus presenting a serious challenge to classical models for electron acceleration in astrophysical environments. In this paper, we apply our recently developed time-dependent self-similar analytical model describing electrostatic acceleration in the explosive reconnection region around the pulsar wind termination shock to the 2007 September flare. This event was unique in that it displayed both long-duration "wave" and short-duration "sub-flare" features. The unusual temporal variation makes this flare an especially interesting test for our model. We demonstrate that our model can reproduce the time-dependent gamma-ray spectrum for this event, as well as the associated gamma-ray light curve, obtained by integrating the spectrum for photon energies >100 MeV. This establishes that our time-dependent electrostatic acceleration model can explain both wave and sub-flare transients, which lends further support to the theoretical framework we have developed. We also further examine the validity of the self-similar electric and magnetic field evolution implied by our model. We conclude that strong electrostatic acceleration driven by shock-induced magnetic reconnection is able to power the Crab Nebula gamma-ray flares by energizing the electrons on sub-Larmor timescales.

The light emitted by all galaxies over the history of the Universe produces the extragalactic background light (EBL) at ultraviolet, optical, and infrared wavelengths. The EBL is a source of opacity for $\gamma$ rays via photon-photon interactions, leaving an imprint in the spectra of distant $\gamma$-ray sources. We measure this attenuation using {739} active galaxies and one gamma-ray burst detected by the {\it Fermi} Large Area Telescope. This allows us to reconstruct the evolution of the EBL and determine the star-formation history of the Universe over 90\% of cosmic time. Our star-formation history is consistent with independent measurements from galaxy surveys, peaking at redshift $z\sim2$. Upper limits of the EBL at the epoch of re-ionization suggest a turnover in the abundance of faint galaxies at $z\sim 6$.

I show that by assuming a standard Blandford-Konigl jet, it is possible to determine the bulk Lorentz factor and angle to the line of sight of self-similar parsec-scale blazar jets by using five measured quantities: redshift, core radio flux, extended radio flux, the magnitude of the core shift between two frequencies, and apparent jet opening angle. From the bulk Lorentz factor and angle computed with this method, one can compute other jet properties such as the Doppler factor, magnetic field strength, and intrinsic jet opening angle. I use data taken from the literature and marginalize over nuisance parameters associated with the electron distribution and equipartition to compute these quantities, although the errors are large. Results are generally consistent with constraints from other methods. Primary sources of uncertainty are the errors on the core shift measurements and the uncertainty in the electron spectral index.

The very-high-energy (VHE, $\gtrsim 100$ GeV) $\gamma$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. We investigate the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and put it in context with possible emission scenarios. We collect photometric data in radio, optical, X-ray, and $\gamma$ ray. We study both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy $\gamma$ ray\ (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100$^\circ$. (...) The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars, i.e. inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to flat spectrum radio quasars (FSRQs) and low-frequency peaked BL Lac objects (LBL).

The $\gamma$-ray flares from the Crab nebula observed by {\it AGILE} and {\it Fermi}-LAT between 2007-2013 reached GeV photon energies and lasted several days. The strongest emission, observed during the 2011 April "super-flare," exceeded the quiescent level by more than an order of magnitude. These observations challenge the standard models for particle acceleration in pulsar wind nebulae, because the radiating electrons have energies exceeding the classical radiation-reaction limit for synchrotron. Particle-in-cell simulations have suggested that the classical synchrotron limit can be exceeded if the electrons also experience electrostatic acceleration due to shock-driven magnetic reconnection. In this paper, we revisit the problem using an analytic approach based on solving a fully time-dependent electron transport equation describing the electrostatic acceleration, synchrotron losses, and escape experienced by electrons in a magnetically confined plasma "blob" as it encounters and passes through the pulsar-wind termination shock. We show that our model can reproduce the $\gamma$-ray spectra observed during the rising and decaying phases of each of the two sub-flare components of the 2011 April super-flare. We integrate the spectrum for photon energies $\ge 100\,$MeV to obtain the light curve for the event, which agrees with the observations. We find that strong electrostatic acceleration occurs on both sides of the termination shock, driven by magnetic reconnection. We also find that the dominant mode of particle escape changes from diffusive escape to advective escape as the blob passes through the shock.

We introduce a new theoretical model to describe the emitting region in a blazar jet. We assume a one-zone leptonic picture, and construct the particle transport equation for a plasma blob experiencing low-energy, monoenergetic particle injection, energy dependent particle escape, shock acceleration, adiabatic expansion, stochastic acceleration, synchrotron radiation, and external Compton radiation from the dust torus and broad line region. We demonstrate that a one-zone leptonic model is able to explain the IR though {\gamma}-ray spectrum for 3C 279 in 2008-2009. We determine that the broad-line region seed photons cannot be adequately described by a single average distribution, but rather we find that a stratified broad line region provides an improvement in the estimation of the distance of the emitting region from the black hole. We calculate that the jet is not always in equipartition between the particles and magnetic field, and find that stochastic acceleration provides more energy to the particles than does shock acceleration, where the latter is also overshadowed by adiabatic losses. We further introduce a novel technique to implement numerical boundary conditions and determine the global normalization for the electron distribution, based on analysis of stiff ordinary differential equations. Our astrophysical results are compared with those obtained by previous authors.

The extragalactic background light (EBL), from ultraviolet to infrared wavelengths, is predominantly due to emission from stars, accreting black holes and reprocessed light due to Galactic dust. The EBL can be studied through the imprint it leaves, via $\gamma$-$\gamma$ absorption of high-energy photons, in the spectra of distant $\gamma$-ray sources. The EBL has been probed through the search for the attenuation it produces in the spectra of BL Lacertae (BL Lac) objects and individual $\gamma$-ray bursts (GRBs). GRBs have significant advantages over blazars for the study of the EBL especially at high redshifts. Here we analyze a combined sample of twenty-two GRBs, detected by the Fermi Large Area Telescope between 65 MeV and 500 GeV. We report a marginal detection (at the ~2.8 $\sigma$ level) of the EBL attenuation in the stacked spectra of the source sample. This measurement represents a first constraint of the EBL at an effective redshift of ~1.8. We combine our results with prior EBL constraints and conclude that Fermi-LAT is instrumental to constrain the UV component of the EBL. We discuss the implications on existing empirical models of EBL evolution.

Ultraluminous X-ray sources (ULXs) exceed the Eddington luminosity for a $\approx 10M_\odot$ black hole. The recent detection of black hole mergers by the gravitational wave detector ALIGO indicates that black holes with masses $> 10 M_\odot$ do indeed exist. Motivated by this, we explore a scenario where ULXs consist of black holes formed by the collapse of high-mass, low-metallicity stars, and that these ULXs become binary black holes (BBHs) that eventually merge. We use empirical relations between the number of ULXs and the star formation rate and host galaxy metallicity to estimate the ULX formation rate and the BBH merger rate at all redshifts. This assumes the ULX rate is directly proportional to the star formation rate for a given metallicity, and that the black hole accretion rate is distributed as a log-normal distribution. We include an enhancement in the ULX formation rate at earlier epochs due to lower mean metallicities. With simplified assumptions, our model is able to reproduce both the rate and mass distribution of BBH mergers in the nearby universe inferred from the detection of GW 150914, LVT 151012, GW 151226, and GW 170104 by ALIGO if the peak accretion rate of ULXs is a factor $\approx$1 --- 300 greater than the Eddington rate. Our predictions of the BBH merger rate, mass distribution, and redshift evolution can be tested by ALIGO in the near future, which in turn can be used to explore connections between the ULX formation and BBH merger rates over cosmic time.

Quasi-simultaneous observations of the Flat Spectrum Radio Quasar PKS 2326-502 were carried out in the gamma-ray, X-ray, UV, optical, near-infrared, and radio bands. Thanks to these observations we are able to characterize the spectral energy distribution of the source during two flaring and one quiescent gamma-ray states. These data were used to constrain one-zone leptonic models of the spectral energy distributions of each flare and investigate the physical conditions giving rise to them. While modeling one flare only required changes to the electron spectrum, the other flare needed changes in both the electron spectrum and the size of the emitting region with respect to the quiescent state. These results are consistent with an emerging pattern of two broad classes of flaring states seen in blazars.

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

After a long low-activity period, a gamma-ray flare from the narrow-line Seyfert 1 PKS 1502+036 (z=0.4089) was detected by the Large Area Telescope (LAT) on board Fermi in 2015. On 2015 December 20 the source reached a daily peak flux, in the 0.1-300 GeV band, of (93 $\pm$ 19) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, attaining a flux of (237 $\pm$ 71) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 3-hr time-scales, which corresponds to an isotropic luminosity of (7.3 $\pm$ 2.1) $\times$10$^{47}$ erg/s. The gamma-ray flare was not accompanied by significant spectral changes. We report on multi-wavelength radio-to-gamma-ray observations of PKS 1502+036 during 2008 August-2016 March by Fermi-LAT, Swift, XMM-Newton, Catalina Real-Time Transient Survey, and the Owens Valley Radio Observatory (OVRO). An increase in activity was observed on 2015 December 22 by Swift in optical, UV, and X-rays. The OVRO 15 GHz light curve reached the highest flux density observed from this source on 2016 January 12, indicating a delay of about three weeks between the gamma-ray and 15 GHz emission peaks. This suggests that the gamma-ray emitting region is located beyond the broad line region. We compared the spectral energy distribution (SED) of an average activity state with that of the flaring state. The two SED, with the high-energy bump modelled as an external Compton component with seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. The fit of the disc emission during the average state constrains the black hole mass to values lower than 10$^8$ solar masses. The SED, high-energy emission mechanisms, and gamma-ray properties of the source resemble those of a flat spectrum radio quasar.

The discovery by the Large Area Telescope on board Fermi of variable gamma-ray emission from radio-loud narrow-line Seyfert 1 (NLSy1) galaxies revealed the presence of a possible third class of Active Galactic Nuclei (AGN) with relativistic jets in addition to blazars and radio galaxies. Considering that NLSy1 are usually hosted in spiral galaxies, this finding poses intriguing questions about the nature of these objects and the formation of relativistic jets. We report on a systematic investigation of the gamma-ray properties of a sample of radio-loud NLSy1, including the detection of new objects, using 7 years of Fermi-LAT data with the new Pass 8 event-level analysis. In addition we discuss the radio-to-very-high-energy properties of the gamma-ray emitting NLSy1, their host galaxy, and black hole mass in the context of the blazar scenario and the unification of relativistic jets at different scales.

I study the location of the $\gamma$-ray emission in blazar jets by creating a Compton-scattering approximation valid for all anisotropic radiation fields in the Thomson through Klein-Nishina regimes, which is highly accurate and can speed up numerical calculations by up to a factor $\sim10$. I apply this approximation to synchrotron self-Compton, and external Compton-scattering of photons from the accretion disk, broad-line region (BLR), and dust torus. I use a stratified BLR model and include detailed Compton-scattering calculations of a spherical and flattened BLR. I create two dust torus models, one where the torus is an annulus, and one where it is an extended disk. I present detailed calculations of the photoabsorption optical depth using my detailed BLR and dust torus models, including the full angle dependence. I apply these calculations to the emission from a relativistically moving blob traveling through these radiation fields. The ratio of $\gamma$-ray to optical flux produces a predictable pattern that could help locate the $\gamma$-ray emission region. I show that the bright flare from 3C 454.3 in 2010 November detected by the Fermi Large Area Telescope is unlikely to originate from a single blob inside the BLR since it moves outside the BLR in a time shorter than the flare duration, although emission by multiple blobs inside the BLR is possible; and $\gamma$-rays are unlikely to originate from outside the BLR from scattering of photons from an extended dust torus, since then the cooling timescale would be too long to explain the observed short variability.

We report on Fermi Large Area Telescope (LAT) and multi-wavelength results on the recently-discovered very-high-energy (VHE, $E>$ 100 GeV) blazar S4 0954+65 ($z=0.368$) during an exceptionally bright optical flare in 2015 February. During the time period (2015 February, 13/14, or MJD 57067) when the MAGIC telescope detected VHE $\gamma$-ray emission from the source, the Fermi-LAT data indicated a significant spectral hardening at GeV energies, with a power-law photon index of $1.8 \pm 0.1$---compared with the 3FGL value (averaged over four years of observation) of $2.34 \pm 0.04$. In contrast, Swift/XRT data showed a softening of the X-ray spectrum, with a photon index of $1.72 \pm 0.08$ (compared with $1.38 \pm 0.03$ averaged during the flare from MJD 57066 to 57077), possibly indicating a modest contribution of synchrotron photons by the highest-energy electrons superposed on the inverse Compton component. Fitting of the quasi-simultaneous ($<1$ day) broadband spectrum with a one-zone synchrotron plus inverse-Compton model revealed that GeV/TeV emission could be produced by inverse-Compton scattering of external photons from the dust torus. We emphasize that a flaring blazar showing high flux of $\gtrsim 1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ ($E>$ 100 MeV) and a hard spectral index of $\Gamma_{\rm GeV} < 2.0$ detected by Fermi-LAT on daily time scales is a promising target for TeV follow-up by ground-based Cherenkov telescopes to discover high-redshift blazars, investigate their temporal variability and spectral features in the VHE band, and also constrain the intensity of the extragalactic background light.

Electromagnetic radiation from blazar jets often displays strong variability, extending from radio to $\gamma$-ray frequencies. In a few cases, this variability has been characterized using Fourier time lags, such as those detected in the X-rays from Mrk~421 using BeppoSAX. The lack of a theoretical framework to interpret the data has motivated us to develop a new model for the formation of the X-ray spectrum and the time lags in blazar jets based on a transport equation including terms describing stochastic Fermi acceleration, synchrotron losses, shock acceleration, adiabatic expansion, and spatial diffusion. We derive the exact solution for the Fourier transform of the electron distribution, and use it to compute the Fourier transform of the synchrotron radiation spectrum and the associated X-ray time lags. The same theoretical framework is also used to compute the peak flare X-ray spectrum, assuming that a steady-state electron distribution is achieved during the peak of the flare. The model parameters are constrained by comparing the theoretical predictions with the observational data for Mrk~421. The resulting integrated model yields, for the first time, a complete first-principles physical explanation for both the formation of the observed time lags and the shape of the peak flare X-ray spectrum. It also yields direct estimates of the strength of the shock and the stochastic MHD wave acceleration components in the Mrk~421 jet.

Blazars are active galactic nuclei with relativistic jets pointed at the Earth, making them extremely bright at essentially all wavelengths, from radio to gamma rays. I review the modeling of this broadband spectral energy distributions of these objects, and what we have learned, with a focus on gamma rays.

Distant BL Lacertae objects emit $\gamma$ rays which interact with the extragalactic background light (EBL), creating electron-positron pairs, and reducing the flux measured by ground-based imaging atmospheric Cherenkov telescopes (IACTs) at very-high energies (VHE). These pairs can Compton-scatter the cosmic microwave background, creating a $\gamma$-ray signature at slightly lower energies observable by the \fermi\ Large Area Telescope (LAT). This signal is strongly dependent on the intergalactic magnetic field (IGMF) strength ($B$) and its coherence length ($L_B$). We use IACT spectra taken from the literature for 5 VHE-detected BL Lac objects, and combine it with LAT spectra for these sources to constrain these IGMF parameters. Low $B$ values can be ruled out by the constraint that the cascade flux cannot exceed that observed by the LAT. High values of $B$ can be ruled out from the constraint that the EBL-deabsorbed IACT spectrum cannot be greater than the LAT spectrum extrapolated into the VHE band, unless the cascade spectrum contributes a sizable fraction of the LAT flux. We rule out low $B$ values ($B< 10^{-19}$ G for $L_B\ge1$\ Mpc) at $>5\sigma$ in all trials with different EBL models and data selection, except when using $>1$ GeV spectra and the lowest EBL models. We were not able to constrain high values of $B$.

Following the discovery of high-energy (HE; $E>10\,{\rm MeV}$) and very-high-energy (VHE; $E>100\,{\rm GeV}$) $\gamma$-ray emission from the low-frequency-peaked BL~Lac (LBL) object AP Librae, its electromagnetic spectrum is studied over 60 octaves in energy. Contemporaneous data in radio, optical and UV together with the $\gamma$-ray data are used to construct the most precise spectral energy distribution of this source. The data have been found to be modeled with difficulties with single zone homogeneous leptonic synchrotron self-Compton (SSC) radiative scenarios due to the unprecedented width of the high-energy component when compared to the lower-energy component. The two other LBL objects also detected at VHE appear to have similar modeling difficulties. Nevertheless, VHE $\gamma$ rays produced in the extended jet could account for the VHE flux observed by H.E.S.S.

We report results on multiband observations from radio to gamma-rays of the two radio-loud narrow-line Seyfert 1 (NLSy1) galaxies PKS 2004-447 and J1548+3511. Both sources show a core-jet structure on parsec scale, while they are unresolved at the arcsecond scale. The high core dominance and the high variability brightness temperature make these NLSy1 galaxies good gamma-ray source candidates. Fermi-LAT detected gamma-ray emission only from PKS 2004-447, with a gamma-ray luminosity comparable to that observed in blazars. No gamma-ray emission is observed for J1548+3511. Both sources are variable in X-rays. J1548+3511 shows a hardening of the spectrum during high activity states, while PKS 2004-447 has no spectral variability. A spectral steepening likely related to the soft excess is hinted below 2 keV for J1548+3511, while the X-ray spectra of PKS 2004-447 collected by XMM-Newton in 2012 are described by a single power-law without significant soft excess. No additional absorption above the Galactic column density or the presence of an Fe line is detected in the X-ray spectra of both sources.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10{\deg}), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

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

The strong variability of blazars can be characterized by power spectral densities (PSDs) and Fourier frequency-dependent time lags. In previous work, we created a new theoretical formalism for describing the PSDs and time lags produced via a combination of stochastic particle injection and emission via the synchrotron, synchrotron self-Compton, and external Compton (EC) processes. This formalism used the Thomson cross section and simple $\delta$-function approximations to model the synchrotron and Compton emissivities. Here we expand upon this work, using the full Compton cross section and detailed and accurate emissivities. Our results indicate good agreement between the PSDs computed using the $\delta$-function approximations and those computed using the accurate expressions, provided the observed photons are produced primarily by electrons with energies exceeding the lower limit of the injected particle population. Breaks are found in the PSDs at frequencies corresponding to the cooling timescales of the electrons primarily responsible for the observed emission, and the associated time lags are related to the difference in electron cooling timescales between the two energy channels, as expected. If the electron cooling timescales can be determined from the observed time lags and/or the observed EC PSDs, then one could in principle use the method developed here to determine the energy of the external seed photon source for EC, which is an important unsolved problem in blazar physics.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

We performed a systematic X-ray study of eight nearby $\gamma$-ray bright radio galaxies with {\em Suzaku} for understanding the origin of their X-ray emissions. The {\em Suzaku} spectra for five of those have been presented previously, while the remaining three (M\,87, PKS\,0625$-$354, and 3C\,78) are presented here for the first time. Based on the Fe-K line strength, X-ray variability, and X-ray power-law photon indices, and using additional information on the [O III] line emission, we argue for a jet origin of the observed X-ray emission in these three sources. We also analyzed five years of {\em Fermi} Large Area Telescope (LAT) GeV gamma-ray data on PKS\,0625$-$354 and 3C\,78 to understand these sources within the blazar picture. We found significant $\gamma$-ray variability in the former object. Overall, we note that the {\em Suzaku} spectra for both PKS\,0625$-$354 and 3C\,78 are rather soft, while the LAT spectra are unusually hard when compared with other $\gamma$-ray detected low-power (FR\,I) radio galaxies. We demonstrate that the constructed broad-band spectral energy distributions of PKS\,0625$-$354 and 3C\,78 are well described by a one-zone synchrotron/synchrotron self-Compton model. The results of the modeling indicate lower bulk Lorentz factors compared to those typically found in other BL Lac objects, but consistent with the values inferred from modeling other LAT-detected FR\,I radio galaxies. Interestingly, the modeling also implies very high peak ($\sim 10^{16}$\,Hz) synchrotron frequencies in the two analyzed sources, contrary to previously-suggested scenarios for FR I/BL Lac unification. We discuss the implications of our findings in the context of the FR\,I/BL Lac unification schemes.

Before the launch of the Fermi satellite only two classes of Active Galactic Nuclei (AGN) were known to generate relativistic jets and thus to emit up to the gamma-ray energy range: blazars and radio galaxies, both hosted in giant elliptical galaxies. The first four years of observations by the Large Area Telescope (LAT) on board Fermi confirmed that these two populations represent the most numerous identified sources in the extragalactic gamma-ray sky, but the discovery of variable gamma-ray emission from 5 radio-loud Narrow-Line Seyfert 1 (NLSy1) galaxies revealed the presence of a possible emerging third class of AGN with relativistic jets. Considering that NLSy1 are thought to be hosted in spiral galaxies, this finding poses intriguing questions about the nature of these objects, the knowledge of the development of relativistic jets, and the evolution of radio-loud AGN. In this context, the study of the radio-loud NLSy1 from radio to gamma-rays has received increasing attention. Here we discuss the radio-to-gamma-rays properties of the gamma-ray emitting NLSy1, also in comparison with the blazar scenario.

High redshift blazars are among the most powerful objects in the Universe. Although they represent a significant fraction of the extragalactic hard X-ray sky, they are not commonly detected in gamma-rays. High redshift (z>2) objects represent <10 per cent of the AGN population observed by Fermi so far, and gamma-ray flaring activity from these sources is even more uncommon. The characterization of the radio-to-gamma-ray properties of high redshift blazars represent a powerful tool for the study of both the energetics of such extreme objects and the Extragalactic Background Light. We present results of a multi-band campaign on TXS 0536+145, which is the highest redshift flaring gamma-ray blazar detected so far. At the peak of the flare the source reached an apparent isotropic gamma-ray luminosity of 6.6x10^49 erg/s, which is comparable with the luminosity observed from the most powerful blazars. The physical properties derived from the multi-wavelength observations are then compared with those shown by the high redshift population. In addition preliminary results from the high redshift flaring blazar PKS 2149-306 will be discussed.

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

We report on multifrequency observations performed during 2012 December-2013 August of the first narrow-line Seyfert 1 galaxy detected in gamma rays, PMN J0948+0022 ($z$ = 0.5846). A gamma-ray flare was observed by the Large Area Telescope on board Fermi during 2012 December-2013 January, reaching a daily peak flux in the 0.1-100 GeV energy range of (155 $\pm$ 31) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 2013 January 1, corresponding to an apparent isotropic luminosity of about 1.5$\times$10$^{48}$ erg s$^{-1}$. The gamma-ray flaring period triggered Swift and VERITAS observations in addition to radio and optical monitoring by OVRO, MOJAVE, and CRTS. A strong flare was observed in optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the gamma-ray flare, reaching a record flux for this source from optical to gamma rays. VERITAS observations at very high energy (E > 100 GeV) during 2013 January 6-17 resulted in an upper limit of F (> 0.2 TeV) < 4.0$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$. We compared the spectral energy distribution (SED) of the flaring state in 2013 January with that of an intermediate state observed in 2011. The two SEDs, modelled as synchrotron emission and an external Compton scattering of seed photons from a dust torus, can be modelled by changing both the electron distribution parameters and the magnetic field.

We performed a systematic X-ray study of eight nearby gamma-ray bright radio galaxies with Suzaku for understanding the origin of their X-ray emissions. The Suzaku spectra for five of those have been presented previously, while the remaining three (M87, PKS0625-354, and 3C78) are presented here for the first time. Based on the Fe-K line strength, X-ray variability, and X-ray power-law photon indices, and using additional information on the [O III] line emission, we argue for a jet origin of the observed X-ray emission in these three sources. We also analyzed five years of Fermi Large Area Telescope (LAT) GeV gamma-ray data on PKS0625-354 and 3C78 to understand these sources within the blazar picture. We found significant gamma-ray variability in the former object. Overall, we note that the Suzaku spectra for both PKS0625-354 and 3C78 are rather soft, while the LAT spectra are unusually hard when compared with other gamma-ray detected low-power (FR I) radio galaxies. We demonstrate that the constructed broad-band spectral energy distributions of PKS0625-354 and 3C78 are well described by a one-zone synchrotron/synchrotron self-Compton model. The results of the modeling indicate lower bulk Lorentz factors compared to those typically found in other BL Lac objects, but consistent with the values inferred from modeling other LAT-detected FR\,I radio galaxies. Interestingly, the modeling also implies very high peak ($\sim 10^{16}$ Hz) synchrotron frequencies in the two analyzed sources, contrary to previously-suggested scenarios for FR I/BL Lac unification. We discuss the implications of our findings in the context of the FR I/BL Lac unification schemes.

We report results of a multiband monitoring campaign of the flat spectrum radio quasar TXS 0536+145 at redshift 2.69. This source was detected during a very high gamma-ray activity state in 2012 March by the Large Area Telescope on board Fermi, becoming the gamma-ray flaring blazar at the highest redshift detected so far. At the peak of the flare the source reached an apparent isotropic gamma-ray luminosity of 6.6 x 10^49 erg/s which is comparable to the values achieved by the most luminous blazars. This activity triggered radio-to-X-rays monitoring observations by Swift, Very Long Baseline Array, European VLBI Network, and Medicina single-dish telescope. Significant variability was observed from radio to X-rays supporting the identification of the gamma-ray source with TXS 0536+145. Both the radio and gamma-ray light curves show a similar behaviour, with the gamma-rays leading the radio variability with a time lag of about 4-6 months. The luminosity increase is associated with a flattening of the radio spectrum. No new superluminal component associated with the flare was detected in high resolution parsec-scale radio images. During the flare the gamma-ray spectrum seems to deviate from a power law, showing a curvature that was not present during the average activity state. The gamma-ray properties of TXS 0536+145 are consistent with those shown by the high-redshift gamma-ray blazar population.

Context. Blazars are a subset of active galactic nuclei (AGN) with jets that are oriented along our line of sight. Variability and spectral energy distribution (SED) studies are crucial tools for understanding the physical processes responsible for observed AGN emission. Aims. We report peculiar behaviour in the bright gamma-ray blazar PKS 1424-418 and use its strong variability to reveal information about the particle acceleration and interactions in the jet. Methods. Correlation analysis of the extensive optical coverage by the ATOM telescope and nearly continuous gamma-ray coverage by the Fermi Large Area Telescope is combined with broadband, time-dependent modeling of the SED incorporating supplemental information from radio and X-ray observations of this blazar. Results. We analyse in detail four bright phases at optical-GeV energies. These flares of PKS 1424-418 show high correlation between these energy ranges, with the exception of one large optical flare that coincides with relatively low gamma-ray activity. Although the optical/gamma-ray behaviour of PKS 1424-418 shows variety, the multiwavelength modeling indicates that these differences can largely be explained by changes in the flux and energy spectrum of the electrons in the jet that are radiating. We find that for all flares the SED is adequately represented by a leptonic model that includes inverse Compton emission from external radiation fields with similar parameters. Conclusions. Detailed studies of individual blazars like PKS 1424-418 during periods of enhanced activity in different wavebands are helping us identify underlying patterns in the physical parameters in this class of AGN.

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

Blazars display strong variability on multiple timescales and in multiple radiation bands. Their variability is often characterized by power spectral densities (PSDs) and time lags plotted as functions of the Fourier frequency. We develop a new theoretical model based on the analysis of the electron transport (continuity) equation, carried out in the Fourier domain. The continuity equation includes electron cooling and escape, and a derivation of the emission properties includes light travel time effects associated with a radiating blob in a relativistic jet. The model successfully reproduces the general shapes of the observed PSDs and predicts specific PSD and time lag behaviors associated with variability in the synchrotron, synchrotron self-Compton (SSC), and external Compton (EC) emission components, from sub-mm to gamma-rays. We discuss applications to BL Lacertae objects and to flat-spectrum radio quasars (FSRQs), where there are hints that some of the predicted features have already been observed. We also find that FSRQs should have steeper PSD power-law indices than BL Lac objects at Fourier frequencies < 10^{-4} Hz, in qualitative agreement with previously reported observations by the Fermi Large Area Telescope.

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

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

The discovery of gamma-ray emission from 5 radio-loud narrow-line Seyfert 1 galaxies revealed the presence of a possible emerging third class of AGNs with relativistic jets, in addition to blazars and radio galaxies. The existence of relativistic jets also in this subclass of Seyfert galaxies opened an unexplored research space for our knowledge of the radio-loud AGNs. Here, we discuss the radio-to-gamma-rays properties of the gamma-ray emitting narrow-line Seyfert 1 galaxies, also in comparison with the blazar scenario.

PKS2142-75 (a.k.a. 2FGL J2147.4-7534) is a flat-spectrum radio quasar that was observed quasi-simultaneously by a suite of instruments across the elec- tromagnetic spectrum during two flaring states in April 2010 and August 2011 as well as a quiescent state in December 2011 through January 2012. The results of these campaigns and model spectral energy distributions (SEDs) from the active and quiescent states are presented. The SED model parameters of PKS 2142-75 indicate that the two flares of the source are created by unique physical conditions. SED studies of flat spectrum radio quasars are beginning to indicate that there might be two types of flares, those that can be described purely by changes in the electron distribution and those that require changes in other parameters, such as the magnetic field strength or the size of the emitting region.

Fermi has provided the largest sample of gamma-ray selected blazars to date. In this work we use a uniformly selected set of 211 BL Lacertae (BL Lac) objects detected by it Fermi during its first year of operation. We have obtained redshift constraints for 206 out of the 211 BL Lacs in our sample making it the largest and most complete sample of BL Lacs available in the literature. We use this sample to determine the luminosity function of BL Lacs and its evolution with cosmic time. We find that for most BL Lac classes, the evolution is positive with a space density peaking at modest redshift (z~1.2). The low-luminosity, high-synchrotron peaked (HSP) BL Lacs are an exception, showing strong negative evolution, with number density increasing for z$\lesssim$0.5. Since this rise corresponds to a drop-off in the density of flat-spectrum radio quasars (FSRQs), a possible interpretation is that these HSPs represent an accretion-starved end-state of an earlier merger-driven gas-rich phase. We additionally find that the known BL Lac correlation between luminosity and photon spectral index persists after correction for the substantial observational selection effects with implications for the so called `blazar sequence'. Finally, estimating the beaming corrections to the luminosity function, we find that BL Lacs have an average Lorentz factor of $\gamma=6.1^{+1.1}_{-0.8}$, and that most are seen within 10$^{\circ}$ of the jet axis.

Before the launch of the Fermi satellite only two classes of AGNs were known to produce relativistic jets and thus emit up to the gamma-ray energy range: blazars and radio galaxies, both hosted in giant elliptical galaxies. The first four years of observations by the Large Area Telescope on board Fermi confirmed that these two are the most numerous classes of identified sources in the extragalactic gamma-ray sky, but the discovery of gamma-ray emission from 5 radio-loud narrow-line Seyfert 1 galaxies revealed the presence of a possible emerging third class of AGNs with relativistic jets. Considering that narrow-line Seyfert 1 galaxies seem to be typically hosted in spiral galaxy, this finding poses intriguing questions about the nature of these objects, the onset of production of relativistic jets, and the cosmological evolution of radio-loud AGN. Here, we discuss the radio-to-gamma-rays properties of the gamma-ray emitting narrow-line Seyfert 1 galaxies, also in comparison with the blazar scenario.

The narrow-line Seyfert 1 galaxy SBS 0846+513 was first detected by the Large Area Telescope (LAT) on-board Fermi in 2011 June-July when it underwent a period of flaring activity. Since then, as Fermi continues to accumulate data on this source, its flux has been monitored on a daily basis. Two further gamma-ray flaring episodes from SBS 0846+513 were observed in 2012 May and August, reaching a daily peak flux integrated above 100 MeV of (50+/-12)x10^-8 ph/cm^2/s, and (73+/-14)x10^-8 ph/cm^2/s on May 24 and August 7, respectively. Three outbursts were detected at 15 GHz by the Owens Valley Radio Observatory 40-m telescope in 2012 May, 2012 October, and 2013 January, suggesting a complex connection with the gamma-ray activity. The most likely scenario suggests that the 2012 May gamma-ray flare may not be directly related to the radio activity observed over the same period, while the two gamma-ray flaring episodes may be related to the radio activity observed at 15 GHz in 2012 October and 2013 January. The gamma-ray flare in 2012 May triggered Swift observations that confirmed that SBS 0846+513 was also exhibiting high activity in the optical, UV and X-ray bands, thus providing a firm identification between the gamma-ray source and the lower-energy counterpart. We compared the spectral energy distribution (SED) of the flaring state in 2012 May with that of a quiescent state. The two SEDs, modelled as an external Compton component of seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. No significant evidence of thermal emission from the accretion disc has been observed. Interestingly, in the 5 GHz radio luminosity vs. synchrotron peak frequency plot SBS 0846+513 seems to lie in the flat spectrum radio quasar part of the so-called `blazar sequence'.

The first statistically significant detection of the cosmic \gamma-ray horizon (CGRH) that is independent of any extragalactic background light (EBL) model is presented. The CGRH is a fundamental quantity in cosmology. It gives an estimate of the opacity of the Universe to very high energy (VHE) \gamma-ray photons due to photon-photon pair production with the EBL. The only estimations of the CGRH to date are predictions from EBL models and lower limits from \gamma-ray observations of cosmological blazars and \gamma-ray bursts. Here, we present homogeneous synchrotron/synchrotron self-Compton (SSC) models of the spectral energy distributions of 15 blazars based on (almost) simultaneous observations from radio up to the highest energy \gamma-rays taken with the Fermi satellite. These synchrotron/SSC models predict the unattenuated VHE fluxes, which are compared with the observations by imaging atmospheric Cherenkov telescopes. This comparison provides an estimate of the optical depth of the EBL, which allows a derivation of the CGRH through a maximum likelihood analysis that is EBL-model independent. We find that the observed CGRH is compatible with the current knowledge of the EBL.

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

Gamma rays from distant blazars interact with the extragalactic background light, creating electron-positron pairs, and reducing the gamma-ray flux measured by ground-based atmospheric Cherenkov gamma-ray telescopes. These pairs can Compton-scatter the cosmic microwave background, creating a gamma-ray signature observable by the Fermi Large Area Telesope (LAT). The signature is also dependent on the intergalactic magnetic field (IGMF), since it can deflect the pairs from our line of sight, reducing the gamma-ray emission. We present preliminary constraints on the IGMF using Fermi-LAT and Cherenkov telescope observations, ruling out both very large and very small values of the IGMF strength.

Blazars are the most plentiful gamma-ray source at GeV energies, and despite detailed study, there is much that is not known about these sources. In this review I explore some recent results on blazars, including the controversy of the "blazar sequence", the curvature in the LAT spectra, and the location along the jet of the gamma-ray emitting region. I conclude with a discussion of alternative modeling possibilities.

We report on multiwavelength observations of the blazar PKS 0537-441 (z = 0.896) obtained from microwaves through gamma rays by SMA, REM, ATOM, Swift and Fermi during 2008 August-2010 April. Strong variability has been observed in gamma rays, with two major flaring episodes (2009 July and 2010 March) and a harder-when-brighter behaviour, quite common for FSRQs and low-synchrotron-peaked BL Lacs, in 2010 March. In the same way the SED of the source cannot be modelled by a simple synchrotron self-Compton model, as opposed to many BL Lacs, but the addition of an external Compton component of seed photons from a dust torus is needed. The 230 GHz light curve showed an increase simultaneous with the gamma-ray one, indicating co-spatiality of the mm and gamma-ray emission region likely at large distance from the central engine. The low, average, and high activity SED of the source could be fit changing only the electron distribution parameters, but two breaks in the electron distribution are necessary. The ensuing extra spectral break, located at NIR-optical frequencies, together with that in gamma rays seem to indicate a common origin, most likely due to an intrinsic feature in the underlying electron distribution. An overall correlation between the gamma-ray band with the R-band and K-band has been observed with no significant time lag. On the other hand, when inspecting the light curves on short time scales some differences are evident. In particular, flaring activity has been detected in NIR and optical bands with no evident gamma-ray counterparts in 2009 September and November. Moderate variability has been observed in X-rays with no correlation between flux and photon index. An increase of the detected X-ray flux with no counter part at the other wavelengths has been observed in 2008 October, suggesting once more a complex correlation between the emission at different energy bands.

Before the launch of the Fermi satellite only two classes of AGN were known to generate relativistic jets and thus emit up to the gamma-ray energy range: blazars and radio galaxies, both hosted in giant elliptical galaxies. The first two years of observations by the Large Area Telescope (LAT) on board Fermi confirmed that these two are the most numerous classes of identified sources in the extragalactic gamma-ray sky, but the discovery of variable gamma-ray emission from 5 radio-loud Narrow-Line Seyfert 1 galaxies (NLSy1s) revealed the presence of a possible emerging third class of AGN with relativistic jets. Considering also that NLSy1s are typically hosted in spiral galaxy, this finding poses intriguing questions about the nature of these objects, the onset of production of relativistic jets, and the cosmological evolution of radio-loud AGN. Here, we report on a preliminary investigation of the properties of this sample of radio-loud NLSy1 at MeV-GeV photon energies, utilizing the four-year accumulation of Fermi LAT data. In addition we briefly discuss some radio-to-gamma-rays properties of the gamma-ray emitting NLSy1 in the context of the blazar scenario.

The blazar sequence is reflected in a correlation of the peak luminosity versus peak frequency of the synchrotron component of blazars. This correlation has been considered one of the fundamental pieces of evidence for the existence of a continuous sequence that includes low-power BL Lacertae objects through high-power flat spectrum radio quasars. A correlation between the Compton domi- nance, the ratio of the Compton to synchrotron luminosities, and the peak synchrotron frequency is another piece of evidence for the existence of the blazar sequence explored by Fossati et al. (1998). Since that time, however, it has essentially been ignored. We explore this correlation with a sample based on the second LAT AGN catalog (2LAC), and show that is is particularly important, since it is independent of redshift. We reproduce the trends in our sample with a simple model that in- cludes synchrotron and Compton cooling in the slow- and fast-cooling regimes, and angle-dependent effects.

Does the "blazar sequence" exist, or is it a result of a selection effect, due to the difficulty in measuring the redshifts of blazars with both high synchrotron peak frequencies (\gtrsim 10^{15} Hz) and luminosities (\gtrsim 10^{46} erg s^{-1})? We explore this question with a sample of blazars from the Second Catalog of Active Galactic Nuclei (AGN) from the Fermi Large Area Telescope (LAT). The Compton dominance, the ratio of the peak of the Compton to the synchrotron peak luminosities, is essentially a redshift-independent quantity, and thus crucial to answering this question. We find that a correlation exists between Compton dominance and the peak frequency of the synchrotron component for all blazars in the sample, including ones with unknown redshift. We then construct a simple model to explain the blazar properties in our sample, where the difference between sources is due to only the magnetic field of the blazar jet emitting region, the external radiation field energy density, and the jet angle to the line of sight, with the magnetic field strength and external energy density being correlated. This model can reproduce the trends of the blazars in the sample, and predicts blazars may be discovered in the future with high synchrotron peak frequencies and luminosities. At the same time the simple model reproduces the lack of high-synchrotron peaked blazars with high Compton dominances (\gtrsim 1).

The presence of a relativistic jet in some radio-loud Narrow-Line Seyfert 1s (NLSy1) galaxies, first suggested by their variable radio emission and the flat radio spectra, is now confirmed by the Fermi-LAT detection of five NLSy1s in gamma rays. In particular, a strong gamma-ray flare from SBS 0846+513 was observed in 2011 June by Fermi-LAT reaching a gamma-ray luminosity (0.1-300 GeV) of about 10^48 erg/s, comparable to that of bright flat spectrum radio quasars. Apparent superluminal velocity in the jet was inferred from 2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. Both the power released by this object during the flaring activity and the apparent superluminal velocity are strong indicators of the presence of a relativistic jet as powerful as those in blazars. In addition, variability and spectral properties in radio and gamma-ray bands indicate a blazar-like behaviour, suggesting that, except for some distinct optical characteristics, SBS 0846+513 could be considered as a young blazar at the low end of the blazar's black hole mass distribution.

The flat spectrum radio quasar (FSRQ) PKS 2123-463 was associated in the First Fermi-LAT source catalog with the gamma-ray source 1FGL J2126.1-4603, but when considering the full first two years of Fermi observations, no gamma-ray source at a position consistent with this FSRQ was detected, and thus PKS 2123-463 was not reported in the Second Fermi-LAT source catalog. On 2011 December 14 a gamma-ray source positionally consistent with PKS 2123-463 was detected in flaring activity by Fermi-LAT. This activity triggered radio-to-X-ray observations by the Swift, GROND, ATCA, Ceduna, and KAT-7 observatories. Results of the localization of the gamma-ray source over 41 months of Fermi-LAT operation are reported here in conjunction with the results of the analysis of radio, optical, UV and X-ray data collected soon after the gamma-ray flare. The strict spatial association with the lower energy counterpart together with a simultaneous increase of the activity in optical, UV, X-ray and gamma-ray bands led to a firm identification of the gamma-ray source with PKS 2123-463. A new photometric redshift has been estimated as z = 1.46+/-0.05 using GROND and Swift/UVOT observations, in rough agreement with the disputed spectroscopic redshift of z = 1.67. We fit the broadband spectral energy distribution with a synchrotron/external Compton model. We find that a thermal disk component is necessary to explain the optical/UV emission detected by Swift/UVOT. This disk has a luminosity of about 1.8x10^46 erg/s, and a fit to the disk emission assuming a Schwarzschild (i.e., nonrotating) black hole gives a mass of about 2x10^9 solar masses. This is the first black hole mass estimate for this source.

We report Fermi-LAT observations of the radio-loud AGN SBS 0846+513 (z=0.5835), optically classified as a Narrow-Line Seyfert 1 galaxy, together with new and archival radio-to-X-ray data. The source was not active at gamma-ray energies during the first two years of Fermi operation. A significant increase in activity was observed during 2010 October-2011 August. In particular a strong gamma-ray flare was observed in 2011 June reaching an isotropic gamma-ray luminosity (0.1-300 GeV) of 1.0x10^48 erg/s, comparable to that of the brightest flat spectrum radio quasars, and showing spectral evolution in gamma rays. An apparent superluminal velocity of (8.2+/-1.5)c in the jet was inferred from 2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. Both the power released by this object during the flaring activity and the apparent superluminal velocity are strong indications of the presence of a relativistic jet as powerful as those of blazars. In addition, variability and spectral properties in radio and gamma-ray bands indicate blazar-like behaviour, suggesting that, except for some distinct optical characteristics, SBS 0846+513 could be considered as a young blazar at the low end of the blazar's black hole mass distribution.

PKS 2142-758 is a flat spectrum radio quasar which emits few, weak but significant gamma-ray flares in the MeV through GeV energy range. The first flare occured on April 4th, 2010, when the source reached a daily flux of (1.1 \pm 0.3) * 10^-6 ph cm^-2 s^-1 (ATEL #2539) in the 100 MeV to 300 GeV range. This flux represented more than an order of magnitude increase over its quiescent flux. Since the initial flare, this source has been detected in an elevated state within the same energy range from October to November of 2010 and another period ranging from July to August of 2011. During the latest flaring period in 2011 a multi wavelength observing campaign was carried out using the Ceduna radio telescope, the Australian Telescope Compact Array (ATCA), the TANAMI VLBI Array, Swift, the Rapid Eye Mount Telescope (REM), and the Large Area Telescope (LAT) on board Fermi. These quasi-simultaneous data were used to construct a broadband SED of this object in its rare active state. We present these observations and the resulting SED and some preliminary analysis of the constraints they place on the high energy emission from this object.

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

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

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

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

We report on a detailed investigation of the gamma-ray emission from 18 broad line radio galaxies (BLRGs) based on two years of Fermi Large Area Telescope (LAT) data. We confirm the previously reported detections of 3C 120 and 3C 111 in the GeV photon energy range; a detailed look at the temporal characteristics of the observed gamma-ray emission reveals in addition possible flux variability in both sources. No statistically significant gamma-ray detection of the other BLRGs was however found in the considered dataset. Though the sample size studied is small, what appears to differentiate 3C 111 and 3C 120 from the BLRGs not yet detected in gamma-rays is the particularly strong nuclear radio flux. This finding, together with the indications of the gamma-ray flux variability and a number of other arguments presented, indicate that the GeV emission of BLRGs is most likely dominated by the beamed radiation of relativistic jets observed at intermediate viewing angles. In this paper we also analyzed a comparison sample of high accretion-rate Seyfert 1 galaxies, which can be considered radio-quiet counterparts of BLRGs, and found none were detected in gamma-rays. A simple phenomenological hybrid model applied for the broad-band emission of the discussed radio-loud and radio-quiet type 1 active galaxies suggests that the relative contribution of the nuclear jets to the accreting matter is > 1 percent on average for BLRGs, whilst <0.1 percent for Seyfert 1 galaxies.

We report Fermi/LAT observations and broad-band spectral modeling of the radio-loud active galaxy 4C +55.17 (z=0.896), formally classified as a flat-spectrum radio quasar. Using 19 months of all-sky survey Fermi/LAT data, we detect a gamma-ray continuum extending up to an observed energy of 145 GeV, and furthermore we find no evidence of gamma-ray variability in the source over its observed history. We illustrate the implications of these results in two different domains. First, we investigate the origin of the steady gamma-ray emission, where we re-examine the common classification of 4C +55.17 as a quasar-hosted blazar and consider instead its possible nature as a young radio source. We analyze and compare constraints on the source physical parameters in both blazar and young radio source scenarios by means of a detailed multiwavelength analysis and theoretical modeling of its broad-band spectrum. Secondly, we show that the gamma-ray spectrum may be formally extrapolated into the very-high energy (VHE; >= 100 GeV) range at a flux level detectable by the current generation of ground-based Cherenkov telescopes. This enables us to place constraints on models of extragalactic background light (EBL) within LAT energies and features the source as a promising candidate for VHE studies of the Universe at an unprecedented redshift of z=0.896.

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

In this paper we report on the two-year-long Fermi-LAT observation of the peculiar blazar 4C +21.35 (PKS 1222+216). This source was in a quiescent state from the start of science operations of the Fermi Gamma-ray Space Telescope in 2008 August until 2009 September, and then became more active, with gradually increasing flux and some moderately-bright flares. In 2010 April and June, 4C +21.35 underwent a very strong GeV outburst composed of several major flares characterized by rise and decay timescales of the order of a day. During the outburst, the GeV spectra of 4C +21.35 displayed a broken power-law form with spectral breaks observed near 1-3 GeV photon energies. We demonstrate that, at least during the major flares, the jet in 4C +21.35 carried a total kinetic luminosity comparable to the total accretion power available to feed the outflow. We also discuss the origin of the break observed in the flaring spectra of 4C +21.35. We show that, in principle, a model involving annihilation of the GeV photons on the He II Lyman recombination continuum and line emission of "broad line region" clouds may account for such. However, we also discuss the additional constraint provided by the detection of 4C +21.35 at 0.07-0.4 TeV energies by the MAGIC telescope, which coincided with one of the GeV flares of the source. We argue that there are reasons to believe that the $lesssim$,TeV emission of 4C +21.35 (as well as the GeV emission of the source, if co-spatial), is not likely to be produced inside the broad line region zone of highest ionization ($sim 10^{17}$,cm from the nucleus), but instead originates further away from the active center, namely around the characteristic scale of the hot dusty torus surrounding the 4C +21.35 nucleus ($sim 10^{19}$,cm).

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

Single-zone synchrotron self-Compton and external Compton models are widely used to explain broad-band Spectral Energy Distributions (SEDs) of blazars from infrared to gamma-rays. These models bear obvious similarities to the homogeneous synchrotron cloud model which is often applied to explain radio emission from individual components of parsec-scale radio jets. The parsec-scale core, typically the brightest and most compact feature of blazar radio jet, could be the source of high-energy emission. We report on ongoing work to test this hypothesis by deriving the physical properties of parsec-scale radio emitting regions of twenty bright Fermi blazars using dedicated 5-43 GHz VLBA observations and comparing these parameters to results of SED modeling.

The Extragalactic Background Light (EBL) from the infrared (IR) through the ultraviolet (UV) is dominated by emission from stars, either directly or through absorption and reradiation by dust. It can thus give information on the star formation history of the universe. However, it is difficult to measure directly due to foreground radiation fields from the Galaxy and solar system. Gamma-rays from extragalactic sources at cosmological distances (blazars and gamma-ray bursts) interact with EBL photons creating electron-positron pairs, absorbing the gamma-rays. Given the intrinsic gamma-ray spectrum of a source and its redshift, the EBL can in principle be measured. However, the intrinsic gamma-ray spectra of blazars and GRBs can vary considerably from source to source and the from the same source over short timescales. A maximum intrinsic spectrum can be assumed from theoretical grounds, to give upper limits on the EBL absorption from blazars at low redshift with very high energy (VHE) gamma-ray observations with ground-based Atmospheric Cherenkov telescopes. The Fermi-LAT observations of blazars and GRBs can probe EBL absorption at higher redshifts. The lower energy portion of the LAT spectrum of these sources is unattenuated by the EBL, so that extrapolating this to higher energies can give the maximum intrinsic spectrum for a source. Comparing this to the observed higher energy LAT spectrum will then give upper limits on the EBL absorption. For blazars which have been detected by both the Fermi-LAT and at higher energies by Cherenkov telescopes, combined LAT-VHE observations can put more stringent constraints on the low redshift EBL. These procedures above can also be reversed: for sources with an unknown redshift, a given EBL model and gamma-ray spectrum can lead to an upper limit on the source's redshift.

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

We propose a method for setting upper limits to the extragalactic background light (EBL). Our method uses simultaneous {\em Fermi}-LAT and ground-based TeV observations of blazars and is based on the assumption that the intrinsic spectral energy distribution (SED) of TeV blazars lies below the extrapolation of the {\em Fermi}-LAT SED from GeV to TeV energies. By extrapolating the {\em Fermi}-LAT spectrum, which for TeV blazars is practically unattenuated by photon-photon pair production with EBL photons, a firm upper limit on the intrinsic SED at TeV energies is provided. The ratio of the extrapolated spectrum to the observed TeV spectrum provides upper limits to the optical depth for the propagation of the TeV photons due to pair production on the EBL, which in turn sets firm upper limits to EBL models. We demonstrate our method using simultaneous observations from {\em Fermi}-LAT and ground-based TeV telescopes of the blazars \object{PKS 2155-304} and \object{1ES 1218+304}, and show that high EBL density models are disfavored. We also discuss how our method can be optimized and how {\em Fermi} and X-ray monitoring observations of TeV blazars can guide future TeV campaigns, leading to potentially much stronger constraints on EBL models.

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

The extragalactic background light (EBL) from the infrared to the ultraviolet is difficult to measure directly, but can be constrained with a variety of methods. EBL photons absorb gamma-rays from distant blazars, allowing one to use blazar spectra from atmospheric Cherenkov telescopes (ACTs) to put upper limits on the EBL by assuming a blazar source spectrum. Here we apply a simple technique, similar to the one developed by Schroedter (2005), to the most recent very-high energy (VHE) gamma-ray observations of blazars to put upper limits on the EBL energy density. This technique is independent of the EBL model and has well-defined errors on the constraints. Our results are consistent with EBL measurements and constraints but marginally inconsistent with several EBL models.

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

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

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

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.

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

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