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

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

We present multi-epoch optical observations of the blazar 3C 454.3 (z = 0.859) from 2008 August through 2011 December, using the SMARTS Consortium 1.5m+RCSpectrograph and 1.3m+ANDICAM in Cerro Tololo, Chile. The spectra reveal that the broad optical emission lines Mg II, H-beta and H-gamma are far less variable than the optical or gamma-ray continuum. Although, the gamma-rays varied by a factor of 100 above the EGRET era flux, the lines generally vary by a factor of 2 or less. Smaller variations in the gamma-ray flux did not produce significant variation in any of the observed emission lines. Therefore, to first order, the ionizing flux from the disk changes only slowly during large variations of the jet. However, two exceptions in the response of the broad emission lines are reported during the largest gamma-ray flares in 2009 December and 2010 November, when significant deviations from the mean line flux in H-gamma and Mg II were observed. H-gamma showed a maximum 3-sigma and 4-sigma deviation in each flare, respectively, corresponding to a factor of 1.7 and 2.5 increase in flux. Mg II showed a 2-sigma deviation in both flares; no variation was detected in H-beta during either flare. These significant deviations from the mean line flux also coincide with 7mm core ejections reported previously (Jorstad et al. 2012). The correlation of the increased emission line flux with mm core ejections, and gamma-ray, optical and UV flares suggests that the broad line region extends beyond the gamma-emitting region during the 2009 and 2010 flares.

Recent multiwavelength observations of PKS 0208-512 by SMARTS, Fermi, and Swift revealed that gamma-ray and optical light curves of this flat spectrum radio quasars are highly correlated, but with an exception of one large optical flare having no corresponding gamma-ray activity or even detection. On the other hand, recent advances in SNRs observations and plasma simulations both reveal that magnetic field downstream of astrophysical shocks can be largely amplified beyond simple shock compression. These amplifications, along with their associated particle acceleration, might contribute to blazar flares, including the peculiar flare of PKS 0208-512. Using our time dependent multizone blazar emission code, we evaluate several scenarios that may represent such phenomena. This code combines Monte Carlo method that tracks the radiative processes including inverse Compton scattering, and Fokker-Planck equation that follows the cooling and acceleration of particles. It is a comprehensive time dependent code that fully takes into account the light travel time effects. In this study, both the changes of the magnetic field and acceleration efficiency are explored as the cause of blazar flares. Under these assumption, synchrotron self-Compton and external Compton scenarios produce distinct features that favor the external Compton scenario. The optical flares with/without gamma-ray counterparts can be explained by different allocations of energy between the magnetization and particle acceleration, which in turn can be affected by the relative orientation between the magnetic field and the shock flow. We compare the details of the observations and simulation, and highlight what implications this study has on our understanding of relativistic jets.

The Yale/SMARTS optical-near-IR monitoring program has followed the variations in emission of the Fermi-LAT monitored blazars in the southern sky with closely spaced observations since 2008. We report the discovery of an optical-near-IR (OIR) outburst with no accompanying gamma-rays in the blazar PKS 0208-512, one of the targets of this program. While the source undergoes three outbursts of 1 mag or more at OIR wavelengths lasting for longer than 3 months during 2008-2011, only interval 1 and 3 have corresponding bright phases in GeV energies lasting longer than 1 month. The OIR outburst during interval 2 is comparable in brightness and temporal extent to the OIR flares during intervals 1 and 3 which do have gamma-ray counterparts. Gamma-ray and OIR variability are very well-correlated in most cases in the Fermi blazars and the lack of correlation in this case is anomalous. By analyzing the gamma-ray, OIR, and supporting multi-wavelength variability data in details, we speculate that the location of the outburst in the jet during interval 2 was closer to the black hole where the jet is more compact and the magnetic field strength is higher, and the bulk Lorentz factor of the material in the jet is smaller. These result in a much lower Compton dominance and no observable gamma-ray outburst during interval 2.

We report the discovery of an anomalous flare in a bright blazar, namely, PKS 0208-512, one of the targets of the Yale/SMARTS optical-near-infrared (OIR) monitoring program of Fermi blazars. We identify three intervals during which PKS 0208-512 undergoes outbursts at OIR wavelengths lasting for longer than 3 months. Its brightness increases and then decreases again by at least 1 magnitude in these intervals. In contrast, the source undergoes bright phases in GeV energies lasting for longer than 1 month during intervals 1 and 3 only. The OIR outburst during interval 2 is comparable in brightness and temporal extent to the OIR flares during intervals 1 and 3 which do have gamma-ray counterparts. By analyzing the gamma-ray, OIR, and supporting multi-wavelength variability data in details, we speculate that the OIR outburst during interval 2 was caused by a change in the magnetic field without any change in the total number of emitting electrons or Doppler factor of the emitting region. Alternatively, it is possible that the location of the outburst in the jet during interval 2 was closer to the black hole where the jet is more compact and the bulk Lorentz factor of the material in the jet is smaller. We also discuss the complex OIR spectral behavior during these three intervals.

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

