The most powerful persistent accelerators in the Universe are jetted active galaxies. Galaxies whose jets are directed towards Earth, so called blazars, dominate the extragalactic gamma-ray sky. Still, most of the highest-energy accelerators likely elude detection. These extreme blazars, whose radiated energy can peak beyond 10 TeV, are ideal targets to study particle acceleration and radiative processes, and may provide links to cosmic rays and astrophysical neutrinos. The growing number of extreme blazars observed at TeV energies has been critical for the emergence of gamma-ray cosmology, including measurements of the extragalactic background light, tight bounds on the intergalactic magnetic field, and constraints on exotic physics at energies inaccessible with human-made accelerators. Tremendous progress is expected in the decade to come, particularly with the deployment of the Cherenkov Telescope Array.

BL Lac objects can be extreme in two ways: with their synchrotron emission, peaking beyond 1 keV in their spectral energy distribution, or with their gamma-ray emission, peaking at multi-TeV energies up to and beyond 10-20 TeV, like 1ES 0229+200. This second type of extreme BL Lacs --which we can name "TeV-peaked BL Lacs" is not well explained by the usual synchrotron self-Compton scenarios for BL Lacs. These sources are also important as probes for the intergalactic diffuse infrared background and cosmic magnetic fields, as well as possible sites of production of ultra-high-energy cosmic rays and neutrinos. However, all these studies are hindered by their still very limited number. Here I propose a new, simple criterium to select the best candidates for TeV observations, specifically aimed at this peculiar type of BL Lac objects by combining X-ray, gamma-ray and infrared data. It is based on the observation of a clustering towards a high X-ray to GeV gamma-ray flux ratio, and it does not rely on the radio flux or X-ray spectrum. This makes it suitable to find TeV-peaked sources also with very faint radio emission. Taking advantage of the Fermi all-sky gamma-ray survey applied to the ROMA-BZCAT and Sedentary Survey samples, I produce an initial list of 47 TeV-peaked candidates for observations with present and future air-Cherenkov telescopes.

The gamma-ray emission in broad-line blazars is generally explained as inverse Compton (IC) radiation of relativistic electrons in the jet scattering optical-UV photons from the Broad Line Region (BLR), the so-called BLR External Compton scenario. We test this scenario on the Fermi gamma-ray spectra of 106 broad-line blazars detected with the highest significance or largest BLR, by looking for cut-off signatures at high energies compatible with gamma-gamma interactions with BLR photons. We do not find evidence for the expected BLR absorption. For 2/3 of the sources, we can exclude any significant absorption ($\tau_{max}<1$), while for the remaining 1/3 the possible absorption is constrained to be 1.5-2 orders of magnitude lower than expected. This result holds also dividing the spectra in high and low-flux states, and for powerful blazars with large BLR. Only 1 object out of 10 seems compatible with substantial attenuation ($\tau_{max}>5$). We conclude that for 9 out of 10 objects, the jet does not interact with BLR photons. Gamma-rays seem either produced outside the BLR most of the time, or the BLR is ~100x larger than given by reverberation mapping. This means that i) External Compton on BLR photons is disfavoured as the main gamma-ray mechanism, vs IC on IR photons from the torus or synchrotron self-Compton; ii) the Fermi gamma-ray spectrum is mostly intrinsic, determined by the interaction of the particle distribution with the seed-photons spectrum; iii) without suppression by the BLR, broad-line blazars can become copious emitters above 100 GeV, as demonstrated by 3C454.3. We expect the CTA sky to be much richer of broad-line blazars than previously thought.

Hard-TeV BL Lacs are a new type of blazars characterized by a hard intrinsic TeV spectrum, locating the peak of their gamma-ray emission in the spectral energy distribution (SED) above 2-10 TeV. Such high energies are problematic for the Compton emission, using a standard one-zone leptonic model. We study six examples of this new type of BL Lacs in the hard X-ray band with the NuSTAR satellite. Together with simultaneous observations with the SWIFT satellite, we fully constrain the peak of the synchrotron emission in their SED, and test the leptonic synchrotron self-Compton (SSC) model. We confirm the extreme nature of 5 objects also in the synchrotron emission. We do not find evidence of additional emission components in the hard X-ray band. We find that a one-zone SSC model can in principle reproduce the extreme properties of both peaks in the SED, from X-ray up to TeV energies, but at the cost of i) extreme electron energies with very low radiative efficiency, ii) conditions heavily out of equipartition (by 3 to 5 orders of magnitude), and iii) not accounting for the simultaneous UV data, which then should belong to a different emission component, possibly the same as the far-IR (WISE) data. We find evidence of this separation of the UV and X-ray emission in at least two objects. In any case, the TeV electrons must not "see" the UV or lower-energy photons, even if coming from different zones/populations, or the increased radiative cooling would steepen the VHE spectrum.

The Fermi-LAT revealed that the census of the gamma-ray sky is dominated by blazars. Looking for a possible connection between radio and gamma-ray emission is a central issue for understanding the blazar physics, and various works were dedicated to this topic. However, while a strong and significant correlation was found between radio and gamma-ray emission in the 0.1-100 GeV energy range, the connection between radio and very high energy (VHE, E>0.1 TeV) emission is still elusive. The main reason is the lack of a homogeneous VHE sky coverage, due to the operational mode of the imaging atmospheric Cherenkov telescopes. With the present work we aim to quantify and assess the significance of the possible connection between high-resolution radio emission, on milliarcsecond scale, and GeV-TeV gamma-ray emission in blazars. For achieving our goal we extract two large and unbiased blazar samples from the 1FHL and 2FHL Fermi catalogs, above 10 GeV and 50 GeV, respectively. To investigate how the correlation evolves as the gamma-ray energy increases, we perform the same analysis by using the 0.1-300 GeV 3FGL gamma-ray energy fluxes. When we consider the 0.1-300 GeV gamma-ray energy range, we find a strong and significant correlation for all of the blazar sub-classes. Conversely, when we consider the gamma-ray emission above 10 GeV the correlation with the radio emission vanishes, with the exception of the blazar sub-class of high synchrotron peaked objects.

