© ESO, 2013

1. Introduction

The luminous high-redshift (z = 2.218; Stickel & Kühr 1993) quasar of the blazar sub-class 0836+710 (also known as S5 0836+71 or 4C +71.07) is characterized by a flat radio spectrum (α = −0.33; Kühr 1981). It hosts a powerful one-sided radio jet emerging from the core and extending up to kiloparsec scales (Hummel et al. 1992). Very Long Baseline Interferometry (VLBI) images of the source show a complex motion pattern, with jet components moving from apparent subluminal to superluminal velocities (Otterbein et al. 1998). Monitoring with Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE)¹ has shown apparent jet speeds up to β_app = 25 Lister et al. (2009). Internal structure of the jet in 0836+710 has been investigated at 1.6 and 5 GHz using observations with the VLBI Space Observatory Programme (VSOP; Lobanov et al. 1998), suggesting a helical structure for the jet. A more recent study of this source with VLBI data confirmed that the helical structures observed in the jet of 0836+710 are real and not generated artificially by the observing arrays Perucho et al. (2012). VLBA observations also suggest a spine-sheath structure for the jet of this source (Asada et al. 2010).

Blazar 0836+710 has been the subject of several X-ray studies and multiwavelength modeling to understand how X-ray emission is produced and the relationship between X-ray emission and other bands (Fang et al. 2001; Foschini et al. 2006; Sambruna et al. 2007; Gianní et al. 2011). These studies conclude that 0836+710 is a typical blazar, ultimately powered by a black hole with mass 3 × 10⁹  M_⊙ (Ghisellini et al. 2010), requiring an external source of seed photons for the high-energy Compton part of the spectrum.

In the 1990’s this blazar was detected and shown to be variable at >100 MeV γ-ray energies by EGRET on board the Compton Gamma Ray Observatory (CGRO; Thompson et al. 1993), with the name 3EG J0845+7049 in the Third EGRET Catalog (Hartman et al. 1999). Data from other CGRO instruments (BATSE, OSSE, and COMPTEL) at lower γ-ray energies showed that the peak of the spectral energy distribution (SED) falls in the MeV energy range (Malizia et al. 2000; Collmar 2006).

A new era in γ-ray astrophysics began with the launch of the Fermi Gamma-ray Space Telescope (Fermi) in 2008 June. This source was not bright enough in the first three months of normal mission operations to be included in the Fermi Large Area Telescope (LAT) Bright Source List (Abdo et al. 2009). It was, however, associated with 1FGL J0842.2+7054 in the First LAT Catalog (1FGL; Abdo et al. 2010a) and with 2FGL J0841.6+7052 in the Second LAT Catalog (2FGL; Nolan et al. 2012). Regular γ-ray monitoring by Fermi-LAT showed that the source was not active until mid-2011. An outburst from the source on 2011 April 3 was noted by the Fermi-LAT (Ciprini 2011) with an average daily flux nearly a factor of 20 greater than the average flux in the 2FGL catalog. Its γ-ray spectrum is steep, with a photon power-law index of 2.95 ± 0.07 in 2FGL. Even larger flares were seen in late 2011 (Ciprini & Dutka 2011), when follow-up observations with Swift were performed (D’Ammando et al. 2011). The active period for 0836+710 in γ rays ended after 2012 January.

Motivated by the γ-ray flaring activity in 2011, we have assembled a long-term multiwavelength study of this object, with the goal of comparing its temporal and spectral properties during quiescent and flaring intervals and determining whether this high-redshift flat spectrum radio quasar (FSRQ) differs significantly from more nearby γ-ray blazars. Although this source is monitored regularly by radio telescopes and the Fermi-LAT, its coverage at other wavelengths is limited. In this paper, we present the temporal and spectral evolution of 0836+710 as observed by Fermi-LAT and contemporaneous observations at other wavelengths since the time of the Fermi launch. In Sect. 2, observations and data analysis are presented. Results are summarized in Sect. 3, and discussion is given in Sect. 4. Throughout the paper we use cosmological parameters H₀ = 71 km s^-1 Mpc^-1, Ω_m = 0.27, and Ω_Λ = 0.73, consistent with the WMAP results (Komatsu et al. 2009).

2. Observations