We present the first collective evidence that Fermi-detected jets of high kinetic power (L_kin) are dominated by inverse Compton emission from upscattered external photons. Using a sample with a broad range in orientation angle, including radio galaxies and blazars, we find that very high power sources (L_kin > 10^45.5 erg s^{-1}) show a significant increase in the ratio of inverse Compton to synchrotron power (Compton dominance) with decreasing orientation angle, as measured by the radio core dominance and confirmed by the distribution of superluminal speeds. This increase is consistent with beaming expectations for external Compton (EC) emission, but not for synchrotron self Compton (SSC) emission. For the lowest power jets (L_kin < 10^43.5 erg s^{-1}), no trend between Compton and radio core dominance is found, consistent with SSC. Importantly, the EC trend is not seen for moderately high power flat spectrum radio quasars with strong external photon fields. Coupled with the evidence that jet power is linked to the jet speed (Kharb et al. 2010), this finding suggests that external photon fields become the dominant source of seed photons in the jet comoving frame only for the faster and therefore more powerful jets.

In recent work, we have identified two sub-populations of radio-loud AGN which appear to be distinguished by jet structure, where low-efficiency accreting systems produce `weak' jets which decelerate more rapidly than the `strong' jets of black holes accreting near the Eddington limit. The two classes are comprised of: (1) The weak jet sources, corresponding to FR I radio galaxies, having a decelerating or spine-sheath jet with velocity gradients, and (2) The strong jet sources, having fast, collimated jets, and typically displaying strong emission lines. The dichotomy in the \nu_peak-L_peak plane can be understood as a `broken power sequence' in which jets exist on one branch or the other based on the particular accretion mode. We suggest that the intrinsic kinetic power (as measured by low-frequency, isotropic radio emission), the orientation, and the accretion rate of the SMBH system are the the fundamental axes needed for unification of radio-loud AGN by studying a well-characterized sample of several hundred Fermi-detected jets. Finally, we present very recent findings that the most powerful strong jets produce gamma-rays by external Compton rather than SSC emission, placing the dissipation region in these strong jets at a radius inside the BLR and/or molecular torus.

We present the time variability properties of a sample of six blazars, AO 0235+164, 3C 273, 3C 279, PKS 1510-089, PKS 2155-304, and 3C 454.3, at optical-IR as well as gamma-ray energies. These observations were carried out as a part of the Yale/SMARTS program during 2008-2010 that has followed the variations in emission of the bright Fermi-LAT-monitored blazars in the southern sky with closely-spaced observations at BVRJK bands. We find the optical/IR time variability properties of these blazars to be remarkably similar to those at the gamma-ray energies. The power spectral density (PSD) functions of the R-band variability of all six blazars are fit well by simple power-law functions with negative slope such that there is higher amplitude variability on longer timescales. No clear break is identified in the PSD of any of the sources. The average slope of the PSD of R-band variability of these blazars is similar to what was found by the Fermi team for the gamma-ray variability of a larger sample of bright blazars. This is consistent with leptonic models where the optical-IR and gamma-ray emission is generated by the same population of electrons through synchrotron and inverse-Compton processes, respectively. The prominent flares present in the optical-IR as well as the gamma-ray light curves of these blazars are predominantly symmetric, i.e., have similar rise and decay timescales, indicating that the long-term variability is dominated by the crossing time of radiation or a disturbance through the emission region rather than by the acceleration or energy-loss timescales of the radiating electrons. In the blazar 3C 454.3, which has the highest-quality light curves, the location of a large gamma-ray outburst during 2009 December is consistent with being in the jet at ~18 pc from the central engine. This poses strong constraints on the models of high energy emission in the jets of blazars.

The location of the Gamma-ray emission of powerful blazars is a matter of active debate. Is the location within the UV emitting sub-pc scale broad line region, or farther out at pc scales where the molecular torus IR emission dominates? We present a diagnostic that connects three observables, the synchrotron and external Compton peak frequencies and the Compton dominance (the ratio of Compton to synchrotron luminosity) to the seed photon energy and energy density. We discuss encouraging preliminary results and discuss how to use our diagnostic to understand the location of the Gamma-ray emission as a function of source power through the use of multiwavelength observations.

In examining a select sample of over 200 blazars of known jet kinetic power (L_kin) and well-characterized SEDs, we found (Meyer et al., 2011) that Intermediate synchrotron-peaking (ISP) blazars may have lower gamma-ray output than high synchrotron-peaking (HSP) blazars of similar L_kin, consistent with our hypothesis that ISP blazars are less-beamed versions of HSP blazars, rather than a distinct population. Further, by using the radio core dominance as a measure of relative beaming, we find that gamma-ray luminosity depends on beaming in a consistent way for blazars ranging over all jet kinetic powers (10^42 - 10^46 ergs/s). We re-examine the gamma-ray properties of this core sample of blazars using the 1-year LAT catalog (Abdo et al., 2010). We find that for weak jets, the ratio of inverse Compton to synchrotron emission remains constant with increased beaming, consistent with an SSC model for the jet emission, while the most powerful jets show a strong increase in Compton dominance with orientation, consistent with an external Compton (EC) emission model.

We present multiwavelength data for twelve blazars observed from 2008-2010 as part of an ongoing optical-infrared photometric monitoring project. Sources were selected to be bright, southern (dec < 20 deg) blazars observed by the Fermi Gamma-Ray Space Telescope, with daily and weekly gamma-ray fluxes made available from the start of the Fermi mission. Light curves are presented for the twelve blazars in BVRJK at near-daily cadence. We find that optical and infrared fluxes are well correlated in all sources. Gamma-ray bright flat spectrum radio quasars (FSRQs) in our sample have optical/infrared emission correlated with gamma-rays consistent with inverse Compton-scattering models for GeV emission. In FSRQs, the variability amplitude decreases towards optical/IR wavelengths, consistent with the presence of a thermal emission component from the accretion disk varying on significantly longer timescales than the jet synchrotron emission. In BL Lac objects, variability is mainly constant across wavelengths, consistent with a weak or radiatively inefficient disk. FSRQs have redder optical-infrared colors when they are brighter, while BL Lac objects show no such trend. Several objects show complicated color-magnitude behavior: AO 0235+164 appears in two different states depending on whether it is gamma-ray bright or not. OJ 287 and 3C 279 show some hysteresis tracks in their color-magnitude diagrams. Individual flares may be achromatic or otherwise depart from the trend, suggesting different jet components becoming important at different times. We present a time-dependent spectral energy distribution of the bright FSRQ 3C 454.3 during its December 2009 flare, which is well fit by an external Compton model in the bright state, although day to day changes pose challenges to a simple one-zone model. All data from the SMARTS monitoring program are publicly available on our website.