The Fermi Large Area Telescope (LAT) revealed that blazars, representing the most extreme radio-loud active galactic nuclei (AGN) population, dominate the census of the gamma-ray sky, and a significant correlation was found between radio and gamma-ray emission in the 0.1-100 GeV energy range. However, the possible connection between radio and very high energy (VHE, E>0.1 TeV) emission still remains elusive, owing to the lack of a homogeneous coverage of the VHE sky. The main goal of this work is to quantify and assess the significance of a possible connection between the radio emission on parsec scale measured by the very long baseline interferometry (VLBI) and GeV-TeV gamma-ray emission in blazars, which is a central issue for understanding the blazar physics and the emission processes. We investigate the radio VLBI and high energy gamma-ray emission by using two large and unbiased AGN samples extracted from the first and second Fermi-LAT catalogs of hard gamma-ray sources detected above 10 GeV (1FHL) and 50 GeV (2FHL). For comparison, we perform the same correlation analysis by using the 0.1-300 GeV gamma-ray energy flux provided by the third Fermi-LAT source catalog. We find that the correlation strength and significance depend on the gamma-ray energy range with a different behavior among the blazar sub-classes. Overall, the radio and gamma-ray emission above 10 GeV turns out to be uncorrelated for the full samples and for all of the blazar sub-classes with the exception of high synchrotron peaked (HSP) objects, which show a strong and significant correlation. On the contrary, when 0.1-300 GeV gamma-ray energies are considered, a strong and significant correlation is found for the full blazar sample as well as for all of the blazar sub-classes. We interpret and explain this correlation behavior within the framework of the blazar spectral energy distribution properties.

We revisit the blazar sequence exploiting the complete, flux limited sample of blazars with known redshift detected by the Fermi satellite after 4 years of operations (the 3LAC sample). We divide the sources into gamma-ray luminosity bins, collect all the archival data for all blazars, and construct their spectral energy distribution (SED). We describe the average SED of blazars in the same luminosity bin through a very simple, completely phenomenological function consisting of two broken power laws connecting with a power law of fixed slope describing the radio emission. We do that separately for BL Lacs and for flat spectrum radio quasars (FSRQs) and also for all blazars together. The main results are: i) FSRQs display approximately the same SED as the luminosity increases, except for the fact that the relative importance of the high energy peak increases; ii) as a consequence, X-ray spectra of FSRQs become harder for larger luminosities; iii) BL Lacs form indeed a sequence: they become redder (i.e. the peak frequencies becomes smaller) for increasing luminosities, with a steeper gamma-ray slope and a larger dominance of the high energy peak; iv) for all blazars (BL Lacs+FSRQs) these properties becomes more prominent, as the highest luminosity bin is populated mostly by FSRQs and the lowest luminosity bin mostly by BL Lacs. This agrees with the original blazar sequence, although BL Lacs never have an average gamma-ray slope as hard as found in the original sequence. v) At high luminosities, a large fraction of FSRQs shows signs of thermal emission from the accretion disc, contributing in the optical-UV.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of radio-loud Active Galactic Nuclei. For a summary, we refer to the paper.

Blazars are highly variable sources over timescales that can be as low as minutes. This is the case of the High Energy Peaked BL Lac (HBL) objects showing strong variability in X-rays, which highly correlate with that of the TeV emission. The degree of this correlation is still debated, particularly when the flaring activity is followed down to very short time scales. This correlation could challenge the synchrotron-self-Compton scenario in which one relativistic electron population dominates the entire radiative output. We argue that the LOFT Large Area Detector (10 m$^2$, LAD), thanks to its unprecedented timing capability, will allow us to detect the X-ray counterpart (2-50 keV) of the very fast variability observed at TeV energies, sheding light on the nature of X-TeV connection. We will discuss the test case of PKS 2155-304, showing as it would be possible to look for any X-ray variability occurring at very short timescales, never explored so far. This will put strong constraints on the size and the location of any additional electron population in the multi-zone scenario. Under this perspective, LOFT and the CTA observatories, planned to operate in the same time frame, will allow us to investigate in depth the connection between X-ray and TeV emissions. We also discuss the potentialities of LOFT in measuring the change in spectral curvature of the synchrotron spectra in HBLs which will make possible to directly study the mechanism of acceleration of highly energetic electrons. LOFT timing capability will be also promising in the study of Flat Spectrum Radio Quasars (FSRQs) with flux $\ge 1$ mCrab. Constraints to the location of the high energy emission will be given by temporal investigation on second timescale and spectral trend analysis on minute timescales. This represents a further link with CTA because of the rapid (unexpected) TeV emission recently detected in some FSRQs.

