W. Max-Moerbeck

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.

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 Beq ≥ 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. The optical/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.

Radio-to-[gamma]-ray monitoring of the narrow-line Seyfert 1 galaxy PMN J0948+0022 from 2008 to 2011

We present more than three years of observations at different frequencies, from radio to high-energy γ-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z = 0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of γ-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, which is able to develop an isotropic luminosity at γ-rays of the order of 10 48 erg s −1 , at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets with correlated multiwavelength variability (γ-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of γ-ray spectra before and including 2011 data suggested that there was a softening of the highenergy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, which are characterised by an outburst at γ-rays or very low/high flux at other wavelengths. The observed variability can largely be explained by changes in the injected power, the bulk Lorentz factor of the jet, or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at γ-rays is 2.3 ± 0.5 days. These small values underline the need of highly sampled multiwavelength campaigns to better understand the physics of these sources.

A method for the estimation of the significance of cross-correlations in unevenly sampled red-noise time series

We present a practical implementation of a Monte Carlo method to estimate the significance of cross-correlations in unevenly sampled time series of data, whose statistical properties are modelled with a simple power-law power spectral density. This implementation builds on published methods; we introduce a number of improvements in the normalization of the cross-correlation function estimate and a bootstrap method for estimating the significance of the cross-correlations. A closely related matter is the estimation of a model for the light curves, which is critical for the significance estimates. We present a graphical and quantitative demonstration that uses simulations to show how common it is to get high cross-correlations for unrelated light curves with steep power spectral densities. This demonstration highlights the dangers of interpreting them as signs of a physical connection. We show that by using interpolation and the Hanning sampling window function we are able to reduce the effects of red-noise leakage and to recover steep simple power-law power spectral densities. We also introduce the use of a Neyman construction for the estimation of the errors in the power-law index of the power spectral density. This method provides a consistent way to estimate the significance of cross-correlations in unevenly sampled time series of data.

The most powerful flaring activity from the NLSy1 PMN J0948+0022

We report on multifrequency observations performed during 2012 December–2013 August of the first narrow-line Seyfert 1 galaxy detected in γ-rays, PMN J0948+0022 (z = 0.5846). A γ-ray flare was observed by the Large Area Telescope on board Fermi during 2012 December–2013 January, reaching a daily peak flux in the 0.1–100 GeV energy range of (155 ± 31) × 10−8 ph cm−2 s−1 on 2013 January 1, corresponding to an apparent isotropic luminosity of ∼1.5 × 1048u2009ergu2009s−1. The γ-ray flaring period triggered Swift and Very Energetic Radiation Imaging Telescope Array System (VERITAS) observations in addition to radio and optical monitoring by Owens Valley Radio Observatory, Monitoring Of Jets in Active galactic nuclei with VLBA Experiments, and Catalina Real-time Transient Survey. A strong flare was observed in optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the γ-ray flare, reaching a record flux for this source from optical to γ-rays. VERITAS observations at very high energy (E > 100 GeV) during 2013 January 6–17 resulted in an upper limit of F>0.2u2009TeV < 4.0 × 10−12u2009phu2009cm−2u2009s−1. We compared the spectral energy distribution (SED) of the flaring state in 2013 January with that of an intermediate state observed in 2011. The two SEDs, modelled as synchrotron emission and an external Compton scattering of seed photons from a dust torus, can be modelled by changing both the electron distribution parameters and the magnetic field.

The F-GAMMA programme: multi-frequency study of active galactic nuclei in the Fermi era - Programme description and the first 2.5 years of monitoring

Context. To fully exploit the scientific potential of the Fermi mission for the physics of active galactic nuclei (AGN), we initiated the F-GAMMA programme. Between 2007 and 2015 the F-GAMMA was the prime provider of complementary multi-frequency monitoring in the radio regime. n nAims. We quantify the radio variability of γ-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the γ-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. n nMethods. The F-GAMMA performed monthly monitoring of a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-yr data set to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and γ-ray fluxes. n nResults. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100 m and IRAM 30 m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60–80 GHz as expected by shock-in-jet models. Compared to flat-spectrum radio quasars (FSRQs), BL Lacertae objects (BL Lacs) show systematically lower variability amplitudes, brightness temperatures, and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BLu2009Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes, brightness temperatures, and Doppler factors than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3σ, implying that γ rays are produced very close to the mm-band emission region.

37 GHz observations of narrow-line Seyfert 1 galaxies

Observations performed at Metsahovi Radio Observatory at 37 GHz are presented for a sample of 78 radio-loud and radio-quiet narrow-line Seyfert 1 (NLS1) galaxies, together with additional lower and higher frequency radio data from RATAN-600, Owens Valley Radio Observatory, and the Planck satellite. Most of the data have been gathered between February 2012 and April 2015 but for some sources even longer light curves exist. The detection rate at 37 GHz is around 19%, which is comparable to other populations of active galactic nuclei presumed to be faint at radio frequencies, such as BL Lac objects. Variability and spectral indices are determined for sources with enough detections. Based on the radio data, many NLS1 galaxies show a blazar-like radio spectra exhibiting significant variability. The spectra at a given time are often inverted or convex. The source of the high-frequency radio emission in NLS1 galaxies, detected at 37 GHz, is most probably a relativistic jet rather than star formation. Jets in NLS1 galaxies are therefore expected to be a much more common phenomenon than earlier assumed.

The awakening of the γ-ray narrow-line Seyfert 1 galaxy PKS 1502+036

After a long low-activity period, a γ-ray flare from the narrow-line Seyfert 1 PKS 1502+036 (z = 0.4089) was detected by the Large Area Telescope (LAT) on board Fermi in 2015. On 2015 December 20, the source reached a daily peak flux, in the 0.1–300 GeV band, of (93 ± 19) × 10^(−8) ph cm^(−2) s^(−1), attaining a flux of (237 ± 71) × 10^(−8) ph cm^(−2) s^(−1) on 3-h time-scales, which corresponds to an isotropic luminosity of (7.3 ± 2.1) × 10^(47) erg s^(−1). The γ-ray flare was not accompanied by significant spectral changes. We report on multiwavelength radio-to-γ-ray observations of PKS 1502+036 during 2008 August–2016 March by Fermi-LAT, Swift, XMM–Newton, Catalina Real-Time Transient Survey and the Owens Valley Radio Observatory (OVRO). An increase in activity was observed on 2015 December 22 by Swift in optical, UV and X-rays. The OVRO 15 GHz light curve reached the highest flux density observed from this source on 2016 January 12, indicating a delay of about three weeks between the γ-ray and 15 GHz emission peaks. This suggests that the γ-ray-emitting region is located beyond the broad-line region. We compared the spectral energy distribution (SED) of an average activity state with that of the flaring state. The two SED, with the high-energy bump modelled as an external Compton component with seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. The fit of the disc emission during the average state constrains the black hole mass to values lower than 10^8 M_⊙. The SED, high-energy emission mechanisms and γ-ray properties of the source resemble those of a flat spectrum radio quasar.

Optical and radio variability of the northern VHE gamma-ray emitting BL Lacertae objects

We compare the variability properties of very high energy gamma-ray emitting BL Lac objects in the optical and radio bands. We use the variability information to distinguish multiple emission components in the jet, to be used as a guidance for spectral energy distribution modelling. Our sample includes 32 objects in the Northern sky that have data for at least 2 years in both bands. We use optical R-band data from the Tuorla blazar monitoring program and 15 GHz radio data from the Owens Valley Radio Observatory blazer monitoring program. We estimate the variability amplitudes using the intrinsic modulation index, and study the time-domain connection by cross-correlating the optical and radio light curves assuming power law power spectral density. Our sample objects are in general more variable in the optical than radio. We find correlated flares in about half of the objects, and correlated long-term trends in more than 40% of the objects. In these objects we estimate that at least 10%-50% of the optical emission originates in the same emission region as the radio, while the other half is due to faster variations not seen in the radio. This implies that simple single-zone spectral energy distribution models are not adequate for many of these objects.

Radio follow-up of the γ-ray flaring gravitational lens JVAS B0218+357

We present results on multifrequency Very Long Baseline Array (VLBA) monitoring observations of the double-image gravitationally lensed blazar JVAS B0218+357. Multi-epoch observations started less than one month after the γ-ray flare detected in 2012 by the Large Area Telescope on board Fermi, and spanned a 2-month interval. The radio light curves did not reveal any significant flux density variability, suggesting that no clear correlation between the high-energy and low-energy emission is present. This behaviour was confirmed also by the long-term Owens Valley Radio Observatory monitoring data at 15 GHz. The milliarcsecond-scale resolution provided by the VLBA observations allowed us to resolve the two images of the lensed blazar, which have a core-jet structure. No significant morphological variation is found by the analysis of the multi-epoch data, suggesting that the region responsible for the γ-ray variability is located in the core of the active galactic nuclei, which is opaque up to the highest observing frequency of 22 GHz.