We recently argued (Meyer 2011) that the collective properties of radio loud active galactic nuclei point to the existence of two families of sources, one of powerful sources with single velocity jets and one of weaker jets with significant velocity gradients in the radiating plasma. These families also correspond to different accretion modes and therefore different thermal and emission line intrinsic properties: powerful sources have radiatively efficient accretion disks, while in weak sources accretion must be radiatively inefficient. Here, after we briefly review of our recent work, we present the following findings that support our unification scheme: (i) along the broken sequence of aligned objects, the jet kinetic power increases. (ii) in the powerful branch of the sequence of aligned objects the fraction of BLLs decreases with increasing jet power. (iii) for powerful sources, the fraction of BLLs increases for more un-aligned objects, as measured by the core to extended radio emission. Our results are also compatible with the possibility that a given accretion power produces jets of comparable kinetic power.

We revisit the concept of a blazar sequence that relates the synchrotron peak frequency ({\nu}peak) in blazars with synchrotron peak luminosity (Lpeak, in {\nu}L{\nu}) using a large sample of radio-loud AGN. We present observational evidence that the blazar sequence is formed from two populations in the synchrotron {\nu}peak - Lpeak plane, each forming an upper edge to an envelope of progressively misaligned blazars, and connecting to an adjacent group of radio galaxies having jets viewed at much larger angles to the line of sight. When binned by jet kinetic power (Lkin; as measured through a scaling relationship with extended radio power), we find that radio core dominance decreases with decreasing synchrotron Lpeak, revealing that sources in the envelope are generally more misaligned. We find population-based evidence of velocity gradients in jets at low kinetic powers (~ 10^42-10^44.5 erg/s), corresponding to FR I radio galaxies and most BL Lacs. These low jet power 'weak jet' sources, thought to exhibit radiatively inefficient accretion, are distinguished from the population of non-decelerating, low synchrotron-peaking (LSP) blazars and FR II radio galaxies ('strong' jets) which are thought to exhibit radiatively efficient accretion. The two-population interpretation explains the apparent contradiction of the existence of highly core-dominated, low-power blazars at both low and high synchrotron peak frequencies, and further implies that most intermediate synchrotron peak (ISP) sources are not intermediate in intrinsic jet power between LSP and high synchrotron-peaking (HSP) sources, but are more misaligned versions of HSP sources with similar jet powers.

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

We use Monte Carlo/Fokker-Planck simulations to study the X-ray time lags. Our results show that soft lags will be observed as long as the decay of the flare is dominated by radiative cooling, even when acceleration and cooling timescales are similar. Hard lags can be produced in presence of a competitive achromatic particle energy loss mechanism if the acceleration process operates on a timescale such that particles are slowly moved towards higher energy while the flare evolves. In this type of scenario, the {\gamma} -ray/X-ray quadratic relation is also reproduced.

The blazar 3C 454.3 was revealed by the Fermi Gamma-ray Space Telescope to be in an exceptionally high flux state in July 2008. Accordingly, we performed a multi-wavelength monitoring campaign on this blazar using IR and optical observations from the SMARTS telescopes, optical, UV and X-ray data from the Swift satellite, and public-release gamma-ray data from Fermi. We find an excellent correlation between the IR, optical, UV and gamma-ray light curves, with a time lag of less than one day. The amplitude of the infrared variability is comparable to that in gamma-rays, and larger than at optical or UV wavelengths. The X-ray flux is not strongly correlated with either the gamma-rays or longer wavelength data. These variability characteristics find a natural explanation in the external Compton model, in which electrons with Lorentz factor gamma~10^(3-4) radiate synchrotron emission in the infrared-optical and also scatter accretion disk or emission line photons to gamma-ray energies, while much cooler electrons (gamma~10^(1-2)) produce X-rays by scattering synchrotron or other ambient photons.

(Abridged) We present a detailed analysis of week-long simultaneous observations of the blazar Mrk421 at 2-60 keV X-rays (RXTE) and TeV gamma-rays (Whipple and HEGRA) in 2001. The unprecedented quality of this dataset enables us to establish firmly the existence of the correlation between the TeV and X-ray luminosities, and to start unveiling some of its more detailed characteristics, in particular its energy dependence, and time variability. The source shows strong, highly correlated variations in X-ray and gamma-ray. No evidence of X-ray/gamma-ray interband lag is found on the full week dataset (<3 ks). However, a detailed analysis of the March 19 flare reveals that data are not consistent with the peak of the outburst in the 2-4 keV X-ray and TeV band being simultaneous. We estimate a 2.1+/-0.7 ks TeV lag. The amplitudes of the X-ray and gamma-ray variations are also highly correlated, and the TeV luminosity increases more than linearly w.r.t. the X-ray one. The strong correlation supports the standard model in which a unique electrons population produces the X-rays by synchrotron radiation and the gamma-ray component by inverse Compton scattering. However, for the individual best observed flares the gamma-ray flux scales approximately quadratically w.r.t. the X-ray flux, posing a serious challenge to emission models for TeV blazars. Rather special conditions and/or fine tuning of the temporal evolution of the physical parameters of the emission region are required in order to reproduce the quadratic correlation.

