A. Reimer

Simultaneous Planck, Swift, and Fermi observations of X-ray and γ-ray selected blazars

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and -ray bands, and we compare our results to those of a companion paper presenting simultaneous Planck and multi-frequency observations of 104 radio-loud northern active galactic nuclei selected at radio frequencies. While we confirm several previous results, our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi Large Area Telescope (LAT), whereas 30 to 40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the -ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-millimetre spectral slope of blazars is quite flat, withh i 0 up to about 70 GHz, above which it steepens toh i 0:65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency ( S ) in the spectral energy distribution (SED) of FSRQs is the same in all the blazar samples withh S i = 10 13:1 0:1 Hz, while the mean inverse-Compton peak frequency,h IC i, ranges from 10 21 to 10 22 Hz. The distributions of S and of IC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs; their shapes strongly depend on the selection method. The Compton dominance of blazars ranges from less than 0.2 to nearly 100, with only FSRQs reaching values larger than about 3. Its distribution is broad and depends strongly on the selection method, with -ray selected blazars peaking at 7 or more, and radio-selected blazars at values close to 1, thus implying that the common assumption that the blazar power budget is largely dominated by high-energy emission is a selection e ect. A comparison of our multi-frequency data with theoretical predictions shows that simple homogeneous SSC models cannot explain the simultaneous SEDs of most of the -ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and S predicted by the blazar sequence.

The Luminosity Function of Fermi-detected Flat-Spectrum Radio Quasars

Fermi has provided the largest sample of γ-ray-selected blazars to date. In this work we use a complete sample of flat spectrum radio quasars (FSRQs) detected during the first year of operation to determine the luminosity function (LF) and its evolution with cosmic time. The number density of FSRQs grows dramatically up to redshift ~0.5-2.0 and declines thereafter. The redshift of the peak in the density is luminosity dependent, with more luminous sources peaking at earlier times; thus the LF of γ-ray FSRQs follows a luminosity-dependent density evolution similar to that of radio-quiet active galactic nuclei. Also, using data from the Swift Burst Alert Telescope we derive the average spectral energy distribution (SED) of FSRQs in the 10 keV-300 GeV band and show that there is no correlation between the luminosity at the peak of the γ-ray emission component and its peak frequency. Using this luminosity-independent SED with the derived LF allows us to predict that the contribution of FSRQs to the Fermi isotropic γ-ray background is 9.3^(+1.6)_(–1.0%)(±3% systematic uncertainty) in the 0.1-100 GeV band. Finally we determine the LF of unbeamed FSRQs, finding that FSRQs have an average Lorentz factor of γ = 11.7^(+3.3)_(– 2.2), that most are seen within 5° of the jet axis, and that they represent only ~0.1% of the parent population.

NONTHERMAL HIGH-ENERGY EMISSION FROM COLLIDING WINDS OF MASSIVE STARS

Colliding winds of massive star binary systems are considered as potential sites of nonthermal high-energy photon production. Motivated by the detection of synchrotron radio emission from the colliding wind location, we here investigate the properties of high-energy photon production in colliding winds of long-period WR+OB systems. Analytical formulae for the steady state proton- and electron-particle spectra are derived assuming diffusive particle accelerationoutofapoolofthermalwindparticles,takinginto account adiabaticand allrelevantradiative losses,and includeadvection/convection out of the wind collision zone. This includes analytical approximations for the electron energy losses in the Klein-Nishina transition regime. For the first time in the context of CWB systems, our calculations use the full Klein-Nishina cross section and account for the anisotropy of the inverse Compton scattering process. This leads to orbitalfluxvariations byupto several orders of magnitude, which may, however,beblurredby the system’s geometry. Both anisotropy and Klein-Nishina effects may yield characteristic spectral and variability signatures in the � -ray domain. Since propagation effects lead to a deficit of low-energy particles in the convectiondominated zone, one expects imprints in the radiation spectra. If protons are accelerated to at least several GeV, � 0 -decay � -rays might be observable, depending on the injected electron-to-proton ratio. We show that photon-photon pair production is generally not negligible, potentially affecting the emitted spectrum above � 50 GeV, depending on orbital phase and system inclination. The calculations are applied to the archetypal WR+OB systems WR 140 and WR 147 to predict their expected spectral and temporal characteristics and to assess their detectability with current and upcoming � -ray experiments. Subject headingg binaries: general — gamma rays: theory — radiation mechanisms: nonthermal — stars: early-type — stars: winds, outflows

