SSC radiation in BL Lacertae sources, the end of the tether^{}
A. Paggi^{1}  F. Massaro^{2}  V. Vittorini^{3}  A. Cavaliere^{1}  F. D'Ammando^{1,3}  F. Vagnetti^{1}  M. Tavani^{1,3}
1  Dipartimento di Fisica, Università di Roma Tor Vergata, Via della Ricerca Scientifica 1, 00133 Roma, Italy
2 
Harvard, Smithsonian Astrophysical Observatory, 60 Garden Street, Cambridge, MA 02138, USA
3 
INAF, via Fosso del Cavaliere 1, 00100 Roma, Italy
Received 31 March 2009 / Accepted 12 July 2009
Abstract
Context. The synchrotronself Compton (SSC) radiation process is widely held to provide a close representation of the double peaked spectral energy distributions from BL Lac Objects (BL Lacs). This subclass of Active Galactic Nuclei is marked by nonthermal beamed radiations, highly variable on timescales of days or less. Their outbursts in the
rays relative to the optical/X rays might be surmised to be enhanced in BL Lacs as these photons are upscattered via the inverse Compton (IC) process.
Aims. From the observed correlations among the spectral parameters (peak frequencies, fluxes and curvature) during optical/Xray variations we aim at predicting corresponding correlations in the ray band, and the actual relations between the ray and the Xray variability consistent with the SSC emission process.
Methods. We start from the homogeneous singlezone SSC source model, with logparabolic energies distributions of emitting electron as required by the Xray data of many sources. We find relations among spectral parameters of the IC radiation in both the Thomson (for Low energy BL Lacs) and the KleinNishina regimes (mainly for High energy BL Lacs); whence we compute how variability is driven by a smooth increase of key source parameters, primarily the root mean square electron energy.
Results. In the KleinNishina regime the model predicts for HBLs lower inverse Compton fluxes relative to synchrotron, and milder ray relative to Xray variations. Stronger ray flares observed in some HBLs like Mrk 501 are understood in terms of additional, smooth increases also of the emitting electron density. However, episodes of rapid flares as recently reported at TeV energies are beyond the reach of the single component SSC model with one dominant varying parameter. Furthermore, spectral correlations at variance with our predictions, as well as TeV emissions in LBL objects (like BL Lacertae itself) cannot be explained in terms of the simple HSZ SSC model, and in these cases the source may require additional electron populations in more elaborate structures like decelerated relativistic outflows or subjet scenarios.
Conclusions. We provide a comprehensive benchmark to straightforwardly gauge the capabilities and effectiveness of the SSC radiation process. The single component SSC source model in the Thomson regime turns out to be adequate for many LBL sources. In the mild KleinNishina regime it covers HBL sources undergoing variations driven by smooth increase of a number of source parameters. However, the simple model meets its limits with the fast/strong flares recently reported for a few sources in the TeV range; these clearly require sudden accelerations of emitting electrons in a second source component.
Key words: galaxies: active  galaxies: jets  BL Lacertae objects: general  radiation mechanisms: nonthermal
1 Introduction
Blazars are among the brightest Active Galactic Nuclei (AGNs) observed, with inferred (isotropic) luminosities up to . Actually these sources are widely held to be AGNs radiating from a relativistic jet closely aligned to the observer line of sight; this emits nonthermal radiations with observed fluxes enhanced by relativistic effects so as to often overwhelm other, thermal or reprocessed emissions.
BL Lacs in particular are Blazars showing no or just weak emission lines. Their spectra may be represented as a continuous spectral energy distribution (SED) marked by two peaks: a lower frequency one, widely interpreted as synchrotron emission by highly relativistic electrons; and a higher frequency counterpart, believed to be inverse Compton upscattering by the same electrons of seed photons either emitted by an external source (external Compton, EC, see e.g. Sikora et al. 1994), or provided by the photons of the very synchrotron radiation (synchrotronself Compton, SSC; see, e.g., Marscher & Gear 1985; Maraschi et al. 1992).
The BL Lacs are usually classified in terms of the frequency of their synchrotron peak (see Padovani & Giommi 1995): for the low frequency BL Lacs (LBLs) the peak lies in the infraredoptical bands, while the inverse Compton (IC) component peaks at MeV energies; instead, high frequency BL Lacs (HBLs) feature a first peak in the Xray band and the second one at energies or beyond; intermediate frequency BL Lacs (IBLs) stand in the middle, with the first peak at optical frequencies and the second peak at 1 GeV.
BL Lac Objects also show rapid variability on timescales of days or shorter, during which the sources undergo strong flux variations often named ``flares''.
The homogeneous single zone (HSZ) SSC model (based, in particular, on a single electron population) constitutes an attractively simple source structure worth to be extensively tested and pushed to its very limits, as we discuss below; following Katarzynski et al. (2005), Tramacere et al. (2007), we focus on deriving the simultaneous changes of the synchrotron and IC fluxes of BL Lacs to be expected from an updated version of SSC. We use a realistic logparabolic energy distribution of the electrons radiating the synchrotron photons and upscattering them by IC in either the Thomson or the Klein Nishina (KN) regime.
2 Background relations
2.1 Logparabolic electron distribution and spectra
The SEDs of BL Lacs exhibit two distinct peaks or humps. To describe these in detail a number of spectral models have been used over the years: powerlaws with an exponential cutoff, broken powerlaws, or spectra curved in the shape of a logparabola. The latter are being increasingly used, as they generally provide more accurate fits than broken powerlaws in the optical Xray bands, with residuals uniformly low throughout a wide energy range (see discussion by Landau et al. 1986; Tanihata et al. 2004; Massaro et al. 2004a; Tramacere et al. 2007; Nieppola et al. 2006; Donnarumma et al. 2009)^{}.
So we adopt logparabolae in view of their optimal performance as a precise fitting tool. Moreover, they are also effective in yielding simple, analytical expressions (given in their exact form in Appendix A, following Eqs. (A.6), (A.9) and (A.15)) for peak frequencies and fluxes in both in the inverse Compton regimes, Thomson and KleinNishina. Finally, the physically interesting feature of logparabolic spectra is their direct link to the acceleration processes of the emitting particles; in fact, as shown in Appendix B, logparabolic emitted spectra simply relate to logparabolic particle energy distributions, that in turn naturally arise from systematic and stochastic acceleration processes as dealt with by the FokkerPlanck equation.
Thus for the synchrotron SED we write
here gives the constant component of the spectral index for the energy flux and is the spectral curvature. In the following we shall denote with the frequency where the SED peaks and by S the related flux. Spectral curvature values are observed to be around 0.2 (Massaro et al. 2004b) and decreasing with time and/or frequency as given by Eqs. (B.11) and (B.13) of Appendix A; since the decrease is mild during flares (which are essentially dominated by systematic acceleration processes), we can consider a nearly constant spectral curvature to a first approximation. Then one finds the scaling relations and , where n is the particle density, R is the emitting region radius, B is the magnetic field, is the beaming factor, and is the square root of the mean quadratic energy (thereafter rms energy) of the particle distribution.
An analogous expression holds for the inverse Compton radiation
which in the Thomson regime the IC has less curved spectra with . At nearly constant curvature the parameters to focus on are and for the SED peak frequency and flux, respectively. In the extreme KleinNishina regime we obtain instead more curved spectra with , and and hold (see Appendix A for more details, and massaro2008a,b).
Such spectra are radiated by electron populations with a logparabolic energy distribution, namely
The result is easily inferred and straightforward to compute on using the deltafunction approximations for the single particle emissions (see Appendix A); it is found to hold also for the exact shape of the latter with the numerical simulations outlined in Fig. 1, and reported in full detail in Figs. B.1 and B.2.
Figure 1: Two examples of exact, steady SEDs provided by numerical simulations (code by Massaro 2007) after the HSZ SSC model from relativistic electrons with a logparabolic energy distribution (see Appendix B). The two characteristic humps correspond to synchrotron emission ( left) and to inverse Compton scattering ( right), either in the Thomson regime (blue line) or in the KN regime (red line). Standard sources that are so fitted include the low state of 0716+714 (an IBL, blue crosses) observations of November 1996 by BeppoSAX, AIT and EGRET (Tagliaferri et al. 2003); and 1ES 1553+113 (a HBL, red crosses) observations of April 2005 by Swift and AprilMay 2005 by MAGIC (Massaro 2007). 