We discuss the correlation between the evolution of the TeV emission and X-ray radiation observed in high-energy peaked BL Lac objects. We describe such a correlation by a simple power law. In the first part of this work we present correlations obtained for the activity of Mrk~501 observed in 1997 April and for the activity of Mrk~421 observed in 2000 February. We also discuss results of similar correlation obtained for this source by other authors. They report an almost quadratic correlations observed between the evolution of the TeV and X-ray emission. In the second part of this paper we present a phenomenological model which describes the evolution of the synchrotron and inverse Compton emission of a simple spherical homogeneous source. Neglecting the radiative cooling of the particles we derive analytical expressions that describe the evolution. Then we use a numerical code to investigate the impact of radiative cooling on the evolution. The quadratic correlation observed during the decay phase of the flare observed in Mrk~421 on 2001 March 19 appears problematic for this basic modeling. The quadratic correlation can be explained only for specific choices of the spectral bands used for the calculation. Therefore, looking for more robust solutions, we investigate the evolution of the emission generated by a cylindrical source. However this model does not provide robust solutions for the problem of a quadratic correlation. factor. Finally we briefly discuss the possible influence of the light travel time effect on our results.

Following the detection of strong TeV gamma-ray flares from the BL Lac object 1ES 1959+650 with the Whipple 10 m Cherenkov telescope on May 16 and 17, 2002, we performed intensive Target of Opportunity (ToO) radio, optical, X-ray and TeV gamma-ray observations from May 18, 2002 to August 14, 2002. Observations with the X-ray telescope RXTE and the Whipple and HEGRA gamma-ray telescopes revealed several strong flares, enabling us to sensitively test the X-ray/gamma-ray flux correlation properties. Although the X-ray and gamma-ray fluxes seemed to be correlated in general, we found an ``orphan'' gamma-ray flare that was not accompanied by an X-ray flare. After describing in detail the radio (UMRAO, VLA), optical (Boltwood, Abastumani), X-ray (RXTE) and gamma-ray (Whipple, HEGRA) light curves and Spectral Energy Distributions (SEDs) we present initial modeling of the SED with a simple Synchrotron Self-Compton (SSC) model. With the addition of another TeV blazar with good broadband data, we consider the set of all TeV blazars to begin to look for a connection of the jet properties to the properties of the central accreting black hole thought to drive the jet. Remarkably, the temporal and spectral X-ray and gamma-ray emission characteristics of TeV blazars are very similar, even though the masses estimates of their central black holes differ by up to one order of magnitude.

We report the BeppoSAX observations of 6 Flat Spectrum Radio Quasars. Three of them have a clear detection up to 100 keV with the PDS instrument. For 4 objects the X-ray spectrum is satisfactorily fitted by a power-law continuum with Galactic absorption. 2251+158 show the presence of absorption higher than the galactic value, while the spectrum of the source 0208-512 shows a complex structure, with evidence of absorption at low energy. We construct the Spectral Energy Distributions adding historical data to the broad band X-ray spectra obtained with $Beppo$SAX and reproduce them with a one-zone Synchrotron-Inverse Compton model (including both SSC and External Compton). The implications are briefly discussed.

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.

PKS 2155--304 is one of the brightest blazars in the X-ray band. It was repeatedly monitored with BeppoSax during three long campaigns of about 2 days each in November of 1996, 1997 and 1999. The source underwent different states of intensity and was clearly variable with successive flares detected. This paper presents temporal and spectral analysis to study the X-ray variability trends for a blazar. The rich phenomenology is interpreted in the context of a model where relativistic electrons are accelerated through internal shocks taking place in the jets. The most important parameter turns out to be the initial time interval between the two shells ejected from the central engine to produce the flare, which may determine the structure of the shock and in turn the physical quantities of the emitting region to reproduce the observed trends of the X-ray variability (abridged).

We present the BeppoSAX observations of 3C 279 performed in 1997 January 13-23 simultaneously with gamma-ray pointings with EGRET on board CGRO and optical observations at the Torino Observatory (1997 January 11-18). ISO data close to this epoch are also available. We compare the derived Spectral Energy Distribution with those obtained at all other epochs with adequate multiwavelength coverage. Simple spectral models fitted to the multiepoch SEDs suggest that changes in intensity can be ascribed mainly to the variation of the bulk Lorentz factor of the plasma in the jet.