Recent observations of blazars at high energy (HE, 0.1-100 GeV) and very high energy (VHE, >0.1 TeV) have provided important constraints on the intensity and spectrum of the diffuse Extragalactic Background Light (EBL), shedding light on its main origin. Several issues remain open, however, in particular in the mid- and far-infrared bands and in the blazar emission at multi-TeV energies. This review summarizes the observational and theoretical progress in the study of the EBL with gamma-rays and the most promising future improvements, which are mainly expected from spectra in the multi-TeV range.

The high-energy-peaked BL Lac H 2356-309 (z=0.165) was detected by HESS at very high energies (VHE, >100 GeV) with relatively high significance in the years 2004-2007, allowing a good determination of its gamma-ray spectrum. After correction for the interaction with the diffuse extragalactic background light (EBL), the VHE spectrum is flat (Gamma~1.9-2) over a decade in energy, locating the gamma-ray peak around or above 0.6-1 TeV. This is remarkably at odds with the interpretation and modeling provided by HESS, which do not correspond to the source properties and can be excluded with high confidence. The overall GeV-to-TeV characteristics of H 2356-309 seem intermediate between the "TeV-peaked" (Fermi-faint) and "100 GeV-peaked" (Fermi-bright) BL Lac objects, and difficult to reconcile with the shape of the synchrotron emission in a single-zone SSC scenario.

Fermi-LAT spectra at high energies (HE, 0.1-100 GeV) are often extrapolated to very high energies (VHE, >100 GeV) and considered either a good estimate or an upper limit for the blazars intrinsic VHE spectrum. This assumption seems not well justified, neither theoretically nor observationally. Besides being often softer, observations do indicate that spectra at VHE could be also harder than at HE, even when adopting the limit of Gamma=1.5. Results based on such straightforward GeV-TeV extrapolations are in general not reliable, and should be considered with caution.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

Very high energy (VHE) gamma-rays from extragalactic sources, interacting by gamma-gamma collisions with diffuse intergalactic radiation fields, provide an alternative way to constrain the diffuse background light, completely independent of direct measurements. The limits depend however on our knowledge of the physics of the gamma-ray sources. After clarifying the interplay between background light and VHE spectra, I summarize the extent and validity of the obtainable limits, and where future improvements can be expected.

1ES 0414+009 (z = 0.287) is a distant high-frequency-peaked BL Lac object, and has long been considered a likely emitter of very-high energy (VHE, E>100 GeV) gamma-rays due to its high X-ray and radio flux. Observations in the VHE gamma-ray band and across the electromagnetic spectrum can provide insights into the origin of highly energetic particles present in the source and the radiation processes at work. Because of the distance of the source, the gamma-ray spectrum might provide further limits on the level of the Extragalactic Background Light (EBL). We report observations made between October 2005 and December 2009 with H.E.S.S., an array of four imaging atmospheric Cherenkov telescopes. Observations at high energies (HE, 100 MeV - 100 GeV) with the Fermi-LAT instrument in the first 20 months of its operation are also reported. To complete the multi-wavelength picture, archival UV and X-ray observations with the Swift satellite and optical observations with the ATOM telescope are also used. Based on the observations with H.E.S.S., 1ES 0414+009 is detected for the first time in the VHE band. An excess of 224 events is measured, corresponding to a significance of 7.8 sigma. The photon spectrum of the source is well described by a power law, with photon index of 3.45 \pm 0.25stat \pm 0.20syst. The integral flux above 200 GeV is (1.88 \pm 0.20stat \pm 0.38syst) \times10-12 cm-2 s-1. Observations with the Fermi-LAT in the first 20 months of operation show a flux between 200 MeV and 100 GeV of (2.3 \pm 0.2stat) \times 10-9 erg cm-2 s-1, and a spectrum well described by a power-law function with a photon index 1.85 \pm 0.18. Swift/XRT observations show an X-ray flux between 2 and 10 keV of (0.8 - 1) \times 10-11 erg cm-2 s-1, and a steep spectrum (2.2 - 2.3). Combining X-ray with optical-UV data, a fit with a log-parabolic function locates the synchrotron peak around 0.1 keV. ...

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

We study the applicability of the idea of internal absorption of gamma-rays produced through synchrotron radiation of ultrarelativistic protons in highly magnetized blobs to 1ES 0229+200 and 3C 66A, the two TeV blazars which show unusually hard intrinsic gamma-ray spectra after being corrected for the intergalactic absorption. We show that for certain combinations of reasonable model parameters, even with quite modest energy requirements, the scenario allows a self-consistent explanation of the non-thermal emission of these objects in the keV, GeV, and TeV energy bands.

The high energy peaked BL Lac (HBL) object 1ES 0414+009 (z=0.287) is a distant very high-energy (VHE, E > 100 GeV) blazars with well-determined redshift. This source was detected with the High Energy Stereoscopic System (H.E.S.S.) between October 2005 and September 2009. It was also detected with the Fermi Large Area Telescope (LAT) in 21 months of data. The combined high energy (HE) and VHE spectra, once corrected for gamma-gamma absorption on the extragalactic background light (EBL), indicate a Compton peak located above few TeV, among the highest in the BL Lac class.

