Source author record

A. C. S. Readhead

A. C. S. Readhead appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

82works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

82 published item(s)

preprint2025arXiv

Pinpointing the location of the gamma-ray emitting region in the FSRQ 4C+01.28

The FSRQ 4C+01.28 is a bright and highly variable radio and $γ$-ray emitter. We aim to pinpoint the location of the $γ$-ray emitting region within its jet in order to derive strong constraints on $γ$-ray emission models for blazar jets. We use radio and $γ$-ray data obtained with ALMA, OVRO, SMA and Fermi/LAT to study the cross-correlation between $γ$-ray and multi-frequency radio light curves. Moreover, we employ VLBA observations at 43 GHz over a period of around nine years to study the parsec-scale jet kinematics. To pinpoint the location of the $γ$-ray emitting region, we use a model in which outbursts shown in the $γ$-ray and radio light curves are produced when moving jet components pass through the $γ$-ray emitting and the radio core regions. We find two bright and compact newly ejected jet components that are likely associated with a high activity period visible in the $γ$-ray and radio light curves. The kinematic analysis of the VLBA observations leads to a maximum apparent jet speed of $β_{app}=19\pm10$ and an upper limit on the viewing angle of $ϕ$ < 4 deg. We determine the power law indices that are characterizing the jet geometry, brightness temperature distribution, and core shift to be $l=0.974\pm0.098$, $s=-3.31\pm0.31$, and $k_r=1.09\pm0.17$, which are in agreement with a conical jet in equipartition. A cross-correlation analysis shows that the radio light curves follow the $γ$-ray light curve. We pinpoint the location of the $γ$-ray emitting region with respect to the jet base to the range of $2.6\,\mathrm{pc}\leq d_γ\leq20\,\mathrm{pc}$. Our derived observational limits places the location of $γ$-ray production in 4C+01.28 beyond the expected extent of the broad-line region (BLR) and therefore challenges blazar-emission models that rely on inverse Compton up-scattering of seed photons from the BLR.

preprint2022arXiv

A multiband study and exploration of the radio wave - $γ$-ray connection in 3C 84

Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilise the availability of radio and $γ$-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculate the multiband time lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process regression. We find that the jet particle and magnetic field energy densities are in equipartition ($k_\textrm{r} = 1.08\pm0.18$). The jet apex is located $z_\textrm{91.5 GHz}= 22 - 645$ $R_\textrm{s}$ ($2 - 20 \times 10^{-3}$ pc) upstream of the 3 mm radio core; at that position, the magnetic field amplitude is $B_\textrm{core}^\textrm{91.5 GHz}= 3 - 10$ G. Our results are in good agreement with earlier studies, which utilised very-long-baseline interferometry. Furthermore, we investigate the temporal relation between the ejection of radio and $γ$-ray flares. Our results are in favour of the $γ$-ray emission being associated with the radio emission. We are able to tentatively connect the ejection of features identified at 43 and 86 GHz to prominent $γ$-ray flares. Finally, we compute the multiplicity parameter $λ$ and the Michel magnetisation $σ_\textrm{M}$ and find that they are consistent with a jet launched by the Blandford & Znajek 1977 mechanism.

preprint2022arXiv

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.

preprint2022arXiv

Jet kinematics in the transversely stratified jet of 3C 84 A two-decade overview

3C84 (NGC1275) is one of the brightest radio sources in the mm radio-bands, which led to a plethora of VLBI observations at numerous frequencies over the years. They reveal a two-sided jet structure, with an expanding but not well-collimated parsec-scale jet, pointing southward. High resolution mm-VLBI observations allow the study and imaging of the jet base on the sub-parsec scale. This could facilitate the investigation of the nature of the jet origin, also in view of the previously detected two-railed jet structure and east-west oriented core region seen with RadioAstron at 22 GHz. We produce VLBI images of this core and inner jet region, observed during the past twenty years at 15, 43, and 86 GHz. We determine the kinematics of the inner jet and ejected features at 43 and 86 GHz and compare their ejection times with radio and $γ$-ray variability. For the moving jet features, we find an average velocity of $β^\textrm{avg}_\textrm{app} = 0.055-0.22$c ($μ^\textrm{avg} = 0.04-0.18\,$mas/yr). From the time-averaged VLBI images at the three frequencies, we measure the transverse jet width along the bulk flow. On the $\leq 1.5$ parsec scale, we find a clear trend of the jet width being frequency dependent, with the jet being narrower at higher frequencies. This stratification is discussed in the context of a spine-sheath scenario, and is compared to other possible interpretations. From quasi-simultaneous observations at 43 and 86\,GHz we obtain spectral index maps, revealing a time-variable orientation of the spectral index gradient, due to structural variability of the inner jet.

preprint2022arXiv

New Tests of Millilensing in the Blazar PKS 1413+135

Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational millilensing, PKS 1413+135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the millilens itself. We discuss two possible candidates for the putative millilens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989-day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in August 2022, May 2025, and February 2028.

preprint2022arXiv

The C-Band All-Sky Survey (C-BASS): Template Fitting of Diffuse Galactic Microwave Emission in the Northern Sky

The C-Band All-Sky Survey (C-BASS) has observed the Galaxy at 4.76GHz with an angular resolution of $0.73^\circ$ full-width half-maximum, and detected Galactic synchrotron emission with high signal-to-noise ratio over the entire northern sky ($δ> -15^{\circ}$). We present the results of a spatial correlation analysis of Galactic foregrounds at mid-to-high ($b > 10^\circ$) Galactic latitudes using a preliminary version of the C-BASS intensity map. We jointly fit for synchrotron, dust, and free-free components between $20$ and $1000$GHz and look for differences in the Galactic synchrotron spectrum, and the emissivity of anomalous microwave emission (AME) when using either the C-BASS map or the 408MHz all-sky map to trace synchrotron emission. We find marginal evidence for a steepening ($\left<Δβ\right> = -0.06\pm0.02$) of the Galactic synchrotron spectrum at high frequencies resulting in a mean spectral index of $\left<β\right> = -3.10\pm0.02$ over $4.76-22.8$GHz. Further, we find that the synchrotron emission can be well modelled by a single power-law up to a few tens of GHz. Due to this, we find that the AME emissivity is not sensitive to changing the synchrotron tracer from the 408MHz map to the 4.76GHz map. We interpret this as strong evidence for the origin of AME being spinning dust emission.

preprint2022arXiv

The Unanticipated Phenomenology of the Blazar PKS~2131$-$021: A Unique Super-Massive Black Hole Binary Candidate

Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within $<2\%$ and $\sim 10\%$, respectively, straddling a 20-year period when this variation was absent. Our simulated light curves accurately reproduce the ``red noise'' of this object, and Lomb-Scargle, weighted wavelet Z-transform, and least-squares sine wave analyses demonstrate conclusively, at the $4.6σ$ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to $2.082\pm 0.003$ years in the rest frame at the $z=1.285$ redshift of PKS~2131$-$021. The periodic variation in PKS~2131$-$021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily fits the observations. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS~2131$-$021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field.

preprint2021arXiv

Detection of Spectral Variations of Anomalous Microwave Emission with QUIJOTE and C-BASS

Anomalous Microwave Emission (AME) is a significant component of Galactic diffuse emission in the frequency range $10$-$60\,$GHz and a new window into the properties of sub-nanometre-sized grains in the interstellar medium. We investigate the morphology of AME in the $\approx10^{\circ}$ diameter $λ$ Orionis ring by combining intensity data from the QUIJOTE experiment at $11$, $13$, $17$ and $19\,$GHz and the C-Band All Sky Survey (C-BASS) at $4.76\,$GHz, together with 19 ancillary datasets between $1.42$ and $3000\,$GHz. Maps of physical parameters at $1^{\circ}$ resolution are produced through Markov Chain Monte Carlo (MCMC) fits of spectral energy distributions (SEDs), approximating the AME component with a log-normal distribution. AME is detected in excess of $20\,σ$ at degree-scales around the entirety of the ring along photodissociation regions (PDRs), with three primary bright regions containing dark clouds. A radial decrease is observed in the AME peak frequency from $\approx35\,$GHz near the free-free region to $\approx21\,$GHz in the outer regions of the ring, which is the first detection of AME spectral variations across a single region. A strong correlation between AME peak frequency, emission measure and dust temperature is an indication for the dependence of the AME peak frequency on the local radiation field. The AME amplitude normalised by the optical depth is also strongly correlated with the radiation field, giving an overall picture consistent with spinning dust where the local radiation field plays a key role.

preprint2021arXiv

Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We constrain the VHE duty cycle of the source to be < 16 2014-like flares per year (95% confidence). For the first time for this source, a broadband low-state SED is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 +- 0.93 stat ) x 10^21 cm^-2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modeled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet.

preprint2020arXiv

A decade of multi-wavelength observations of the TeV blazar 1ES 1215+303: Extreme shift of the synchrotron peak frequency and long-term optical-gamma-ray flux increase

Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed.

preprint2020arXiv

An intermittent extreme BL Lac: MWL study of 1ES 2344+514 in an enhanced state

Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $ν_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $γ$-ray telescope FACT during a high $γ$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $γ$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $γ$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $Γ$ index of the intrinsic spectrum in the VHE $γ$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz

preprint2020arXiv

Testing two-component models on very-high-energy gamma-ray emitting BL Lac objects

Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

preprint2020arXiv

The C-Band All-Sky Survey (C-BASS): Total intensity point-source detection over the northern sky

We present a point-source detection algorithm that employs the second order Spherical Mexican Hat wavelet filter (SMHW2), and use it on C-BASS northern intensity data to produce a catalogue of point-sources. This catalogue allows us to cross-check the C-BASS flux-density scale against existing source surveys, and provides the basis for a source mask which will be used in subsequent C-BASS and cosmic microwave background (CMB) analyses. The SMHW2 allows us to filter the entire sky at once, avoiding complications from edge effects arising when filtering small sky patches. The algorithm is validated against a set of Monte Carlo simulations, consisting of diffuse emission, instrumental noise, and various point-source populations. The simulated source populations are successfully recovered. The SMHW2 detection algorithm is used to produce a $4.76\,\mathrm{GHz}$ northern sky source catalogue in total intensity, containing 1784 sources and covering declinations $δ\geq-10^{\circ}$. The C-BASS catalogue is matched with the Green Bank 6\,cm (GB6) and Parkes-MIT-NRAO (PMN) catalogues over their areas of common sky coverage. From this we estimate the $90$ per cent completeness level to be approximately $610\,\mathrm{mJy}$, with a corresponding reliability of $98$ per cent, when masking the brightest $30$ per cent of the diffuse emission in the C-BASS northern sky map. We find the C-BASS and GB6 flux-density scales to be consistent with one another to within approximately $4$ per cent.

preprint2020arXiv

The Great Markarian 421 Flare of February 2010: Multiwavelength variability and correlation studies

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $γ$-ray (VERITAS, MAGIC), high-energy (HE) $γ$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Metsähovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($δ\gtrsim 33$) and the size of the emission region ($ δ^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

preprint2020arXiv

The Relativistic Jet Orientation and Host Galaxy of the Peculiar Blazar PKS 1413+135

PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and is one of the rare objects in this class with a visible counterjet. BL Lac objects have jet axes aligned close to the line of sight. It has also been classified as a compact symmetric object, which have jet axes not aligned close to the line of sight. Intensive efforts to understand this blazar have hitherto failed to resolve even the questions of the orientation of the relativistic jet, and the host galaxy. Answering these two questions is important as they challenge our understanding of jets in active galactic nuclei and the classification schemes we use to describe them. We show that the jet axis is aligned close to the line of sight and PKS 1413+135 is almost certainly not located in the apparent host galaxy, but is a background object in the redshift range $0.247 < z < 0.5$. The intervening spiral galaxy at $z = 0.247$ provides a natural host for the putative lens responsible for symmetric achromatic variability and is shown to be a Seyfert 2 galaxy. We also show that, as for the radio emission, a "multizone" model is needed to account for the high-energy emission.

preprint2018arXiv

Multi-wavelength characterization of the blazar S5~0716+714 during an unprecedented outburst phase

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $γ$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

preprint2018arXiv

The extreme HBL behaviour of Markarian 501 during 2012

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

preprint2016arXiv

Insights into the emission of the blazar 1ES 1011+496 through unprecedented broadband observations during 2011 and 2012

1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $γ$-rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to a rather incomplete characterization of the broadband spectral energy distribution (SED). We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron-self-Compton (SSC) scenario is able to explain the observed broadband spectrum. We analyzed VHE to radio data from 2011 and 2012 collected by MAGIC, $Fermi$-LAT, $Swift$, KVA, OVRO, and Metsähovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE. The VHE spectrum was fit with a simple power law with a photon index of $3.69\pm0.22$ and a flux above 150 GeV of $(1.46\pm0.16)\times10^{-11}$ ph cm$^{-2}$ s$^{-1}$. 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of $1.8\pm0.4\,c$, which is the highest speed statistically significantly measured so far in a high-frequency-peaked BL Lac. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field to the optical emission.

preprint2016arXiv

Investigating the peculiar emission from the new VHE gamma-ray source H1722+119

The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, $E > 100$ GeV) $γ$-ray band with a statistical significance of 5.9 $σ$. The integral flux above 150 GeV is estimated to be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous high energy (HE, 100 MeV $ < E < 100$ GeV) $γ$-ray observations from Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution (SED) of H1722+119 shows surprising behaviour in the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model.

preprint2016arXiv

Optical polarization map of the Polaris Flare with RoboPol

The stages before the formation of stars in molecular clouds are poorly understood. Insights can be gained by studying the properties of quiescent clouds, such as their magnetic field structure. The plane-of-the-sky orientation of the field can be traced by polarized starlight. We present the first extended, wide-field ($\sim$10 $\rm deg^2$) map of the Polaris Flare cloud in dust-absorption induced optical polarization of background stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first application of the wide-field imaging capabilities of RoboPol. The data were taken in the R-band and analysed with the automated reduction pipeline of the instrument. We present in detail optimizations in the reduction pipeline specific to wide-field observations. Our analysis resulted in reliable measurements of 641 stars with median fractional linear polarization 1.3%. The projected magnetic field shows a large scale ordered pattern. At high longitudes it appears to align with faint striations seen in the Herschel-SPIRE map of dust emission (250 $μm$), while in the central 4-5 deg$^2$ it shows an eddy-like feature. The overall polarization pattern we obtain is in good agreement with large scale measurements by Planck of the dust emission polarization in the same area of the sky.

preprint2016arXiv

Radio follow-up of the gamma-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 gamma-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 gamma-ray variability is located in the core of the AGN, which is opaque up to the highest observing frequency of 22 GHz.

preprint2016arXiv

RoboPol: Do optical polarization rotations occur in all blazars?

We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.

preprint2016arXiv

RoboPol: First season rotations of optical polarization plane in blazars

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.

preprint2016arXiv

RoboPol: optical polarization-plane rotations and flaring activity in blazars

We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.

preprint2016arXiv

RoboPol: The optical polarization of gamma-ray--loud and gamma-ray--quiet blazars

We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3σ level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.

preprint2016arXiv

The awakening of the gamma-ray narrow-Line Seyfert 1 galaxy PKS 1502+036

After a long low-activity period, a gamma-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 $\pm$ 19) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, attaining a flux of (237 $\pm$ 71) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 3-hr time-scales, which corresponds to an isotropic luminosity of (7.3 $\pm$ 2.1) $\times$10$^{47}$ erg/s. The gamma-ray flare was not accompanied by significant spectral changes. We report on multi-wavelength radio-to-gamma-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 gamma-ray and 15 GHz emission peaks. This suggests that the gamma-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$ solar masses. The SED, high-energy emission mechanisms, and gamma-ray properties of the source resemble those of a flat spectrum radio quasar.

preprint2016arXiv

The F-GAMMA program: Multi-frequency study of Active Galactic Nuclei in the Fermi era. Program description and the first 2.5 years of monitoring

To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly 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-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. 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 FSRQs, 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 BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes as well as 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 3sigma, implying that gamma rays are produced very close to the mm-band emission region.

preprint2015arXiv

A combined radio and GeV gamma-ray view of the 2012 and 2013 flares of Mrk 421

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio frequencies. In the present study, we start exploring this unique event and compare it to a less extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Observatory 40-m telescope, 95 GHz millimeter data from the Combined Array for Research in Millimeter-Wave Astronomy, and GeV gamma-ray data from the Fermi Gamma-ray Space Telescope. The radio light curves during the flaring periods in 2012 and 2013 have very different appearances, both in shape and peak flux density. Assuming that the radio and gamma-ray flares are physically connected, we attempt to model the most prominent sub-flares of the 2012 and 2013 activity periods by using the simplest possible theoretical framework. We first fit a one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate parameters describing the emission region. We then model the major gamma-ray and radio flares of 2012 using the same framework. The 2012 gamma-ray flare shows two distinct spikes of similar amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can accommodate this by adding plausible time variations to the Doppler beaming factor. In the second scenario, a varying Doppler factor is not needed, but the SSC model parameters require fine tuning. Both alternatives indicate that the sharp radio flare, if physically connected to the preceding gamma-ray flares, can be reproduced only for a very specific choice of parameters.

preprint2015arXiv

C-Band All-Sky Survey: A First Look at the Galaxy

We present an analysis of the diffuse emission at 5 GHz in the first quadrant of the Galactic plane using two months of preliminary intensity data taken with the C-Band All Sky Survey (C-BASS) northern instrument at the Owens Valley Radio Observatory, California. Combining C-BASS maps with ancillary data to make temperature-temperature plots we find synchrotron spectral indices of $β= -2.65 \pm 0.05$ between 0.408 GHz and 5 GHz and $ β= -2.72 \pm 0.09$ between 1.420 GHz and 5 GHz for $-10^{\circ} < |b| < -4^{\circ}$, $20^{\circ} < l < 40^{\circ}$. Through the subtraction of a radio recombination line (RRL) free-free template we determine the synchrotron spectral index in the Galactic plane ($ |b| < 4^{\circ}$) to be $β= -2.56 \pm 0.07$ between 0.408 GHz and 5 GHz, with a contribution of $53 \pm 8$ per cent from free-free emission at 5\,GHz. These results are consistent with previous low frequency measurements in the Galactic plane. By including C-BASS data in spectral fits we demonstrate the presence of anomalous microwave emission (AME) associated with the HII complexes W43, W44 and W47 near 30 GHz, at 4.4 sigma, 3.1 sigma and 2.5 sigma respectively. The CORNISH VLA 5 GHz source catalogue rules out the possibility that the excess emission detected around 30\;GHz may be due to ultra-compact HII regions. Diffuse AME was also identified at a 4 sigma level within $30^{\circ} < l < 40^{\circ}$, $-2^{\circ} < b < 2^{\circ}$ between 5 GHz and 22.8 GHz.

