Scheduling and calibration strategy for continuous radio monitoring of 1700 sources every three days

The Owens Valley Radio Observatory 40 meter telescope is currently monitoring a sample of about 1700 blazars every three days at 15 GHz, with the main scientific goal of determining the relation between the variability of blazars at radio and gamma-rays as observed with the Fermi Gamma-ray Space Telescope. The time domain relation between radio and gamma-ray emission, in particular its correlation and time lag, can help us determine the location of the high-energy emission site in blazars, a current open question in blazar research. To achieve this goal, continuous observation of a large sample of blazars in a time scale of less than a week is indispensable. Since we only look at bright targets, the time available for target observations is mostly limited by source observability, calibration requirements and slewing of the telescope. Here I describe the implementation of a practical solution to this scheduling, calibration, and slewing time minimization problem. This solution combines ideas from optimization, in particular the traveling salesman problem, with astronomical and instrumental constraints. A heuristic solution using well stablished optimization techniques and astronomical insights particular to this situation, allow us to observe all the sources in the required three days cadence while obtaining reliable calibration of the radio flux densities. Problems of this nature will only be more common in the future and the ideas presented here can be relevant for other observing programs.Comment: Published in Proc. SPIE. 9149, Observatory Operations: Strategies, Processes, and Systems V, 91492

Identification of γ-ray emission from 3C 345 and NRAO 512

For more than 15 years, since the days of the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray Observatory (CGRO; 1991−2000), it has remained an open question why the prominent blazar 3C 345 was not reliably detected at γ-ray energies ≥ 20 MeV. Recently a bright γ-ray source (0FGL J1641.4+3939/1FGL J1642.5+3947), potentially associated with 3C 345, was detected by the Large Area Telescope (LAT) on Fermi. Multiwavelength observations from radio bands to X-rays (mainly GASP-WEBT and Swift) of possible counterparts (3C 345, NRAO 512, B3 1640 + 396) were combined with 20 months of Fermi-LAT monitoring data (August 2008 − April 2010) to associate and identify the dominating γ-ray emitting counterpart of 1FGL J1642.5+3947. The source 3C 345 is identified as the main contributor for this γ-ray emitting region. However, after November 2009 (15 months), a significant excess of photons from the nearby quasar NRAO 512 started to contribute and thereafter was detected with increasing γ-ray activity, possibly adding flux to 1FGL J1642.5+3947. For the same time period and during the summer of 2010, an increase of radio, optical and X-ray activity of NRAO 512 was observed. No γ-ray emission from B3   1640 + 396 was detected

The Spectral Energy Distribution of Fermi Bright Blazars

We have conducted a detailed investigation of the broadband spectral properties of the γ-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi γ-ray spectra with Swift, radio, infra-red, optical, and other hard X-ray/γ-ray data, collected within 3 months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous spectral energy distributions (SED) for 48 LBAS blazars. The SED of these γ-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual log ν-log ν F _ν representation, the typical broadband spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SED to characterize the peak intensity of both the low- and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broadband colors (i.e., the radio to optical, α_(ro), and optical to X-ray, α_(ox), spectral slopes) and from the γ-ray spectral index. Our data show that the synchrotron peak frequency (ν^S _(peak)) is positioned between 10^(12.5) and 10^(14.5) Hz in broad-lined flat spectrum radio quasars (FSRQs) and between 10^(13) and 10^(17) Hz in featureless BL Lacertae objects. We find that the γ-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron-inverse Compton scenarios. However, simple homogeneous, one-zone, synchrotron self-Compton (SSC) models cannot explain most of our SED, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. More complex models involving external Compton radiation or multiple SSC components are required to reproduce the overall SED and the observed spectral variability. While more than 50% of known radio bright high energy peaked (HBL) BL Lacs are detected in the LBAS sample, only less than 13% of known bright FSRQs and LBL BL Lacs are included. This suggests that the latter sources, as a class, may be much fainter γ-ray emitters than LBAS blazars, and could in fact radiate close to the expectations of simple SSC models. We categorized all our sources according to a new physical classification scheme based on the generally accepted paradigm for Active Galactic Nuclei and on the results of this SED study. Since the LAT detector is more sensitive to flat spectrum γ-ray sources, the correlation between ν ^S _(peak) and γ-ray spectral index strongly favors the detection of high energy peaked blazars, thus explaining the Fermi overabundance of this type of sources compared to radio and EGRET samples. This selection effect is similar to that experienced in the soft X-ray band where HBL BL Lacs are the dominant type of blazars

Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010

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

Radio to gamma-ray variability study of blazar S5 0716+714

We present the results of a series of radio, optical, X-ray, and γ-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multifrequency observations were obtained using several ground- and space-based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend on a time scale of ~350 days. Episodes of fast variability recur on time scales of ~60−70 days. The intense and simultaneous activity at optical and γ-ray frequencies favors the synchrotron self-Compton mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/γ-ray activity period. The radio flares are characterized by a rising and a decaying stage and agrees with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield robust and self-consistent lower limits of δ ≥ 20 and equipartition magnetic field B_eq ≥ 0.36 G. Causality arguments constrain the size of emission region θ ≤ 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and γ-rays. Theoptical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods

Rapid TeV Gamma-Ray Flaring of BL Lacertae

We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6 minute exposure, when the integral flux above 200 GeV reached (3.4 ± 0.6) × 10^(–6) photons m^(–2) s^(–1), roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be 13 ± 4 minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of 3.6 ± 0.4, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results

Constraining the limiting brightness temperature and Doppler factors for the largest sample of radio bright blazars

Relativistic effects dominate the emission of blazar jets complicating our understanding of their intrinsic properties. Although many methods have been proposed to account for them, the variability Doppler factor method has been shown to describe the blazar populations best. We use a Bayesian hierarchical code called {\it Magnetron} to model the light curves of 1029 sources observed by the Owens Valley Radio Observatory's 40-m telescope as a series of flares with an exponential rise and decay, and estimate their variability brightness temperature. Our analysis allows us to place the most stringent constraints on the equipartition brightness temperature i.e., the maximum achieved intrinsic brightness temperature in beamed sources which we found to be $\rm \langle T_{eq}\rangle=2.78\times10^{11}K\pm26\%$. Using our findings we estimated the variability Doppler factor for the largest sample of blazars increasing the number of available estimates in the literature by almost an order of magnitude. Our results clearly show that $\gamma$-ray loud sources have faster and higher amplitude flares than $\gamma$-ray quiet sources. As a consequence they show higher variability brightness temperatures and thus are more relativistically beamed, with all of the above suggesting a strong connection between the radio flaring properties of the jet and $\gamma$-ray emission.Comment: 14 pages, 8 figures, accepted for publication in AP