We have considered all blazars observed in the X-ray band and for which the slope of the X-ray spectrum is available. We have collected 421 spectra of 268 blazars, including 12 archival unpublished ASCA spectra of 7 blazars whose analysisis presented here. The X-ray spectra of blazars show trends as a function of their power, confirming that the blazar overall energy distribution can be parameterized on the basis of one parameter only, i.e. the bolometric luminosity. This is confirmed by the relatively new hard (2-10 keV) X-ray data. Our results confirm the idea that in low power objects the X-ray emission mechanism is the synchrotron process, dominating both the soft and the hard X-ray emissions. Low energy peaked BL Lac objects are intermediate, often showing harder spectra in the hard X-ray band, suggesting that the synchrotron process dominates in the soft band, with the inverse Compton process dominating at high energies. The most powerful objects have X-ray spectra that are flat both in the soft and in the hard band, consistent with a dominating inverse Compton component.

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 discuss the preliminary results of an extensive effort to address the fundamental, and yet un-answered, question that can be trivialized as: ``are there more blue or red blazars?''. This problematic is tightly connected with the much debated issue of the unified picture(s) of radio-loud AGNs, which in turn revolves around the existence, and the properties of relativistic jets. We address this question by comparing -- simultaneously -- the properties of the collection of heterogeneously selected samples that are available now, with the predictions of a set of plausible unifications scenarios. We show that it is already possible to make significant progress even by using only the present samples. The important role of selection effects is discussed. For instance we show that the multiple flux selections typical of available surveys could induce some of the correlations found in color-color diagrams. These latter results should apply to any study of flux limited samples.

We conducted a multi-frequency campaign for the TeV blazar Mrk~421 in 1998 April. The campaign started from a pronounced high amplitude flare recorded by SAX and Whipple; ASCA observation started three days later. In the X-ray data, we detected multiple flares, occuring on time scales of about one day. ASCA data clearly reveal spectral variability. The comparison of the data from ASCA, EUVE and RXTE indicates that the variability amplitudes in the low energy synchrotron component are larger at higher photon energies. In TeV Gamma-rays, large intra-day variations -- which were correlated with the X-ray flux -- were observed when results from three Cherenkov telescopes are combined. The RMS variability of TeV Gamma--rays was similar to that observed in hard X-rays, above 10 keV. The X-ray light curve reveals flares which are almost symmetric for most of cases, implying the dominant time scale is the light crossing time through the emitting region. The structure function analysis based on the continuous X-ray light curve of seven days indicates that the characteristic time scale is ~0.5 day. The analysis of ASCA light curves in various energy bands appears to show both soft (positive) and hard (negative) lags. These may not be real, as systematic effects could also produce these lags, which are all much smaller than an orbit. If the lags of both signs are real, these imply that the particle acceleration and X-ray cooling time scales are similar.

We present the preliminary results of a long BeppoSAX observation of the BL Lac object Mkn501 carried out in June 1999. The source was fainter than found during the BeppoSAX pointings of 1997 and 1998, but is still detected with a good signal-to-noise ratio up to ~40 keV. The X-ray spectrum in the energy range 0.1-40 keV, produced through synchrotron radiation, is steeper than in the previous years, it is clearly curved, and peaks (in nu*F_nu) at ~0.5 keV. This energy is much lower than those at which the synchrotron component was found to peak in 1997 and 1998. Some intraday variability suggests that activity of the source on small time scales accompanies the large long time scale changes of brightness and spectrum.

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 present BeppoSAX observations of three gamma-ray emitting quasars, namely 0836+710, 1510-089 and 2230+114. The three objects have been detected up to 100 keV showing extremely flat power-law spectra above 2 keV (energy index alpha_{2-10}=0.3-0.5). The soft X-ray spectrum of 0836+710 implies either an absorption column density higher than the galactic one or an intrinsically very hard slope (alpha_{0.1-1}=-0.2) below 1 keV. 1510-089 shows a soft excess, with the low energy spectrum steeper (alpha_{0.1-1}=1.6) than the high energy power-law. The results are discussed in the framework of current Inverse Compton models for the high energy emission of Flat Spectrum Radio Quasars and are used to estimate the physical quantities in the jet emitting region and to shed light on the energy transport mechanisms in jets. Finally we discuss the estimates of the jet luminosity in the context of the Blandford & Znajek mechanism for jet production.

Mkn 421 was repeatedly observed with BeppoSAX in 1997-1998. This is the first of two papers where we present the results of a thorough temporal and spectral analysis of all the data available to us, focusing in particular on the flare of April 1998, which was simultaneously observed also at TeV energies. Here we focus on the time analysis, while the spectral analysis and physical interpretation are presented in the companion paper. The detailed study of the flare in different energy bands reveals very important new results: i) hard photons lag the soft ones by 2-3 ks -a behavior opposite to what is normally found in high energy peak BL Lacs X-ray spectra; ii) the flare light curve is symmetric in the softest X-ray band, while it becomes increasingly asymmetric at higher energies, with the decay being progressively slower than the rise; iii) the flux decay of the flare can be intrinsically achromatic if a stationary underlying emission component is present. The temporal and spectral information obtained challenge the simplest models currently adopted for the (synchrotron) emission and most importantly provide clues on the particle acceleration process.

Mkn 421 was repeatedly observed with BeppoSAX in 1997-1998. The source showed a very rich phenomenology, with remarkable spectral variability. This is the second of two papers presenting the results of a thorough temporal and spectral analysis of all the data available to us, focusing in particular on the flare of April 1998, which was simultaneously observed also at TeV energies. The spectral analysis and correlations are presented in this paper, while the data reduction and timing analysis are the content of the companion paper. The spectral evolution during the flare has been followed over few ks intervals, allowing us to detect for the first time the peak of the synchrotron component shifting to higher energies during the rising phase, and then receding. This spectral analysis nicely confirms the delay of the flare at the higher energies, which in Paper I we quantified as a hard lag of a few ks. Furthermore, at the highest energies, evidence is found of variations of the inverse Compton component. The spectral and temporal information obtained challenge the simplest models currently adopted for the (synchrotron) emission and most importantly provide clues on the particle acceleration process. A scenario accounting for all the observational constraints is discussed, where electrons are injected at progressively higher energies during the development of the flare, and the achromatic decay is ascribed to the source light crossing time exceeding the particle cooling timescales.

