• ### Stochastic modeling of multiwavelength variability of the classical BL Lac object OJ 287 on timescales ranging from decades to hours(1709.04457)

July 10, 2018 astro-ph.HE
We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.
• ### Gamma-ray and Optical Oscillations of 0716+714, MRK 421, and BL Lac(1701.04454)

Jan. 29, 2017 astro-ph.HE
We examine the 2008-2016 $\gamma$-ray and optical light curves of three bright BL Lac objects, 0716+714, MRK 421, BL Lac, which exhibit large structured variability. We searched for periodicities by using a fully Bayesian approach. For two out of three sources investigated no significant periodic variability was found. In the case of BL Lac we detected a periodicity of ~ 680 days. Although the signal related to this is modest, the coincidence of the periods in both gamma and optical bands is indicative of a physical relevance. Considering previous literature results, possibly related $\gamma$-ray and optical periodicities of about one year time scale are proposed in 4 bright $\gamma$-ray blazars out of the 10 examined in detail. Comparing with results from periodicity search of optical archives of quasars, the presence of quasi-periodicities in blazars might be more frequent by a large factor. This suggests the intriguing possibility that the basic conditions for their observability are related to the relativistic jet in the observer direction, but the overall picture remains uncertain.
• ### Polarization angle swings in blazars: The case of 3C 279(1603.00249)

June 21, 2016 astro-ph.CO, astro-ph.HE
Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity.
• ### Primary black hole spin in OJ287 as determined by the General Relativity centenary flare(1603.04171)

March 14, 2016 astro-ph.HE
OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.
• We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.
• ### Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113(1509.02063)

Oct. 12, 2015 astro-ph.HE
We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.
• ### Analyzing polarization swings in 3C 279(1311.3126)

Nov. 15, 2013 astro-ph.CO, astro-ph.HE
Quasar 3C 279 is known to exhibit episodes of optical polarization angle rotation. We present new, well-sampled optical polarization data for 3C 279 and introduce a method to distinguish between random and deterministic electric vector position angle (EVPA) variations. We observe EVPA rotations in both directions with different amplitudes and find that the EVPA variation shows characteristics of both random and deterministic cases. Our analysis indicates that the EVPA variation is likely dominated by a random process in the low brightness state of the jet and by a deterministic process in the flaring state.
• ### Redshift constraints for RGB 0136+391 and PKS 0735+178 from deep optical imaging(1209.4755)

Sept. 21, 2012 astro-ph.CO
We present the results of deep I-band imaging of two BL Lacerate objects, RGB 0136+391 and PKS 0735+178, during an epoch when the optical nucleus was in a faint state in both targets. In PKS 0735+178 we find a significant excess over a point source, which, if fitted by the de Vaucouleurs model, corresponds to a galaxy with I = 18.64 +- 0.11 and r_eff = 1.8 +- 0.4 arcsec. Interpreting this galaxy as the host galaxy of PKS 0735+178 we derive z = 0.45 +- 0.06 using the host galaxy as a "standard candle". We also discuss the immediate optical environment of PKS 0735+178 and the identity of the MgII absorber at z = 0.424. Despite of the optimally chosen epoch and deep imaging we find the surface brightness profile of RGB 0136+391 to be consistent with a point source. By determining a lower limit for the host galaxy brightness by simulations, we derive z > 0.40 for this target.
• ### Variability of the blazar 4C 38.41 (B3 1633+382) from GHz frequencies to GeV energies(1207.3979)

July 17, 2012 astro-ph.CO, astro-ph.HE
The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and \gamma-ray flux apparently fades in time, likely because of an increasing optical to \gamma-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor \delta. Under the hypothesis that in the period 2008-2011 all the \gamma-ray and optical variability on a one-week timescale were due to changes in \delta, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle \theta\ only, then \theta\ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model.
• The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.
• ### Large scale environments of z<0.4 active galaxies(1109.1558)

Sept. 29, 2011 astro-ph.CO
Properties of galaxies depend on their large-scale environment. As the influence of active galactic nuclei (AGN) in galaxy evolution is becoming more evident, their large scale environments may help us understand the evolutionary processes leading to activity. The effect of activity can be seen particularly by showing if different types of active galaxies are formed by similar mechanisms. Our aim is to study the supercluster-scale environments of active galaxies up to redshift 0.4. Our data includes quasars, BL Lac objects, Seyfert and radio galaxies. We use a three-dimensional low-resolution luminosity-density field constructed of a sample of luminous red galaxies in the seventh data release of the Sloan Digital Sky Survey. We calculate the average density of this field in a volume of a 3\,$h^{-1}$Mpc sphere around each AGN for estimating the environmental density levels of different types of AGN. This analysis gives us the distribution of AGN in the global environment of superclusters, filaments, and voids. Our results show that while radio-quiet quasars and Seyfert galaxies are mostly located in low-density regions, radio galaxies have higher environmental densities. BL Lac objects usually have low-density environments, but some of them are also in very high density regions. Our results give support to the theory of galaxy evolution where galaxies are affected by two modes of AGN feedback: quasar mode that turns a blue star-forming galaxy into a red and dead one, and radio mode that regulates the growth of massive elliptical galaxies. We conclude that quasars are in lower density environments than radio galaxies most likely because the galaxies in rich environments have evolved to a state suitable for radio-loud activity earlier. Galaxies in poor environment have evolved slower, and are still going through the earlier quasar mode feedback in galaxy evolution.
• ### AGILE detection of extreme gamma-ray activity from the blazar PKS 1510-089 during March 2009. Multifrequency analysis(1103.3647)

March 18, 2011 astro-ph.HE
We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak.
• ### Another look at the BL Lacertae flux and spectral variability(1009.2604)