The time-dependent one-zone hadronic model - First principles

We present a time-dependent approach to the one-zone hadronic model in the case where the photon spectrum is produced by ultrarelativistic protons interacting with soft photons that are produced from protons and low magnetic fields. Assuming that protons are injected at a certain rate in a homogeneous spherical volume containing a magnetic field, the evolution of the system can be described by five coupled kinetic equations, for protons, electrons, photons, neutrons, and neutrinos. Photopair and photopion interactions are modelled using the results of Monte-Carlo simulations and, in particular from the SOPHIA code for the latter. The coupling of energy losses and injection introduces a self-consistency in our approach and allows the study of the comparative relevancy of processes at various conditions, the efficiency of the conversion of proton luminosity to radiation, the resulting neutrino spectra, and the effects of time variability on proton injection, among other topics. We present some characteristic examples of the temporal behaviour of the system and show that this can be very different from the one exhibited by leptonic models. Furthermore, we argue that, contrary to the wide-held belief, there are parameter regimes where the hadronic models can become quite efficient. However, to keep the free parameters at a minimum and facilitate an in-depth study of the system, we have only concentrated on the case where protons are injected; i.e., we did not consider the effects of a co-accelerated leptonic component.

Assessing the Significance of Apparent Correlations between Radio and Gamma-Ray Blazar Fluxes

Whether or not 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: although it may fail to confirm an existing intrinsic correlation in a small sample that cannot be split into many redshift bins, it will 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.

Fermi-LAT upper limits on gamma-ray emission from colliding wind binaries

Context. Colliding wind binaries (CWBs) are thought to give rise to a plethora of physical processes including acceleration and interaction of relativistic particles. Observation of synchrotron radiation in the radio band confirms there is a relativistic electron population in CWBs. Accordingly, CWBs have been suspected sources of high-energy -ray emission since the COS-B era. Theoretical models exist that characterize the underlying physical processes leading to particle acceleration and quantitatively predict the non-thermal energy emission observable at Earth. Aims. We strive to find evidence of -ray emission from a sample of seven CWB systems: WR 11, WR 70, WR 125, WR 137, WR 140, WR 146, and WR 147. Theoretical modelling identified these systems as the most favourable candidates for emitting -rays. We make a comparison with existing -ray flux predictions and investigate possible constraints. Methods. We used 24 months of data from the Large Area Telescope (LAT) on-board the Fermi Gamma Ray Space Telescope to perform a dedicated likelihood analysis of CWBs in the LAT energy range. Results. We find no evidence of -ray emission from any of the studied CWB systems and determine corresponding flux upper limits. For some CWBs the interplay of orbital and stellar parameters renders the Fermi -LAT data not sensitive enough to constrain the parameter space of the emission models. In the cases of WR140 and WR147, the Fermi -LAT upper limits appear to rule out some model predictions entirely and constrain theoretical models over a significant parameter space. A comparison of our findings to the CWB Car is made.

A Hard X-Ray View of Two Distant VHE Blazars: 1ES 1101–232 and 1ES 1553+113