2.1. Fermi-LAT analysis

The Fermi-LAT is a γ-ray telescope operating from 20 MeV to >300 GeV. The instrument is an array of 4 × 4 identical towers, each one consisting of a tracker (where the photon interacts by pair conversion and the resulting particles are tracked) and a calorimeter (where the energy of the pair-converted photon is measured). The entire instrument is covered by an anticoincidence detector to help in rejecting the charged-particle background. Further details on the description of LAT are given by Atwood et al. (2009). The LAT normally operates in sky-survey mode in which the whole sky is observed every 2 orbits (~3 h).

The Fermi-LAT data (100 MeV−300 GeV) considered for this analysis cover the period from 2008 August 4 (MJD 54 682) to 2012 January 31 (MJD 55 957). We restricted the analysis to a region of interest centered on the source and having a radius of 10°. Only events belonging to the Source class were used. To limit contamination from the γ-ray-bright Earth limb, only events with zenith angle <100° were selected. The data analysis was performed with the standard binned maximum-likelihood analysis tool gtlike, and the corresponding P7_V6 instrument response functions (IRF) were used. The source model includes point sources from the 2FGL catalog (Nolan et al. 2012) within 20° of 0836+710. The most important nearby source for this analysis is 2FGL J0721.9+7120, associated with the bright, variable blazar S5 0716+714. The background model also includes a component for the Galactic diffuse emission (gal_2yearp7v6_v0.fits), and an isotropic component (iso_p7v6source.txt) that represents the extragalactic diffuse emission as well as residual background from cosmic-ray particles mis-identified as γ rays. These resources are provided separately and with the Fermi-LAT Science Tools package².

2.1.1. Temporal behavior

We divided the observations into two time periods: a quiescent period from 2008 August through 2011 February, and an active period from 2011 March through 2012 January. Light curves were extracted using likelihood analysis, fitting the source spectrum by a power-law spectral model for these two periods³. The flux was evaluated by integrating the fitted model from 100 MeV−300 GeV. A detection was defined by a Test Statistic TS > 4, where TS = 2(log L₁ − log L₀), where L is the likelihood of the data given the model with (L₁) or without a source (L₀) in the model Mattox et al. (1996). Although the scanning mode of Fermi produces essentially continuous observations, we have chosen to omit time intervals that have only upper limits from the plots for clarity. The weekly light curve for the quiescent period, calculated with a spectral index fixed to the value in the 2FGL catalog, is presented in Fig. 1. Error bars shown are statistical. In addition, all results have systematic uncertainties of 10% at 100 MeV, 5% at 560 MeV, 10% at 10 GeV and above (Ackermann et al. 2012). For most time periods 0836+710 was too faint in the LAT data for detailed studies on short time scales.

Fig. 1 Fermi-LAT light curve for 0836+710 during the quiescent period, MJD 54 682−55 621 (2008 August 4−2011 March 1). The source flux is shown with 7-day time binning for E > 100 MeV. Points with TS < 4 are omitted from this figure as well as from Figs. 2 and 3.

In Fig. 2, the weekly light curve for the active period is presented, again with the same fixed spectral index as in Fig. 1. Within this active period, the interval from 2011 October 25 (MJD 55 859) to 2011 December 25 (MJD 55 920) showed a significant flare, noted in the figure. Figure 3 shows the E > 100 MeV flux with one-day binning during the active period, again with the same fixed spectral index as in Fig. 1. The flare interval shows evidence of sub-flares within the time period of enhanced emission.

Fig. 2 Fermi-LAT light curve for 0836+710 during the active period, MJD 55 621−55 957 (2011 March 1−2012 January 31). The source flux is shown with 7-day time binning for E > 100 MeV.

Fig. 3 Fermi-LAT light curve for 0836+710 during the active period, MJD 55 621−55957 (2011 March 1−2012 January 31). The source flux is shown with 1-day time binning for E > 100 MeV.

Although its high redshift and steep energy spectrum might suggest that 0836+710 would produce too few γ rays to enable detailed variability studies, its 2011 episode (c) flaring activity was as bright as that of many lower-redshift blazars. In particular we note the following:

• 1.

  Even during the quiescent state, the blazar is regularly detectedon one-week time scales, and almost always on one-month timescales (see Sect. 3). There are few time intervals inwhich the γ-ray emission is “off”. The particle acceleration and interaction processes that produce γ rays must therefore be operating most of the time, not just during flaring activity.

• 2.

  As shown in Fig. 1, there is some indication of variability during this relatively quiet period.

• 3.

  Like many of the blazars seen by the LAT (Abdo et al. 2010b), variability is seen on a variety of time scales, with much of the variation taking place on times extending over weeks or months rather than for shorter durations.

• 4.

  In addition to the general variation, short flares are seen in 0836+710, as seen in Fig. 3. On 2011 November 1 (MJD 55 866) 0836+710 flared for one day in γ rays, reaching a flux F(>100 MeV) = (3.0 ± 0.3) × 10^-6 ph cm^-2 s^-1, nearly 40 times brighter than the average flux in the 2FGL catalog. The rise was a factor of 8 in a week and a factor of 4 in 24 h or less. A similar one-day spike was seen six days later. To estimate the rise times of the flare, we tried smaller time bins (1 day, 12 h, and 6 h) and used an exponential fit (Ackermann et al. 2010), producing an e-folding rise time of ~0.3 days. The counting statistics limit the accuracy of the fit. These flares are comparable in rise time to the fastest flaring activity observed by the LAT for other blazars such as 3C 454.3 (Ackermann et al. 2010) and PKS 1502+106 (Abdo et al. 2010c). Additional flares lasting several days appear around MJD 55 885 and MJD 55 920. Because the flaring time scale sets a limit on the size of the emitting region due to light-crossing constraints, the different flaring times seen for 0836+710 suggest that the γ rays are produced either in one small region with varying external inputs (energetic particles or seed photons) or in emitting regions of different sizes.

The peak daily flux for E > 100 MeV represents an apparent isotropic luminosity of 8 × 10⁴⁷ erg s^-1, using a luminosity distance of 18 Gpc. This luminosity, while not the most extreme seen among γ-ray blazars, puts 0836+710 among the most powerful such objects.

Fig. 4 Energy spectrum of 0836+710 during the quiescent period. The red shaded region is the likelihood power-law fit to the entire energy range, showing the statistical uncertainty. The blue shaded region is the likelihood log-parabola fit to the entire energy range, showing the statistical uncertainty. The reference, or pivot, energy for each of the shaded regions is 0.32 GeV.

2.1.2. Gamma-ray spectra

We performed the spectral analysis separately for the quiescent and active periods. In addition to the likelihood fits for the full energy range, we computed the flux values for smaller energy intervals independently. This approach provides flux measurements in each band that can be plotted along with the maximum likelihood model fit to the whole energy range, providing a sense of how well that model describes the data.

As a first step, we analyzed the γ-ray spectral shape of 0836+710 during the quiescent and active periods using a power-law spectral model, as this was the functional form reported in the 2FGL catalog. In the analysis for the quiescent state, the power-law index Γ for the spectrum is (2.91 ± 0.07) and the integral flux F_{E    >    100MeV} = (5.22 ± 0.43) × 10^-8 ph cm^-2s^-1. The pivot energy, for which the correlation between the photon flux and the spectral index is minimal, is 0.32 GeV. Following the approach of the 2FGL catalog, we also tested a log-parabola spectral model $\log \left(\frac{\mathrm{d}\mathit{N}}{\mathrm{d}\mathit{E}}\right)\mathrm{=}\log \mathrm{\left(}\mathit{K}\mathrm{\right)}\mathrm{-}\mathit{\alpha }\log \left(\frac{\mathit{E}}{{\mathit{E}}_{\mathrm{0}}}\right)\mathrm{-}\mathit{\beta }{\log }^{\mathrm{2}}\left(\frac{\mathit{E}}{{\mathit{E}}_{\mathrm{0}}}\right)\mathrm{·}$(1)The parameters are K (normalization), α (spectral slope at E₀) and the curvature β, with E₀ an arbitrary reference energy, chosen to be the same value used in the 2FGL catalog, E₀ = 0.32 GeV. To compare the two spectral fits, we calculated TS_Curv as twice the difference of the log(likelihood) for the power law and the log parabola. As shown in Fig. 4, the log-parabola form is compatible within uncertainties with a spectrum with no curvature, β = (0.07 ± 0.06), and provides no better fit than the power-law model, with TS_Curv = 1.3.

