Blazars are characterized by largely aperiodic variability on timescales ranging from minutes to decades across the electromagnetic spectrum. The Transiting Exoplanet Survey Satellite (TESS) mission provides continuous sampling of blazar variability on timescales ranging from tens of minutes to 27 days for a single sector observation. Proper removal of the background, thermal ramping, and onboard systematic effects are crucial to the extraction of a reliable blazar light curve. Multiple publicly available procedures have been created to correct for these effects. Using ground-based observations from the Zwicky Transit Facility (ZTF) and the Asteroid Terrestrial-impact Last Alert System (ATLAS) as “ground truth” observations, we compare six different methods (regression, cotrending basis vectors, pixel-level decorrelation, eleanor , quaver , and simple differential photometry) to each other, and to our “ground truth” observations, to identify which methods properly correct light curves of a sample of 11 bright blazars, including the suspected neutrino source TXS 0506+056. In addition to comparing the resulting light curves, we compare the slopes of the power spectral densities, perform least-square fitting to simultaneous ZTF data, and quantify other statistical qualities generated from the light curves of each method. We find that only three of the six methods compared (simple differential photometry, eleanor , and quaver ) produce TESS light curves consistent with the ground-based ZTF and ATLAS observations.

We present second observations by K2 of OJ 287 and seven other γ -ray active galactic nuclei (AGNs) obtained in 2017–2018, second and third observations of the lobe-dominated, steep-spectrum quasar 3C 207, and observations of nine additional blazars not previously observed with K2. The AGNs were observed simultaneously with K2 and the Fermi Large Area Telescope for 51–81 days. Our full sample, observed in 2014–2018, contained 16 BL Lacertae objects (BL Lacs), nine flat-spectrum radio quasars (FSRQs), and four other γ -ray AGNs. Twelve BL Lacs and seven FSRQs exhibited fast, jagged light curves, while four BL Lacs and two FSRQs had slow, smooth light curves. Some objects changed their red-noise character significantly between repeated K2 observations. The optical characteristics of OJ 287 derived from the short-cadence K2 light curves changed between observations made before and after the predicted passage of the suspected secondary supermassive black hole through the accretion disk of the primary supermassive black hole. The average slopes of the periodogram power spectral densities of the BL Lacs’ and FSRQs’ light curves differed significantly, by ≈12%, with the BL Lac slopes being steeper, and a Kolmogorov–Smirnov test with a p -value of 0.039 indicates that these samples probably come from different populations; however, this result is not as strongly supported by power spectral response (PSRESP) analyses. Differences in the origin of the jets from the ergosphere or accretion disk in these two classes could produce such a disparity, as could different sizes or locations of emission regions within the jets.

We characterize the variability in nearly continuous optical observations of the bright radio-loud quasar 3C 273 and nine additional active galactic nuclei detected with the Fermi Gamma-ray Space Telescope Large Area Telescope (Fermi-LAT). Optical observations were obtained during the K2 mission with the Kepler spacecraft for periods of 49 to 83 days conducted with ≃1 minute (short) or ≃30 minutes (long) cadences in 2015–2017. 3C 273 was quiescent during the course of the observations, varying by only a factor of 1.02. Three objects, PKS 0047+023, PKS 1216-10, and PKS B2320-035, were active, varying by factors of 1.8–3.4. Six other objects were comparatively quiet, varying by factors of less than 1.4. Power spectral densities (PSDs) were calculated for each object. Overall, the slopes of most PSDs, as well as those we reported in an earlier paper, were in the range −2.0 to −2.7 and are consistent with those produced by turbulence in the relativistic jet, and not by “hot spots” in the disk emission. Mechanisms operating in the jet other than turbulence, such as “mini-jets” or “jet-in-jets”, may also produce the observed range of PSD slopes. Both accretion disk and jet models are plausible origins for the 3C 273 optical variability during the K2 observations.