Open with DEXTER 
Table 1: Source parameters for some BL Lac sources.
2.2 Source parameters
Using the previous scalings, with the relation
between the observed variation time
and the size R of the emitting region (assuming spherical symmetry in the rest frame), we obtain relations of source parameters with the observables. In the Thomson regime we find
For the extreme KN regime we find instead:
Note that a delay of the order of may elapse between synchrotron and IC spectral variations.
Numerical values for the prefactors are given in Appendix C, where normalizations are chosen referring to a typical LBL source with , , , ; and to a typical HBL source with , , , . Values of physical parameters for some sources are given in Table 1. Note that the simple Eqs. (4)(8) provide parameter values in close agreement with the ones needed to fit the observations of LBL or IBL sources in Thomson regime; on the other hand the extreme KN limit underlying Eqs. (9)(13) would yield disagreement for HBLs, indicating that these sources just border the KN regime.
We evaluate the position
of the synchrotron peak when the transition between the two IC regimes occurs, to find
or equivalently
So we conclude that in LBLs, which feature the synchrotron peak at optical/IR frequencies with , the majority of photons of the IC peak are upscattered in the Thomson regime; whereas in extreme HBLs, which show the synchrotron peak in the Xray band with , the majority of IC peak photons are upscattered in the KleinNishina regime.
2.3 Spectral correlations during Xray spectral variations
Correlations in the synchrotron emission between S and with can be used to pinpoint the main driver of the spectral changes during Xray variations (see also Tramacere et al. 2007). In fact, the synchrotron SED peak scales (for nearly constant spectral curvature as anticipated in Sect. 2.1) as at the peak frequency itself scaling as .
Thus on writing the dependence of S on in the form of a powerlaw , we expect to apply when the spectral changes are driven mainly by variations of the electrons rms energy; for dominant changes of the magnetic field; if changes in the beaming factor dominate; formally applies for changes only in the number density of the emitting particles. Results by the latter authors focus on and B as dominant drivers. Starting from such correlations we aim at predicting the expected correlations and flux variations in the ray band.
Particle rms energy variations are of particular interest because they naturally arise as a consequence of systematic plus stochastic acceleration of the electron population (details in Appendix B).
3 ray spectra
3.1 Spectral correlations
The above results lead us to expect specific correlations between SED peak frequencies and fluxes as follows. In the case of dominant rms particle energy variations we expect for the synchrotron emission
while for the IC we have
(17) 
for LBL objects, and
(18) 
for extreme HBLs.
In the case of magnetic field variations we have for the synchrotron emission
(19) 
while for the IC we have
(20) 
for LBL objects, and
formally with for extreme HBLs.
Figure 2: To illustrate the differences of the IC SEDs in the Thomson and in the extreme KN regime, we show the results of numerical simulations (code by Massaro 2007) after the HSZ SSC model, for different sources with larger and larger rms energy of the radiating particles and other parameters kept constant (including the curvature of the particle energy distribution). The vertical line indicates the frequency of the synchrotron peak where the IC scattering changes over the Thomson to the KN regime (see Eq. (14)). Actual sources in their evolution (represented in Figs. B.1 and B.2) vary their spectral curvature somewhat and only span a limited range in frequency and flux to the left (LBLs) or to the right (HBLs) of the line, with only the IBLs likely to cross it. 