For the active time period (2011 March through 2012 January), the fit results in Γ = (2.75 ± 0.03) for a power-law model and gives a flux value F_E > 100MeV = (2.59 ± 0.09) × 10^-7 ph cm^-2s^-1. Figure 5 shows the energy spectrum of 0836+710 during this interval. The γ-ray energy spectrum exhibited a significant deviation from a power law, as shown in the figure, with TS_Curv = 43.8. The log-parabola spectral parameters are α = (2.62 ± 0.04) and β = (0.24 ± 0.04).

We also performed a spectral analysis for the bright γ-ray flare (see Figs. 1 or 2) separately. For this flare state, the power-law fit gives Γ = (2.65 ± 0.04) and flux value F_{E    >    100MeV} = (6.05    ±    0.20) × 10^-7 ph cm^-2 s^-1. During this flare, the γ-ray energy spectrum also exhibited a significant deviation from a power law, as shown in Fig. 6, with TS_Curv = 58.9. The log-parabola spectral parameters are α = (2.50 ± 0.05) and β = (0.32 ± 0.05).

The power-law representation of the spectrum in the lower-flux states suggests that the LAT energy range lies well above the energy of the peak in the SED, while one possibility for the appearance of curvature in the flaring spectrum might be that the SED peak has shifted to higher energies (although still below the LAT energy range, making any explanation for the appearance of curvature somewhat speculative).

Fig. 5 Energy spectrum of 0836+710 during the full active period (2011 March through 2012 January). The red shaded region is the likelihood power-law fit to the entire energy range, showing the statistical uncertainty. The blue shaded region is the likelihood log-parabola fit to the entire energy range, showing the statistical uncertainty. The reference, or pivot, energy for each of the shaded regions is 0.32 GeV.

Fig. 6 Energy spectrum of 0836+710 during the flare (see Figs. 1 or 2). The red shaded region is the likelihood power-law fit to the entire energy range, showing the statistical uncertainty. The blue shaded region is the likelihood log-parabola fit to the entire energy range, showing the statistical uncertainty. The reference, or pivot, energy for each of the shaded regions is 0.32 GeV.

Some Fermi-LAT blazar spectra have shown a trend toward a “harder-when-brighter” pattern (Abdo et al. 2010d), often with only a moderate spectral evolution, such as 3C 273 (Abdo et al. 2010e) or PKS 1510−089 (Abdo et al. 2010f), but in some cases dramatically, such as GB 1310+487 (Sokolovsky et al. 2013). There is an indication of a “harder-when-brighter” pattern seen for the observations reported here, although the primary feature is the appearance of curvature in the spectrum during the active phase.

2.2. Swift analysis

2.2.1. XRT