Sept. 14, 2010 astro-ph.CO
The GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) monitored BL Lacertae in 2008-2009 at radio, near-IR, and optical frequencies. During this period, high-energy observations were performed by XMM-Newton, Swift, and Fermi. We analyse these data with particular attention to the calibration of Swift UV data, and apply a helical jet model to interpret the source broad-band variability. The GASP-WEBT observations show an optical flare in 2008 February-March, and oscillations of several tenths of mag on a few-day time scale afterwards. The radio flux is only mildly variable. The UV data from both XMM-Newton and Swift seem to confirm a UV excess that is likely caused by thermal emission from the accretion disc. The X-ray data from XMM-Newton indicate a strongly concave spectrum, as well as moderate flux variability on an hour time scale. The Swift X-ray data reveal fast (interday) flux changes, not correlated with those observed at lower energies. We compare the spectral energy distribution (SED) corresponding to the 2008 low-brightness state, which was characterised by a synchrotron dominance, to the 1997 outburst state, where the inverse-Compton emission was prevailing. A fit with an inhomogeneous helical jet model suggests that two synchrotron components are at work with their self inverse-Compton emission. Most likely, they represent the radiation from two distinct emitting regions in the jet. We show that the difference between the source SEDs in 2008 and 1997 can be explained in terms of pure geometrical variations. The outburst state occurred when the jet-emitting regions were better aligned with the line of sight, producing an increase of the Doppler beaming factor. Our analysis demonstrates that the jet geometry can play an extremely important role in the BL Lacertae flux and spectral variability.
• ### The December 2009 gamma-ray flare of 3C 454.3: the multifrequency campaign(1005.3263)

May 18, 2010 astro-ph.CO
During the month of December, 2009 the blazar 3C 454.3 became the brightest gamma-ray source in the sky, reaching a peak flux F ~2000E-8 ph/cm2/s for E > 100 MeV. Starting in November, 2009 intensive multifrequency campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results of a 2-month campaign involving AGILE, INTEGRAL, Swift/XRT, Swift/BAT, RossiXTE for the high-energy observations, and Swift/UVOT, KANATA, GRT, REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and additional optical data. We detected a long-term active emission phase lasting ~1 month at all wavelengths: in the gamma-ray band, peak emission was reached on December 2-3, 2009. Remarkably, this gamma-ray super-flare was not accompanied by correspondingly intense emission in the optical/UV band that reached a level substantially lower than the previous observations in 2007-2008. The lack of strong simultaneous optical brightening during the super-flare and the determination of the broad-band spectral evolution severely constrain the theoretical modelling. We find that the pre- and post-flare broad-band behavior can be explained by a one-zone model involving SSC plus external Compton emission from an accretion disk and a broad-line region. However, the spectra of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 2009 cannot be satisfactorily modelled by a simple one-zone model. An additional particle component is most likely active during these states.
• ### H-alpha monitoring of OJ 287 in 2005-08(1004.2617)

April 15, 2010 astro-ph.CO
• ### Host galaxy subtraction of TeV candidate BL Lacertae objects(0709.2533)

Oct. 10, 2007 astro-ph
We have derived R-band host galaxy fluxes of 20 TeV candidate BL Lacertae objects as a function of aperture radius and FWHM. The results are given as correction tables, that list the fluxes (in mJy) of all contaminating'' sources (host galaxy + significant nearby objects) as a function of aperture radius and FWHM. We found that the derived fluxes depend strongly on aperture radius, but the FWHM has only a minor effect (a few percent). We also discuss the implications of our findings to optical monitoring programs and potential sources of error in our derived fluxes. During this work we have also constructed new calibration star sequences for 9 objects and present the finding charts and calibrated magnitudes.
• ### Day-Scale Variability of 3C 279 and Searches for Correlations in Gamma-Ray, X-Ray, and Optical Bands(astro-ph/0105247)

May 15, 2001 astro-ph
Light curves of 3C 279 are presented in optical (R-band), X-rays (RXTE/PCA), and gamma rays (CGRO/EGRET) for 1999 Jan-Feb and 2000 Jan-Mar. During both of those epochs the gamma-ray levels were high, and all three observed bands demonstrated substantial variation, on time scales as short as one day. Correlation analyses provided no consistent pattern, although a rather significant optical/gamma-ray correlation was seen in 1999, with a gamma-ray lag of ~2.5 days, and there are other suggestions of correlations in the light curves. For comparison, correlation analysis is also presented for the gamma-ray and X-ray light curves during the large gamma ray flare in 1996 Feb and the two gamma-bright weeks leading up to it; the correlation at that time was strong, with a gamma-ray/X-ray offset of no more than 1 day.
• ### Broadband Spectral Analysis of PKS 0528+134: A Report on Six Years of EGRET Observations(astro-ph/9901106)

June 25, 1999 astro-ph
The multiwavelength spectra of PKS 0528+134 during six years of observations by EGRET have been analyzed using synchrotron self-Compton (SSC) and external radiation Compton (ERC) models. We find that a two-component model, in which the target photons are produced externally to the gamma-ray emitting region, but also including an SSC component, is required to suitably reproduce the spectral energy distributions of the source. Our analysis indicates that there is a trend in the observed properties of PKS 0528+134, as the source goes from a gamma-ray low state to a flaring state. We observe that during the higher gamma-ray states, the bulk Lorentz factor of the jet increases and the ERC component dominates the high-energy emission. Our model calculations indicate the trend that the energies of the electrons giving rise to the synchrotron peak decreases, and the power-ratio of the gamma-ray and low energy spectral components increases, as the source goes from a low to a high gamma-ray state.