Open with DEXTER 
To complete the above picture, we consider how saturation affects correlations. Denoting the total number of particles with , in the jet frame we have for the synchrotron luminosity , while for IC luminosity we have and in Thomson and KN regime, respectively. Consider the case where the total power available to the jet is limited, e.g., in the BZ mechanism (Blandford & Znajek 1977) of power extraction from a maximally rotating black hole with an extreme, radiation pressure dominated disk that holds the pressure of the magnetic field facing the black hole horizon (see discussion by Cavaliere & D'Elia 2002, and references therein). Then the power has an upper bound ; as a consequence, we expect an effect of saturation on the SED peaks as the emitted power, which represents a substantial fraction of the total jet power (see Ghisellini et al. 2008), approaches . For the total power, including radiative, kinetic and magnetic components, we envisage , that is, the particles effectively accelerated to the increasing decrease in number on approaching the BZ limit.
The data analyses by Tramacere et al. (2007, 2009) for the variations of Mrk 421 are consistent with the beginning of such a saturation effect, as they show a declining value of when the peak moves toward the highest energies and approaches . Partial saturation effects can also be explained in terms of systematic and stochastic acceleration processes; in particular, when the stochastic prevails over the systematic acceleration (as is bound to occur at high energies, see Appendix A) we expect to decrease toward 0.6, i.e. .
3.2 Related variabilities
We now focus on related variabilities between IC and synchrotron fluxes; here we are mainly interested on increases of peak frequencies and fluxes, as driven by acceleration processes of the emitting particles or by a growing magnetic field (see also Katarzynski et al. 2005).
For LBLs, that radiate mostly in the Thomson regime, we predict from Eqs. (16)(21) a quadratic or linear variation of C with respect to S, depending on the parameter that mainly drives source variations. That is to say, we expect
in the case of dominant particle rms energy variations, or
for dominant magnetic field variations (see also Figs. 2 and 4).
Instead, for HBLs that radiate closer to the KleinNishina regime, we obtain a weaker ray variability due to the decreasing KN cross section, leading in the extreme to
in case of dominant particle rms energy variations, or
for dominant magnetic field variations. Expected correlations are summarized in Table 2.
From Figs. 2 and 4 it is seen that the transition to the KN regime has three effects on the IC spectrum: first, it brakes the peak frequency increase, because the energy the photon can gain is limited to the total electron energy ; second, it reduces the flux increase at the SED peak, as a consequence of reduced cross section; third, it increases the spectral curvature near the peak as a consequence of frequency compression.
Viceversa, variations observed in IC section of the spectrum are related via Eqs. (22)(25) to variations expected in the synchrotron emission; in particular, in HBLs ray variations ought to have enhanced counterparts in X rays.
Up to now Figs. 2 and 4 may be interpreted as a collection of different sources or of unrelated states of a single source; henceforth we will focus on the interpretation in terms of evolution in a flaring source where the parameters vary in a continuous fashion; in particular, increases under the drive of systematic and stochastic electron accelerations, in keeping with the description in terms of a continuity or kinetic equation of the FokkerPlanck type as given in detail in Appendix B. Signatures of such an evolution are provided not only by continuous growth of the peak frequencies, but even more definitely by the irreversible decrease of the spectral curvature (in terms of time or frequency) under the drive of the stochastic acceleration, see Eq. (B.13) and the observations by Tramacere et al. (2007, 2009). In fact, for a single electron population the curvature is to slowly decrease or stay nearly constant on the timescale of the systematic acceleration. So a sudden increase of curvature signals the injection of a new population/component.
In the following we will neglect radiative cooling, which does not strongly affect the spectral shape around the peaks; such a process will be particularly relevant when the source, after flaring up, relaxes back to a lower state upon radiating away its excess energies from high frequencies downwards (see Appendix B).
Table 2: Spectral correlations (for ). We have denoted with or (B) the variations driven by increases of rms electron energy or magnetic field, respectively.
Figure 3: Computed changes of spectral parameters corresponding to extended increase of rms particle energy with the other source parameters kept constant. As functions of (where is the synchrotron SED peak frequency, see Eq. (14)), the full line represents the ratio , the dashed line the ratio and the dotdashed line the curvature ratio at the peaks. Circles represent the parameters of 0716+714 as observed on November 1996 and September 2007 (Tagliaferri et al. 2003; Chen et al. 2008), and diamonds represent Mrk 501 observations of April 1997 (DjannatiAtai et al. 1999; Massaro et al. 2006), each with the same colours as the lines. Changes of n only would yield at all frequencies, represented by the horizontal dotted line; added to a pure rms energy increase, a moderate increase of n canaccount for the upward deviations of observed values. 

Open with DEXTER 
Figure 4: As in Fig. 2, we show SEDs obtained from numerical simulations (code by Massaro 2007) of HSZ SSC radiations, for different sources with larger and larger magnetic field B and other parameters kept constant. The vertical line again indicates the transition frequency from the Thomson to the KN regime (see Eq. (15)); as B increases, the IC scattering changes its regime, and reduced cross section yields similar compression effects as in Fig. 2. 

Open with DEXTER 
Figure 5: Computed changes of spectral parameters as driven by an extended increase of the magnetic field with the other source parameters kept constant. The black full line represents the value of the ratio , the red dashed line represents the value of the ratio and the blue dotdashed line represents the ratio , for different values of (where is the synchrotron SED peak frequency, see Eq. (15)), data points are the same as in Fig. 3; compared with the latter considerably larger changes of n are needed to recover the observed values. 

Open with DEXTER 
3.3 Specific sources
We have applied the above simple expectations to two sources, 0716+714 (IBL) and Mrk 501 (HBL), among the few to date to provide extended simultaneous coverage of two different source states in X rays and rays.
The first is widely considered to be an IBL. In the low state observed by Automatic Imaging Telescope (AIT) in November 14 1996, the synchrotron peak at about at about falls quite below the threshold frequency expressed by Eqs. (14) or (15), so we expect for this source the IC scattering to occur mostly in the Thomson regime; simultaneous ray observations were carried out by EGRET. The high activity state of the source is described by data relative to GASP project of the WEBT and AGILEGRID observations of 2007 September 712 (see Fig. 6). The source has shown similarly low flux levels during EGRET observations (Tagliaferri et al. 2003); so for flares with short duty cycle, it is not unreasonable to evaluate variations between the two distant states in the absence of closer observations; curvature variation is not required to fit the data, as we expect at lower energies where systematic dominates the stochastic acceleration (see Appendix A). As shown in Figs. 3 and 5, peak flux and frequency variations of 0716+714 appear to be in between the lines representing increasing rms particle energy and increasing magnetic field, and so they may be described in terms of HSZ SSC with simultaneous variations of these two parameters.
For the HBL source Mrk 501 we consider the two states described by Massaro et al. (2006), with simultaneous BeppoSAX and CAT observations of April 7 and 16 1997. To understand the behavior of this source, from Figs. 3 and 5 it is seen that variations of B alone are ineffective and require complementary n increase by a factor of about 30; variations of are adequate in the KN regime, and only require a complementary doubling of n (which by itself would yield in both scattering regimes) to yield higher ray flux increase. Note that in going to higher energies the spectral curvature decreases somewhat (from to for the low and high state, respectively) in keeping with model predictions; such a behavior is interesting as it marks a smooth growth, if anything, in the number of emitting particles as contrasted with sudden, substantial reinjection of nearly monoenergetic electrons that would suddenly reverse the otherwise irreversible decrease of r and b. Such a smooth growth obtains in a scenario of an expanding blastwave that progressively involves more electrons (see, e.g., Ostriker & McKee 1988; Lapi et al. 2005; Vietri 2006).
We stress that in HBL sources even small flux variations in rays are expected to have enhanced and observable counterparts in X rays. We suggest that such variations should be checked upon ray ``alarms'', the inverse triggering relative to usual. The absence of such lower energies counterparts will indicate, for example, a flare driven dominantly by a particle number density increase associated with a magnetic field decrease , causing but a decreasing synchrotron peak frequency; such a kind of flare may easily drown into the primary synchrotron component.
Figure 6: Spectral fit of two different states of the IBL source 0716+714. The low state (blue line) data are relative to 1996 November 14 simultaneous observation by AIT, BeppoSAX and EGRET; the high state (red line) data are relative to 2007 September 712 simultaneous observations by GASPWEBT and AGILEGRID (Tagliaferri et al. 2003; Chen et al. 2008). 

Open with DEXTER 
3.4 Limiting timescales
As shown in Appendix B, systematic and stochastic acceleration processes occur in the sources on timescales and , respectively. Limits to the variability timescales relate to the apparent size of the emission region by (see Sect. 2.2); rapid variations require small emitting sources, but these could be optically thick to the pair production process. In fact, ray photons collide with less energetic ones to produce e pairs, and the cross section of this process tops at a value around (where is the Thomson cross section) when the ray photon frequency and the target photons frequency satisfy the relation .
For the
rays to escape from the emitting region the source ``compactness'' intervenes and can set lower bounds on the beaming factor (Cavaliere & Morrison 1980; Massaro 2007; Begelman et al. 2008). In fact the optical depth for this process is expressed in terms of observational quantities (primed quantities refer to the jet rest frame) as
(26) 
where D is the source distance expressed in , with the shorthand . Escaping ray radiation requires the source to be optically thin, that is, , resulting in
(27) 
Thus rapid flux variations set a lower bound on the beaming factor; the highest bound between the previous relation and is to be relevant.
The exponent contains a weak logarithmic dependence on that may be neglected to a first approximation, and has a magnitude which depends on the source spectral properties; for definiteness, we use for b the typical value 0.2, and for . For a LBL with we obtain , while for an extreme HBL with , we obtain . A variation timescale min for a LBL with implies , while in an extreme HBL with we find .
In the specific case of PKS 2155304 discussed by Begelman et al. (2008), the emitted power corresponds to , with , to give and for min. We note that the source, even if widely considered a HBL, shows peculiar behaviors with respect to other TeV HBLs (see Massaro et al. 2008b); moreover, it does not satisfy the condition and so it cannot be considered an extreme HBL, as discussed below.
Extreme values of the beaming factor have been proposed to account for large peak separations (Konopelko et al. 2003) or to explain very rapid spectral variations (Begelman et al. 2008), formally with min for the ray flares of PKS 2155304 in 2006 July 28 (Aharonian et al. 2007) and Mrk 501 in June 30 and July 9, 2007 (Albert et al. 2007a). Note that with BL Lac spectra realistically curved, the observed flux variations may be enhanced due to a slope effect best legible on the differential flux ; that is, when fluxes are measured at frequencies where the spectral slope is steep (as may be the case for PKS 2155304), a strong observed flux variation implies only a mild variation of the peak flux and requires smaller values of . Otherwise, a peak variation on a scale min would require in the SSC model for LBLs where (see Eq. (7)), and for an extreme HBL where (see Eq. (12)).
As stated under Sect. 3.3, in HBL sources even small flux variations in rays should have enhanced counterparts in X rays. These may become hard to observe for example in a flare dominantly driven by particle number density increase associated with a magnetic field decrease (i.g., with ), leading to a synchrotron emission easily drowned into other components. On the other hand, components with appearing in an HBL spectrum go beyond the HSZ SSC model with variations of one dominant parameter, and therefore require a more elaborate source structure.
In such cases the next natural scenarios are provided by decelerated relativistic outflows (Georganopoulos & Kazanas 2003), or by nested spinelayer jets (Tavecchio & Ghisellini 2008) and jets in a jet (Giannios et al. 2009).
4 Beyond the SSC
In Table 2 we have provided a comprehensive benchmark to gauge the performance of the HSZ SSC model for flaring BL Lac objects. This allows us to easily recognize events that may be accounted for within the model with physical variations of key source parameters, or that instead require more complex source structures.
Toward that purpose, we have used realistic logparabolic spectral shapes produced by logparabolic energy distributions of the emitting electrons, and studied IC radiation both in the Thomson and the KleinNishina regimes.
In the model we expect S to dominate C fluxes; in fact, moving to higher energies in a collection of different sources or in a prolonged evolution of a given source, we expect the emission to drift out of the pure Thomson and approach the KN regime, where the cross section decreases and limits the IC fluxes enforcing . Such a relation is found to hold in a number of sources, and in particular in HBLs (see for example Tagliaferri et al. 2008).
From Table 2 we stress here the source spectral variations predicted in rays. For LBLs, where the IC scattering mostly occurs in the Thomson regime, we recover the standard quadratic increase in ray fluxes with respect to IRoptical ones, expressed by for dominant particle rms energy variations (see Eq. (22)). Instead for HBLs, in which the IC scattering approaches the KleinNishina regime, we expect smaller or even vanishing increases of the ray flux relative to X rays, that is, for dominant particle rms energy variations (see Eq. (24)).
Comparing with specific sources, we find the HSZ SSC model to be adequate for most LBLs at the present observational levels; whilst for example the HBL source Mrk 501 in April 1997 showed a ray flux increase appreciably stronger than expected if it were dominated by just one driving parameter. We explain this behavior within the model in terms of an additional, smooth and moderate increase (by a factor around 2) in number density of the electrons responsible for the emission (see Sect. 3.3). Such conditions can still be provided by a single electron population; this should be marked by a continuously decreasing spectral curvature b, as indicated by current data (Massaro et al. 2006), a feature providing a potentially powerful signature to closely check on further data. Instead, for dominant magnetic field variations we would expect for LBL sources (see Eq. (23)), and for HBL sources (see Eq. 25)); but in that case the required increase in particle number density should be much higher by a factor of about 30, hard to interpret in terms of a single electron population.
On the other hand, on using the previous relations inversely, we see that even small ray variations in HBL sources ought to have enhanced counterparts in the X rays, unless flare activity were driven by an additional jet component, for example one with higher magnetic field but lower particle density; thus the corresponding emission is easily overwhelmed by, or drowned into the main synchrotron, but then the rest frame acceleration times must be short enough to involve a significant fraction of the emitting region. We suggest this as a critical test for the simple model.
Another limitation arises from rapid peak flux variations requiring large values of as shown in Sect. 3.4. In the extreme, the few sources with particularly fast ray increases observed so far, like PKS 2155304 in 2006 July 28 and Mrk 501 in June 30 and July 9, 2007 (see Sect. 3.4), require additional components with very high beaming factors that go definitely beyond the simple SSC model.
Finally, another problem for the model arises from LBL objects showing substantial emissions in the GeVTeV range that cannot be explained it terms of the simple HSZ SSC model (considering that second order IC scattering is negligible, see Massaro 2007) because the IC peak falls for these sources around ; this may be the case for BL Lacertae itself (Albert et al. 2007b) and possibly similar sources like M87 (Georganopoulos et al. 2005) and Cen A (Lenian et al. 2008; Aharonian et al. 2009).
In all these cases the source may require more elaborate structure, like decelerated relativistic outflows or subjet scenarios (see Georganopoulos & Kazanas 2003; Tavecchio & Ghisellini 2008; Giannios et al. 2009). We stress that the injection of a second, monoenergetic electron population is expected to be marked by a sudden increase of the spectral curvature.
5 Discussion and conclusions
We propose that the sources of BL Lac type may be conveniently ordered in a succession spanning from smooth variations of one or a few dominant SSC parameter, up to the appearance of truly different components that ultimately break through the limits of the simple HSZ SSC model. Even more so at increasing energies and frequencies, that would imply weaker and weaker ray fluxes relative to Xray, owing to KN cross section effects. This picture is supported by the following two lines of evidence.
First, a recent statistical study of Third EGRET Catalogue data (Mukherjee et al. 2001; Casula 2008; Vagnetti et al., in prep.) shows for BL Lac objects observed at a weak ray variability on average, compared to the FSRQ Blazars. Also the timestructure functions for these two classes of objects indicate a similar trend, within the limitations of the sample.
Second, the first AGILEGRID Catalogue of high confidence gammaray sources (Pittori et al. 2009) apparently shows few if any new BL Lac sources other than those already detected by EGRET. This circumstance suggests the sources observed by AGILE to be more powerful than the average; it leads to expect for the majority of the sources either reduced flare activity in this band, or weak average fluxes. We relate these features to reduced KleinNishina cross section that tends to limit both the average fluxes and the flares.
With Fermi, the first three months of observations Abdo et al. (2009) have yielded a substantial number of new BL Lac sources, and in particular a larger fraction of HBLs close to the ratio known from radio surveys (Nieppola et al. 2006). This is consistent with our expectations, considering the better instrumental sensitivity, especially at high energies generally favorable to the HBLs; for these sources we suggest to test the model's outcomes: C lower than S, and generally small (though possibly fast) flux variations compared to the LBLs
We plan to check our theoretical predictions against the public availability of the first year data from Fermi observations, and with simultaneous, multiwavelength observations comparing the two basic subclasses of LBL and HBL of the BL Lac sources.
Acknowledgements
We thank our referee for useful comments and helpful suggestions. F. Massaro acknowledges the Foundation BLANCEFLOR BoncompagniLudovisi, née Bildt, for the grant awarded him in 2009 to support his research.
References
 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, [arXiv:0902.1559] (In the text)
 Aharonian, F., Akhperjanian, A. G., BazerBachi, A. R., et al. 2007, ApJ, 664, L71 [NASA ADS] [CrossRef] (In the text)
 Aharonian, F., Akhperjanian, A. G., BazerBachi, A. R., et al. 2009,[arXiv:0903.1582] (In the text)
 Albert, J., Aliu, E., Anderhub, H., et al. 2007a, ApJ, 669, 862 [NASA ADS] [CrossRef] (In the text)
 Albert, J., Aliu, E., Anderhub, H., et al. 2007b. ApJ, 666, L17 (In the text)
 Begelman, M. C., Fabian, A., & Rees, M. J. 2008, MNRAS, 384, L19 [NASA ADS] (In the text)
 Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433 [NASA ADS] (In the text)
 Bloom, S. D., & Marscher, A. P. 1996, ApJ, 461, 657 [NASA ADS] [CrossRef]
 Blumenthal, G. R., & Gould, R. J. 1970, Rev. Mod. Phys., 42, 237 [NASA ADS] [CrossRef]
 Casula, V. 2008, Thesis (In the text)
 Cavaliere, A., & D'Elia, V. 2002, ApJ, 571, 226 [NASA ADS] [CrossRef] (In the text)
 Cavaliere, A., & Morrison, P. 1980, ApJ, 238, L63 [NASA ADS] [CrossRef] (In the text)
 Chen, A. W., D'Ammando, F., Villata, M., et al. 2008, A&A, 489, L37 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 DjannatiAtai, A., Piron, F., Barrau, A., et al. 1999, A&A, 350, 17 [NASA ADS] (In the text)
 Donnarumma, I., Vittorini, V., Vercellone. S., et al. 2009, ApJ, 691, L13 [NASA ADS] [CrossRef] (In the text)
 Fossati, G., Buckley, J. H., Bond, I. H., et al. 2008, ApJ, 677, 906 [NASA ADS] [CrossRef]
 Georganopoulos, M., & Kazanas, D. 2003, ApJ, 594, 27 [NASA ADS] [CrossRef] (In the text)
 Georganopoulos, M., Perlman, E. S., & Kazanas, D. 2005, ApJ, 634, L33 [NASA ADS] [CrossRef] (In the text)
 Giannios, D., Uzdensky, D. A., & Begelman, M. C. 2009, [arXiv:0901.1877]
 Giommi, P., Colafrancesco, S., Cutini, S., et al 2008, A&A, 487, L49 [NASA ADS] [CrossRef] [EDP Sciences]
 Inoue, S., & Takahara, F. 1996, ApJ, 463, 555 [NASA ADS] [CrossRef]
 Jones, F. C. 1968, Phys. Rev., 167, 1159 [NASA ADS] [CrossRef] (In the text)
 Kaplan, S. A. 1956, Sov. Phys., 2, 2 (In the text)
 Kardashev, N. S. 1962, SvA, 6, 317 [NASA ADS]
 Kataoka, J., Mattox, J. R., Quinn, J., et al. 1999, ApJ, 514, 138 [NASA ADS] [CrossRef]
 Katarzynski, K., Ghisellini, G., Tavecchio, F., et al. 2005, A&A, 433, 479 [NASA ADS] [CrossRef] [EDP Sciences]
 Konopelko, A., Mastichiadis, A., Kirk, J., et al. 2003, ApJ, 597, 851 [NASA ADS] [CrossRef] (In the text)
 Landau, R., Golisch, B., Jones, T., et al. 1986, ApJ, 308, 78 [NASA ADS] [CrossRef] (In the text)
 Lapi, A., Cavaliere, A., & Menci, N. 2005, ApJ, 619, 90 [NASA ADS] [CrossRef] (In the text)
 Lenain, J.P., Boisson, C., & Sol, H. 2008, [arXiv:0807.2733] (In the text)
 Maraschi, L., Ghisellini, G., & Celotti, A. 1992, ApJ, 397, L5 [NASA ADS] [CrossRef] (In the text)
 Marscher, A. P. 1996, Energy Transport in Radio Galaxies and Quasars, ASP Conf. Ser., 100, 45
 Marscher, A. P., & Gear, W. K. 1985, ApJ, 298, 114 [NASA ADS] [CrossRef] (In the text)
 Massaro, F. 2007, Ph.D. Thesis (In the text)
 Massaro, E., Perri, M., Giommi, P., & Nesci, R. 2004a, A&A, 413, 489 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Massaro, E., Perri, M., Giommi, P., Nesci, R., & Verrecchia, F. 2004b, A&A, 422, 103 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Massaro, E., Tramacere, A., Perri, M., Giommi, P., & Tosti, G. 2006, A&A, 448, 861 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Massaro, F., Tramacere, A., Cavaliere, A., Perri, M., & Giommi, P. 2008a, A&A, 478, 395 [NASA ADS] [CrossRef] [EDP Sciences]
 Massaro, F., Giommi, P., Tosti, G., et al. 2008b, A&A, 489, 1047 [NASA ADS] [CrossRef] [EDP Sciences]
 Mukherjee, R., et al. 2001, in High Energy GammaRay Astronomy, ed. F. A. Aharonian, & H. J. Völk (New York: AIP), AIP Conf. Proc., 558, 324 (In the text)
 Nieppola, E., Tornikoski, M., & Valtaoja, E. 2006, A&A, 445, 441 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Ostriker, J. P., & McKee, C. F. 1988, Rev. Mod. Phys., 61, 1 [NASA ADS] [CrossRef] (In the text)
 Padovani, P., & Giommi, P. 1995, ApJ, 444, 567 [NASA ADS] [CrossRef] (In the text)
 Paggi, A. 2007, Thesis (In the text)
 Pittori, C., Verrecchia, F., Chen, A. W., et al. 2009, [arXiv:0902.2959] (In the text)
 Rybicki, G. B., & Lightman, A. P. 1979, Radiative Processes in Astrophysics (New York: Wiley)
 Sikora, M., Begelman, M. C., & Rees, M. J. 1994, ApJ, 421, 123 [NASA ADS] [CrossRef] (In the text)
 Tagliaferri, G., Ravasio, M., & Ghisellini, G. 2003, A&A, 400, 477 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Tagliaferri, G., Foschini, L., Ghisellini, G., et al. 2008, ApJ, 679, 1029 [NASA ADS] [CrossRef] (In the text)
 Tanihata, C., Kataoka, J., Takahashi, T., et al. 2004, ApJ, 601, 759 [NASA ADS] [CrossRef] (In the text)
 Tavecchio, F., & Ghisellini, G. 2008, [arXiv:0810.0134] (In the text)
 Tavecchio, F., Maraschi, L., & Ghisellini, G. 1998, ApJ, 509, 608 [NASA ADS] [CrossRef]
 Tramacere, A., Massaro, F., & Cavaliere, A. 2007, A&A, 466, 521 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Tramacere, A., Giommi, P., Perri, M., Verrecchia, F., & Tosti, G. 2009, [arXiv:0901.4124v1]
 Vietri, M. 2006, Astrofisica delle alte energie, Torino, Bollati Boringhieri (In the text)