preprint2015arXiv

High Radio Frequency Properties and Variability of Brightest Cluster Galaxies

We consider the high radio frequency (15 GHz - 353 GHz) properties and variability of 35 Brightest Cluster Galaxies (BCGs). These are the most core-dominated sources drawn from a parent sample of more than 700 X-ray selected clusters, thus allowing us to relate our results to the general population. We find that >6.0% of our parent sample (>15.1% if only cool-core clusters are considered) contain a radio-source at 150 GHz of at least 3mJy (~1x10^23 W/Hz at our median redshift of z~0.13). Furthermore, >3.4% of the BCGs in our parent sample contain a peaked component (Gigahertz Peaked Spectrum, GPS) in their spectra that peaks above 2 GHz, increasing to >8.5% if only cool-core clusters are considered. We see little evidence for strong variability at 15 GHz on short (week-month) time-scales although we see variations greater than 20% at 150 GHz over 6-month times-frames for 4 of the 23 sources with multi-epoch observations. Much more prevalent is long-term (year-decade time-scale) variability, with average annual amplitude variations greater than 1% at 15 GHz being commonplace. There is a weak trend towards higher variability as the peak of the GPS-like component occurs at higher frequency. We demonstrate the complexity that is seen in the radio spectra of BCGs and discuss the potentially significant implications of these high-peaking components for Sunyaev-Zel'dovich cluster searches.

preprint2015arXiv

The 2009 multiwavelength campaign on Mrk 421: Variability and correlation studies

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

preprint2015arXiv

The Q/U Imaging ExperimenT: Polarization Measurements of the Galactic Plane at 43 and 95 GHz

We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations. The resulting maps represent the deepest large-area Galactic polarization observations published to date at the relevant frequencies with instrumental rms noise varying between 1.8 and 2.8 uK deg, 2.3-6 times deeper than corresponding WMAP and Planck maps. The angular resolution is 27.3' and 12.8' FWHM at Q- and W-band, respectively. We find excellent agreement between the QUIET and WMAP maps over the entire fields, and no compelling evidence for significant residual instrumental systematic errors in either experiment, whereas the Planck 44 GHz map deviates from these in a manner consistent with reported systematic uncertainties for this channel. We combine QUIET and WMAP data to compute inverse-variance-weighted average maps, effectively retaining small angular scales from QUIET and large angular scales from WMAP. From these combined maps, we derive constraints on several important astrophysical quantities, including a robust detection of polarized synchrotron spectral index steepening of ~0.2 off the plane, as well as the Faraday rotation measure toward the Galactic center (RM=-4000 +/- 200 rad m^-2), all of which are consistent with previously published results. Both the raw QUIET and the co-added QUIET+WMAP maps are made publicly available together with all necessary ancillary information.

preprint2015arXiv

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

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. 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 one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC 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 to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to 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 the variation of the parameters related to the emitting region itself ($δ$, $B$ and $R$), in addition to the parameters related to the particle population.

preprint2015arXiv

Why Haven't Many of the Brightest Radio Loud Blazars Been Detected by Fermi ?

We use the complete MOJAVE 1.5 Jy sample of active galactic nuclei (AGN) to examine the gamma-ray detection statistics of the brightest radio-loud blazars in the northern sky. We find that 23% of these AGN were not detected above 0.1 GeV by the Fermi LAT during the 4-year 3FGL catalog period partly because of an instrumental selection effect, and partly due to their lower Doppler boosting factors. Blazars with synchrotron peaks in their spectral energy distributions located below $10^{13.4}$ Hz also tend to have high-energy peaks that lie below the 0.1 GeV threshold of the LAT, and are thus less likely to be detected by Fermi. The non-detected AGN in the 1.5 Jy sample also have significantly lower 15 GHz radio modulation indices and apparent jet speeds, indicating that they have lower than average Doppler factors. Since the effective amount of relativistic Doppler boosting is enhanced in gamma-rays (particularly in the case of external inverse-Compton scattering), this makes them less likely to appear in the 3FGL catalog. Based on their observed properties, we have identified several bright radio-selected blazars that are strong candidates for future detection by Fermi.

preprint2014arXiv

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 modeled 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.

preprint2014arXiv

Connection between optical and gamma-ray variability in blazars

We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.

preprint2014arXiv

Detection of significant cm to sub-mm band radio and gamma-ray correlated variability in Fermi bright blazars

The exact location of the gamma-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA program) and gamma-ray (0.1-300 GeV) ~ 3.5 year light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. Our results reveal: (i) the first highly significant multi-band radio and gamma-ray correlations (radio lagging gamma rays) when averaging over the whole sample, (ii) average time delays (source frame: 76+/-23 to 7+/-9 days), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence tau_r,gamma (nu) ~ nu^-1, in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk gamma-ray production region typically located within/upstream of the 3 mm core region (tau_3mm,gamma=12+/-8 days), (iv) mean distances between the region of gamma-ray peak emission and the radio "tau=1 photosphere" decreasing from 9.8+/-3.0 pc (11 cm) to 0.9+/-1.1 pc (2 mm) and 1.4+/-0.8 pc (0.8 mm), (v) 3 mm/gamma-ray correlations in 9 individual sources at a significance level where one is expected by chance (probability: 4 times 10^-6), (vi) opacity and "time lag core shift" estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk gamma-ray production region from the supermassive black hole (SMBH) of ~ 0.8-1.6 pc, i.e. at the outer edge of the Broad Line Region (BLR) or beyond. A 3 mm tau=1 surface at ~ 2-3 pc from the jet-base (i.e. well outside the "canonical BLR") finally suggests that BLR material extends to several pc distances from the SMBH.

preprint2014arXiv

Early-time polarized optical light curve of GRB 131030A

We report the polarized optical light curve of a gamma-ray burst afterglow obtained using the RoboPol instrument. Observations began 655 seconds after the initial burst of gamma-rays from GRB131030A, and continued uninterrupted for 2 hours. The afterglow displayed a low, constant fractional linear polarization of $p = (2.1 \pm 1.6)\,\%$ throughout, which is similar to the interstellar polarization measured on nearby stars. The optical brightness decay is consistent with a forward-shock propagating in a medium of constant density, and the low polarization fraction indicates a disordered magnetic field in the shock front. This supports the idea that the magnetic field is amplified by plasma instabilities on the shock front. These plasma instabilities produce strong magnetic fields with random directions on scales much smaller than the total observable region of the shock, and the resulting randomly-oriented polarization vectors sum to produce a low net polarization over the total observable region of the shock.

preprint2014arXiv

First broadband characterization and redshift determination of the VHE blazar MAGIC J2001+439

We aim to characterize the broadband emission from 2FGL J2001.1+4352, which has been associated with the unknown-redshift blazar MG4 J200112+4352. Based on its gamma-ray spectral properties, it was identified as a potential very high energy (VHE; E > 100 GeV) gamma-ray emitter. The source was observed with MAGIC first in 2009 and later in 2010 within a multi-instrument observation campaign. The MAGIC observations yielded 14.8 hours of good quality stereoscopic data. The object was monitored at radio, optical and gamma-ray energies during the years 2010 and 2011. The source, named MAGIC J2001+439, is detected for the first time at VHE with MAGIC at a statistical significance of 6.3 σ (E > 70 GeV) during a 1.3-hour long observation on 2010 July 16. The multi-instrument observations show variability in all energy bands with the highest amplitude of variability in the X-ray and VHE bands. We also organized deep imaging optical observations with the Nordic Optical Telescope in 2013 to determine the source redshift. We determine for the first time the redshift of this BL Lac object through the measurement of its host galaxy during low blazar activity. Using the observational evidence that the luminosities of BL Lac host galaxies are confined to a relatively narrow range, we obtain z = 0.18 +/- 0.04. Additionally, we use the Fermi-LAT and MAGIC gamma-ray spectra to provide an independent redshift estimation, z = 0.17 +/- 0.10. Using the former (more accurate) redshift value, we adequately describe the broadband emission with a one-zone SSC model for different activity states and interpret the few-day timescale variability as produced by changes in the high-energy component of the electron energy distribution.

preprint2014arXiv

MAGIC gamma-ray and multifrequency observations of flat spectrum radio quasar PKS 1510-089 in early 2012

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

preprint2014arXiv

The connection between the 15 GHz radio and gamma-ray emission in blazars