The Swift X-Ray Telescope (XRT; Burrows et al. 2005) has regularly provided X-ray coverage of blazars. We ran the “Swift-XRT data products generator” (Evans et al. 2009), available at the University of Leicester website⁴, to create 3 light curves, binned by observation: one in the 0.3−10 keV range and two in the 0.3−1.5 keV (soft) and 1.5−10 keV (hard) sub-ranges. The Swift archival data cover a long period of time, from 2006 April 1 to 2012 April 3, and the light curves show significant variations. Indeed, the 0.3−10 keV count rate (corrected for dead or hot pixels, the finite aperture of the extraction region, and possible pile-up) varies by a factor of 3.5 (from minimum (0.27    ±    0.09) ct s^-1 to maximum (0.90 ± 0.03) ct s^-1), suggesting we are observing a flaring X-ray source. This activity is not coupled with a spectral evolution. The uniformity of the spectrum is confirmed by the fact that the hardness ratio, defined as the ratio between the soft and the hard count rate, does not significantly vary over the course of 6 years, fluctuating around 1.1 and 1.2 (with errors on the order of 0.1). This constant spectral behavior allows us to use all the observations simultaneously in the spectral analysis to increase the signal to noise ratio. We used the xrtgrblc tools available within the HEASoft package released by HEASARC⁵ to generate 3 sets of source and background spectral files and ancillary files. The sets correspond to the two observations in 2006, the single pointing in 2007, and all the observations afterward, respectively. Indeed, in these 3 periods the response matrix is different and the changes reflect different operating modes of the XRT. The 3 source PHA files were binned using a grouping with a minimum of 20 counts per energy bin. This allowed us to use the χ² statistic within XSPEC. We fitted the 0.3−10 keV spectral data with a power law and two layers of absorption: one fixed at the Galactic value, nH_Gal = 2.8 × 10²⁰ cm^-2 (Kalberla et al. 2005), and one free to vary. The additional absorption, nH_add = (2.1 ± 0.7) × 10²⁰ cm^-2, is needed. We ran the F-test within XSPEC and we found that the latter component is significant, with a probability under the null hypothesis of 5 × 10^-8. We also fitted the spectra assuming that the absorber is intrinsic (i.e., located at the source distance, z = 2.218). In this case, the level of absorption is higher but with a much larger uncertainty: nH_int = (20.2 ± 6.0) × 10²⁰ cm^-2. Based on the fact that typically blazars do not show any significant intrinsic absorption along the line of sight, we tend to believe that the radio estimate of the column density is simply underestimating the amount of Galactic hydrogen column density by a factor of less than 2. Assuming that the absorption is Galactic, the best fit is obtained with a spectral slope of Γ = 1.43 ± 0.03. The observed and unabsorbed fluxes in the 0.3−10 keV range are (2.11 ± 0.26) × 10^-11 erg cm^-2 s^-1 and (2.23 ± 0.30) × 10^-11 erg cm^-2 s^-1, respectively. The conversion factor from count rate to observed flux is 4.8 × 10^-11 erg cm^-2 ct^-1. Results after 2007, which are contemporaneous with the observations at other wavelengths reported here, are consistent with the results accumulated over the longer time interval.

2.2.2. UVOT

The Ultraviolet/Optical Telescope (UVOT; Roming et al. 2005) provides coverage simultaneous to the XRT for Swift observations. Querying the Swift UVOT archive, one can notice that the source was observed with all the filters during each pointing only in the first years of the mission (2006 and 2007) and during the 2011 and 2012 target of opportunity observations. In the period 2008−2010, only the “filter of the day” was used. This gives us a partial, random coverage in the U, W1, M2, and W2 filters. We analyzed the public data using the standard tools from the Swift analysis Web page of HEASARC⁶. For each filter, we combined all the exposures within a single observation to estimate the monochromatic flux (corrected for the finite aperture of the extraction region, coincidence loss, and large scale sensitivity). We summed both the sky images and the exposure maps using uvotimsum. The photometry has been obtained running uvotsource and using a circular source extraction region (with radius varying between 3′′ and 5′′, depending on the source intensity) and an annular background centered on the source (with an inner radius not less than 15′′). No significant evolution is seen, although the optical and UV light curves are sparsely populated. The flux densities across the 6-year span vary by less than 50%. Based on the paucity of the UVOT coverage, no conclusive result can be obtained about time variability.

2.2.3. BAT

The Swift Burst Alert Telescope (BAT; Barthelmy et al. 2005) is a wide-field hard X-ray telescope operating in the energy range from 14−195 keV. The BAT results used here are from the 70-month catalog⁷. The average flux in the BAT energy range is (6.8 ± 0.4) × 10^-11 erg cm^-2 s^-1. The source has shown no significant variability within this time frame, which ended in 2010 September, before the start of the γ-ray flaring activity. The power-law spectrum in the BAT catalog has photon index Γ = 1.59    ±    0.08, thus marginally steeper than the XRT spectrum.

2.3. Optical observations

R-band optical data come from the blazar monitoring program of the Abastumani Observatory (Kurtanidze et al. 2007), reduced using their standard software. The R-band photometry for 0836+710 was obtained by comparison to two nearby reference stars recommended for photometry by Villata et al. (1997)⁸.

2.4. Radio observations

0836+710 is one of the blazars studied as part of the Fermi-GST AGN Multi-frequency Monitoring Alliance (F-GAMMA) radio monitoring program (Fuhrmann et al. 2007), with an overall frequency range spanning 2.64 GHz to 142 GHz. The lower frequency radio data (2.6, 4.85, 8.35, 10.45, 14.6, 23, 32 and 43 GHz) are obtained with the Effelsberg 100-m telescope, while the millimeter IRAM 30-m observations (86 and 142 GHz) are closely coordinated with the more general flux monitoring conducted by IRAM. Data from both programs are included here.