Online Material
Appendix A: Logparabolic spectra
In this Appendix we show how electron populations with a logparabolic energy distribution of the form expressed by Eq. (3), that is,
emit logparabolic spectra via the SSC process. The related particle synchrotron emissivity
(A.2) 
is easily computed on using the close approximation to the single particle emission in the shape of a deltafunction (see Rybicki & Lightmann 1979), that is, with and (with B measured in Gauss) is the synchrotron critical frequency. This leads to (Massaro et al. 2004a) a logparabolic differential flux
(A.3) 
and to a SED again of logparabolic shape
(A.4) 
its slope at the synchrotron reference frequency is given in terms of s, by
(A.5) 
the spectral curvature by
and the peak value occurs at a frequency ( , n, B and R are defined in Sect. 2.1 of the main text).
For IC radiation in the Thomson regime we may write to a fair approximation
(see Rybicky & Lightmann 1979) where
is the power radiated by a singleparticle IC scattering in the Thomson regime, having denoted with
the synchrotron radiation density^{}
. We obtain once again a logparabolic SED
(A.7) 
where the slope at the IC reference frequency is given by
(A.8) 
and the spectral curvature reads
The peak value is attained at a frequency .
For the KleinNishina (KN) regime instead it necessary to consider the convolution
(A.10) 
where and are the electron frequencies before and after the scattering, respectively, is the number spectrum of seed photons, and is the full Compton kernel (Jones 1968). Only in the extreme KN regime one may again approximate
(A.11) 
on approximating with its mean value for a homogeneous, spherical optically thin source
(A.12) 
one obtains again a SED with a logparabolic shape
(A.13) 
Now the slope at
(A.14) 
is steeper, and the spectral curvature
is larger than in the Thomson regime. The peak value occurs at a frequency .
The transition between the two regimes occurs when
,
where
,
that is, when
(A.16) 
holds, or equivalently
(A.17) 
We close these calculations with two remarks. First, the (primed) quantities used here refer to the rest frame of the emitting region, while the observed (unprimed) quantities must be multiplied by powers of the beaming factor
(A.18) 
being the bulk Lorentz factor of the relativistic electron flow in the emitting region, the angle between its velocity and the line of sight, and ; so we obtain , , and , . Second, we note that the actual specific powers radiated by a single particle differ from deltafunction shape, slightly for the synchrotron emission and considerably for the IC radiation in the Thomson regime (see discussion in Rybicky & Lightmann 1979). However, the above spectral shapes still approximatively apply, as confirmed by numerical simulations (Massaro 2007); in detail, the convolution with a broader single particle power yields a less curved spectrum, so that Eqs. (A.6) and (A.9) become (Massaro et al. 2006) and , respectively.
Appendix B: Particle acceleration processes
Figure B.1: Example of time evolution of the distribution (with peak energy ) due to stochastic and systematic accelerations ( upper panel) with , from the initial value , and of the related evolution for the SSC SEDs ( lower panel). In terms of the stochastic acceleration time , the time interval between each pair of lines is , corresponding to an observed time interval of about two days. 