Since mid-2007 we have carried out a dedicated long-term monitoring programme at 15 GHz using the Owens Valley Radio Observatory 40 meter telescope. One of the main goals of this programme is to study the relation between the radio and gamma-ray emission in blazars and to use it as a tool to locate the site of high energy emission. Using this large sample of objects we are able to characterize the radio variability, and study the significance of correlations between the radio and gamma-ray bands. We find that the radio variability of many sources can be described using a simple power law power spectral density, and that when taking into account the red-noise characteristics of the light curves, cases with significant correlation are rare. We note that while significant correlations are found in few individual objects, radio variations are most often delayed with respect to the gamma-ray variations. This suggests that the gamma-ray emission originates upstream of the radio emission. Because strong flares in most known gamma-ray-loud blazars are infrequent, longer light curves are required to settle the issue of the strength of radio-gamma cross-correlations and establish confidently possible delays between the two. For this reason continuous multiwavelength monitoring over a longer time period is essential for statistical tests of jet emission models.

preprint2014arXiv

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 gamma rays, PMN J0948+0022 ($z$ = 0.5846). A gamma-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 $\pm$ 31) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 2013 January 1, corresponding to an apparent isotropic luminosity of about 1.5$\times$10$^{48}$ erg s$^{-1}$. The gamma-ray flaring period triggered Swift and VERITAS observations in addition to radio and optical monitoring by OVRO, MOJAVE, and CRTS. A strong flare was observed in optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the gamma-ray flare, reaching a record flux for this source from optical to gamma rays. VERITAS observations at very high energy (E > 100 GeV) during 2013 January 6-17 resulted in an upper limit of F (> 0.2 TeV) < 4.0$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-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.

preprint2014arXiv

The Parsec-scale Structure, Kinematics, and Polarization of Radio-Loud Narrow-Line Seyfert 1 Galaxies

Several narrow-line Seyfert 1 galaxies (NLS1s) have now been detected in gamma rays, providing firm evidence that at least some of this class of active galactic nuclei (AGN) produce relativistic jets. The presence of jets in NLS1s is surprising, as these sources are typified by comparatively small black hole masses and near- or super-Eddington accretion rates. This challenges the current understanding of the conditions necessary for jet production. Comparing the properties of the jets in NLS1s with those in more familiar jetted systems is thus essential to improve jet production models. We present early results from our campaign to monitor the kinematics and polarization of the parsec-scale jets in a sample of 15 NLS1s through multifrequency observations with the Very Long Baseline Array. These observations are complemented by fast-cadence 15 GHz monitoring with the Owens Valley Radio Observatory 40m telescope and optical spectroscopic monitoring with with the 2m class telescope at the Guillermo Haro Astrophysics Observatory in Cananea, Mexico.

preprint2014arXiv

The Q/U Imaging Experiment: Polarization Measurements of Radio Sources at 43 and 95 GHz

We present polarization measurements of extragalactic radio sources observed during the Cosmic Microwave Background polarization survey of the Q/U Imaging Experiment (QUIET), operating at 43 GHz (Q-band) and 95 GHz (W-band). We examine sources selected at 20 GHz from the public, $>$40 mJy catalog of the Australia Telescope (AT20G) survey. There are $\sim$480 such sources within QUIET's four low-foreground survey patches, including the nearby radio galaxies Centaurus A and Pictor A. The median error on our polarized flux density measurements is 30--40 mJy per Stokes parameter. At S/N $> 3$ significance, we detect linear polarization for seven sources in Q-band and six in W-band; only $1.3 \pm 1.1$ detections per frequency band are expected by chance. For sources without a detection of polarized emission, we find that half of the sources have polarization amplitudes below 90 mJy (Q-band) and 106 mJy (W-band), at 95% confidence. Finally, we compare our polarization measurements to intensity and polarization measurements of the same sources from the literature. For the four sources with WMAP and Planck intensity measurements $>1$ Jy, the polarization fraction are above 1% in both QUIET bands. At high significance, we compute polarization fractions as much as 10--20% for some sources, but the effects of source variability may cut that level in half for contemporaneous comparisons. Our results indicate that simple models---ones that scale a fixed polarization fraction with frequency---are inadequate to model the behavior of these sources and their contributions to polarization maps.

preprint2014arXiv

The RoboPol optical polarization survey of gamma-ray - loud blazars

We present first results from RoboPol, a novel-design optical polarimeter operating at the Skinakas Observatory in Crete. The data, taken during the May - June 2013 commissioning of the instrument, constitute a single-epoch linear polarization survey of a sample of gamma-ray - loud blazars, defined according to unbiased and objective selection criteria, easily reproducible in simulations, as well as a comparison sample of, otherwise similar, gamma-ray - quiet blazars. As such, the results of this survey are appropriate for both phenomenological population studies and for tests of theoretical population models. We have measured polarization fractions as low as $0.015$ down to $R$ magnitude of 17 and as low as $0.035$ down to 18 magnitude. The hypothesis that the polarization fractions of gamma-ray - loud and gamma-ray - quiet blazars are drawn from the same distribution is rejected at the $10^{-3}$ level. We therefore conclude that gamma-ray - loud and gamma-ray - quiet sources have different optical polarization properties. This is the first time this statistical difference is demonstrated in optical wavelengths. The polarization fraction distributions of both samples are well-described by exponential distributions with averages of $\langle p \rangle =6.4 ^{+0.9}_{-0.8}\times 10^{-2}$ for gamma-ray--loud blazars, and $\langle p \rangle =3.2 ^{+2.0}_{-1.1}\times 10^{-2}$ for gamma-ray--quiet blazars. The most probable value for the difference of the means is $3.4^{+1.5}_{-2.0}\times 10^{-2}$. The distribution of polarization angles is statistically consistent with being uniform.

preprint2014arXiv

The RoboPol Pipeline and Control System

We describe the data reduction pipeline and control system for the RoboPol project. The RoboPol project is monitoring the optical $R$-band magnitude and linear polarization of a large sample of active galactic nuclei that is dominated by blazars. The pipeline calibrates and reduces each exposure frame, producing a measurement of the magnitude and linear polarization of every source in the $13'\times 13'$ field of view. The control system combines a dynamic scheduler, real-time data reduction, and telescope automation to allow high-efficiency unassisted observations.

preprint2014arXiv

Time-correlation between the radio and gamma-ray activity in blazars and the production site of the gamma-ray emission

In order to determine the location of the gamma-ray emission site in blazars, we investigate the time-domain relationship between their radio and gamma-ray emission. Light-curves for the brightest detected blazars from the first 3 years of the mission of the Fermi Gamma-ray Space Telescope are cross-correlated with 4 years of 15GHz observations from the OVRO 40-m monitoring program. The large sample and long light-curve duration enable us to carry out a statistically robust analysis of the significance of the cross-correlations, which is investigated using Monte Carlo simulations including the uneven sampling and noise properties of the light-curves. Modeling the light-curves as red noise processes with power-law power spectral densities, we find that only one of 41 sources with high quality data in both bands shows correlations with significance larger than 3-sigma (AO 0235+164), with only two more larger than even 2.25-sigma (PKS 1502+106 and B2 2308+34). Additionally, we find correlated variability in Mrk 421 when including a strong flare that occurred in July-September 2012. These results demonstrate very clearly the difficulty of measuring statistically robust multiwavelength correlations and the care needed when comparing light-curves even when many years of data are used. This should be a caution. In all four sources the radio variations lag the gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. Continuous simultaneous monitoring over a longer time period is required to obtain high significance levels in cross-correlations between gamma-ray and radio variability in most blazars.

preprint2014arXiv

Two active states of the narrow-line gamma-ray-loud AGN GB 1310+487

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

preprint2013arXiv

A quasi-periodic oscillation in the blazar J1359+4011

The OVRO 40-m telescope has been monitoring the 15 GHz radio flux density of over 1200 blazars since 2008. The 15 GHz light curve of the flat spectrum radio quasar J1359+4011 shows a strong and persistent quasi-periodic oscillation. The time-scale of the oscillation varies between 120 and 150 days over a $\sim4$ year time span. We interpret this as the active galactic nucleus mass-scaled analog of low-frequency quasi-periodic oscillations from Galactic microquasars, or as evidence of modulation of the accretion flow by thermal instabilites in the "inner" accretion disc.

preprint2013arXiv

Intrinsic brightness temperatures of blazar jets at 15 GHz

We have developed a new Bayesian Markov Chain Monte Carlo method to deconvolve light curves of blazars into individual flares, including proper estimation of the fit errors. We use the method to fit 15\,GHz light curves obtained within the OVRO 40-m blazar monitoring program where a large number of AGN have been monitored since 2008 in support of the Fermi Gamma-Ray Space Telescope mission. The time scales obtained from the fitted models are used to calculate the variability brightness temperature of the sources. Additionally, we have calculated brightness temperatures of a sample of these objects using Very Long Baseline Array data from the MOJAVE survey. Combining these two data sets enables us to study the intrinsic brightness temperature distribution in these blazars at 15 GHz. Our preliminary results indicate that the mean intrinsic brightness temperature in a sample of 14 sources is near the equipartition brightness temperature of $\sim 10^{11}$K.

preprint2013arXiv

Multifrequency studies of the narrow-line Seyfert 1 galaxy SBS 0846+513