Observations of very high energy gamma-rays from blazars provide information about acceleration mechanisms occurring in their innermost regions. Studies of variability in these objects allow a better understanding of the mechanisms at play. To investigate the spectral and temporal variability of VHE (>100 GeV) gamma-rays of the well-known high-frequency-peaked BL Lac object PKS 2155-304 with the H.E.S.S. imaging atmospheric Cherenkov telescopes over a wide range of flux states. Data collected from 2005 to 2007 are analyzed. Spectra are derived on time scales ranging from 3 years to 4 minutes. Light curve variability is studied through doubling timescales and structure functions, and is compared with red noise process simulations. The source is found to be in a low state from 2005 to 2007, except for a set of exceptional flares which occurred in July 2006. The quiescent state of the source is characterized by an associated mean flux level of 4.32 +/-0.09 x 10^-11 cm^-2 s^-1 above 200 GeV, or approximately 15% of the Crab Nebula, and a power law photon index of 3.53 +/-0.06. During the flares of July 2006, doubling timescales of ~2 min are found. The spectral index variation is examined over two orders of magnitude in flux, yielding different behaviour at low and high fluxes,which is a new phenomenon in VHE gamma-ray emitting blazars. The variability amplitude characterized by the fractional r.m.s. is strongly energy-dependent and is proportional to E^(0.19 +/- 0.01). The light curve r.m.s. correlates with the flux. This is the signature of a multiplicative process which can be accounted for as a red noise with a Fourier index of ~2. This unique data set shows evidence for a low level gamma-ray emission state from PKS 2155-304, which possibly has a different origin than the outbursts. The discovery of the light curve lognormal behaviour might be an indicator ..

AIMS: The properties of the broad-band emission from the high-frequency peaked BL Lac H 2356-309 (z=0.165) are investigated. METHODS: Very High Energy (VHE; E > 100 GeV) observations of H 2356-309 were performed with the High Energy Stereoscopic System (HESS) from 2004 through 2007. Simultaneous optical/UV and X-ray observations were made with the XMM-Newton satellite on June 12/13 and June 14/15, 2005. NRT radio observations were also contemporaneously performed in 2005. ATOM optical monitoring observations were also made in 2007. RESULTS: A strong VHE signal, ~13 sigma total, was detected by HESS after the four years HESS observations (116.8 hrs live time). The integral flux above 240 GeV is I(>240 GeV) = (3.06 +- 0.26 {stat} +- 0.61 {syst}) x 10^{-12} cm^{-2} s^{-1}, corresponding to ~1.6% of the flux observed from the Crab Nebula. A time-averaged energy spectrum is measured from 200 GeV to 2 TeV and is characterized by a power law (photon index of Gamma = 3.06 +- 0.15 {stat} +- 0.10 {syst}). Significant small-amplitude variations in the VHE flux from H 2356-309 are seen on time scales of months and years, but not on shorter time scales. No evidence for any variations in the VHE spectral slope are found within these data. The XMM-Newton X-ray measurements show a historically low X-ray state, characterized by a hard, broken-power-law spectrum on both nights. CONCLUSIONS: The broad-band spectral energy distribution (SED) of the blazar can be adequately fit using a simple one-zone synchrotron self-Compton (SSC) model. In the SSC scenario, higher VHE fluxes could be expected in the future since the observed X-ray flux is at a historically low level.

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

After a brief critical overview of the main properties of blazars and their classification, some significant results from recent multiwavelength observations are summarized, in the context of the jet physics.

Recent VHE observations have unveiled a surprising behaviour in two well-known blazars at opposite sides of the blazar sequence. PKS 2155-304 have shown for the first time in an HBL a large Compton dominance, high gamma-ray luminosities and a cubic relation between X-ray and VHE fluxes. 3C 279 is the first FSRQ detected at VHE. The high luminosity required to overcome the significant absorption caused by the BLR emission cannot be easily reconciled with the historical and quasi-simultaneous SED properties. Both cases shed a new light on the structure and ambient fields of blazars. Contrary to previous claims, it is also shown that 3C 279 --as any FSRQ-- cannot in general provide robust constraints on the EBL.

The BL Lac object 1ES 0229+200 (z=0.14) has been detected by HESS during observations taking place in 2005-2006. The TeV spectrum, when corrected for the absorption of gamma-ray photons through the interaction with the extragalactic background light, is extremely hard, even if the most conservative level for the background is considered. The case of 1ES 0229+200 is very similar to that of 1ES 1101-232, for which a possible explanation, in the framework of the standard one-zone synchrotron-self Compton model, is that the high-energy emission is synchrotron-self Compton radiation of electrons distributed as a power law with a large value of the minimum energy. In this scenario the hard TeV spectrum is accompanied by a very hard synchrotron continuum below the soft X-ray band. We will show that recent Swift observations of 1ES 0229+200 in the critical UV-X-ray band strongly support this model, showing the presence of the expected spectral break and hard continuum between the UV and the X-ray bands.

We use 3 years of data from the Swift/BAT survey to select a complete sample of X-ray blazars above 15 keV. This sample comprises 26 Flat-Spectrum Radio Quasars (FSRQs) and 12 BL Lac objects detected over a redshift range of 0.03<z<4.0. We use this sample to determine, for the first time in the 15--55 keV band, the evolution of blazars. We find that, contrary to the Seyfert-like AGNs detected by BAT, the population of blazars shows strong positive evolution. This evolution is comparable to the evolution of luminous optical QSOs and luminous X-ray selected AGNs. We also find evidence for an epoch-dependence of the evolution as determined previously for radio-quiet AGNs. We interpret both these findings as a strong link between accretion and jet activity. In our sample, the FSRQs evolve strongly, while our best-fit shows that BL Lacs might not evolve at all. The blazar population accounts for 10-20 % (depending on the evolution of the BL Lacs) of the Cosmic X-ray background (CXB) in the 15--55 keV band. We find that FSRQs can explain the entire CXB emission for energies above 500 keV solving the mystery of the generation of the MeV background. The evolution of luminous FSRQs shows a peak in redshift ($z_c$=4.3$\pm0.5$) which is larger than the one observed in QSOs and X-ray selected AGNs. We argue that FSRQs can be used as tracers of massive elliptical galaxies in the early Universe.