Open with DEXTER 
Figure B.2: Example of time evolution of the distribution (with peak energy ) due to stochastic and systematic accelerations ( upper panel) with , from the initial value , and of the related evolution for the SSC SEDs ( lower panel). In terms of the stochastic acceleration time , the time interval between each pair of lines is , corresponding to an observed time interval of about two days. 

Open with DEXTER 
Here we derive a logparabolic electron energy distribution
from a kinetic continuity equation of the FokkerPlanck type; following Kardashev (1962), in the jet rest frame this reads
where is the particle energy, t denotes time, and and describe systematic and stochastic acceleration rates, occurring on timescales and , respectively. For example, in the picture of Fermi accelerations (e.g., Vietri 2006), the accelerations rates and can be expressed in terms of physical quantities related to processes occurring in shocks; in this framework a plane shock front of thickness moves with speed and gas clouds of average size move downstream of the shock with speed V; it is found that and hold (see Kaplan 1956; Paggi 2007), . On the other hand, numerical values of and are directly derived from the emitted spectrum as shown later.
The FokkerPlanck Eq. (B.1) describes the evolution of the electron distribution function; with an initially monoenergetic distribution in the form of a deltafunction
(n is the initial particle number density) the solution at subsequent times t takes the form of the logparabolic energy distribution assumed in Eq. (3) of the main text (see also A.1), reading
Here the time depending slope at is given by
(B.3) 
meanwhile, the curvature of driven by the diffusive (stochastic) term in Eq. (B.1), irreversibly decreases in time after
(B.4) 
from the large initial values corresponding to the initially monoenergetic distribution. Correspondingly, the timedependent height at follows
(B.5) 
To wit, Eq. (B.2) describes the evolution of the electron distribution, growing broader and broader under the effect of stochastic acceleration, while its peak moves from to the current position
(B.6) 
under the contrasting actions of the systematic and stochastic accelerations. An important quantity to focus on for the emission properties is the rms energy
=  (B.7) 
which is also the position for the peak of the distribution .
We have already derived in Appendix A the shapes of the spectra (synchrotron, and IC in both the Thomson and KN regimes) emitted by the distribution given in Eq. (3); here we stress the time dependence of their main spectral features. We can write for the synchrotron emission^{}
for IC emission we have in the Thomson regime
(B.9) 
and in the extreme KN regime
Note that during flares, since , we expect the curvature to vary so little that we can therefore approximatively write (see main text).
Now we focus the above relations for the simple case of time independent
and ,
when
;
then we have
while for the rms energy we have
so for the peak frequency and the spectral curvature we have
where is a normalization frequency. It is seen that soon after the injection the curvature (proportional to r) drops rapidly, then progressively decreases more and more gently, while the peak frequency still increases.
The value of
can be evaluated from observing the synchrotron spectral curvatures b_{2} and b_{1} at two times t_{2} and t_{1}, respectively (recall that
); denoting with t_{2}  t_{1} this time interval, we have
(B.14) 
on the other hand, form observing the related synchrotron peaks and , the value of can be evaluated as
(B.15) 
For example, in the case of Mrk 501 in the states of 7 and 16 April 1997, we obtain (on assuming ) and , corresponding to acceleration times and . Note that with the current data the evaluations of and turn out to be affected by uncertainties considerably larger than the single curvatures and .
If the total energy available to the jet is limited (e.g., by the BZ limit, see text) we expect that and cannot grow indefinitely, but are to attain a limiting value. At low energies where holds, we have and as before; at higher energies when reaches its limit, holds, leading to and . Eventually also reaches its limit, and both the fluxes and the peak frequencies cannot grow any more.
Appendix C: Prefactors for Eqs. (4)(13)
The HSZ SSC model is, as stated before, characterized by five parameters: the rms particle energy
,
the particle density n, the magnetic field B, the size of the emitting region R and the beaming factor .
So the model may be constrained by five observables that we denote with
(C.1) 
(C.2) 
where the indexes i,j,k,h express the normalizations demonstrated below; in addition, we denote with the time in days for the source variations, and with D the distance of the source in Gpc.
In the Thomson regime we find
B  =  
(C.3)  
=  13.5  
(C.4)  
R  =  
(C.5)  
n  =  
(C.6)  
=  (C.7) 
For the extreme KN regime we obtain
B  =  
(C.8)  
=  
(C.9)  
R  =  
(C.10)  
n  =  
(C.11)  
=  
(C.12) 
Footnotes
 ... tether^{}
 Appendices are only available in electronic form at http://www.aanda.org
 ...Donnarumma et al. 2009)^{}
 The logparabolic particle energy distribution could be made asymmetric by adding a powerlaw tail at low energies (see Massaro et al. 2006) to represent also low frequency data, e.g., radio in addition to optical/Xray data; however, here we will focus on the SEDs around their peaks and so we do not use such additions.
 ... density^{}
 Note that in Thomson regime no direct relation appears between electrons emitting peak synchrotron photons and those responsible of IC peak radiation (see Tavecchio et al. 1998).
 ... emission^{}
 Note that singular behaviors of S and C for t=0 are a consequence of the (differential) definition of the SED as , and relate to the initial singularity of the particle energy distribution. Integrated quantities like behave regularly.