The narrow-line Seyfert 1 galaxy SBS 0846+513 was first detected by the Large Area Telescope (LAT) on-board Fermi in 2011 June-July when it underwent a period of flaring activity. Since then, as Fermi continues to accumulate data on this source, its flux has been monitored on a daily basis. Two further gamma-ray flaring episodes from SBS 0846+513 were observed in 2012 May and August, reaching a daily peak flux integrated above 100 MeV of (50+/-12)x10^-8 ph/cm^2/s, and (73+/-14)x10^-8 ph/cm^2/s on May 24 and August 7, respectively. Three outbursts were detected at 15 GHz by the Owens Valley Radio Observatory 40-m telescope in 2012 May, 2012 October, and 2013 January, suggesting a complex connection with the gamma-ray activity. The most likely scenario suggests that the 2012 May gamma-ray flare may not be directly related to the radio activity observed over the same period, while the two gamma-ray flaring episodes may be related to the radio activity observed at 15 GHz in 2012 October and 2013 January. The gamma-ray flare in 2012 May triggered Swift observations that confirmed that SBS 0846+513 was also exhibiting high activity in the optical, UV and X-ray bands, thus providing a firm identification between the gamma-ray source and the lower-energy counterpart. We compared the spectral energy distribution (SED) of the flaring state in 2012 May with that of a quiescent state. The two SEDs, modelled as an external Compton component of seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. No significant evidence of thermal emission from the accretion disc has been observed. Interestingly, in the 5 GHz radio luminosity vs. synchrotron peak frequency plot SBS 0846+513 seems to lie in the flat spectrum radio quasar part of the so-called `blazar sequence'.

preprint2013arXiv

Multiwavelength observations of the gamma-ray emitting narrow-line Seyfert 1 PMN J0948+0022 in 2011

We report on radio-to-gamma-ray observations during 2011 May-September of PMN J0948+0022, the first narrow-line Seyfert 1 (NLSy1) galaxy detected in gamma-rays by Fermi-LAT. Strong variability was observed in gamma-rays, with two flaring periods peaking on 2011 June 20 and July 28. The variability observed in optical and near-infrared seems to have no counterpart in gamma-rays. This different behaviour could be related to a bending and inhomogeneous jet or a turbulent extreme multi-cell scenario. The radio spectra showed a variability pattern typical of relativistic jets. The XMM spectrum shows that the emission from the jet dominates above 2 keV, while a soft X-ray excess is evident in the low-energy part of the X-ray spectrum. Models where the soft emission is partly produced by blurred reflection or Comptonisation of the thermal disc emission provide good fits to the data. The X-ray spectral slope is similar to that found in radio-quiet NLSy1, suggesting that a standard accretion disc is present, as expected from the high accretion rate. Except for the soft X-ray excess, unusual in jet-dominated AGNs, PMN J0948+0022 shows all characteristics of the blazar class.

preprint2013arXiv

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

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency 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 at 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 gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement 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 a 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 gamma-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.

preprint2013arXiv

Second Season QUIET Observations: Measurements of the CMB Polarization Power Spectrum at 95 GHz

The Q/U Imaging ExperimenT (QUIET) has observed the cosmic microwave background (CMB) at 43 and 95GHz. The 43-GHz results have been published in QUIET Collaboration et al. (2011), and here we report the measurement of CMB polarization power spectra using the 95-GHz data. This data set comprises 5337 hours of observations recorded by an array of 84 polarized coherent receivers with a total array sensitivity of 87 uK sqrt(s). Four low-foreground fields were observed, covering a total of ~1000 square degrees with an effective angular resolution of 12.8', allowing for constraints on primordial gravitational waves and high-signal-to-noise measurements of the E-modes across three acoustic peaks. The data reduction was performed using two independent analysis pipelines, one based on a pseudo-Cl (PCL) cross-correlation approach, and the other on a maximum-likelihood (ML) approach. All data selection criteria and filters were modified until a predefined set of null tests had been satisfied before inspecting any non-null power spectrum. The results derived by the two pipelines are in good agreement. We characterize the EE, EB and BB power spectra between l=25 and 975 and find that the EE spectrum is consistent with LCDM, while the BB power spectrum is consistent with zero. Based on these measurements, we constrain the tensor-to-scalar ratio to r=1.1+0.9-0.8 (r<2.8 at 95% C.L.) as derived by the ML pipeline, and r=1.2+0.9-0.8 (r<2.7 at 95% C.L.) as derived by the PCL pipeline. In one of the fields, we find a correlation with the dust component of the Planck Sky Model, though the corresponding excess power is small compared to statistical errors. Finally, we derive limits on all known systematic errors, and demonstrate that these correspond to a tensor-to-scalar ratio smaller than r=0.01, the lowest level yet reported in the literature.

preprint2013arXiv

The Cosmic Evolution of Fermi BL Lacertae Objects

Fermi has provided the largest sample of gamma-ray selected blazars to date. In this work we use a uniformly selected set of 211 BL Lacertae (BL Lac) objects detected by it Fermi during its first year of operation. We have obtained redshift constraints for 206 out of the 211 BL Lacs in our sample making it the largest and most complete sample of BL Lacs available in the literature. We use this sample to determine the luminosity function of BL Lacs and its evolution with cosmic time. We find that for most BL Lac classes, the evolution is positive with a space density peaking at modest redshift (z~1.2). The low-luminosity, high-synchrotron peaked (HSP) BL Lacs are an exception, showing strong negative evolution, with number density increasing for z$\lesssim$0.5. Since this rise corresponds to a drop-off in the density of flat-spectrum radio quasars (FSRQs), a possible interpretation is that these HSPs represent an accretion-starved end-state of an earlier merger-driven gas-rich phase. We additionally find that the known BL Lac correlation between luminosity and photon spectral index persists after correction for the substantial observational selection effects with implications for the so called `blazar sequence'. Finally, estimating the beaming corrections to the luminosity function, we find that BL Lacs have an average Lorentz factor of $γ=6.1^{+1.1}_{-0.8}$, and that most are seen within 10$^{\circ}$ of the jet axis.

preprint2013arXiv

The ordinary life of the gamma-ray emitting narrow-line Seyfert 1 galaxy PKS 1502+036

We report on multifrequency observations of the gamma-ray emitting narrow-line Seyfert 1 galaxy PKS 1502+036 performed from radio to gamma-rays during 2008 August-2012 November by Fermi-LAT, Swift (XRT and UVOT), OVRO, VLBA, and VLA. No significant variability has been observed in gamma-rays, with 0.1-100 GeV flux that ranged between (3-7)x10^-8 ph/cm^2/s using 3-month time bins. The photon index of the LAT spectrum (Gamma=2.60+/-0.06) and the apparent isotropic gamma-ray luminosity, L(0.1-100 GeV)= 7.8x10^45 erg/s, over 51 months are typical of a flat spectrum radio quasar. The radio spectral variability and the one-sided structure, in addition to the observed gamma-ray luminosity, suggest a relativistic jet with a high Doppler factor. In contrast to SBS 0846+513, the VLBA at 15 GHz did not observe superluminal motion for PKS 1502+036. Despite having the optical characteristics typical of a narrow-line Seyfert 1 galaxy, radio and gamma-ray properties of PKS 1502+036 are found to be similar to those of a blazar at the low end of the black hole mass distribution for blazars. This is in agreement with what has been found in the case of the other gamma-ray emitting narrow-line Seyfert 1 SBS 0846+513.

preprint2013arXiv

VLBA observations of a rare multiple quasar imaging event caused by refraction in the interstellar medium

We report on the first detection of the theoretically-predicted rare phenomenon of multiple parsec-scale imaging of an active galactic nucleus induced by refractive effects due to localized foreground electron density enhancements, e.g., in an AU-scale plasma lens(es) in the ionized component of the Galactic interstellar medium. We detected multiple imaging in the low galactic latitude (b=-2 deg) quasar 2023+335 from the 15.4 GHz MOJAVE observations when the source was undergoing an ESE. While the parsec-scale jet of the source normally extends along PA -20 deg, in the 28 May 2009 and 23 July 2009 images a highly significant multi-component pattern of secondary images is stretched out nearly along the constant galactic latitude line with a local PA 40 deg, indicating that the direction of relative motion of the plasma lens is close to orbital. Weaker but still detectable imaging patterns at similar position angles are sporadically manifest at several other epochs. Modeling the ESE that occurred in early 2009 and lasted ~0.14 yr, we determined that the foreground screen has a double-lens structure, with proper motion (~6.8 mas/yr), and angular size (~0.27 mas). We also found that the angular separation between the two brightest sub-images roughly follows a wavelength-squared dependence expected from plasma scattering. Furthermore, by analyzing archival non-simultaneous VLBA observations covering a wide frequency range from 1.4 to 86 GHz, we found that the scattered angular size of the VLBI core follows a $ν^{-1.89}$ dependence, implying the presence of a turbulent, refractive dominated scattering screen that has a confined structure or is truncated transverse to the line of sight toward 2023+335.

preprint2012arXiv

A Circularly Symmetric Antenna Design With High Polarization Purity and Low Spillover