2.5. Other recent observations

The Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010) mapped the sky at wavelengths of 3.4, 4.6, 12, and 22 μm between 2010 January and August, during the quiescent phase for 0836+710. The WISE all-sky catalog includes photometric information about 0836+710 in all four wavelength bands⁹.

The Planck early release data included observations of 0836+710 in 2009 October and 2010 March, both also during the quiescent period. The SED for the source using these data has been included in Planck Collaboration (2011) and Giommi et al. (2012). The frequency range of Planck extends higher than that included in the F-GAMMA program.

Optical polarization and 43 GHz VLBA observations of 0836+710 between 2011 April and December by Jorstad and Marscher (2013) as part of the Boston University blazar program¹⁰ have shown that the onset of the γ-ray active period in 2011 April coincided with the emission of a new radio component with β_app = 20 ± 2, together with a significant change of the optical polarization. They argue that the γ-ray emission is produced more than 20 pc from the central black hole.

3. Multiwavelength results

3.1. Long term multiwavelength light curve

Although 0836+710 is monitored regularly by radio telescopes, its coverage at other wavelengths has been fairly sparse in recent years. Figure 7 summarizes the long-term flux history of the source, anchored in the top panel by the results from Fermi-LAT and in the bottom two panels by the F-GAMMA program.

Fig. 7 Long term multiwavelength light curve of 0836+710. The LAT data are accumulated in 28-day bins. See the text for descriptions of the other data sets.

A quantitative analysis comparing the various light curves has been carried out using the Discrete Cross Correlation Function (DCCF; Edelson & Krolik 1988). Although the X-ray and UV data are too sparse for any meaningful correlation analysis, the DCCF has been applied by comparing the γ-ray data with the optical and the various radio bands.

Fig. 8 Weekly Fermi-LAT γ-ray fluxes (black) compared with Abastumani R-band optical fluxes (magenta). The data points are connected by lines to guide the eye, although only the points shown by the + symbol are actual measurements.

In the radio bands the source is highly active and is flaring during the γ-ray activity, but without a simple 1 to 1 correlation with the γ rays. The medium/high radio frequency light curves show prominent activity throughout the whole period (starting at MJD ~ 54 500), with flares at MJD ~ 54 700, ~55 250 and ~55 750 and with up to a factor of 4 change in amplitude. The strongest radio flares lead the γ-ray flares by hundreds of days, making any physical connection improbable. The lower-frequency radio bands (bottom panel of Fig. 7) do not exhibit any significant short-term flaring activity. Except at the lowest (2.64 GHz) frequency, all the lower-frequency radio bands show a generally rising trend during most of the flaring Fermi observations, a pattern that has been seen since the EGRET era, e.g. Valtaoja & Terasranta (1995). Quantitatively, the radio/γ-ray DCCFs are all consistent with “chance correlation”, which means similar to the typical DCCF levels obtained when the 0836+710 radio light curve is correlated with the LAT light curves of 130 other AGNs for comparison (see Larsson 2011). The radio/γ-ray situation can probably best be described as complex and ambiguous. Statistically, active galactic nuclei seen by the LAT have correlated radio and γ-ray fluxes (Ackermann et al. 2011). Radio/γ-ray flux comparisons for individual blazars often show little correlation, however (e.g. Abdo et al. 2010g). 0836+710 appears to be a similar case. The appearance of high-frequency radio flares without corresponding γ-ray activity suggests that the radio emission for 0836+710 must contain one or more components not related to the source of the γ rays. Such a condition could occur, for example, if a batch of synchrotron-radiating electrons lost energy before encountering a source of photons to upscatter to γ-ray energies. Another possibility is that the radio emission comes from a different region than the γ-ray emission, at a different distance from the central black hole.