All Tables
Table 1: Source parameters for some BL Lac sources.
Table 2: Spectral correlations (for ). We have denoted with or (B) the variations driven by increases of rms electron energy or magnetic field, respectively.
All Figures
Figure 1: Two examples of exact, steady SEDs provided by numerical simulations (code by Massaro 2007) after the HSZ SSC model from relativistic electrons with a logparabolic energy distribution (see Appendix B). The two characteristic humps correspond to synchrotron emission ( left) and to inverse Compton scattering ( right), either in the Thomson regime (blue line) or in the KN regime (red line). Standard sources that are so fitted include the low state of 0716+714 (an IBL, blue crosses) observations of November 1996 by BeppoSAX, AIT and EGRET (Tagliaferri et al. 2003); and 1ES 1553+113 (a HBL, red crosses) observations of April 2005 by Swift and AprilMay 2005 by MAGIC (Massaro 2007). 

Open with DEXTER  
In the text 
Figure 2: To illustrate the differences of the IC SEDs in the Thomson and in the extreme KN regime, we show the results of numerical simulations (code by Massaro 2007) after the HSZ SSC model, for different sources with larger and larger rms energy of the radiating particles and other parameters kept constant (including the curvature of the particle energy distribution). The vertical line indicates the frequency of the synchrotron peak where the IC scattering changes over the Thomson to the KN regime (see Eq. (14)). Actual sources in their evolution (represented in Figs. B.1 and B.2) vary their spectral curvature somewhat and only span a limited range in frequency and flux to the left (LBLs) or to the right (HBLs) of the line, with only the IBLs likely to cross it. 