We present $Beppo$SAX observations of the $\gamma $-ray emitting quasars 0836+710, 1510-089 and 2230+114. All the objects have been detected in the PDS up to 100 keV and have extremely flat power-law spectra above 2 keV ($\alpha _x$=0.3--0.5). 0836+710 shows absorption higher than the galactic value and marginal evidence for the presence of the redshifted 6.4 keV Iron line. 1510-089 shows a spectral break around 1 keV, with the low energy spectrum steeper ($\alpha_l$=1.6) than the high energy power-law ($\alpha_h$=0.3). The data are discussed in the light of current Inverse Compton models for the high energy emission.

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.

ON 231 was observed with BeppoSAX in May and June 1998, following an exceptional optical outburst which occured in April--May. We measured the X-ray spectrum from 0.1 up to 100 keV. In both occasions the spectrum had a concave shape, with a break detected at about 4 and 2.5 keV, respectively. We interpret the steeper component at energies below the break as due to synchrotron emission and the extremely flat component at energies above the break as due to inverse Compton emission. This is so far the best example in which both the synchrotron and the Inverse Compton component are detected simultaneously and with the same instruments in the X-ray spectrum of a blazar. We detect a decrease of the X-ray flux of about 30% below the break between the first and the second observation, and smaller variability amplitude between 4 and 10 keV. During the May observation we also detected a fast variability event with the flux below 4 keV increasing by about a factor of three in 5 hours. Above 4 keV no variability was detected. We discuss these results in the framework of synchrotron self-Compton models.

Mkn 421 was repeatedly observed with BeppoSAX in 1997-1998. We present highlights of the results of the thorough temporal and spectral analysis discussed by Fossati et al. (1999) and Maraschi et al. (1999), focusing on the flare of April 1998, which was simultaneously observed also at TeV energies. The detailed study of the flare in different energy bands reveals a few very important new results: (a) hard photons lag the soft ones by 2-3 ks *a behavior opposite to what is normally found in High energy peak BL Lacs X-ray spectra*; (b) the flux decay of the flare can be intrinsically achromatic if a stationary underlying emission component is present. Moreover the spectral evolution during the flare has been followed by extracting X-ray spectra on few ks intervals, allowing to detect for the first time the peak of the synchrotron component shifting to higher energies during the rising phase, and then receding. The spectral analysis confirms the delay in the flare at the higher energies, as above a few keV the spectrum changes only after the peak of the outburst has occurred. The spectral and temporal information obtained challenge the simplest models currently adopted for the (synchrotron) emission and most importantly provide clues on the particle acceleration process. A theoretical picture accounting for all the observational constraints is discussed, where electrons are injected at low energies and then progressively accelerated during the development of the flare.

We present a detailed Cross Correlation analysis of the light curves of the X-ray bright BL Lac object Mkn 421 observed with $Beppo$SAX period 1997-1998. Monte Carlo simulations are performed in order to estimate the uncertainties. We found that in one occasion the lag is negative.

Markarian 421 was observed for about four days with BeppoSAX and the Whipple Observatory gamma-ray telescope in April 1998. A pronounced, well-defined, flare with hour-scale variability was observed simultaneously in X-rays and very high energy gamma-rays. These data provide the first evidence that the X-ray and TeV intensities are well correlated on time-scales of hours. While the rise of the flare occurred on a similar time-scale in the two wavebands, the decay of the flare was much more rapid in gamma rays, providing the first clear indication that the X-ray and gamma-ray emission may not be completely correlated in Markarian 421.

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 present the results of BeppoSAX observations of PKS 2155-304 during an intense gamma-ray flare. The source was in a high X-ray state. A temporal analysis of the data reveals a tendency of the amplitude of variations to increase with energy, and the presence of a soft lag with a timescale of the order 10^3 s. A curved continuum spectrum, with no evidence of spectral features, extends up to ~50 keV, while there is indication of a flatter component emerging at higher energies, consistent with the interpretation of the broad band spectral energy distribution (SED) as due to synchrotron self-Compton (SSC) emission from a single region. Notably, the fitting of the SED with such a model is consistent with an interpretation of the detected soft lag as due to radiative cooling, supporting the idea that radiation losses play an important role in variability. The observed shifts of the SED peaks between the lowest and highest flux levels can be accounted for by an increase of the break energy in the relativistic particle spectrum. The model predicts emission at TeV energies in good agreement with the recently reported detection.

Mkn 421 was observed for about two days with BeppoSAX, prior to and partly overlapping the start of a 1 week continuous exposure with ASCA in April 1998, as part of a world-wide multiwavelength campaign. A pronounced, well defined, flare observed in X-rays was also observed simultaneously at TeV energies by the Whipple Observatory's 10 m gamma-ray telescope. These data provide the first evidence that the X-ray and TeV intensities are well correlated on time-scales of hours.

Results obtained with BeppoSAX observations of blazars within various collaborative programs are presented. The spectral similarity "paradigm", whereby the spectral energy distributions of blazars follow a sequence, leading to a unified view of the whole population, is briefly illustrated. We concentrate on recent observations of flares and associated spectral variability for three objects at the "blue" end of the spectral sequence, namely PKS 2155-304, Mkn 421 and Mkn 501. The results are discussed in terms of a general analytic synchrotron self-Compton interpretation of the overall spectrum. The physical parameters of the quasi-stationary emission region can be derived with some confidence, while the variability mechanism(s) must be complex.