We characterize the variability in nearly continuous optical observations of nine active galactic nuclei (AGNs) detected with the Fermi Gamma Ray Space Telescope Large Area Telescope (Fermi-LAT), obtained during the K2 mission with the Kepler spacecraft with 1-minute or 30-minute cadences. The binary black hole candidate OJ 287, six other γ-ray blazars, the lobe-dominated radio quasar 3C 207, and the X-ray-bright quasar RBS 1273 were observed. The most dramatic variability was seen in PKS 1352−104, which varied by a factor of six in 4 days. OJ 287 varied by a factor of 2.8, six other blazars varied by factors of 1.14–2.23, and RBS 1273 varied by a factor of only 1.03. Power spectral densities (PSDs) were calculated for each object. Six of the nine AGNs were detected by the Fermi-LAT during the K2 observations. The slopes of the PSDs for these six contemporaneously detected AGNs were not significantly different from the slopes of the three contemporaneously undetected AGNs. The PSD slopes of the three BL Lac objects and the four flat-spectrum radio quasars were not significantly different, consistent with the relativistic jet acting as the dominant source of the emission in both subclasses of blazars. Overall, the slopes of most PSDs are consistent with turbulence in the relativistic jet, but the steepest PSDs require the addition of another mechanism that produces additional flux variations on timescales of several days during major flaring episodes.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the 1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

ABSTRACT BL Lac, the eponymous blazar, flared to historically high levels at millimeter, infrared, X-ray, and gamma-ray wavelengths in 2012. We present observations made with Herschel, Swift, NuSTAR, Fermi , the Submillimeter Array, CARMA, and the VLBA in 2012–2013, including three months with nearly daily sampling at several wavebands. We have also conducted an intensive campaign of 30 hr with every-orbit observations by Swift and NuSTAR , accompanied by Herschel , and Fermi observations. The source was highly variable at all bands. Time lags, correlations between bands, and the changing shapes of the spectral energy distributions can be explained by synchrotron radiation and inverse Compton emission from nonthermal seed photons originating from within the jet. The passage of four new superluminal very long baseline interferometry knots through the core and two stationary knots about 4 pc downstream accompanied the high flaring in 2012–2013. The seed photons for inverse Compton scattering may arise from the stationary knots and from a Mach disk near the core where relatively slow-moving plasma generates intense nonthermal radiation. The 95 spectral energy distributions obtained on consecutive days form the most densely sampled, broad wavelength coverage for any blazar. The observed spectral energy distributions and multi-waveband light curves are similar to simulated spectral energy distributions and light curves generated with a model in which turbulent plasma crosses a conical shock with a Mach disk.

Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The radio spectra peak at high frequencies and exhibit a variety of shapes. For a small set of low-z sources, we find a spectral upturn at high frequencies, indicating the presence of intrinsic cold dust. Variability can generally be approximated by achromatic variations, while sources with clear signatures of evolving shocks appear to be limited to the strongest outbursts.

Context. Because of its proximity, Mrk 421 is one of the best sources on which to study the nature of BL Lac objects. Its proximity allows us to characterize its broadband spectral energy distribution (SED). Aims. The goal is to better understand the mechanisms responsible for the broadband emission and the temporal evolution of Mrk 421. These mechanisms may also apply to more distant blazars that cannot be studied with the same level of detail. Methods. A flare occurring in March 2010 was observed for 13 consecutive days (from MJD 55 265 to MJD 55 277) with unprecedented wavelength coverage from radio to very high energy (VHE; E> 100 GeV) γ-rays with MAGIC, VERITAS, Whipple, Fermi-LAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We modeled the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigated the physical parameters, and evaluated whether the observed broadband SED variability can be associated with variations in the relativistic particle population. Results. The activity of Mrk 421 initially was high and then slowly decreased during the 13-day period. The flux variability was remarkable at the X-ray and VHE bands, but it was minor or not significant at the other bands. The variability in optical polarization was also minor. These observations revealed an almost linear correlation between the X-ray flux at the 2–10 keV band and the VHE γ-ray flux above 200 GeV, consistent with the γ-rays being produced by inverse-Compton scattering in the Klein-Nishina regime in the framework of SSC models. The one-zone SSC model can describe the SED of each day for the 13 consecutive days reasonably well, which once more shows the success of this standard theoretical scenario to describe the SEDs of VHE BL Lacs such as Mrk 421. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission, while the other smaller zone, which is spatially separated from the first, contributes to the daily variable emission occurring at X-rays and VHE γ-rays. The second blob is assumed to have a smaller volume and a narrow electron energy distribution with 3 × 10^4<γ< 6 × 10^5, where γ is the Lorentz factor of the electrons. Such a two-zone scenario would naturally lead to the correlated variability at the X-ray and VHE bands without variability at the optical/UV band, as well as to shorter timescales for the variability at the X-ray and VHE bands with respect to the variability at the other bands. Conclusions. Both the one-zone and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly with the underlying particle population. This shows that the particle acceleration and cooling mechanism that produces the radiating particles might be the main mechanism responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement with the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by varying the parameters related to the emitting region itself (δ, B and R), in addition to the parameters related to the particle population.