Figure 8 shows an overlay of the Abastumani R-band optical data with the weekly Fermi-LAT γ-ray fluxes during the active period. There is a trend for the optical flux to be high during γ-ray enhancements, but no clear correlation of variability on shorter time scales. To evaluate the possible longer-term correlation, we constructed the DCCF using the R-band data and the 28-day Fermi-LAT γ-ray fluxes shown in Fig. 7. In this case, shown in Fig. 9, the optical/γ-ray DCCF indicates a correlation with time lag ~0−100 days (γ rays leading optical). The 90% confidence level in the figure is estimated independently for each lag bin and is based on the probability of “chance correlations” in the same way as for the radio − γ-ray correlation analysis. Other FSRQs, e.g. PKS1510−089 (Abdo et al. 2010f) and 3C 279 (Hayashida et al. 2012), have shown optical − γ-ray correlations with an optical emission delayed with respect to the γ rays, strengthening the conclusion that the correlation is real, even if the value of the lag is uncertain. The correlation also provides strong support to the identification of the γ-ray source with 0836+710.

Fig. 9 Solid curve: Discrete Cross Correlation Function comparing the available Abastumani R-band optical data points and the 28-day Fermi-LAT γ-ray flux values as a function of the lag between the optical and γ-ray data. Positive lag indicates that the γ rays lead the optical. Dashed curve: DCCF for a calculated correlation probability of 90%.

Fig. 10 The 0836+710 SED shows contemporaneous data for the quiescent (open symbols) and active (filled symbols) periods. For the LAT data, the filled symbols show the spectrum during the flare, 2011 October 25 (MJD 55 859) to 2011 December 25 (MJD 55 920).

3.2. Spectral energy distribution

Figure 10 shows the broadband SED of the source. The γ-ray part of the spectrum shows the largest variability, a common occurrence for blazars.

The combination of the WISE results with the optical and Swift UVOT observations confirms the presence of a “blue bump” feature, likely associated with the accretion disk, although the feature seems to have a rather narrow shape. Such a feature has been included in SED models (Sambruna et al. 2007; Ghisellini et al. 2010) when the blazar was in a quiescent state. The homogeneous, one-zone Synchrotron-Self-Compton (SSC) plus External Compton (EC) model with an additional accretion disk component matches the current infrared through γ-ray data reasonably well for the quiescent state, although it fails to match the radio data, indicating that much of the radio emission does not come from the compact jet that dominates most other parts of the spectrum. Sambruna et al. (2007) suggest that a higher γ-ray flux could be explained by an injected electron spectrum that extends to higher energies than the modeled one. Such a change would presumably push both the synchrotron and Compton peaks to higher energies, possibly resulting in a curved γ-ray spectrum during the flaring period. The lack of MeV or far-IR observations during the flaring state makes it hard to evaluate this possibility in quantitative terms. We have not, therefore, attempted a new SED model for this work.

We did use the comparison of optical to γ-ray variability to obtain a clue to the primary origin of the γ-ray emission. If the γ rays were produced primarily by the SSC process, the γ-ray variability should go approximately as the square of the optical variability, while an EC source would predict a linear variation (Maraschi et al. 1994). The complication of the accretion disk optical component makes this analysis only a rough approximation, obtained as follows:

• 1.

  Extrapolating the synchrotron peak through the WISE points,the contribution of synchrotron in the R filter is a factor of 10 lower than the observed values.

• 2.

  The lowest flux in the R filter is 0.55 mJy; therefore the synchrotron part of this should account for 0.55/10 = ~ 0.05 mJy.

• 3.

  If we assume that the variation in the R filter (0.9 mJy at the highest peak) is solely due to synchrotron (and not to the accretion disk), than the synchrotron flux changed by a factor of (0.9−0.55)/0.05 ~ 7.

• 4.

  The ratio between the high and low γ-ray flux levels is about 14, certainly smaller than would be expected from an SSC model.

• 5.

  The variation in the synchrotron is comparable (linear) to the variation of the Compton, pointing to a model dominated by EC.

4. Conclusions