We propose a new physical model for the unification of all Blazars. The blazar phenomenology can be accounted for by a sequence in the source power and intensity of the diffuse radiation field surrounding the relativistic jet. This regulates the equilibrium electron distribution and hence the shape of the spectral energy distribution (SED), and in turn the classification into a blazar "flavour". The processes inside and beside the emitting jet are the same throughout the blazar class. Objects along this sequence would be observationally classified as High-frequency-peaked BL Lac objects (HBL), Low-frequency BL Lac objects (LBL), Highly Polarized Quasars (HPQ) and Lowly Polarized Quasars (LPQ).

The BL Lac object Mkn 501 was observed with the BeppoSAX satellite at three epochs in April-May 1998, simultaneously with the Whipple and HEGRA Cherenkov telescopes. The X-ray spectrum is well detected up to 70 keV and it exhibits, at all epochs, a continuous curvature, which is here modeled with three power-laws of increasingly steeper index at larger energies. In the nu*f_nu representation the spectrum exhibits a peak at ~20 keV, which is interpreted as the maximum of the synchrotron emission. This implies that the synchrotron peak energy has lowered by an order of magnitude with respect to the powerful X-ray outburst observed in April 1997. The simultaneous TeV flux was comparable to the lowest levels observed for Mkn 501, possibly suggesting that the peak of the inverse Compton radiation had also shifted toward lower energies.

PKS 2155-304 was observed with BeppoSAX in November 1997 for 64 ksec (total elapsed time 33.5 hours) and, for the first time, simultaneously in gamma-rays with EGRET on board the Compton Gamma Ray Observatory and with the ground based TeV telescope CANGAROO, during a phase of high brightness in the X-ray band. The LECS and MECS light curves show a pronounced flare (with an excursion of a factor 3.5 between min and max), with evidence of spectral hardening at maximum intensity. The source is weakly detected by the PDS in the 12-100 keV band with no significant evidence of variability. The broad band X-ray data from Beppo SAX are compared with the gamma-ray results and discussed in the framework of homogenous synchrotron self Compton models.

The phenomenology of gamma-ray bright blazars can be accounted for by a sequence in the source power and intensity of the diffuse radiation field surrounding the relativistic jet. Correspondingly, the equilibrium particle distribution peaks at different energies. This leads to a trend in the observed properties: an increase of the observed power corresponds to: 1) a decrease in the frequencies of the synchrotron and inverse Compton peaks; 2) an increase in the ratio of the powers of the high and low energy spectral components. Objects along this sequence would be observationally classified respectively as high frequency BL Lac objects, low frequency BL Lac objects, highly polarized quasars and lowly polarized quasars. The proposed scheme is based on the correlations among the physical parameters derived in the present paper by applying to 51 gamma ray loud blazars two of the most accepted scenarios for the broad band emission of blazars, namely the synchrotron self--Compton and external Compton models, and explains the observational trends presented by Fossati et al. (1998) in a companion paper, dealing with the spectral energy distributions of all blazars. This gives us confidence that our scheme applies to all blazars as a class.

The apparency of low energy absorption features in Blazar LECS spectra has been reported during the early phase of the BeppoSAX mission (e.g. by Grandi et al. 1997 for 3C~273 SVP data; Giommi et al. 1997 for PKS 2155--304 SVP data), with data reduction and analysis tools released prior to July 1997. We set up a robust method to investigate the reliability of absorption features detections in LECS data of bright blazars. The analysis shows that similar features are common to all datasets (with possibly the exception of 3C~273 SVP data). This suggests that most of the features observed from 0.2 to 1 keV are possibly due to miscalibration of the instrument response. The reported problems have been overcome with the latest release of calibrations (September 1997).

The blazar Mkn 421 has been observed, as part of the AO1 Core Program, five times from 2 to 7 May 1997. In the LECS+MECS energy band the spectrum shows convex curvature, well represented by a broken power--law. Flux variability (more than a factor 2) has been detected over the entire 0.1--10 keV range, accompanying which the spectrum steepens with the decrease in intensity. Mkn 421 has also been detected with the PDS instrument. Our preliminary analysis indicates that the PDS spectrum lies significantly above the extrapolation from the MECS, suggesting a contribution from a flatter high energy component.

We collect data from the radio to the gamma-ray range for three complete samples of blazars: the Slew Survey and the 1Jy samples of BL Lacs and the 2Jy sample of Flat Spectrum Radio-Loud Quasars (FSRQs). The fraction of objects detected in gamma-rays (E > 100 MeV) is 17%, 26% and 40% in the three samples respectively. Except for the Slew Survey sample, gamma-ray detected sources do not differ either from other sources in each sample, nor from all the gamma-ray detected sources, in terms of the distributions of redshift, radio and X-ray luminosities and of the broad band spectral indices (radio to optical and radio to X-ray). We compute average Spectral Energy Distributions (SEDs) from radio to gamma rays for groups of blazars binned according to radio luminosity, irrespective of the original classification as BL Lac or FSRQ. The resulting SEDs show a remarkable continuity in that: i) the first peak occurs in different frequency ranges for different luminosity classes, with most luminous sources peaking at lower frequencies; ii) the peak frequency of the gamma-ray component correlates with the peak frequency of the lower energy one; iii) the luminosity ratio between the high and low frequency components increases with bolometric luminosity. The continuity of properties among different classes of blazars and the systematic trends of the SEDs as a function of luminosity favor a unified view of the blazar phenomenon: a single parameter, related to luminosity, seems to govern the physical properties and radiation mechanisms in the relativistic jets present in BL Lac objects as well as in FSRQ. The general implications of this unified scheme are discussed.

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.