We present an analysis of the multiwavelength behaviour of the blazar OJ 248 at z = 0.939 in the period 2006-2013. We use low-energy data (optical, near-infrared, and radio) obtained by 21 observatories participating in the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), as well as data from the Swift (optical-UV and X-rays) and Fermi (gamma-rays) satellites, to study flux and spectral variability and correlations among emissions in different bands. We take into account the effect of absorption by the Damped Lyman Alpha intervening system at z = 0.525. Two major outbursts were observed in 2006-2007 and in 2012-2013 at optical and near-IR wavelengths, while in the high-frequency radio light curves prominent radio outbursts are visible peaking at the end of 2010 and beginning of 2013, revealing a complex radio-optical correlation. Cross-correlation analysis suggests a delay of the optical variations after the gamma-ray ones of about a month, which is a peculiar behaviour in blazars. We also analyse optical polarimetric and spectroscopic data. The average polarization percentage P is less than 3 per cent, but it reaches about 19 per cent during the early stage of the 2012-2013 outburst. A vague correlation of P with brightness is observed. There is no preferred electric vector polarisation angle and during the outburst the linear polarization vector shows wide rotations in both directions, suggesting a complex behaviour or structure of the jet and possible turbulence. The analysis of 140 optical spectra acquired at the Steward Observatory reveals a strong Mg II broad emission line with an essentially stable flux of 6.2 e-15 erg cm-2 s-1 and a full width at half-maximum of 2053 km s-1.

We present an analysis of the multiwavelength behaviour of the blazar OJ 248 at z = 0.939 in the period 2006–2013. We use low-energy data (optical, near-infrared, and radio) obtained by 21 observatories participating in the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), as well as data from the Swift (optical–UV and X-rays) and Fermi (γ-rays) satellites, to study flux and spectral variability and correlations among emissions in different bands. We take into account the effect of absorption by the damped Lyman α intervening system at z = 0.525. Two major outbursts were observed in 2006–2007 and in 2012–2013 at optical and near-IR wavelengths, while in the high-frequency radio light curves prominent radio outbursts are visible peaking at the end of 2010 and beginning of 2013, revealing a complex radio-optical correlation. Cross-correlation analysis suggests a delay of the optical variations after the γ-ray ones of about a month, which is a peculiar behaviour in blazars. We also analyse optical polarimetric and spectroscopic data. The average polarization percentage P is less than 3%, but it reaches ∼ 19% during the early stage of the 2012–2013 outburst. A vague correlation of P with brightness is observed. There is no preferred electric vector polarisation angle and during the outburst the linear polarization vector shows wide rotations in both directions, suggesting a complex behaviour/structure of the jet and possible turbulence. The analysis of 140 optical spectra acquired at the Steward Observatory reveals a strong Mg II broad emission line with an essentially stable flux of 6.2 × 10 erg cm s and a full width at half-maximum of 2053 km s.

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