We describe the development of two circularly symmetric antennas with high polarization purity and low spillover. Both were designed to be used in an all-sky polarization and intensity survey at 5 GHz (the C-Band All-Sky Survey, C-BASS). The survey requirements call for very low levels of cross-polar leakage and far-out sidelobes. Two different existing antennas, with 6.1-m and 7.6-m diameter primaries, were adapted by replacing the feed and secondary optics, resulting in identical beam performances of 0.73deg FWHM, cross-polarization better than -50 dB, and far-out sidelobes below -70 dB. The polarization purity was realized by using a symmetric low-loss dielectric foam support structure for the secondary mirror, avoiding the need for secondary support struts. Ground spill-over was largely reduced by using absorbing baffles around the primary and secondary mirrors, and by the use of a low-sidelobe profiled corrugated feedhorn. The 6.1-m antenna and receiver have been completed and test results show that the optics meet their design goals.

preprint2012arXiv

A giant radio flare from Cygnus X-3 with associated Gamma-ray emission

With frequent flaring activity of its relativistic jets, Cygnus X-3 is one of the most active microquasars and is the only Galactic black hole candidate with confirmed high energy Gamma-ray emission, thanks to detections by Fermi/LAT and AGILE. In 2011, Cygnus X-3 was observed to transit to a soft X-ray state, which is known to be associated with high-energy Gamma-ray emission. We present the results of a multi-wavelength campaign covering a quenched state, when radio emission from Cygnus X-3 is at its weakest and the X-ray spectrum is very soft. A giant (~ 20 Jy) optically thin radio flare marks the end of the quenched state, accompanied by rising non-thermal hard X-rays. Fermi/LAT observations (E >100 MeV) reveal renewed Gamma-ray activity associated with this giant radio flare, suggesting a common origin for all non-thermal components. In addition, current observations unambiguously show that the Gamma-ray emission is not exclusively related to the rare giant radio flares. A 3-week period of Gamma-ray emission is also detected when Cygnus X-3 was weakly flaring in radio, right before transition to the radio quenched state. No Gamma rays are observed during the ~ one-month long quenched state, when the radio flux is weakest. Our results suggest transitions into and out of the ultrasoft X-ray (radio quenched) state trigger Gamma-ray emission, implying a connection to the accretion process, and also that the Gamma-ray activity is related to the level of radio flux (and possibly shock formation), strengthening the connection to the relativistic jets.

preprint2012arXiv

Assessing the Significance of Apparent Correlations Between Radio and Gamma-ray Blazar Fluxes

Whether a correlation exists between the radio and gamma-ray flux densities of blazars is a long-standing question, and one that is difficult to answer confidently because of various observational biases which may either dilute or apparently enhance any intrinsic correlation between radio and gamma-ray luminosities. We introduce a novel method of data randomization to evaluate quantitatively the effect of these biases and to assess the intrinsic significance of an apparent correlation between radio and gamma-ray flux densities of blazars. The novelty of the method lies in a combination of data randomization in luminosity space (to ensure that the randomized data are intrinsically, and not just apparently, uncorrelated) and significance assessment in flux space (to explicitly avoid Malmquist bias and automatically account for the limited dynamical range in both frequencies). The method is applicable even to small samples that are not selected with strict statistical criteria. For larger samples we describe a variation of the method in which the sample is split in redshift bins, and the randomization is applied in each bin individually; this variation is designed to yield the equivalent to luminosity-function sampling of the underlying population in the limit of very large, statistically complete samples. We show that for a smaller number of redshift bins, the method yields a worse significance, and in this way it is conservative in that it does not assign a stronger, artificially enhanced significance. We demonstrate how our test performs as a function of number of sources, strength of correlation, and number of redshift bins used, and we show that while our test is robust against common-distance biases and associated false positives for uncorrelated data, it retains the power of other methods in rejecting the null hypothesis of no correlation for correlated data.

preprint2012arXiv

First Season QUIET Observations: Measurements of CMB Polarization Power Spectra at 43 GHz in the Multipole Range 25 <= ell <= 475

The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43GHz and 95GHz, operating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the polarization of the CMB. QUIET primarily targets the B modes from primordial gravitational waves. The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galactic synchrotron radiation. Between 2008 October and 2010 December, >10,000hours of data were collected, first with the 19-element 43GHz array (3458hours) and then with the 90-element 95GHz array. Each array observes the same four fields, selected for low foregrounds, together covering ~1000deg^2. This paper reports initial results from the 43GHz receiver which has an array sensitivity to CMB fluctuations of 69uK sqrt(s). The data were extensively studied with a large suite of null tests before the power spectra, determined with two independent pipelines, were examined. Analysis choices, including data selection, were modified until the null tests passed. Cross correlating maps with different telescope pointings is used to eliminate a bias. This paper reports the EE, BB and EB power spectra in the multipole range ell=25-475. With the exception of the lowest multipole bin for one of the fields, where a polarized foreground, consistent with Galactic synchrotron radiation, is detected with 3sigma significance, the E-mode spectrum is consistent with the LCDM model, confirming the only previous detection of the first acoustic peak. The B-mode spectrum is consistent with zero, leading to a measurement of the tensor-to-scalar ratio of r=0.35+1.06-0.87. The combination of a new time-stream double-demodulation technique, Mizuguchi-Dragone optics, natural sky rotation, and frequent boresight rotation leads to the lowest level of systematic contamination in the B-mode power so far reported, below the level of r=0.1

preprint2012arXiv

OVRO 40m blazar monitoring program: Understanding the relationship between 15 GHz radio variability properties and gamma-ray activity in blazars

A large sample of known and likely gamma-ray blazars has been monitored twice per week since late 2007 at 15 GHz with the Owens Valley Radio Observatory (OVRO) 40-meter Telescope. The sample contains about 1700 sources, including the initial sample of 1158 sources above declination -20 degrees from the Candidate Gamma-Ray Blazar Survey (CGRaBS) plus all the blazars associated with Fermi-LAT detections as released in the Fermi AGN catalogs. Using statistical likelihood analyses, we compare the variability amplitude for various sub-populations within our sample. These include comparisons of gamma-ray-loud versus quiet objects, BL Lac objects versus flat-spectrum radio quasars, and a study of the variability amplitude trend with redshift. To learn about the location of the gamma-ray emission region we study the significance of peaks in the radio/gamma-ray cross-correlation using Monte Carlo simulations. First results for 52 sources with data from both the high-confidence Fermi Large Area Telescope Bright AGN Sample and the first 2 years of our monitoring program are presented. We find that assuming a power spectral density with power law slope of -2 at 15 GHz and -1.5 at gamma-ray energies, 7 of our objects show cross-correlations at the 3sigma level. We are now studying the physical significance of these correlations by further exploring the range of power law slopes that are consistent with the data. An extension of this to a larger sample and longer light curves is underway and preliminary results are presented. We also describe KuPol, the new digital Ku-band receiver being constructed for the 40-meter telescope. This new receiver will provide total intensity and linear polarization measurements over the 12-18 GHz band, with 16 MHz spectral resolution. The polarization data will provide important clues about the magnetic field configuration in the radio emission region.

preprint2012arXiv

Parsec-Scale Radio Properties of Gamma-Ray Emitting Blazars

The parsec-scale radio properties of blazars detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been investigated using observations with the Very Long Baseline Array (VLBA). Comparisons between LAT and non-LAT detected samples were made using contemporaneous data. In total, 232 sources were used in the LAT-detected sample. This very large, radio flux-limited sample of active galactic nuclei (AGN) provides insights into the mechanism that produces strong gamma-ray emission. It has been found that LAT-detected BL Lac objects are very similar to the non-LAT BL Lac objects in most properties, although LAT BL Lac objects may have longer jets. The LAT flat spectrum radio quasars (FSRQs) are significantly different from non-LAT FSRQs and are likely extreme members of the FSRQ population. Contemporaneous observations showed a strong correlation, whereas no correlation is found using archival radio data. Most of the differences between the LAT and non-LAT populations are related to the cores of the sources, indicating that the gamma-ray emission may originate near the base of the jets (i.e., within a few pc of the central engine). There is some indication that LAT-detected sources may have larger jet opening angles than the non-LAT sources. Strong core polarization is significantly more common among the LAT sources, suggesting that gamma-ray emission is related to strong, uniform magnetic fields at the base of the jets of the blazars. Observations of sources in two epochs indicate that core fractional polarization was higher when the objects were detected by the LAT. Included in our sample are several non-blazar AGN such as 3C84, M82, and NGC 6251.

preprint2012arXiv

Radio and gamma-ray follow-up of the exceptionally high activity state of PKS 1510-089 in 2011

We investigate the radio and gamma-ray variability of the flat spectrum radio quasar PKS 1510-089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major gamma-ray flares detected in 2011 July and October. During the latter episode both the gamma-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380 deg of the optical polarization angle close in time with the rapid and strong gamma-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands is observed about three weeks later, close in time with another gamma-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimeter regime followed with some time delay at centimeter down to decimeter wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the gamma-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the gamma-ray variability is not within the broad line, but a few parsecs downstream along the jet.