This long-term multiwavelength study of 0836+710 shows no dramatic differences between this high-redshift blazar and other more nearby γ-ray blazars. Variability is seen on a wide range of time scales at multiple wavelengths, with the strongest variability seen in the γ-ray band. Although the optical and γ-ray light curves show some evidence of long-term correlation, the data do not provide a clear measurement of the time offset. Comparing the optical and γ-ray variability shows that an EC process is favored as the main emission mechanism of the γ-ray activity. The lack of correlation between the radio and γ-ray fluxes indicates that the radio contains components in addition to those that produce the γ-ray emission. The most unusual feature, the change in γ-ray spectral shape during the flare, might be explained by a model in which the particles in the jet are accelerated to higher energies during the high-flux state. Continued regular multiwavelength monitoring of this source will be important for resolving the ambiguities resulting from sparse coverage and for characterizing in more detail the broadband SED.

Acknowledgments

The Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat à l’Énergie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucléaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d’Études Spatiales in France. This research is partly based on observations with the 100-m telescope of the MPIfR (Max-Planck-Institut für Radioastronomie) at Effelsberg. This work has made use of observations with the IRAM 30-m telescope. Part of this work was supported by the German Deutsche Forschungsgemeinschaft, DFG project number Ts 17/2–1. This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. The Abastumani Observatory team acknowledges financial support by the Georgian National Science Foundation through grant GNSF/ST09/521 4-320. Part of this work was supported by the COST Action MP0905, “Black Holes in a Violent Universe”. A. Akyuz acknowledges financial support by The Council of Higher Education in Turkey. The authors extend special thanks to F. D’Ammando for many helpful suggestions during the preparation of this manuscript.

References

All Figures

	Fig. 1 Fermi-LAT light curve for 0836+710 during the quiescent period, MJD 54 682−55 621 (2008 August 4−2011 March 1). The source flux is shown with 7-day time binning for E > 100 MeV. Points with TS < 4 are omitted from this figure as well as from Figs. 2 and 3.
In the text

	Fig. 2 Fermi-LAT light curve for 0836+710 during the active period, MJD 55 621−55 957 (2011 March 1−2012 January 31). The source flux is shown with 7-day time binning for E > 100 MeV.
In the text

	Fig. 3 Fermi-LAT light curve for 0836+710 during the active period, MJD 55 621−55957 (2011 March 1−2012 January 31). The source flux is shown with 1-day time binning for E > 100 MeV.
In the text

	Fig. 4 Energy spectrum of 0836+710 during the quiescent period. The red shaded region is the likelihood power-law fit to the entire energy range, showing the statistical uncertainty. The blue shaded region is the likelihood log-parabola fit to the entire energy range, showing the statistical uncertainty. The reference, or pivot, energy for each of the shaded regions is 0.32 GeV.
In the text

Fig. 5 Energy spectrum of 0836+710 during the full active period (2011 March through 2012 January). The red shaded region is the likelihood power-law fit to the entire energy range, showing the statistical uncertainty. The blue shaded region is the likelihood log-parabola fit to the entire energy range, showing the statistical uncertainty. The reference, or pivot, energy for each of the shaded regions is 0.32 GeV.

In the text

	Fig. 6 Energy spectrum of 0836+710 during the flare (see Figs. 1 or 2). The red shaded region is the likelihood power-law fit to the entire energy range, showing the statistical uncertainty. The blue shaded region is the likelihood log-parabola fit to the entire energy range, showing the statistical uncertainty. The reference, or pivot, energy for each of the shaded regions is 0.32 GeV.
In the text

	Fig. 7 Long term multiwavelength light curve of 0836+710. The LAT data are accumulated in 28-day bins. See the text for descriptions of the other data sets.
In the text

	Fig. 8 Weekly Fermi-LAT γ-ray fluxes (black) compared with Abastumani R-band optical fluxes (magenta). The data points are connected by lines to guide the eye, although only the points shown by the + symbol are actual measurements.
In the text

	Fig. 9 Solid curve: Discrete Cross Correlation Function comparing the available Abastumani R-band optical data points and the 28-day Fermi-LAT γ-ray flux values as a function of the lag between the optical and γ-ray data. Positive lag indicates that the γ rays lead the optical. Dashed curve: DCCF for a calculated correlation probability of 90%.
In the text

	Fig. 10 The 0836+710 SED shows contemporaneous data for the quiescent (open symbols) and active (filled symbols) periods. For the LAT data, the filled symbols show the spectrum during the flare, 2011 October 25 (MJD 55 859) to 2011 December 25 (MJD 55 920).
In the text