As the largest mass concentrations in the local Universe, nearby clusters of galaxies and their central galaxies are prime targets in searching for indirect signatures of dark matter annihilation (DMA). We seek to constrain the dark matter annihilation emission component from multi-frequency observations of the central galaxy of the Virgo cluster. The annihilation emission component is modeled by the prompt and inverse-Compton gamma rays from the hadronization of annihilation products from generic weakly interacting dark matter particles. This component is fitted to the excess of the observed data above the spectral energy distribution (SED) of the jet in M87, described with a best-fit synchrotron-self-Compton (SSC) spectrum. While this result is not sufficiently significant to claim a detection, we emphasize that a dark matter "double hump signature" can be used to unambiguously discriminate the dark matter emission component from the variable jet-related emission of M87 in future, more extended observation campaigns.

Clusters of galaxies, such as the Virgo cluster, host enormous quantities of dark matter, making them prime targets for efforts in indirect dark matter detection via potential radiative signatures from annihilation of dark matter particles and subsequent radiative losses of annihilation products. However, a careful study of ubiquitous astrophysical backgrounds is mandatory to single out potential evidence for dark matter annihilation. Here, we construct a multiwavelength spectral energy distribution for the central radio galaxy in the Virgo cluster, M87, using a state-of-the-art numerical Synchrotron Self Compton approach. Fitting recent Chandra, Fermi-LAT and Cherenkov observations, we probe different dark matter annihilation scenarios including a full treatment of the inverse Compton losses from electrons and positrons produced in the annihilation. It is shown that such a template can substantially improve upon existing dark matter detection limits.

Context. 3C 279, the first quasar discovered to emit VHE gamma-rays by the MAGIC telescope in 2006, was reobserved by MAGIC in January 2007 during a major optical flare and from December 2008 to April 2009 following an alert from the Fermi space telescope on an exceptionally high gamma -ray state. Aims. The January 2007 observations resulted in a detection on January 16 with significance 5.2 sigma, corresponding to a F(> 150 GeV) (3.8 \pm 0.8) \cdot 10^-11 ph cm^-2 s^-1 while the overall data sample does not show significant signal. The December 2008 - April 2009 observations did not detect the source. We study the multiwavelength behavior of the source at the epochs of MAGIC observations, collecting quasi-simultaneous data at optical and X-ray frequencies and for 2009 also gamma-ray data from Fermi. Methods. We study the light curves and spectral energy distribution of the source. The spectral energy distributions of three observing epochs (including the February 2006, which has been previously published in Albert et al. 2008a) are modeled with one-zone inverse Compton models and the emission on January 16, 2007 also with two zone model and with a lepto-hadronic model. Results. We find that the VHE gamma-ray emission detected in 2006 and 2007 challenges standard one-zone model, based on relativistic electrons in a jet scattering broad line region photons, while the other studied models fit the observed spectral energy distribution more satisfactorily.

We discuss simplified models for photo-meson production in cosmic accelerators, such as Active Galactic Nuclei and Gamma-Ray Bursts. Our self-consistent models are directly based on the underlying physics used in the SOPHIA software, and can be easily adapted if new data are included. They allow for the efficient computation of neutrino and photon spectra (from pi^0 decays), as a major requirement of modern time-dependent simulations of the astrophysical sources and parameter studies. In addition, the secondaries (pions and muons) are explicitely generated, a necessity if cooling processes are to be included. For the neutrino production, we include the helicity dependence of the muon decays which in fact leads to larger corrections than the details of the interaction model. The separate computation of the pi^0, pi^+, and pi^- fluxes allows, for instance, for flavor ratio predictions of the neutrinos at the source, which are a requirement of many tests of neutrino properties using astrophysical sources. We confirm that for charged pion generation, the often used production by the Delta(1232)-resonance is typically not the dominant process in Active Galactic Nuclei and Gamma-Ray Bursts, and we show, for arbitrary input spectra, that the number of neutrinos are underestimated by at least a factor of two if they are obtained from the neutral to charged pion ratio. We compare our results for several levels of simplification using isotropic synchrotron and thermal spectra, and we demonstrate that they are sufficiently close to the SOPHIA software.

In this paper we present a fully selfconsistent SSC model with particle acceleration due to shock and stochastic acceleration (Fermi-I and Fermi-II-Processes respectively) to model the quiescent spectral energy distribution (SED) observed from PKS 2155. The simultaneous August/September 2008 multiwavelength data of H.E.S.S., Fermi, RXTE, SWIFT and ATOM give new constraints to the high-energy peak in the SED concerning its curvature. We find that, in our model, a monoenergetic injection of electrons at $\gamma_0 = 910$ into the model region, which are accelerated by Fermi-I- and Fermi-II-processes while suffering synchrotron and inverse Compton losses, finally leads to the observed SED of PKS 2155-30.4 shown in H.E.S.S. and Fermi-LAT collaborations (2009). In contrast to other SSC models our parameters arise from the jet's microphysics and the spectrum is evolving selfconsistently from diffusion and acceleration. The $\gamma_0$-factor can be interpreted as two counterstreaming plasmas due to the motion of the blob at a bulk factor of $\Gamma = 58$ and opposed moving upstream electrons at moderate Lorentz factors with an average of $\gamma_u \approx 8$.

We employ a time-dependent synchrotron-self-Compton code for modeling contemporaneous multiwavelength data of the blazar 1 ES 1218+30.4 The input parameters of the model are used to infer physical parameters of the emitting region. An acceptable fit to the data is obtained by taking into account a stellar emission component in the optical regime due to the host galaxy. The physical parameters inferred from the fit are in line with particle acceleration due to the Fermi mechanism providing s = 2.1 spectra. From the properties of the host galaxy in the optical, we estimate the central black hole mass and thus confirm that the jet power injected into the emission region is in the sub-Eddington regime, as expected for BL Lacertae objects.

We report the results of a multi-wavelength campaign on the blazar Mrk 421 during outburst. We observed four strong flares at X-ray energies that were not seen at other wavelengths (partially because of missing data). From the fastest rise in the X-rays, an upper limit could be derived on the extension of the emission region. A time lag between high-energy and low-energy X-rays was observed, which allowed an estimation of the magnetic-field strength. The spectral analysis of the X-rays revealed a slight spectral hardening of the low-energy (3 - 43 keV) spectral index. The hardness-ratio analysis of the Swift-XRT (0.2 - 10 keV) data indicated a small correlation with the intensity; i. e., a hard-to-soft evolution was observed. At the energies of IBIS/ISGRI (20 - 150 keV), such correlations are less obvious. A multiwavelength spectrum was composed and the X-ray and bolometric luminosities are calculated.