We report the discovery of faint very high energy (VHE, E > 100 GeV) gamma-ray emission from the radio galaxy Centaurus A in deep observations performed with the H.E.S.S. experiment. A signal with a statistical significance of 5.0 sigma is detected from the region including the radio core and the inner kpc jets. The integral flux above an energy threshold of ~250 GeV is measured to be 0.8% of the flux of the Crab Nebula and the spectrum can be described by a power law with a photon index of 2.7 +/- 0.5_stat +/- 0.2_sys. No significant flux variability is detected in the data set. The discovery of VHE gamma-ray emission from Centaurus A reveals particle acceleration in the source to >TeV energies and, together with M 87, establishes radio galaxies as a class of VHE emitters.

The time evolution of the electromagnetic emission from blazars, in particular high frequency peaked sources (HBLs), displays irregular activity not yet understood. In this work we report a methodology capable of characterizing the time behavior of these variable objects. The Maximum Likelihood Blocks (MLBs) is a model-independent estimator which sub-divides the light curve into time blocks, whose length and amplitude are compatible with states of constant emission rate of the observed source. The MLBs yields the statistical significance in the rate variations and strongly suppresses the noise fluctuations in the light curves. We apply the MLBs for the first time on the long term X-ray light curves (RXTE/ASM) of Mkn~421,Mkn~501, 1ES 1959+650 and 1ES 2155-304, which consist of more than 10 years of observational data (1996-2007). Using the MLBs interpretation of RXTE/ASM data, the integrated time flux distribution is determined for each single source considered. We identify in these distributions the characteristic level as well as the flaring states of the blazars. All the distributions show a significant component at negative flux values, most probably caused by an uncertainty in the background subtraction and by intrinsic fluctuations of RXTE/ASM. This effect interests in particular short time observations. In order to quantify the probability that the intrinsic fluctuations give rise to a false identification of a flare, we study a population of very faint sources and their integrated time flux distribution. We determine duty cycle or fraction of time a source spent in the flaring state of the source Mkn~421, Mkn~501, 1ES 1959+650 and 1ES 2155-304. Moreover, we study the random coincidences between flares and generic sporadic events such as high energy neutrinos or flares in other wavelengths.

TeV blazars are known as prominent nonthermal emitters across the entire electromagnetic spectrum with their photon power peaking in the X-ray and TeV bands. If distant, absorption of gamma-ray photons by the extragalactic background light (EBL) alters the intrinsic TeV spectral shape, thereby affecting the overall interpretation. Suzaku observations for two of the more distant TeV blazars known to date, 1ES 1101-232 and 1ES 1553+113, were carried out in 2006 May and July, respectively, including a quasi-simultaneous coverage with the state-of-the-art Cerenkov telescope facilities. We report on the resulting data sets with emphasis on the X-ray band and set in context to their historical behavior. During our campaign, we did not detect any significant X-ray or gamma-ray variability. 1ES 1101-232 was found in a quiescent state with the lowest X-ray flux ever measured. The combined XIS and HXD PIN data for 1ES 1101-232 and 1ES 1553+113 clearly indicate spectral curvature up to the highest hard X-ray data point (~30 keV), manifesting as softening with increasing energy. We describe this spectral shape by either a broken power law or a log-parabolic fit with equal statistical goodness of fits. The combined 1ES 1553+113 very high energy spectrum (90-500 GeV) did not show any significant changes with respect to earlier observations. The resulting contemporaneous broadband spectral energy distributions of both TeV blazars are discussed in view of implications for intrinsic blazar parameter values, taking into account the gamma-ray absorption in the EBL.

The energy spectra of TeV gamma-rays from blazars, after being corrected for intergalatic absorption in the Extragalactic Background Light (EBL), appear unusually hard, a fact that poses challenges to the conventional models of particle acceleration in TeV blazars and/or to the EBL models. In this paper we show that the internal absorption of gamma-rays caused by interactions with dense narrow-band radiation fields in the vicinity of compact gamma-ray production regions can lead to the formation of gamma-ray spectra of an almost arbitrary hardness. This allows significant relaxation of the current tight constraints on particle acceleration and radiation models, although at the expense of enhanced requirements to the available nonthermal energy budget. The latter, however, is not a critical issue, as long as it can be largely compensated by the Doppler boosting, assuming very large ($\geq 30$) Doppler factors of the relativistically moving gamma-ray production regions. The suggested scenario of formation of hard gamma-ray spectra predicts detectable synchrotron radiation of secondary electron-positron pairs which might require a revision of the current ``standard paradigm'' of spectral energy distributions of gamma-ray blazars. If the primary gamma-rays are of hadronic origin related to $pp$ or $p \gamma$ interactions, the ``internal gamma-ray absorption'' model predicts neutrino fluxes close to the detection threshold of the next generation high energy neutrino detectors.

