Vasiliki Pavlidou

BLAZARS IN THE FERMI ERA: THE OVRO 40 m TELESCOPE MONITORING PROGRAM

The Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope provides an unprecedented opportunity to study gamma-ray blazars. To capitalize on this opportunity, beginning in late 2007, about a year before the start of LAT science operations, we began a large-scale, fast-cadence 15 GHz radio monitoring program with the 40 m telescope at the Owens Valley Radio Observatory. This program began with the 1158 northern (δ > –20°) sources from the Candidate Gamma-ray Blazar Survey and now encompasses over 1500 sources, each observed twice per week with about 4 mJy (minimum) and 3% (typical) uncertainty. Here, we describe this monitoring program and our methods, and present radio light curves from the first two years (2008 and 2009). As a first application, we combine these data with a novel measure of light curve variability amplitude, the intrinsic modulation index, through a likelihood analysis to examine the variability properties of subpopulations of our sample. We demonstrate that, with high significance (6σ), gamma-ray-loud blazars detected by the LAT during its first 11 months of operation vary with almost a factor of two greater amplitude than do the gamma-ray-quiet blazars in our sample. We also find a significant (3σ) difference between variability amplitude in BL Lacertae objects and flat-spectrum radio quasars (FSRQs), with the former exhibiting larger variability amplitudes. Finally, low-redshift (z < 1) FSRQs are found to vary more strongly than high-redshift FSRQs, with 3σ significance. These findings represent an important step toward understanding why some blazars emit gamma-rays while others, with apparently similar properties, remain silent.

Cosmic Gamma-Ray Background from Star-Forming Galaxies

The origin of the extragalactic gamma-ray background is a pressing cosmological mystery. The Fermi Gamma-Ray Space Telescope has recently measured the intensity and spectrum of this background; both are substantially different from previous measurements. We present a novel calculation of the gamma-ray background from normal star-forming galaxies. Contrary to long-standing expectations, we find that numerous but individually faint normal galaxies may comprise the bulk of the Fermi signal, rather than rare but intrinsically bright active galaxies. This result has wide-ranging implications, including: the possibility to probe the cosmic star formation history with gamma rays; the ability to infer the cosmological evolution of cosmic rays and galactic magnetic fields; and an increased likelihood of identifying subdominant components from rare sources (e.g., dark matter clumps) through their large anisotropy.

Diffuse Gamma Rays from Local Group Galaxies

Diffuse γ-ray radiation in galaxies is produced by cosmic-ray interactions with the interstellar medium. With the completion of EGRET observations, the only extragalactic object from which there has been a positive detection of diffuse γ-ray emission is the Large Magellanic Cloud. We systematically estimate the expected diffuse γ-ray flux from Local Group galaxies and determine their detectability by a new generation of γ-ray observatories such as the Gamma-Ray Large Area Space Telescope (GLAST). For each galaxy, the expected γ-ray flux depends only on its total gas content and its cosmic-ray flux. We present a method for calculating cosmic-ray flux in these galaxies in terms of the observed rate of supernova explosions, in which cosmic-ray acceleration is believed to take place. The difficulty in deriving accurate supernova rates from observational data is a dominant uncertainty in our calculations. We estimate the γ-ray flux for Local Group galaxies and find that our predictions are consistent with the observations for the LMC and with the observational upper limits for the SMC and M31. Both the Andromeda galaxy, with a flux of ~1.0 × 10-8 photons s-1 cm-2 above 100 MeV, and the SMC, with a flux of ~1.7 × 10-8 photons s-1 cm-2 above 100 MeV, are expected to be observable by GLAST. M33 is at the limit of detectability with a flux of ~0.11 × 10-8 s-1 cm-2. Other Local Group galaxies are at least 2 orders of magnitude below GLAST sensitivity.

Imprint of galaxy clustering in the cosmic gamma-ray background

Star-forming galaxies are predicted to contribute considerably to the cosmic gamma-ray background (CGB) as they are confirmed γ-ray emitters and are the most numerous population of γ-ray sources, although individually faint. Even though the Fermi Gamma-ray Space Telescope will be able to resolve few star-forming galaxies individually, their fractional contribution to the CGB should become far more significant than it was for past measurements of the CGB as many of the brighter, formerly unresolved sources will be resolved out. Thus, the clustering feature of galaxies imprinted on the CGB might be detectable by Fermi. In anticipation of such measurements, we calculate the predicted angular auto- and cross-power spectra of the CGB from normal galaxies. We find that the amplitude of the auto-power spectrum is smaller than that for other sources, such as blazars and dark matter annihilation; the shape is also characteristic. We also show that the cross-power spectrum with galaxy surveys features larger amplitude. Fermi should be able to detect the correlation signature in both the auto- and cross-power spectra at angular scales of ~1°–10° after 5 years of operation. Such a detection would be valuable in confirming the level of the star-forming galaxy contribution to the CGB and, more importantly, in serving as a tool in the effort to discriminate between possible origins of the CGB.

Anisotropies in the gamma‐ray sky from millisecond pulsars

Pulsars emerge in the Fermi era as a sizable population of gamma-ray sources. Millisecond pulsars (MSPs) constitute an older subpopulation whose sky distribution extends to high Galactic latitudes, and it has been suggested that unresolved members of this class may contribute a significant fraction of the measured large-scale isotropic gamma-ray background (IGRB). We investigate the possible energy-dependent contribution of unresolved MSPs to the anisotropy of the Fermi-measured IGRB. For observationally motivated MSP population models, we show that the preliminary Fermi anisotropy measurement places an interesting constraint on the abundance of MSPs in the Galaxy and the typical MSP flux, about an order of magnitude stronger than constraints on this population derived from the intensity of the IGRB alone. We also examine the possibility of an MSP component in the IGRB mimicking a dark matter signal in anisotropy-based searches, and conclude that the energy dependence of an anisotropy signature would distinguish MSPs from all but very light dark matter candidates.

Robust identification of isotropic diffuse gamma rays from Galactic dark matter

Dark matter annihilation in Galactic substructure will produce diffuse gamma-ray emission of remarkably constant intensity across the sky, making it difficult to disentangle this Galactic dark matter signal from the extragalactic gamma-ray background. We show that if Galactic dark matter contributes a modest fraction of the measured emission in an energy range accessible to the Fermi Gamma-ray Space Telescope, the energy dependence of the angular power spectrum of the total measured emission could be used to confidently identify gamma rays from Galactic dark matter substructure.

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 realized in nature.

Connection between optical and γ-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 γ-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 γ-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 γ-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 γ-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.

THE DETECTABILITY OF DARK MATTER ANNIHILATION WITH FERMI USING THE ANISOTROPY ENERGY SPECTRUM OF THE GAMMA-RAY BACKGROUND

The energy dependence of the anisotropy (the anisotropy energy spectrum) of the large-scale diffuse gamma-ray background can reveal the presence of multiple source populations. Annihilating dark matter in the substructure of the Milky Way halo could give rise to a modulation in the anisotropy energy spectrum of the diffuse gamma-ray emission measured by Fermi, enabling the detection of a dark matter signal. We determine the detectability of a dark-matter-induced modulation for scenarios in which unresolved blazars are the primary contributor to the measured emission above ~1 GeV and find that in some scenarios pair-annihilation cross sections on the order of the value expected for thermal relic dark matter can produce a detectable feature. We anticipate that the sensitivity of this technique to specific dark matter models could be improved by tailored likelihood analysis methods.

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{\sigma} 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.