TeV blazars are known as prominent nonthermal emitters across the entire electromagnetic spectrum with their photon power peaking in the X-ray and TeV bands. If distant, absorption of ?-ray photons by the extragalactic background light (EBL) alters the intrinsic TeV spectral shape, thereby affecting the overall interpretation. Suzaku observations for two of the more distant TeV blazars known to date, 1ES 1101?232 and 1ES 1553+113, were carried out in 2006 May and July, respectively, including a quasi-simultaneous coverage with the state-of-the-art Cerenkov telescope facilities. We report on the resulting data sets with emphasis on the X-ray band and set in context to their historical behavior. During our campaign, we did not detect any significant X-ray or ?-ray variability. 1ES 1101?232 was found in a quiescent state with the lowest X-ray flux ever measured. The combined XIS and HXD PIN data for 1ES 1101?232 and 1ES 1553+113 clearly indicate spectral curvature up to the highest hard X-ray data point (~30 keV), manifesting as softening with increasing energy. We describe this spectral shape by either a broken power law or a log-parabolic fit with equal statistical goodness of fits. The combined 1ES 1553+113 very high energy spectrum (90-500 GeV) did not show any significant changes with respect to earlier observations. The resulting contemporaneous broadband spectral energy distributions of both TeV blazars are discussed in view of implications for intrinsic blazar parameter values, taking into account the ?-ray absorption in the EBL.

The first full orbit of η Carinae seen by Fermi

Aims. The binary system η Carinae has completed its first 5.54 y orbit since the beginning of science operation of the Fermi Large Area Telescope (LAT). We are now able to investigate the high-energy γ-ray source at the position of η Carinae over its full orbital period. By this, we can address and confirm earlier predictions for temporal and spectral variability. Methods. Newer versions of the LAT datasets, instrument response functions and background models allow for a more accurate analysis. Therefore it is important to re-evaluate the previously analyzed time period along with the new data to further constrain location, spectral shape, and flux time history of the γ-ray source. Results. We confirm earlier predictions of increasing flux values above 10 GeV toward the next periastron passage. For the most recent part of the data sample, flux values as high as those before the first periastron passage in 2008 are recorded. A comparison of spectral energy distributions around periastron and apastron passages reveals strong variation in the high-energy band. This is due to a second spectral component that is present only around periastron. Conclusions. Improved spatial consistency with the γ-ray source at the position of η Carinae along with the confirmation of temporal variability above 10 GeV in conjunction with the orbital period strengthens the argument for unambiguous source identification. Spectral variability provides additional constraints for future modeling of the particle acceleration and γ-ray emission in colliding-wind binary systems.

Implications of the VHE gamma-ray detection of the Quasar 3C279

The MAGIC collaboration recently reported the detection of the quasar 3C279 at > 100 GeV gamma-ray energies. Here we present simultaneous optical (BVRI) and X-ray (RXTE PCA) data from the day of the VHE detection and discuss the implications of the snap-shot spectral energy distribution for jet models of blazars. A one-zone synchrotron-self-Compton origin of the entire SED, including the VHE gamma-ray emission can be ruled out. The VHE emission could, in principle, be interpreted as Compton upscattering of external radiation (e.g., from the broad-line regions). However, such an interpretation would require either an unusually low magnetic field of B ~ 0.03 G or an unrealistically high Doppler factor of Gamma ~ 140. In addition, such a model fails to reproduce the observed X-ray flux. This as well as the lack of correlated variability in the optical with the VHE gamma-ray emission and the substantial gamma-gamma opacity of the BLR radiation field to VHE gamma-rays suggests a multi-zone model. In particular, an SSC model with an emission region far outside the BLR reproduces the simultaneous X-ray -- VHE gamma-ray spectrum of 3C279. Alternatively, a hadronic model is capable of reproducing the observed SED of 3C279 reasonably well. However, the hadronic model requires a rather extreme jet power of L_j ~ 10^{49} erg s^{-1}, compared to a requirement of L_j ~ 2 X 10^{47} erg s^{-1} for a multi-zone leptonic model.

On the Physics of Hadronic Blazar Emission Models

The material of relativistic jets of active galactic nuclei is composed of leptons and hadrons. If a significant fraction of the overall jet power is converted into energizing protons to relativistic energies, along with the leptons, hadronic interactions play a role in converting particle energy into radiation and neutrino production. I will discuss the physics and some properties of hadronic blazar emission models in the framework of recent developments in this field.