Open with DEXTER  
In the text 
Figure 3: Computed changes of spectral parameters corresponding to extended increase of rms particle energy with the other source parameters kept constant. As functions of (where is the synchrotron SED peak frequency, see Eq. (14)), the full line represents the ratio , the dashed line the ratio and the dotdashed line the curvature ratio at the peaks. Circles represent the parameters of 0716+714 as observed on November 1996 and September 2007 (Tagliaferri et al. 2003; Chen et al. 2008), and diamonds represent Mrk 501 observations of April 1997 (DjannatiAtai et al. 1999; Massaro et al. 2006), each with the same colours as the lines. Changes of n only would yield at all frequencies, represented by the horizontal dotted line; added to a pure rms energy increase, a moderate increase of n canaccount for the upward deviations of observed values. 

Open with DEXTER  
In the text 
Figure 4: As in Fig. 2, we show SEDs obtained from numerical simulations (code by Massaro 2007) of HSZ SSC radiations, for different sources with larger and larger magnetic field B and other parameters kept constant. The vertical line again indicates the transition frequency from the Thomson to the KN regime (see Eq. (15)); as B increases, the IC scattering changes its regime, and reduced cross section yields similar compression effects as in Fig. 2. 

Open with DEXTER  
In the text 
Figure 5: Computed changes of spectral parameters as driven by an extended increase of the magnetic field with the other source parameters kept constant. The black full line represents the value of the ratio , the red dashed line represents the value of the ratio and the blue dotdashed line represents the ratio , for different values of (where is the synchrotron SED peak frequency, see Eq. (15)), data points are the same as in Fig. 3; compared with the latter considerably larger changes of n are needed to recover the observed values. 

Open with DEXTER  
In the text 
Figure 6: Spectral fit of two different states of the IBL source 0716+714. The low state (blue line) data are relative to 1996 November 14 simultaneous observation by AIT, BeppoSAX and EGRET; the high state (red line) data are relative to 2007 September 712 simultaneous observations by GASPWEBT and AGILEGRID (Tagliaferri et al. 2003; Chen et al. 2008). 

Open with DEXTER  
In the text 
Figure B.1: Example of time evolution of the distribution (with peak energy ) due to stochastic and systematic accelerations ( upper panel) with , from the initial value , and of the related evolution for the SSC SEDs ( lower panel). In terms of the stochastic acceleration time , the time interval between each pair of lines is , corresponding to an observed time interval of about two days. 

Open with DEXTER  
In the text 
Figure B.2: Example of time evolution of the distribution (with peak energy ) due to stochastic and systematic accelerations ( upper panel) with , from the initial value , and of the related evolution for the SSC SEDs ( lower panel). In terms of the stochastic acceleration time , the time interval between each pair of lines is , corresponding to an observed time interval of about two days. 

Open with DEXTER  
In the text 
Copyright ESO 2009