The underlying physics of giant and mini radio halos in galaxy clusters is still an open question. We find that mini halos (such as in Perseus and Ophiuchus) can be explained by radio-emitting electrons that are generated in hadronic cosmic ray (CR)
interactions with protons of the intracluster medium. By contrast, the hadronic model either fails to explain the extended emission of giant radio halos (as in Coma at low frequencies) or would require a flat CR profile, which can be realized through outward streaming and diffusion of CRs (in Coma and A2163 at 1.4 GHz). We suggest that a second, leptonic component could be responsible for the missing flux in the outer parts of giant halos within a new hybrid scenario and we describe its possible observational consequences. To study the hadronic emission component of the radio halo population statistically, we use a cosmological mock galaxy cluster catalog built from the MultiDark simulation. Because of the properties of CR streaming and the different scalings of the X-ray luminosity (L_X) and the Sunyaev-Zeldovich flux (Y) with gas density, our model can simultaneously reproduce the observed bimodality of radio-loud and radio-quiet clusters at the same L_X as well as the unimodal distribution of radio-halo luminosity versus Y; thereby suggesting a physical solution to this apparent contradiction. We predict radio halo emission down to the mass scale of galaxy groups, which highlights the unique prospects for low-frequency radio surveys (such as the LOFAR Tier 1 survey) to increase the number of detected radio halos by at least an order of magnitude.
Cosmological hydrodynamical simulations of galaxy clusters are still challenged to produce a model for the intracluster medium that matches all aspects of current X-ray and Sunyaev-Zeldovich observations. To facilitate such comparisons with future si
mulations and to enable realistic cluster population studies for modeling e.g., non-thermal emission processes, we construct a phenomenological model for the intracluster medium that is based on a representative sample of observed X-ray clusters. We create a mock galaxy cluster catalog based on the large collisionless N-body simulation MultiDark, by assigning our gas density model to each dark matter cluster halo. Our clusters are classified as cool-core and non cool-core according to a dynamical disturbance parameter. We demonstrate that our gas model matches the various observed Sunyaev-Zeldovich and X-ray scaling relations as well as the X-ray luminosity function, thus enabling to build a reliable mock catalog for present surveys and forecasts for future experiments. In a companion paper, we apply our catalogs to calculate non-thermal radio and gamma-ray emission of galaxy clusters. We make our cosmologically complete multi-frequency mock catalogs for the (non-)thermal cluster emission at different redshifts publicly and freely available online through the MultiDark database (www.multidark.org).
Fermi has been instrumental in constraining the luminosity function and redshift evolution of gamma-ray bright blazars. This includes limits upon the spectrum and anisotropy of the extragalactic gamma-ray background (EGRB), redshift distribution of n
earby Fermi active galactic nuclei (AGN), and the construction of a log(N)-log(S) relation. Based upon these, it has been argued that the evolution of the gamma-ray bright blazar population must be much less dramatic than that of other AGN. However, critical to such claims is the assumption that inverse Compton cascades reprocess emission above a TeV into the Fermi energy range, substantially enhancing the strength of the observed limits. Here we demonstrate that in the absence of such a process, due, e.g., to the presence of virulent plasma beam instabilities that preempt the cascade, a population of TeV-bright blazars that evolve similarly to quasars is consistent with the population of hard gamma-ray blazars observed by Fermi. Specifically, we show that a simple model for the properties and luminosity function is simultaneously able to reproduce their log(N)-log(S) relation, local redshift distribution, and contribution to the EGRB and its anisotropy without any free parameters. Insofar the naturalness of a picture in which the hard gamma-ray blazar population exhibits the strong redshift evolution observed in other tracers of the cosmological history of accretion onto halos is desirable, this lends support for the absence of the inverse Compton cascades and the existence of the beam plasma instabilities.
In principle, the angular anisotropy in the extragalactic gamma-ray background (EGRB) places severe constraints upon putative populations of unresolved gamma-ray point sources. Existing estimates of the EGRB anisotropy have been constructed by excisi
ng known point sources, e.g., taken from the First or 2 Year Fermi-LAT Source Catalog (1FGL or 2FGL, respectively) and statistically analyzing the residual gamma-ray sky maps. We perform an independent check of the EGRB anisotropy limits by comparing the values obtained from the 1FGL-masked sky maps to the signal implied by sources that lie below the 1FGL detection threshold in the more sensitive 2FGL and 1FHL (First Fermi-LAT catalog of >10 GeV sources). As such, our analysis provides an internal consistency check of implications for source counts and spectral index distributions of gamma-ray bright active galactic nuclei obtained from Fermi-LAT data. Based on this, we find evidence for substantially larger anisotropies than those previously reported at energies above 5 GeV, where BL Lac objects are likely to provide the bulk of their contribution to the EGRB. This uncertainty in the EGRB anisotropy cautions against using it as an independent constraint for the high-redshift gamma-ray universe. Moreover, this would suggest that contrary to previous claims, smooth extensions of the resolved point-source population may be able to simultaneously explain both the isotropic and anisotropic components of the EGRB.
A subset of blazars emit TeV gamma rays which annihilate and pair produce on the extragalactic background light. We have argued in Broderick et al. (2011, Paper I) that plasma beam instabilities can dissipate the pairs energy locally. This heats the
intergalactic medium and dramatically increases its entropy after redshift z~2, with important implications for structure formation: (1) This suggests a scenario for the origin of the cool core (CC)/non-cool core (NCC) bimodality in galaxy clusters and groups. Early forming galaxy groups are unaffected because they can efficiently radiate the additional entropy, developing a CC. However, late forming groups do not have sufficient time to cool before the entropy is gravitationally reprocessed through successive mergers - counteracting cooling and raising the core entropy further. Hence blazar heating works different than feedback by active galactic nuclei, which balances radiative cooling but is unable to transform CC into NCC clusters due to the weak coupling to the cluster gas. (2) We predict a suppression of the Sunyaev-Zeldovich power spectrum on angular scales smaller than 5 due to the globally reduced central pressure of groups and clusters forming after z~1. (3) Our redshift dependent entropy floor increases the characteristic halo mass below which dwarf galaxies cannot form by a factor of ~10 (50) at mean density (in voids) over that found in models that include photoionization alone. This prevents the formation of late forming dwarf galaxies (z<2) with masses ranging from 10^{10} to 10^{11} M_sun for redshifts z~2 to 0, respectively. This may help resolve the missing satellite problem in the Milky Way and the void phenomenon of the low observed abundances of dwarf satellites compared to cold dark matter simulations and may bring the observed early star formation histories into agreement with galaxy formation models. (abridged)
Clusters of galaxies, filled with hot magnetized plasma, are the largest bound objects in existence and an important touchstone in understanding the formation of structures in our Universe. In such clusters, thermal conduction follows field lines, so
magnetic fields strongly shape the clusters thermal history; that some have not since cooled and collapsed is a mystery. In a seemingly unrelated puzzle, recent observations of Virgo cluster spiral galaxies imply ridges of strong, coherent magnetic fields offset from their centre. Here we demonstrate, using three-dimensional magnetohydrodynamical simulations, that such ridges are easily explained by galaxies sweeping up field lines as they orbit inside the cluster. This magnetic drape is then lit up with cosmic rays from the galaxies stars, generating coherent polarized emission at the galaxies leading edges. This immediately presents a technique for probing local orientations and characteristic length scales of cluster magnetic fields. The first application of this technique, mapping the field of the Virgo cluster, gives a startling result: outside a central region, the magnetic field is preferentially oriented radially as predicted by the magnetothermal instability. Our results strongly suggest a mechanism for maintaining some clusters in a non-cooling-core state.
Recently, it has been shown that electrons and positrons from dark matter (DM) annihilations provide an excellent fit to the Fermi, PAMELA, and HESS data. Using this DM model, which requires an enhancement of the annihilation cross section over its s
tandard value to match the observations, we show that it immediately implies an observable level of gamma-ray emission for the Fermi telescope from nearby galaxy clusters such as Virgo and Fornax. We show that this DM model implies a peculiar feature from final state radiation that is a distinctive signature of DM. Using the EGRET upper limit on the gamma-ray emission from Virgo, we constrain the minimum mass of substructures within DM halos to be > 5x10^-3 M_sun -- about four orders of magnitudes larger than the expectation for cold dark matter. This limits the cutoff scale in the linear matter power spectrum to k < 35/kpc which can be explained by e.g., warm dark matter. Very near future Fermi observations will strongly constrain the minimum mass to be > 10^3 M_sun: if the true substructure cutoff is much smaller than this, the DM interpretation of the Fermi/PAMELA/HESS data must be wrong. To address the problem of astrophysical foregrounds, we performed high-resolution, cosmological simulations of galaxy clusters that include realistic cosmic ray (CR) physics. We compute the dominating gamma-ray emission signal resulting from hadronic CR interactions and find that it follows a universal spectrum and spatial distribution. If we neglect the anomalous enhancement factor and assume standard values for the cross section and minimum subhalo mass, the same model of DM predicts comparable levels of the gamma-ray emission from DM annihilations and CR interactions. This suggests that spectral subtraction techniques could be applied to detect the annihilation signal.