preprint2012arXiv

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 gamma-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 gamma-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, able to develop an isotropic luminosity at gamma-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 (gamma-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 gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either 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 gamma-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.

preprint2012arXiv

SBS 0846+513: a new gamma-ray emitting Narrow-Line Seyfert 1 galaxy

We report Fermi-LAT observations of the radio-loud AGN SBS 0846+513 (z=0.5835), optically classified as a Narrow-Line Seyfert 1 galaxy, together with new and archival radio-to-X-ray data. The source was not active at gamma-ray energies during the first two years of Fermi operation. A significant increase in activity was observed during 2010 October-2011 August. In particular a strong gamma-ray flare was observed in 2011 June reaching an isotropic gamma-ray luminosity (0.1-300 GeV) of 1.0x10^48 erg/s, comparable to that of the brightest flat spectrum radio quasars, and showing spectral evolution in gamma rays. An apparent superluminal velocity of (8.2+/-1.5)c in the jet was inferred from 2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. Both the power released by this object during the flaring activity and the apparent superluminal velocity are strong indications of the presence of a relativistic jet as powerful as those of blazars. In addition, variability and spectral properties in radio and gamma-ray bands indicate blazar-like behaviour, suggesting that, except for some distinct optical characteristics, SBS 0846+513 could be considered as a young blazar at the low end of the blazar's black hole mass distribution.

preprint2012arXiv

The QUIET Instrument

The Q/U Imaging ExperimenT (QUIET) is designed to measure polarization in the Cosmic Microwave Background, targeting the imprint of inflationary gravitational waves at large angular scales (~ 1 degree). Between 2008 October and 2010 December, two independent receiver arrays were deployed sequentially on a 1.4 m side-fed Dragonian telescope. The polarimeters which form the focal planes use a highly compact design based on High Electron Mobility Transistors (HEMTs) that provides simultaneous measurements of the Stokes parameters Q, U, and I in a single module. The 17-element Q-band polarimeter array, with a central frequency of 43.1 GHz, has the best sensitivity (69 uK sqrt(s)) and the lowest instrumental systematic errors ever achieved in this band, contributing to the tensor-to-scalar ratio at r < 0.1. The 84-element W-band polarimeter array has a sensitivity of 87 uK sqrt(s) at a central frequency of 94.5 GHz. It has the lowest systematic errors to date, contributing at r < 0.01. The two arrays together cover multipoles in the range l= 25-975. These are the largest HEMT-based arrays deployed to date. This article describes the design, calibration, performance of, and sources of systematic error for the instrument.

preprint2011arXiv

Dust-correlated cm-wavelength continuum emission on translucent clouds ζ Oph and LDN 1780

The diffuse cm-wave IR-correlated signal, the "anomalous" CMB foreground, is thought to arise in the dust in cirrus clouds. We present Cosmic Background Imager (CBI) cm-wave data of two translucent clouds, ζ Oph and LDN 1780 with the aim of characterising the anomalous emission in the translucent cloud environment. In ζ Oph, the measured brightness at 31 GHz is 2.4σ higher than an extrapolation from 5 GHz measurements assuming a free-free spectrum on 8 arcmin scales. The SED of this cloud on angular scales of 1{\odot} is dominated by free-free emission in the cm-range. In LDN 1780 we detected a 3 σ excess in the SED on angular scales of 1{\odot} that can be fitted using a spinning dust model. In this cloud, there is a spatial correlation between the CBI data and IR images, which trace dust. The correlation is better with near-IR templates (IRAS 12 and 25 μm) than with IRAS 100 μm, which suggests a very small grain origin for the emission at 31 GHz. We calculated the 31 GHz emissivities in both clouds. They are similar and have intermediate values between that of cirrus clouds and dark clouds. Nevertheless, we found an indication of an inverse relationship between emissivity and column density, which further supports the VSGs origin for the cm-emission since the proportion of big relative to small grains is smaller in diffuse clouds.

preprint2011arXiv

Gamma-Ray Emission from Two Blazars Behind the Galactic Plane: B2013+370 & B2023+336

B2013+370 and B2023+336 are two blazars at low-galactic latitude that were previously proposed to be the counterparts for the EGRET unidentified sources, 3EG J2016+3657 and 3EG J2027+3429. Gamma-ray emission associated with the EGRET sources has been detected by the Fermi Gamma-ray Space Telescope, and the two sources, 1FGL J2015.7+3708 and 1FGL J2027.6+3335, have been classified as unidentified in the 1-year catalog. This analysis of the Fermi-LAT data collected during 31 months reveals that the 1FGL sources are spatially compatible with the blazars, and are significantly variable, supporting the hypothesis of extragalactic origin for the gamma-ray emission. The gamma-ray light curves are compared with 15 GHz radio light curves from the 40-m telescope at the Owens Valley Radio Observatory (OVRO). Simultaneous variability is seen in both bands for the two blazar candidates. The study is completed with the X-ray analysis of 1FGL J2015.7+3708 using Swift observations that were triggered in August 2010 by a Fermi-detected flare. The resulting spectral energy distribution shows a two-component structure typical of blazars. We also identify a second source in the field of view of 1FGL J2027.6+3335 with similar characteristics to the known LAT pulsars. This study gives solid evidence favoring blazar counterparts for these two unidentified EGRET and Fermi sources, supporting the hypothesis that a number of unidentified gamma-ray sources at low galactic latitudes are indeed of extragalactic origin.

preprint2011arXiv

The Birth of a Relativistic Outflow in the Unusual γ-ray Transient Swift J164449.3+573451

Active galactic nuclei (AGN), powered by long-term accretion onto central supermassive black holes, produce relativistic jets with lifetimes of greater than one million yr that preclude observations at birth. Transient accretion onto a supermassive black hole, for example through the tidal disruption of a stray star, may therefore offer a unique opportunity to observe and study the birth of a relativistic jet. On 2011 March 25, the Swift γ-ray satellite discovered an unusual transient source (Swift J164449.3+573451) potentially representing such an event. Here we present the discovery of a luminous radio transient associated with Swift J164449.3+573451, and an extensive set of observations spanning centimeter to millimeter wavelengths and covering the first month of evolution. These observations lead to a positional coincidence with the nucleus of an inactive galaxy, and provide direct evidence for a newly-formed relativistic outflow, launched by transient accretion onto a million solar mass black hole. While a relativistic outflow was not predicted in this scenario, we show that the tidal disruption of a star naturally explains the high-energy properties, radio luminosity, and the inferred rate of such events. The weaker beaming in the radio compared to γ-rays/X-rays, suggests that radio searches may uncover similar events out to redshifts of z ~ 6.

preprint2011arXiv

The first gamma-ray outburst of a Narrow-Line Seyfert 1 Galaxy: the case of PMN J0948+0022 in July 2010

We report on a multiwavelength campaign on the radio-loud Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.5846) performed in 2010 July-September and triggered by a high-energy gamma-ray outburst observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. The peak flux in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^-6 ph cm^-2 s^-1, corresponding to an observed luminosity of 10^48 erg s^-1. Although the source was too close to the Sun position to organize a densely sampled follow-up, it was possible to gather some multiwavelength data that confirmed the state of high activity across the sampled electromagnetic spectrum. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

preprint2011arXiv

The July 2010 outburst of the NLS1 PMN J0948+0022

We report about the multiwavelength campaign on the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z = 0.5846) performed in 2010 July-September and triggered by high activity as measured by Fermi/LAT. The peak luminosity in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^48 erg/s, a level comparable to the most powerful blazars. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

preprint2010arXiv

Characteristics of Gamma-Ray Loud Blazars in the VLBA Imaging and Polarimetry Survey

The radio properties of blazars detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been observed as part of the VLBA Imaging and Polarimetry Survey (VIPS). This large, flux-limited sample of active galactic nuclei (AGN) provides insights into the mechanism that produces strong gamma-ray emission. At lower flux levels, radio flux density does not directly correlate with gamma-ray flux. We find that the LAT-detected BL Lacs tend to be similar to the non-LAT BL Lacs, but that the LAT-detected FSRQs are often significantly different from the non-LAT FSRQs. The differences between the gamma-ray loud and quiet FSRQs can be explained by Doppler boosting; these objects appear to require larger Doppler factors than those of the BL Lacs. It is possible that the gamma-ray loud FSRQs are fundamentally different from the gamma-ray quiet FSRQs. Strong polarization at the base of the jet appears to be a signature for gamma-ray loud AGN.

preprint2010arXiv

IR-correlated 31 GHz radio emission from Orion East

Lynds dark cloud LDN1622 represents one of the best examples of anomalous dust emission, possibly originating from small spinning dust grains. We present Cosmic Background Imager (CBI) 31 GHz data of LDN1621, a diffuse dark cloud to the north of LDN1622 in a region known as Orion East. A broken ring with diameter g\approx 20 arcmin of diffuse emission is detected at 31 GHz, at \approx 20-30 mJy beam$^{-1}$ with an angular resolution of \approx 5 arcmin. The ring-like structure is highly correlated with Far Infra-Red emission at $12-100 μ$m with correlation coefficients of r \approx 0.7-0.8, significant at $\sim10σ$. Multi-frequency data are used to place constraints on other components of emission that could be contributing to the 31 GHz flux. An analysis of the GB6 survey maps at 4.85 GHz yields a $3σ$ upper limit on free-free emission of 7.2 mJy beam$^{-1}$ ($\la 30 per cent of the observed flux) at the CBI resolution. The bulk of the 31 GHz flux therefore appears to be mostly due to dust radiation. Aperture photometry, at an angular resolution of 13 arcmin and with an aperture of diameter 30 arcmin, allowed the use of IRAS maps and the {\it WMAP} 5-year W-band map at 93.5 GHz. A single modified blackbody model was fitted to the data to estimate the contribution from thermal dust, which amounts to $\sim$ 10 per cent at 31 GHz. In this model, an excess of 1.52\pm 0.66 Jy (2.3σ) is seen at 31 GHz. Future high frequency $\sim$ 100-1000 GHz data, such as those from the {\it Planck} satellite, are required to accurately determine the thermal dust contribution at 31 GHz. Correlations with the IRAS $100 μ$m gave a coupling coefficient of $18.1\pm4.4 μ$K (MJy/sr)$^{-1}$, consistent with the values found for LDN1622.