The average Spectral Energy Distributions (SED) of different classes of blazars from the radio to the gamma-ray band are studied adopting a "complete-sample" approach. The SED are globally similar, despite systematic differences in specific frequency ranges, showing two broad peaks. However, i) the first and second peak occur in different frequency ranges for different objects, with a tendency for the most luminous sources to peak at lower frequencies, ii) the ratio between the two peak frequencies *seems to be constant*, iii) the luminosity ratio between the high and low frequency component increases with the bolometric luminosity.

We discuss alternative interpretations of the differences in the Spectral Energy Distributions (SEDs) of BL Lacs found in complete Radio or X-ray surveys. A large body of observations in different bands suggests that the SEDs of BL Lac objects appearing in X-ray surveys differ from those appearing in radio surveys mainly in having a (synchrotron) spectral cut-off (or break) at much higher frequency. In order to explain the different properties of radio and X-ray selected BL Lacs Giommi and Padovani proposed a model based on a common radio luminosity function. At each radio luminosity, objects with high frequency spectral cut-offs are assumed to be a minority. Nevertheless they dominate the X-ray selected population due to the larger X-ray-to-radio-flux ratio. An alternative model explored here (reminiscent of the orientation models previously proposed) is that the X-ray luminosity function is "primary" and that at each X-ray luminosity a minority of objects has larger radio-to-X-ray flux ratio. The predictions of the two scenarios, computed via a Montecarlo technique, are compared with the observed properties of BL Lacs in the two samples extracted respectively from the 1 Jy radio survey and the Einstein Slew Survey. We show that both models can explain a number but not all the observed features. We then propose a completely new approach, based on the idea that the physical parameter which governs the shape of the SEDs, is (or is associated with) the bolometric luminosity. Assuming an empirical relation between spectral shape and luminosity we show that the observational properties of the two surveys can be reproduced at least with the same accuracy as the two previous models.

We discuss alternative interpretations of the differences in the Spectral Energy Distributions (SEDs) of BL Lacs found in complete Radio or X--ray surveys. In order to explain the different properties of radio and X-ray selected within the assumption of a single population we propose a new approach The model is based on the idea that the BL Lacs constitute a ``parameterized family'' and that the physical parameter which governs the shape of the SEDs, is (or is associated with) the bolometric luminosity. Assuming an empirical relation between spectral shape and luminosity the observational properties of the two surveys, including the redshift distributions, can be reasonably well reproduced.

Average Spectral Energy Distributions (SED) for different subgroups of blazars are derived from available homogeneous (but small) data sets, including the gamma-ray band. Comparing Flat Spectrum Radio Quasars (FSRQ) with BL Lacs extracted from radio (RBL) or X-ray surveys (XBL) remarkable differences and similarities are apparent: i) in all cases the overall SED from radio to gamma-rays shows two peaks; ii) the first and second peak fall in different frequency ranges for different objects, with a tendency for the most luminous objects to peak at lower frequencies; iii) the ratio between the two peak frequencies seems to be constant, while the luminosity ratio between the high and low frequency component increases from XBL to RBL and FSRQ. The variability properties, (amplitude and frequency dependence) are similar in different objects if referred to their respective peak frequencies. Finally, comparing spectral snapshots obtained at different epochs, the intensities of the two components at frequencies close and above their respective peaks seem to be correlated. The relevance of these properties for theoretical models is briefly discussed.

ROSAT observations of a large sample of bright gamma-ray (E > 100 MeV) blazars are presented. Results of a detailed spectral analysis in the soft $\sim$ 0.1-2.0 keV energy range are discussed in relation to the overall energy distribution with particular emphasis on the relation between X-ray and gamma-ray properties. A significant anti-correlation between X-ray and gamma-ray spectral shapes of flat radio spectrum quasars (FSRQ) and BL Lacs has been discovered. A different shape in the overall energy distributions from radio to gamma-ray energies between FSRQ and BL Lacs is also implied by the correlation of their broad-band spectral indices $\alpha_{ro}$ and $\alpha_{x \gamma}$. Both the above correlations can be explained if both the IR to UV emission and the hard X-ray to gamma-ray emission originate from the same electron population, via, respectively, the synchrotron process and the inverse Compton mechanism. We suggest that a key parameter for understanding the overall energy distributions of both classes of objects is the energy at which the synchrotron emission peaks in a $\nu-\nu F(\nu)$ representation.

We have observed with ROSAT six blazars from the sample of 31 polarized objects constructed by Impey and Tapia (1988). Results of the spectral analysis in the ROSAT band of three objects, PKS 0301-243, PKS 1055+018 and PG 1424+240 are presented here and compared with the three previously published objects. Using a simple power law model with free absorption a wide range of spectral slopes is found. Objects with emission lines have flatter X-ray spectra than lineless objects. Indication of spectral curvature within the X-ray band is found in the sense of a concavity (soft excess) for the flattest spectrum and a convexity (soft deficiency) for the steepest one. The X-ray spectral shape is discussed in relation to the overall energy distribution for the six objects. The results are compared with other samples of blazars. The X-ray spectral index, alpha_X, is found to be anticorrelated with the broad band spectral index alpha_RO and with bolometric luminosity and redshift. It is suggested that the different spectra of quasar-like and BL Lac-like blazars represent extremes in a continuous distribution rather than two distinct populations.