Since 2002 the VHE (>100 GeV) gamma-ray flux of the high-frequency peaked BL Lac PKS 2155-304 has been monitored with the High Energy Stereoscopic System (HESS). An extreme gamma-ray outburst was detected in the early hours of July 28, 2006 (MJD 53944). The average flux above 200 GeV observed during this outburst is ~7 times the flux observed from the Crab Nebula above the same threshold. Peak fluxes are measured with one-minute time scale resolution at more than twice this average value. Variability is seen up to ~600 s in the Fourier power spectrum, and well-resolved bursts varying on time scales of ~200 seconds are observed. There are no strong indications for spectral variability within the data. Assuming the emission region has a size comparable to the Schwarzschild radius of a ~10^9 solar mass black hole, Doppler factors greater than 100 are required to accommodate the observed variability time scales.

We have observed 1ES 1426+428 with INTEGRAL detecting it up to $\sim$150 keV. The spectrum is hard, confirming that this source is an extreme BL Lac object, with a synchrotron component peaking, in a $\nu F_\nu$ plot, at or above 100 keV, resembling the hard states of Mkn 501 and 1ES 2344+514. All these three sources are TeV emitters, with 1ES 1426+428 lying at a larger redshift (z=0.129): for this source the absorption of high energy photons by the IR cosmic background is particularly relevant. The observed hard synchrotron tail helps the modeling of its spectral energy distribution, giving information on the expected intrinsic shape and flux in the TeV band. This in turn constrains the amount of the poorly known IR background.

The unexpectedly hard spectra measured by HESS for the BLLacs 1ES 1101-232 and H 2356-309 has allowed an upper limit on the Extragalactic Background Light (EBL) to be derived in the optical/near-infrared range, which is very close to the lower limit given by the resolved galaxy counts. This result seems to exclude a large contribution to the EBL from other sources (e.g. Population III stars) and indicates that the intergalactic space is more transparent to gamma-rays than previously thought. A brief discussion of EBL absorption effects on blazar spectra and further observational tests to check this conclusion are presented, including the selection of new candidates for observations with Cherenkov telescopes.

The understanding of acceleration mechanisms in active galactic nuclei (AGN) jets and the measurement of the extragalactic-background-light (EBL) density are closely linked and require the detection of a large sample of very-high-energy (VHE) emitting extragalactic objects at varying redshifts. We report here on the discovery with the H.E.S.S. (High Energy Stereoscopic System) atmospheric-Cherenkov telescopes of the VHE Gamma-ray emission from H2356 - 309, an extreme BL Lac object located at a redshift of 0.165. The observations of this object, which was previously proposed as a southern-hemisphere VHE candidate source, were performed between June and December 2004. The total exposure is 38.9 hours live time, after data quality selection, which yields the detection of a signal at the level of 9.0$\sigma$ (standard deviations) .