preprint2010arXiv

J16021+3326: New Multi-Frequency Observations of a Complex Source

We present multifrequency Very Long Baseline Array (VLBA) observations of J16021+3326. These observations, along with variability data obtained from the Owens Valley Radio Observatory (OVRO) candidate gamma-ray blazar monitoring program, clearly indicate this source is a blazar. The peculiar characteristic of this blazar, which daunted previous classification attempts, is that we appear to be observing down a precessing jet, the mean orientation of which is aligned with us almost exactly.

preprint2010arXiv

Magic constraints on Gamma-ray emission from Cygnus X-3

Cygnus X-3 is a microquasar consisting of an accreting compact object orbiting around a Wolf-Rayet star. It has been detected at radio frequencies and up to high-energy gamma rays (above 100 MeV). However, many models also predict a very high energy (VHE) emission (above hundreds of GeV) when the source displays relativistic persistent jets or transient ejections. Therefore, detecting such emission would improve the understanding of the jet physics. The imaging atmospheric Cherenkov telescope MAGIC observed Cygnus X-3 for about 70 hours between 2006 March and 2009 August in different X-ray/radio spectral states and also during a period of enhanced gamma-ray emission. MAGIC found no evidence for a VHE signal from the direction of the microquasar. An upper limit to the integral flux for energies higher than 250 GeV has been set to 2.2 x 10-12 photons cm-2 s-1 (95% confidence level). This is the best limit so far to the VHE emission from this source. The non-detection of a VHE signal during the period of activity in the high-energy band sheds light on the location of the possible VHE radiation favoring the emission from the innermost region of the jets, where absorption is significant. The current and future generations of Cherenkov telescopes may detect a signal under precise spectral conditions.

preprint2009arXiv

15 GHz Monitoring of Gamma-ray Blazars with the OVRO 40 Meter Telescope in Support of Fermi

We present results from the first two years of our fast-cadence 15 GHz gamma-ray blazar monitoring program, part of the F-GAMMA radio monitoring project. Our sample includes the 1158 blazars north of -20 degrees declination from the Candidate Gamma-Ray Blazar Survey (CGRaBS), which encompasses a significant fraction of the extragalactic sources detected by the Fermi Gamma-ray Space Telescope. We introduce a novel likelihood analysis for computing a time series variability amplitude statistic that separates intrinsic variability from measurement noise and produces a quantitative error estimate. We use this method to characterize our radio light curves. We also present results indicating a statistically significant correlation between simultaneous average 15 GHz radio flux density and gamma-ray photon flux.

preprint2009arXiv

A 31 GHz Survey of Low-Frequency Selected Radio Sources

The 100-m Robert C. Byrd Green Bank Telescope (GBT) and the Owens Valley Radio Observatory (OVRO) 40-m radio telescope have been used to conduct a survey of 3165 known extragalactic radio sources over 143 square degrees of the sky. Target sources were selected from the NRAO VLA Sky Survey in fields observed by the Cosmic Background Imager (CBI); most are extragalactic active galactic nuclei (AGN) with 1.4 GHz flux densities of 3 to 10 mJy. The resulting 31 GHz catalogs are presented in full online. Using a Maximum-Likelihood analysis to obtain an unbiased estimate of the distribution of the 1.4 to 31 GHz spectral indices of these sources, we find a mean 31 to 1.4 GHz flux ratio of 0.110 +/- 0.003 corresponding to a spectral index of alpha=-0.71 +/- 0.01 (S ~ nu^alpha); 9.0 +/- 0.8 % of sources have alpha > -0.5 and 1.2 +/- 0.2 % have alpha > 0. By combining this spectral index distribution with 1.4 GHz source counts we predict 31 GHz source counts in the range 1 mJy < S_31 < 4 mJy, N(>S_31) = (16.7 +/- 1.7) deg^2 (S_31/1 mJy)^(-0.80 +/- 0.07). We also assess the contribution of mJy-level (S_1.4 < 3.4 mJy) radio sources to the 31 GHz CMB power spectrum, finding a mean power of ell (ell+1) C^src_ell/(2 pi) = 44 +/- 14 micro-Kelvin^2 and a 95% upper limit of 80 micro-Kelvin^2 at ell = 2500. Including an estimated contribution of 12 micro-Kelvin^2 from the population of sources responsible for the turn-up in counts below S_1.4 = 1 mJy this amounts to 21 +/- 7 % of what is needed to explain the CBI high-ell excess signal, 275 +/- 63 micro-Kelvin^2. These results are consistent with other measurements of the 31 GHz point source foreground.

preprint2009arXiv

Multi-frequency measurements of the NVSS foreground sources in the Cosmic Background Imager fields. I. Data release

We present the results of the flux density measurements at 4.85 GHz and 10.45 GHz of a sample of 5998 NVSS radio sources with the Effelsberg 100 m telescope. The initial motivation was the need to identify the NVSS radio sources that could potentially contribute significant contaminating flux in the frequency range at which the Cosmic Background Imager experiment operated. An efficient way to achieve this challenging goal has been to compute the high frequency flux density of those sources by extrapolating their radio spectrum. This is determined by the three-point spectral index measured on the basis of the NVSS entry at 1.4 GHz and the measurements at 4.85 GHz and 10.45 GHz carried out with the 100 m Effelsberg telescope. These measurements are important since the targeted sample probes the weak part of the flux density distribution, hence the decision to make the data available. We present the table with flux density measurements of 3434 sources that showed no confusion allowing reliable measurements, their detection rates, their spectral index distribution and an interpretation which explains satisfactorily the observed uncertainties.

preprint2008arXiv

Anomalous Microwave Emission from the HII region RCW175

We present evidence for anomalous microwave emission in the RCW175 \hii region. Motivated by 33 GHz $13\arcmin$ resolution data from the Very Small Array (VSA), we observed RCW175 at 31 GHz with the Cosmic Background Imager (CBI) at a resolution of $4\arcmin$. The region consists of two distinct components, G29.0-0.6 and G29.1-0.7, which are detected at high signal-to-noise ratio. The integrated flux density is $5.97\pm0.30$ Jy at 31 GHz, in good agreement with the VSA. The 31 GHz flux density is $3.28\pm0.38$ Jy ($8.6σ$) above the expected value from optically thin free-free emission based on lower frequency radio data and thermal dust constrained by IRAS and WMAP data. Conventional emission mechanisms such as optically thick emission from ultracompact \hii regions cannot easily account for this excess. We interpret the excess as evidence for electric dipole emission from small spinning dust grains, which does provide an adequate fit to the data.

preprint2006arXiv

Implications of the Cosmic Background Imager Polarization Data

We present new measurements of the power spectra of the E-mode of CMB polarization, the temperature T, the cross-correlation of E and T, and upper limits on the B-mode from 2.5 years of dedicated Cosmic Background Imager (CBI) observations. Both raw maps and optimal signal images in the uv-plane and real space show strong detections of the E-mode (11.7 sigma for the EE power spectrum overall) and no detection of the B-mode. The power spectra are used to constrain parameters of the flat tilted adiabatic Lambda-CDM models: those determined from EE and TE bandpowers agree with those from TT, a powerful consistency check. There is little tolerance for shifting polarization peaks from the TT-forecast locations, as measured by the angular sound crossing scale theta = 100 ell_s = 1.03 +/- 0.02 from EE and TE cf. 1.044 +/- 0.005 with the TT data included. The scope for extra out-of-phase peaks from subdominant isocurvature modes is also curtailed. The EE and TE measurements of CBI, DASI and BOOMERANG are mutually consistent, and, taken together rather than singly, give enhanced leverage for these tests.

preprint1997arXiv

An Anomalous Component of Galactic Emission

We present results from microwave background observations at the Owens Valley Radio Observatory. These observations, at 14.5 and 32 GHz, are designed to detect intrinsic anisotropy on scales of 7'. After point source removal, we detect significant emission with temperature spectral index beta ~ -2 towards the North Celestial Pole (NCP). Comparison of our data with the IRAS 100 micron map of the same fields reveals a strong correlation between this emission and the infrared dust emission. From the lack of detectable H-alpha emission, we conclude that the signals are consistent either with flat-spectrum synchrotron radiation, or with free-free emission from T_e ~ 10^6 K gas, probably associated with a large HI feature known as the NCP Loop. Assuming beta = -2.2, our data indicate a conversion T_f/I_(100 micron) = 0.075*nu(GHz)^-2.2 K/(MJy/sr). The detection of such a component suggests that we should be cautious in any assumptions made regarding foregrounds when designing experiments to map the microwave background radiation.