We used photometry from the Kepler satellite to characterize the variability of four radio-loud active galactic nuclei (AGNs) on timescales from years to minutes. The Kepler satellite produced nearly continuous high precision data sets which provided better temporal coverage than possible with ground based observations. We have now accumulated 11 quarters of data, eight of which were reported in our previous paper. In addition to constructing power spectral densities (PSDs) and characterizing the variability of the last three quarters, we have linked together the individual quarters using a multiplicative scaling process, providing data sets spanning ∼2.8 yr with >98% coverage at a 30 minute sampling rate. We compute PSDs on these connected data sets that yield power law slopes at low frequencies in the approximate range of −1.5 to −2.0, with white noise seen at higher frequencies. These PSDs are similar to those of both the individual quarters and to those of ground-based optical observations of other AGNs. We also have explored a PSD binning method intended to reduce a bias toward shallow slope fits by evenly distributing the points within the PSDs. This tends to steepen the computed PSD slopes, especially when the low frequencies are relatively poorly fit. We detected flares lasting several days in which the brightness increased by ∼15%–20% in one object, as well a smaller flare in another. Two AGNs showed only small, ∼1%–2%, fluctuations in brightness.

We have used Kepler photometry to characterize variability in four radio-loud active galactic nuclei (AGN; three quasars and one object tentatively identified as a Seyfert 1.5 galaxy) on timescales from minutes to months, comparable to the light crossing time of the accretion disk around the central supermassive black hole or the base of the relativistic jet. Kepler's almost continuous observations provide much better temporal coverage than is possible from ground-based observations. We report the first such data analyzed for quasars. We have constructed power spectral densities using eight Kepler quarters of long-cadence (30-minute) data for three AGN, six quarters for one AGN and two quarters of short-cadence (1-minute) data for all four AGN. On timescales longer than about 0.2-0.6 days, we find red noise with mean power-law slopes ranging from -1.8 to -1.2, consistent with the variability originating in turbulence either behind a shock or within an accretion disk. Each AGN has a range of red noise slopes which vary slightly by month and quarter of observation. No quasi-periodic oscillations of astrophysical origin were detected. We detected flares of several days long when brightness increased by 3%-7% in two objects. No flares on timescales of minutes to hours were detected. Our observations imply that the duty cycle for enhanced activity in these radio-loud AGN is small. These well-sampled AGN light curves provide an impetus to develop more detailed models of turbulence in jets and instabilities in accretion disks.

We analyze the multi-frequency behavior of the quasar 3C 454.3 during three prominent gamma-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each gamma-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the gamma-ray outbursts, although the gamma-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the gamma-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the gamma-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that gamma-ray outbursts with variability timescales as short as similar to 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells.

Since the launch of the Fermi satellite, BL Lacertae has been moderately active at gamma-rays and optical frequencies until May 2011, when the source started a series of strong flares. The exceptional optical sampling achieved by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) in collaboration with the Steward Observatory allows us to perform a detailed comparison with the daily gamma-ray observations by Fermi. Discrete correlation analysis between the optical and gamma-ray emission reveals correlation with a time lag of 0 +- 1 d, which suggests cospatiality of the corresponding jet emitting regions. A better definition of the time lag is hindered by the daily gaps in the sampling of the extremely fast flux variations. In general, optical flares present more structure and develop on longer time scales than corresponding gamma-ray flares. Observations at X-rays and at millimetre wavelengths reveal a common trend, which suggests that the region producing the mm and X-ray radiation is located downstream from the optical and gamma-ray-emitting zone in the jet. The mean optical degree of polarisation slightly decreases over the considered period and in general it is higher when the flux is lower. The optical electric vector polarisation angle (EVPA) shows a preferred orientation of about 15 deg, nearly aligned with the radio core EVPA and mean jet direction. Oscillations around it increase during the 2011-2012 outburst. We investigate the effects of a geometrical interpretation of the long-term flux variability on the polarisation. A helical magnetic field model predicts an evolution of the mean polarisation that is in reasonable agreement with the observations. These can be fully explained by introducing slight variations in the compression factor in a transverse shock waves model.

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad nonthermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV γ-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed forMkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of ≲ 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

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

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164 in 2008. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with ~0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary core and in the superluminal knot, both at >12 parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.