We present new BeppoSAX LECS and MECS observations, covering the energy range 0.1 - 10 keV (observer's frame), of four BL Lacertae objects selected from the 1 Jy sample. All sources display a flat (alpha_x ~ 0.7) X-ray spectrum, which we interpret as inverse Compton emission. One object shows evidence for a low-energy steepening (Delta alpha_x ~ 0.9) which is likely due to the synchrotron component merging into the inverse Compton one around ~ 2 keV. A variable synchrotron tail would explain why the ROSAT spectra of our sources are typically steeper than the BeppoSAX ones (Delta alpha_x} ~ 0.7). The broad-band spectral energy distributions fully confirm this picture and model fits using a synchrotron inverse Compton model allow us to derive the physical parameters (intrinsic power, magnetic field, etc.) of our sources. By combining the results of this paper with those previously obtained on other sources we present a detailed study of the BeppoSAX properties of a well-defined sub-sample of 14 X-ray bright (f_x (0.1 - 10 keV) > 3 x 10^{-12} erg/cm^2/s) 1-Jy BL Lacs. We find a very tight proportionality between nearly simultaneous radio and X-ray powers for the 1-Jy sources in which the X-ray band is dominated by inverse Compton emission, which points to a strong link between X-ray and radio emission components in these objects.

The intrinsic BL Lac spectra above few hundreds GeV can be very different from the observed ones due to the absorption effects by the diffuse Extragalactic Background Light (EBL), at present poorly known. With the recent results, there are now 4 sources with good spectral information: Mkn 421 (z=0.031), Mkn 501 (z=0.034), 1ES 1426+428 (z=0.129) and 1ES 1959+650 (z=0.047). Making simple assumptions on the shape of the intrinsic spectra (according to the present blazar knowledge), we have considered the effects of different EBL spectral energy distributions (SED) for the first time on all 4 objects together, deriving constraints for the EBL fluxes. These resulted significantly lower than many direct estimates.

Multi-wavelength observations with TeV photons are an essential diagnostic tool to study the physics of TeV sources. The complex variability of blazars, however (timescales from years down to minutes, with different patterns and SED behaviours), requires a great effort on simultaneous campaigns, which should be performed possibly over several days. Spectral information is essential, and now with the new TeV and X-ray telescopes it can be obtained on timescales less than one hour. The insights from such observations can be tremendous, since recent results have shown that the X-ray and TeV emissions do not always follow the same behaviour, and flares can have different relations between rise and decay times. Unfortunately, the strong pointing constraints of XMM do not allow the full use of this satellite simultaneously with ground telescopes.

The recent HEGRA detection and spectrum of 1ES 1426+428 at TeV energies, once corrected for absorption using present estimates of the diffuse extragalactic IR background, suggest that the high energy peak of the Spectral Energy Distribution (SED) could be much higher than the synchrotron one ($L_c/L_s>10$), and lie at energies above 8-10 TeV. To see if such an SED could be accounted for, we have applied a "finite injection time" SSC model, and present here some preliminary results. Within this model, we found the need of an external ("ambient") contribution to the energy density of seed photons, in order to account for both the high Compton dominance and the hard spectrum.

We present new BeppoSAX LECS, MECS, and PDS observations of four flat-spectrum radio quasars (FSRQ) having effective spectral indices alpha_ro and alpha_ox typical of high-energy peaked BL Lacs. Our sources have X-ray-to-radio flux ratios on average ~ 70 times larger than ``classical'' FSRQ and lie at the extreme end of the FSRQ X-ray-to-radio flux ratio distribution. The collected data cover the energy range 0.1 - 10 keV (observer's frame), reaching ~ 100 keV for one object. The BeppoSAX band in one of our sources, RGB J1629+4008, is dominated by synchrotron emission peaking at ~ 2 x 10^16 Hz, as also shown by its steep (energy index alpha_x ~ 1.5) spectrum. This makes this object the FIRST known FSRQ whose X-ray emission is not due to inverse Compton radiation. Two other sources display a flat BeppoSAX spectrum (alpha_x ~ 0.7), with weak indications of steepening at low X-ray energies. The combination of BeppoSAX and ROSAT observations, (non-simultaneous) multifrequency data, and a synchrotron inverse Compton model suggest synchrotron peak frequencies ~ 10^15 Hz, although a better coverage of their spectral energy distributions is needed to provide firmer values. If confirmed, these values would be typical of ``intermediate'' BL Lacs for which the synchrotron and inverse Compton components overlap in the BeppoSAX band. Our sources, although firmly in the radio-loud regime, have powers more typical of high-energy peaked BL Lacs than of FSRQ, and indeed their radio powers put them near the low-luminosity end of the FSRQ luminosity function. We discuss this in terms of an anti-correlation between synchrotron peak frequency and total power, based on physical arguments, and also as possibly due to a selection effect.

With its very wide energy band (0.1-100 keV), BeppoSAX has played a fundamental role in the blazar field, allowing in particular a better study of the highest synchrotron peaked objects. Here we summarize the results of an observational program performed with the aim to find and study more extreme BL Lac sources. We discuss the SEDs of the observed objects and their impact on the "blazar sequence" scenario, and consider their relevance as possible TeV emitting sources.

XMM observed the BL Lac PKS 2155-304 for a full orbit (about 150 ksec) on 2000 November 19-21. Preliminary results on the temporal and spectral analysis of data from the EPIC PN camera and Optical Monitor are presented. The variability amplitude depends systematically on energy, however the slopes of the structure functions of the light-curves in different bands do not appear to be significantly different. No evidence of time lags is found by cross correlating the light-curves in different bands.

The spectral properties of blazars seem to follow a phenomenological sequence according to the source luminosity. By inferring the source physical parameters through (necessarily) modeling the blazar spectra, we have previously proposed that the sequence arises because the particles responsible for most of the emission suffer increasing radiative losses as the luminosity increases. Here we extend those results by considering the widest possible range of blazar spectral properties. We find a new important ingredient for shaping the spectra of the lowest power objects, namely the role of a finite timescale for the injection of relativistic particles. Only high energy particles radiatively cool in such timescale leading to a break in the particle distribution: particles with this break energy are those emitting most of the power, and this gives raise to a link between blazar spectra and total energy density inside the source, which controls the cooling timescale. The emerging picture requires two phases for the particle acceleration: a first pre-heating phase in which particles reach a characteristic energy as the result of balancing heating and radiative cooling, and a more rapid acceleration phase which further accelerate these particles to form a power law distribution. While in agreement with standard shock theory, this scenario also agrees with the idea that the luminosity of blazars is produced through internal shocks, which naturally lead to shocks lasting for a finite time.

We report on two BeppoSAX observations of BL Lac (2200+420) performed respectively in June and December 1999, as part of a ToO program to monitor blazars in high states of activity. During both runs the source has been detected up to 100 keV, but it showed quite different spectra: in June it was concave with a very hard component above 5-6 keV (alpha_1 ~ 1.6; alpha_2 ~ 0.15); in December it was well fitted by a single power law (alpha ~ 0.6). During the first BeppoSAX observation BL Lac showed an astonishing variability episode: the 0.3 - 2 keV flux doubled in ~ 20 minutes, while the flux above 4 keV was almost constant. This frequency-dependent event is one of the shortest ever recordered for BL Lac objects and places lower limits on the dimension and the magnetic field of the emitting region and on the energy of the synchrotron radiating electrons. A similar but less extreme behaviour is detected also in optical light curves, that display non-simultaneous, smaller fluctuations of ~ 20% in 20 min. We fit the spectral energy distributions with a homogeneous, one-zone model to constrain the emission region in a very simple but effective SSC + external Compton scenario, highlightning the importance of the location of the emitting region with respect to the Broad Line Region and the relative spectral shape dependence. We compare our data with historical radio to gamma-ray Spectral Energy Distributions.

The TeV emission of low power BL Lac objects has been established by the detection of an handful of them. The knowledge of the level of the TeV emission and its spectrum can shed light on the particle acceleration mechanisms, and it is especially important to assess the still uncertain level of the far infrared background radiation, which can absorb the TeV photons through photon--photon interactions. In view of these implications, it is necessary to enlarge the number of TeV detected sources, and to find them at different redshifts. To this aim, we propose a general and simple criterium to select the best TeV candidates, and produce a list of them with flux estimates above 40 GeV, 300 GeV and 1 TeV.

We present BeppoSAX LECS, MECS, and PDS spectra of eleven X-ray selected BL Lacertae objects. Combining these sources with the ones presented elsewhere we have a sample of 21 BL Lacs from the Einstein Medium Sensitivity and Einstein Slew Survey. The sample shows strong correlations of several physical parameters with the peak frequency of the synchrotron branch of the spectral energy distribution. In particular the peak frequency is correlated to the X-ray spectral shape: objects with the peak near to the X-ray band show harder and straighter X-ray spectra than those of the low energy peaked sources. This work shows that the recently proposed unification scenario for different types of blazars can hold also within the class of high frequency peaked BL Lac objects.

We present new BeppoSAX observations of seven BL Lacertae objects selected from the 1 Jy sample plus one additional source. The collected data cover the energy range 0.1 - 10 keV (observer's frame), reaching ~ 50 keV for one source (BL Lac). All sources characterized by a peak in their multifrequency spectra at infrared/optical energies (i.e., of the LBL type) display a relatively flat (alpha_x ~ 0.9) X-ray spectrum, which we interpret as inverse Compton emission. Four objects (2/3 of the LBL) show some evidence for a low-energy steepening which is likely due to the synchrotron tail merging into the inverse Compton component around ~ 1 - 3 keV. If this were generally the case with LBL, it would explain why the 0.1 - 2.4 keV ROSAT spectra of our sources are systematically steeper than the BeppoSAX ones (by ~ 0.5 in alpha_x). The broad-band spectral energy distributions fully confirm this picture and a synchrotron inverse Compton model allows us to derive the physical parameters (intrinsic power, magnetic field, etc.) of our sources. Combining our results with those obtained by BeppoSAX on BL Lacs covering a wide range of synchrotron peak frequencies, we confirm and clarify the dependence of the X-ray spectral index on synchrotron peak frequency originally found in ROSAT data.

We performed an observational program with the X-ray satellite BeppoSAX to study objects with extreme synchrotron peak frequencies (nu_peak > 1 keV). Of the seven sources observed, four showed peak frequencies in the range 1-5 keV, while one (1ES 1426+428) displayed a flat power law spectrum (alpha= 0.92), locating its synchrotron peak at or above 100 keV. This is the third source of this type ever found, after Mkn 501 and 1ES 2344+514. Our data confirm the large nu_peak variability of this class of sources, compared with lower peaked objects. The high synchrotron peak energies, flagging the presence of high energy electrons, make the extreme BL Lacs also good candidates for TeV emission, and therefore good probes for the IR background.

We present the results of two ToO BeppoSAX observations of BL Lac. During the first observation we detected both the synchrotron and the Compton components. Fast time variability was present, but only for the synchrotron component. During the second observation the spectrum was flatter and only the Compton component was present. Four different SED are presented and can be described by a blob moving along a jet and responsible for the SSC emission. The seed photons for the Compton scattering are the synchrotron photons themself plus possibly external photons coming from the broad line region.

Blazars can be classified on the basis of their overall Spectral Energy Distribution (SED). BL Lac objects are usually divided in LBL or HBL (Low or High energy peaked BL Lacs), according to the peak frequency of the synchrotron emission, if in the optical or UV-soft-X band respectively. FSRQs instead are characterized by synchrotron peaks mainly at IR-optical frequencies, similarly to LBLs. Here we report on recent BeppoSAX observations which are unveiling the high energy branch of the range of synchrotron peak frequencies. Four new "extreme" HBLs have been discovered, one of which (1ES 1426+428) peaks near or above 100 keV, in a quiescent state.

We report on the BeppoSAX observations of four "extreme" BL Lacs, selected to have high synchrotron peak frequencies. All have been detected also in the PDS band. For 1ES 0120+340, PKS 0548-322 and H 2356-309 the spectrum is well fitted by a convex broken power-law, thus locating the synchrotron peak around 1 - 4 keV. 1ES 1426+428 presents a flat energy spectral index (alpha_x= 0.92) up to 100 keV, thus constraining the synchrotron peak to lie near or above that value. For their extreme properties, all sources could be strong TeV emitters.

BeppoSAX observed 8 times the gamma-ray bright blazar PKS 0528+134 in Feb. and Mar. 1997, during a multiwavelength campaign involving EGRET, onboard CGRO, and ground based telescopes. The source was in its faintest X-ray state observed so far, with an unabsorbed [2-10] keV flux of 2.5E-12 erg/cm2/s. The spectrum can be fitted with a power law, with energy index alpha_X=0.49+-0.07 between 0.1 and 10 keV. The high energy instrument PDS detected a flux in the 15-100 keV band which appears disconnected from the lower X-ray energies, most likely due to a contaminating source. We consider our findings in the context of the overall spectral energy distribution and discuss their implications for synchrotron and inverse Compton models.

We report on the BeppoSAX observations of the gamma-ray blazar PKS 0528+134 performed in Feb and Mar 1997, during a multiwavelength campaign involving EGRET and ground based telescopes. The source was in a faint and hard state, with energy spectral index alpha=0.48+-0.04 between 0.1 and 10 keV, and [2-10] keV flux of 2.7E-12 erg/cm2/s. No significant variability was observed. The source was detected in the 20-120 keV band by the PDS, with a flux lying slightly above the extrapolation from lower X-ray energies. Comparing this low state with previous higher states of the source, there is an indication that the X-ray spectrum hardens and the gamma-ray spectrum steepens when the source is fainter.