Showing 101–150 of 168 results for author: Pfrommer, C

The Effect of Nonlinear Landau Damping on Ultrarelativistic Beam Plasma Instabilities

Authors: Philip Chang, Avery E. Broderick, Christoph Pfrommer, Ewald Puchwein, Astrid Lamberts, Mohamad Shalaby

Abstract: Very-high energy gamma-rays from extragalactic sources pair-produce off of the extragalactic background light, yielding an electron-positron pair beam. This pair beam is unstable to various plasma instabilities, especially the "oblique" instability, which can be the dominant cooling mechanism for the beam. However, recently, it has been claimed that nonlinear Landau damping renders it physically i… ▽ More Very-high energy gamma-rays from extragalactic sources pair-produce off of the extragalactic background light, yielding an electron-positron pair beam. This pair beam is unstable to various plasma instabilities, especially the "oblique" instability, which can be the dominant cooling mechanism for the beam. However, recently, it has been claimed that nonlinear Landau damping renders it physically irrelevant by reducing the effective damping rate to a low level. Here, we show with numerical calculations that the effective damping rate is $8\times 10^{-4}$ of the growth rate of the linear instability, which is sufficient for the "oblique" instability to be the dominant cooling mechanism of these pair beams. In particular, we show that previous estimates of this rate ignored the exponential cutoff in the scattering amplitude at large wavenumber and assumed that the damping of scattered waves entirely depends on collisions, ignoring collisionless processes. We find that the total wave energy eventually grows to approximate equipartition with the beam by increasingly depositing energy into long wavelength modes. As we have not included the effect of nonlinear wave-wave interactions on these long wavelength modes, this scenario represents the "worst-case" scenario for the oblique instability. As it continues to drain energy from the beam at a faster rate than other processes, we conclude that the "oblique" instability is sufficiently strong to make it the physically dominant cooling mechanism for high-energy pair beams in the intergalactic medium. △ Less

Submitted 14 October, 2014; v1 submitted 7 October, 2014; originally announced October 2014.

Comments: 6 pages, 3 figures, accepted to ApJ

On the Cluster Physics of Sunyaev-Zel'dovich and X-ray Surveys IV: Characterizing Density and Pressure Clumping due to Infalling Substructures

Authors: N. Battaglia, J. R. Bond, C. Pfrommer, J. L. Sievers

Abstract: Understanding the outskirts of galaxy clusters at the virial radius (R200) and beyond is critical for an accurate determination of cluster masses and to ensure unbiased cosmological parameter estimates from cluster surveys. This problem has drawn renewed interest due to recent determinations of gas mass fractions beyond R200, which appear to be considerably larger than the cosmic mean, and because… ▽ More Understanding the outskirts of galaxy clusters at the virial radius (R200) and beyond is critical for an accurate determination of cluster masses and to ensure unbiased cosmological parameter estimates from cluster surveys. This problem has drawn renewed interest due to recent determinations of gas mass fractions beyond R200, which appear to be considerably larger than the cosmic mean, and because the clusters' total Sunyaev-Zel'dovich flux receives a significant contribution from these regions. Here, we use a large suite of cosmological hydrodynamical simulations to study the clumpiness of density and pressure and employ different variants of simulated physics, including radiative gas physics and thermal feedback by active galactic nuclei. We find that density and pressure clumping closely trace each other as a function of radius, but the bias on density remains on average < 20% within the virial radius R200. At larger radius, clumping increases steeply due to the continuous infall of coherent structures that have not yet passed the accretion shock. Density and pressure clumping increase with cluster mass and redshift, which probes on average dynamically younger objects that are still in the process of assembling. The angular power spectra of gas density and pressure show that the clumping signal is dominated by comparably large substructures with scales >R200/5, signaling the presence of gravitationally-driven "super-clumping". In contrast, the angular power spectrum of the dark matter (DM) shows an almost uniform size distribution due to unimpeded subhalos. The quadrupolar anisotropy dominates the signal and correlates well across different radii as a result of the prolateness of the DM potential. We provide a synopsis of the radial dependence of the clusters' non-equilibrium measures (kinetic pressure support, ellipticity, and clumping) that all increase sharply beyond R200. △ Less

Submitted 13 May, 2014; originally announced May 2014.

Comments: 17 pages, 16 figures, submitted to ApJ, comments welcome

On the Evolution of the Cosmic Ionizing Background

Authors: Fabio Fontanot, Stefano Cristiani, Christoph Pfrommer, Guido Cupani, Eros Vanzella

Abstract: We study the observed cosmic ionizing background as a constraint on the nature of the sources responsible for the reionization of the Universe. In earlier work, we showed that extrapolations of the Ultra-Violet Luminosity Function (LF) of Lyman Break Galaxies (LBGs) at fixed Lyman continuum photon escape fraction are not able to reproduce the redshift evolution of this background. Here, we employ… ▽ More We study the observed cosmic ionizing background as a constraint on the nature of the sources responsible for the reionization of the Universe. In earlier work, we showed that extrapolations of the Ultra-Violet Luminosity Function (LF) of Lyman Break Galaxies (LBGs) at fixed Lyman continuum photon escape fraction are not able to reproduce the redshift evolution of this background. Here, we employ extrapolations of the high-z LFs to describe the contribution of LBGs to the ionizing photon rate, taking into account the smoothing of the baryonic perturbations, due to the background itself (i.e. the filtering mass), as well as a possible sharp increase of the escape fraction in dwarf galaxies. Under the hypothesis of a dominant contribution of LBGs to cosmic reionization, our results suggest that sources fainter than the current observational limits should be characterised by escape fractions of the order of ~0.1-0.3 (larger than the current estimates for bright galaxies) to account for a z>6 reionization and the measured evolution of cosmic ionizing background, at the same time. The contribution to the background from quasars turns out to be relevant at z<3. Overall, our results support the case for dedicated observations of faint galaxies in the rest-frame UV, in order to better determine their physical properties. Observed escape fractions outside our proposed range bear relevant consequences on the nature of the astrophysical sources responsible for cosmic reionization and/or its buildup process. △ Less

Submitted 2 December, 2013; originally announced December 2013.

Comments: 9 pages, 6 figures, MNRAS accepted

On the Physics of Radio Halos in Galaxy Clusters: Scaling Relations and Luminosity Functions

Authors: Fabio Zandanel, Christoph Pfrommer, Francisco Prada

Abstract: 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 hal… ▽ More 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-Zel'dovich 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. △ Less

Submitted 9 December, 2015; v1 submitted 19 November, 2013; originally announced November 2013.

Comments: Accepted for publication in MNRAS. 22 pages, 9 figures. This article supersedes arXiv:1207.6410. Updated to match the published version. Latest version updated to match the Erratum published after the correction of the mock catalogs

Report number: MNRAS Manuscript stt2250

Journal ref: MNRAS 438, 124-144, 2014

A Phenomenological Model for the Intracluster Medium that matches X-ray and Sunyaev-Zel'dovich observations

Authors: Fabio Zandanel, Christoph Pfrommer, Francisco Prada

Abstract: 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-Zel'dovich observations. To facilitate such comparisons with future simulations and to enable realistic cluster population studies for modeling e.g., non-thermal emission processes, we construct a phenomenological model… ▽ More 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-Zel'dovich observations. To facilitate such comparisons with future simulations 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-Zel'dovich 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). △ Less

Submitted 7 December, 2015; v1 submitted 19 November, 2013; originally announced November 2013.

Comments: Accepted for publication in MNRAS. 8 pages, 4 figures. This article draws heavily from arXiv:1207.6410. Updated to match the published version. Latest version updated to match the Erratum published after the correction of the mock catalogs

Journal ref: MNRAS 438, 116-123, 2014

Measuring the Thermal Sunyaev-Zel'dovich Effect Through the Cross Correlation of Planck and WMAP Maps with ROSAT Galaxy Cluster Catalogs

Authors: Amir Hajian, Nicholas Battaglia, David N. Spergel, J. Richard Bond, Christoph Pfrommer, Jonathan L. Sievers

Abstract: We measure a significant correlation between the thermal Sunyaev-Zel'dovich effect in the Planck and WMAP maps and an X-ray cluster map based on ROSAT. We use the 100, 143 and 343 GHz Planck maps and the WMAP 94 GHz map to obtain this cluster cross spectrum. We check our measurements for contamination from dusty galaxies using the cross correlations with the 220, 545 and 843 GHz maps from Planck.… ▽ More We measure a significant correlation between the thermal Sunyaev-Zel'dovich effect in the Planck and WMAP maps and an X-ray cluster map based on ROSAT. We use the 100, 143 and 343 GHz Planck maps and the WMAP 94 GHz map to obtain this cluster cross spectrum. We check our measurements for contamination from dusty galaxies using the cross correlations with the 220, 545 and 843 GHz maps from Planck. Our measurement yields a direct characterization of the cluster power spectrum over a wide range of angular scales that is consistent with large cosmological simulations. The amplitude of this signal depends on cosmological parameters that determine the growth of structure (σ_8 and Ω_M) and scales as σ_8^7.4 and Ω_M^1.9 around the multipole (ell) ~ 1000. We constrain σ_8 and Ω_M from the cross-power spectrum to be σ_8 (Ω_M/0.30)^0.26 = 0.8 +/- 0.02. Since this cross spectrum produces a tight constraint in the σ_8 and Ω_M plane the errors on a σ_8 constraint will be mostly limited by the uncertainties from external constraints. Future cluster catalogs, like those from eRosita and LSST, and pointed multi-wavelength observations of clusters will improve the constraining power of this cross spectrum measurement. In principle this analysis can be extended beyond σ_8 and Ω_M to constrain dark energy or the sum of the neutrino masses. △ Less

Submitted 12 September, 2013; originally announced September 2013.

Comments: 30 pages, 8 figures

Journal ref: Journal of Cosmology and Astroparticle Physics, Issue 11, article id. 064, pp. (2013)

Introduction to extragalactic sources of very high-energy photons

Authors: C. Pfrommer

Abstract: The launch of the Fermi gamma-ray space telescope and the imaging air Cerenkov telescopes H.E.S.S., MAGIC, and VERITAS have substantially transformed our knowledge of gamma-ray sources in the last decade. The extragalactic gamma-ray sky is teeming with blazars, which are active galactic nuclei whose jet is directed at us. Additionally, there are radio galaxies, starburst and spiral galaxies, and g… ▽ More The launch of the Fermi gamma-ray space telescope and the imaging air Cerenkov telescopes H.E.S.S., MAGIC, and VERITAS have substantially transformed our knowledge of gamma-ray sources in the last decade. The extragalactic gamma-ray sky is teeming with blazars, which are active galactic nuclei whose jet is directed at us. Additionally, there are radio galaxies, starburst and spiral galaxies, and gamma-ray bursts, albeit with smaller numbers. Galaxy clusters have not yet been observed in gamma rays. Here, I will introduce the different gamma-ray emission processes and review what they may tell us about these objects and the underlying acceleration mechanisms. Beyond the study of these fascinating objects, TeV gamma rays from blazars probe the integrated star formation history of the universe. Studies of TeV blazar spectra may provide us with insights into intergalactic magnetic fields or alternatively, may lead us to infer the existence of a novel mechanism that heats the intergalactic medium at late times (for redshifts z<3) and impacts the Lyman-alpha forest and late-time structure formation. The TeV gamma-ray emission may also allow us to probe fundamental physics such as the structure of space time. △ Less

Submitted 29 August, 2013; originally announced August 2013.

Comments: 10 pages, 5 multi-panel figures. Introduction also aimed at non-experts. Proceedings of the conference "Rencontres de Moriond 2013: Very High Energy Phenomena in the Universe" (La Thuile, March 9th - 16th, 2013)

The physics and cosmology of TeV blazars in a nutshell

Authors: C. Pfrommer, A. E. Broderick, P. Chang, E. Puchwein, V. Springel

Abstract: The extragalactic gamma-ray sky at TeV energies is dominated by blazars, a subclass of accreting super-massive black holes with powerful relativistic outflows directed at us. Only constituting a small fraction of the total power output of black holes, blazars were thought to have a minor impact on the universe at best. As we argue here, the opposite is true and the gamma-ray emission from TeV blaz… ▽ More The extragalactic gamma-ray sky at TeV energies is dominated by blazars, a subclass of accreting super-massive black holes with powerful relativistic outflows directed at us. Only constituting a small fraction of the total power output of black holes, blazars were thought to have a minor impact on the universe at best. As we argue here, the opposite is true and the gamma-ray emission from TeV blazars can be thermalized via beam-plasma instabilities on cosmological scales with order unity efficiency, resulting in a potentially dramatic heating of the low-density intergalactic medium. Here, we review this novel heating mechanism and explore the consequences for the formation of structure in the universe. In particular, we show how it produces an inverted temperature-density relation of the intergalactic medium that is in agreement with observations of the Lyman-alpha forest. This suggests that blazar heating can potentially explain the paucity of dwarf galaxies in galactic halos and voids, and the bimodality of galaxy clusters. This also transforms our understanding of the evolution of blazars, their contribution to the extra-galactic gamma-ray background, and how their individual spectra can be used in constraining intergalactic magnetic fields. △ Less

Submitted 28 August, 2013; originally announced August 2013.

Comments: 8 pages, 6 multi-panel figures. Proceedings of the conference "Rencontres de Moriond 2013: Very High Energy Phenomena in the Universe" (La Thuile, March 9th - 16th, 2013)

Search for cosmic-ray induced gamma-ray emission in Galaxy Clusters

Authors: M. Ackermann, M. Ajello, A. Albert, A. Allafort, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, E. D. Bloom, E. Bonamente, E. Bottacini, T. J. Brandt, J. Bregeon, M. Brigida, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, R. C. G. Chaves , et al. (101 additional authors not shown)

Abstract: Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $γ$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching fo… ▽ More Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $γ$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $γ$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7σ$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $γ$-ray flux from individual clusters in our sample. △ Less

Submitted 24 March, 2014; v1 submitted 26 August, 2013; originally announced August 2013.

Comments: Manuscript version accepted for publication in ApJ with expanded set of CR models now including thermal ICM and flat CR profiles; 27 pages, 19 figures, 9 tables; Corresponding authors: S. Zimmer, zimmer@fysik.su.se; J. Conrad, conrad@fysik.su.se; C. Pfrommer, christoph.pfrommer@h-its.org; A. Pinzke, apinzke@fysik.su.se; O. Reimer, olr@slac.stanford.edu

Journal ref: M. Ackermann et al. 2014 ApJ 787 18

MAGIC Gamma-ray Observations of the Perseus Galaxy Cluster

Authors: Fabio Zandanel, Pierre Colin, Saverio Lombardi, Michele Doro, Dorit Eisenacher, Dorothee Hildebrand, Francisco Prada, Christoph Pfrommer, Anders Pinzke

Abstract: In order to detect the gamma-ray emission from cosmic ray (CR) interactions with the intra-cluster medium, the ground-based imaging Cherenkov telescope MAGIC conducted the deepest-to-date observational campaign targeting a galaxy cluster at very high-energies (> 100 GeV) and observed the Perseus cluster for a total of 85 hr during 2009-2011. The observations constrain the average CR-to-thermal pre… ▽ More In order to detect the gamma-ray emission from cosmic ray (CR) interactions with the intra-cluster medium, the ground-based imaging Cherenkov telescope MAGIC conducted the deepest-to-date observational campaign targeting a galaxy cluster at very high-energies (> 100 GeV) and observed the Perseus cluster for a total of 85 hr during 2009-2011. The observations constrain the average CR-to-thermal pressure ratio to be 1-2% and the maximum CR acceleration efficiency at structure formation shocks to be <50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Additionally, assuming that the Perseus radio mini-halo is generated by secondaries created in hadronic CR interactions, the central magnetic field is limited to be > 4-9 muG. This range is well below the field strength inferred from Faraday rotation measurements and, therefore, the hadronic model remains a plausible explanation of the Perseus radio mini-halo. Following this successful campaign, MAGIC is continuing collecting data on Perseus. △ Less

Submitted 9 August, 2013; v1 submitted 2 August, 2013; originally announced August 2013.

Comments: 4 pages, 2 figures, Proc. of 33rd ICRC, Rio de Janeiro, Brazil

Implications of Plasma Beam Instabilities for the Statistics of the Fermi Hard Gamma-ray Blazars and the Origin of the Extragalactic Gamma-Ray Background

Authors: Avery E. Broderick, Christoph Pfrommer, Ewald Puchwein, Philip Chang

Abstract: 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 nearby 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… ▽ More 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 nearby 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. △ Less

Submitted 1 August, 2013; originally announced August 2013.

Comments: 14 pages, 13 figures, submitted to ApJ

Lower Limits on the Anisotropy of the Extragalactic Gamma-Ray Background implied by the 2FGL and 1FHL Catalogs

Authors: Avery E. Broderick, Christoph Pfrommer, Ewald Puchwein, Philip Chang, Kendrick M. Smith

Abstract: 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 excising known point sources, e.g., taken from the First or 2 Year Fermi-LAT Source Catalog (1FGL or 2FGL, respectively) and statistically analyzing the res… ▽ More 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 excising 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. △ Less

Submitted 28 February, 2018; v1 submitted 31 July, 2013; originally announced August 2013.

Comments: 18 pages, 14 figures, published version

Toward a comprehensive model for feedback by active galactic nuclei: new insights from M87 observations by LOFAR, Fermi and H.E.S.S

Authors: Christoph Pfrommer

Abstract: Feedback by active galactic nuclei (AGNs) appears to be critical in balancing radiative cooling of the low-entropy gas at the centers of galaxy clusters and in mitigating the star formation of elliptical galaxies. New observations of M87 enable us to put forward a comprehensive model for the physical heating mechanism. Low-frequency radio observations by LOFAR revealed the absence of fossil cosmic… ▽ More Feedback by active galactic nuclei (AGNs) appears to be critical in balancing radiative cooling of the low-entropy gas at the centers of galaxy clusters and in mitigating the star formation of elliptical galaxies. New observations of M87 enable us to put forward a comprehensive model for the physical heating mechanism. Low-frequency radio observations by LOFAR revealed the absence of fossil cosmic-ray (CR) electrons in the radio halo surrounding M87. This puzzle can be resolved by accounting for the CR release from the radio cocoons and the subsequent mixing of CRs with the dense ambient intracluster gas, which thermalizes the electrons on a timescale similar to the radio halo age of 40 Myr. Hadronic interactions of similarly injected CR protons with the ambient gas should produce an observable gamma-ray signal in accordance with the steady emission of the low state of M87 detected by Fermi and H.E.S.S. Hence, we normalize the CR population to the gamma-ray emission, which shows the same spectral slope as the CR injection spectrum probed by LOFAR, thereby supporting a common origin. We show that CRs, which stream at the Alfven velocity with respect to the plasma rest frame, heat the surrounding thermal plasma at a rate that balances that of radiative cooling on average at each radius. However, the resulting global thermal equilibrium is locally unstable and allows for the formation of the observed cooling multi-phase medium through thermal instability. Provided that CR heating balances cooling during the emerging "cooling flow," the collapse of the majority of the gas is halted around 1 keV - in accordance with X-ray data. We show that both the existence of a temperature floor and the similar radial scaling of the heating and cooling rates are generic predictions of the CR heating model. △ Less

Submitted 24 October, 2013; v1 submitted 21 March, 2013; originally announced March 2013.

Comments: 18 pages, 9 figures, ApJ in print, final version that matches the one to be published, improved predictions for the radio emission, discussion added about the impact of a variable Alfven speed on the global and local thermal stability, all conclusions unchanged

Giant radio relics in galaxy clusters: reacceleration of fossil relativistic electrons?

Authors: Anders Pinzke, S. Peng Oh, Christoph Pfrommer

Abstract: Many bright radio relics in the outskirts of galaxy clusters have low inferred Mach numbers, defying expectations from shock acceleration theory and heliospheric observations that the injection efficiency of relativistic particles plummets at low Mach numbers. With a suite of cosmological simulations, we follow the diffusive shock acceleration as well as radiative and Coulomb cooling of cosmic ray… ▽ More Many bright radio relics in the outskirts of galaxy clusters have low inferred Mach numbers, defying expectations from shock acceleration theory and heliospheric observations that the injection efficiency of relativistic particles plummets at low Mach numbers. With a suite of cosmological simulations, we follow the diffusive shock acceleration as well as radiative and Coulomb cooling of cosmic ray electrons during the assembly of a cluster. We find a substantial population of fossil electrons. When reaccelerated at a shock (through diffusive shock acceleration), they are competitive with direct injection at strong shocks and overwhelmingly dominate by many orders of magnitude at weak shocks, Mach < 3, which are the vast majority at the cluster periphery. Their relative importance depends on cooling physics and is robust to the shock acceleration model used. While the abundance of fossils can vary by a factor of ~10, the typical reaccelerated fossil population has radio brightness in excellent agreement with observations. Fossil electrons with 1 < gamma < 100 (10 < gamma < 10^4) provide the main seeds for reacceleration at strong (weak) shocks; we show that these are well-resolved by our simulation. We construct a simple self-similar analytic model which assumes steady recent injection and cooling. It agrees well with our simulations, allowing rapid estimates and physical insight into the shape of the distribution function. We predict that LOFAR should find many more bright steep-spectrum radio relics, which are inconsistent with direct injection. A failure to take fossil cosmic ray electrons into account will lead to erroneous conclusions about the nature of particle acceleration at weak shocks; they arise from well-understood physical processes and cannot be ignored. △ Less

Submitted 2 September, 2013; v1 submitted 23 January, 2013; originally announced January 2013.

Comments: Published in MNRAS, 25 pages, 20 figures. Cluster region of interest changed to (0.8-1.0) R_200, new figure added, figures and references updated. Conclusions unchanged

On the Cluster Physics of Sunyaev-Zel'dovich and X-ray Surveys III: Measurement Biases and Cosmological Evolution of Gas and Stellar Mass Fractions

Authors: N. Battaglia, J. R. Bond, C. Pfrommer, J. L. Sievers

Abstract: Gas masses tightly correlate with the virial masses of galaxy clusters, allowing for a precise determination of cosmological parameters by means of large-scale X-ray surveys. However, according to recent Suzaku X-ray measurements, gas mass fractions, f_gas, appear to be considerably larger than the cosmic mean at the virial radius, R_200, questioning the accuracy of the cosmological parameter esti… ▽ More Gas masses tightly correlate with the virial masses of galaxy clusters, allowing for a precise determination of cosmological parameters by means of large-scale X-ray surveys. However, according to recent Suzaku X-ray measurements, gas mass fractions, f_gas, appear to be considerably larger than the cosmic mean at the virial radius, R_200, questioning the accuracy of the cosmological parameter estimations. Here, we use a large suite of cosmological hydrodynamical simulations to study measurement biases of f_gas. We employ different variants of simulated physics, including radiative gas physics, star formation, and thermal feedback by active galactic nuclei. Computing the mass profiles in 48 angular cones, whose footprints partition the sphere, we find anisotropic gas and total mass distributions that imply an angular variance of f_gas at the level of 30%. This anisotropic distribution originates from the recent formation epoch of clusters and from the strong internal baryon-to-dark-matter density bias. In the most extreme cones, f_gas can be biased high by a factor of two at R_200 in massive clusters, thereby providing a potential explanation for high f_gas measurements by Suzaku. While projection lowers this factor, there are other measurement biases that may (partially) compensate. We find that at R_200, f_gas is biased high by 20% when assuming hydrostatic equilibrium masses, i.e., neglecting the kinetic pressure, and by another ~10-20% due to the presence of density clumping. At larger radii, both measurement biases increase dramatically. While the cluster sample variance of the true f_gas decreases to a level of 5% at R_200, the sample variance that includes both measurement biases remains fairly constant at the level of 10-20%. The constant redshift evolution of f_gas within R_500 for massive clusters is encouraging for using gas masses to derive cosmological parameters. △ Less

Submitted 28 August, 2013; v1 submitted 18 September, 2012; originally announced September 2012.

Comments: Accepted by ApJ, 24 pages, 17 figures, this version additionally includes projection effects, an improved discussion of hydrostatic equilibrium masses, and compares our simulation results to more data

Constraints on Cosmic Rays, Magnetic Fields, and Dark Matter from Gamma-Ray Observations of the Coma Cluster of Galaxies with VERITAS and Fermi

Authors: T. Arlen, T. Aune, M. Beilicke, W. Benbow, A. Bouvier, J. H. Buckley, V. Bugaev, K. Byrum, A. Cannon, A. Cesarini, L. Ciupik, E. Collins-Hughes, M. P. Connolly, W. Cui, R. Dickherber, J. Dumm, A. Falcone, S. Federici, Q. Feng, J. P. Finley, G. Finnegan, L. Fortson, A. Furniss, N. Galante, D. Gall , et al. (59 additional authors not shown)

Abstract: Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated, suggesting that clusters may also be sources of very high-energy (VHE; E>100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cherenkov telescopes, with complementing Fermi-LAT o… ▽ More Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated, suggesting that clusters may also be sources of very high-energy (VHE; E>100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cherenkov telescopes, with complementing Fermi-LAT observations at GeV energies. No significant gamma-ray emission from the Coma cluster was detected. Integral flux upper limits at the 99% confidence level were measured to be on the order of (2-5)*10^-8\ ph. m^-2 s^-1 (VERITAS, >220 GeV} and ~2*10^-6 ph. m^-2 s^-1 (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain CRs and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be < 16% from VERITAS data and < 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be <50%. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of (2 - 5.5) muG, depending on the radial magnetic-field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this {renders} the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark-matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally-averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, <σv>. (abr.) △ Less

Submitted 3 August, 2012; originally announced August 2012.

Comments: ApJ in press

On the Physics of Radio Halos in Galaxy Clusters: Scaling Relations and Luminosity Functions

Authors: Fabio Zandanel, Christoph Pfrommer, Francisco Prada

Abstract: The underlying physics of giant radio halos and mini halos in galaxy clusters is still an open question, which becomes more pressing with the growing number of detections. In this paper, we explore the possibility that radio-emitting electrons are generated in hadronic cosmic ray (CR) proton interactions with ambient thermal protons of the intra-cluster medium. Our CR model derives from cosmologic… ▽ More The underlying physics of giant radio halos and mini halos in galaxy clusters is still an open question, which becomes more pressing with the growing number of detections. In this paper, we explore the possibility that radio-emitting electrons are generated in hadronic cosmic ray (CR) proton interactions with ambient thermal protons of the intra-cluster medium. Our CR model derives from cosmological hydrodynamical simulations of cluster formation and additionally accounts for CR transport in the form of CR streaming and diffusion. This opens the possibility of changing the radio halo luminosity by more than an order of magnitude on a dynamical time scale. We build a mock galaxy cluster catalog from the large MultiDark N-body LCDM simulation by adopting a phenomenological gas density model for each cluster based on X-ray measurements that matches Sunyaev-Zel'dovich (SZ) and X-ray scaling relations and luminosity function. Using magnetic field strength estimates from Faraday rotation measure studies, our model successfully reproduces the observed surface brightness profiles of giant radio halos (Coma, A2163) as well as radio mini-halos (Perseus, Ophiuchus), while obeying upper limits on the gamma-ray emission in these clusters. Our model is also able to 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 the SZ flux Y; thereby suggesting a physical solution to this apparent contradiction. For a plausible fraction of 10% radio-loud clusters, our model matches the NVSS radio-halo luminosity function. Constructing an analytical radio-halo luminosity function, we demonstrate the unique prospects for low-frequency radio surveys (such as the LOFAR Tier 1 survey) to detect ~3500 radio halos back to redshift two and to probe the underlying physics of radio halos. [abridged] △ Less

Submitted 21 November, 2013; v1 submitted 26 July, 2012; originally announced July 2012.

Comments: This article has been superseded by arXiv:1311.4793 and arXiv:1311.4795

Is dark matter with long-range interactions a solution to all small-scale problems of ΛCDM cosmology?

Authors: Laura G. van den Aarssen, Torsten Bringmann, Christoph Pfrommer

Abstract: The cold dark matter (DM) paradigm describes the large-scale structure of the universe remarkably well. However, there exists some tension with the observed abundances and internal density structures of both field dwarf galaxies and galactic satellites. Here, we demonstrate that a simple class of DM models may offer a viable solution to all of these problems simultaneously. Their key phenomenologi… ▽ More The cold dark matter (DM) paradigm describes the large-scale structure of the universe remarkably well. However, there exists some tension with the observed abundances and internal density structures of both field dwarf galaxies and galactic satellites. Here, we demonstrate that a simple class of DM models may offer a viable solution to all of these problems simultaneously. Their key phenomenological properties are velocity-dependent self-interactions mediated by a light vector messenger and thermal production with much later kinetic decoupling than in the standard case. △ Less

Submitted 6 December, 2012; v1 submitted 25 May, 2012; originally announced May 2012.

Comments: revtex4; 6 pages, 3 figures; minor changes to match published version

Journal ref: Phys. Rev. Lett. 109, 231301 (2012)

Galactic winds driven by cosmic-ray streaming

Authors: M. Uhlig, C. Pfrommer, M. Sharma, B. B. Nath, T. A. Ensslin, V. Springel

Abstract: Galactic winds are observed in many spiral galaxies with sizes from dwarfs up to the Milky Way, and they sometimes carry a mass in excess of that of newly formed stars by up to a factor of ten. Multiple driving processes of such winds have been proposed, including thermal pressure due to supernova-heating, UV radiation pressure on dust grains, or cosmic ray (CR) pressure. We here study wind format… ▽ More Galactic winds are observed in many spiral galaxies with sizes from dwarfs up to the Milky Way, and they sometimes carry a mass in excess of that of newly formed stars by up to a factor of ten. Multiple driving processes of such winds have been proposed, including thermal pressure due to supernova-heating, UV radiation pressure on dust grains, or cosmic ray (CR) pressure. We here study wind formation due to CR physics using a numerical model that accounts for CR acceleration by supernovae, CR thermalization, and advective CR transport. In addition, we introduce a novel implementation of CR streaming relative to the rest frame of the gas. We find that CR streaming drives powerful and sustained winds in galaxies with virial masses M_200 < 10^{11} Msun. In dwarf galaxies (M_200 ~ 10^9 Msun) the winds reach a mass loading factor of ~5, expel ~60 per cent of the initial baryonic mass contained inside the halo's virial radius and suppress the star formation rate by a factor of ~5. In dwarfs, the winds are spherically symmetric while in larger galaxies the outflows transition to bi-conical morphologies that are aligned with the disc's angular momentum axis. We show that damping of Alfven waves excited by streaming CRs provides a means of heating the outflows to temperatures that scale with the square of the escape speed. In larger haloes (M_200 > 10^{11} Msun), CR streaming is able to drive fountain flows that excite turbulence. For halo masses M_200 > 10^{10} Msun, we predict an observable level of H-alpha and X-ray emission from the heated halo gas. We conclude that CR-driven winds should be crucial in suppressing and regulating the first epoch of galaxy formation, expelling a large fraction of baryons, and - by extension - aid in shaping the faint end of the galaxy luminosity function. They should then also be responsible for much of the metal enrichment of the intergalactic medium. △ Less

Submitted 2 April, 2012; v1 submitted 5 March, 2012; originally announced March 2012.

Comments: 25 pages, 14 figures, accepted by MNRAS

Detection of very-high energy γ-ray emission from NGC 1275 by the MAGIC telescopes

Authors: MAGIC Collaboration, J. Aleksić, E. A. Alvarez, L. A. Antonelli, P. Antoranz, M. Asensio, M. Backes, U. Barres de Almeida, J. A. Barrio, D. Bastieri, J. Becerra González, W. Bednarek, A. Berdyugin, K. Berger, E. Bernardini, A. Biland, O. Blanch, R. K. Bock, A. Boller, G. Bonnoli, D. Borla Tridon, I. Braun, T. Bretz, A. Cañellas, E. Carmona , et al. (139 additional authors not shown)

Abstract: We report on the detection of very-high energy (VHE, E>100 GeV) gamma-ray emission from NGC 1275, the central radio galaxy of the Perseus cluster of galaxies. The source has been detected by the MAGIC telescopes with a statistical significance of 6.6 sigma above 100 GeV in 46 hr of stereo observations carried out between August 2010 and February 2011. The measured differential energy spectrum betw… ▽ More We report on the detection of very-high energy (VHE, E>100 GeV) gamma-ray emission from NGC 1275, the central radio galaxy of the Perseus cluster of galaxies. The source has been detected by the MAGIC telescopes with a statistical significance of 6.6 sigma above 100 GeV in 46 hr of stereo observations carried out between August 2010 and February 2011. The measured differential energy spectrum between 70 GeV and 500 GeV can be described by a power law with a steep spectral index of Γ=-4.1+/-0.7stat+/-0.3syst, and the average flux above 100 GeV is F_{gamma}=(1.3+/-0.2stat+/-0.3syst) x 10^-11 cm^-2 s^-1. These results, combined with the power-law spectrum measured in the first two years of observations by the Fermi-LAT above 100 MeV, with a spectral index of Gamma ~ -2.1, strongly suggest the presence of a break or cut-off around tens of GeV in the NGC 1275 spectrum. The light curve of the source above 100 GeV does not show hints of variability on a month time scale. Finally, we report on the nondetection in the present data of the radio galaxy IC 310, previously discovered by the Fermi-LAT and MAGIC. The derived flux upper limit F^{U.L.}_{gamma} (>300 GeV)=1.2 x 10^-12 cm^-2 s^-1 is a factor ~ 3 lower than the mean flux measured by MAGIC between October 2009 and February 2010, thus confirming the year time-scale variability of the source at VHE. △ Less

Submitted 5 March, 2012; v1 submitted 16 December, 2011; originally announced December 2011.

Comments: Corresponding authors are S. Lombardi, P. Colin, D. Hildebrand, and F. Zandanel. Published in A&A; 4 pages, 4 figures

Constraining Cosmic Rays and Magnetic Fields in the Perseus Galaxy Cluster with TeV observations by the MAGIC telescopes

Authors: MAGIC Collaboration, J. Aleksić, E. A. Alvarez, L. A. Antonelli, P. Antoranz, M. Asensio, M. Backes, U. Barres de Almeida, J. A. Barrio, D. Bastieri, J. Becerra González, W. Bednarek, A. Berdyugin, K. Berger, E. Bernardini, A. Biland, O. Blanch, R. K. Bock, A. Boller, G. Bonnoli, D. Borla Tridon, I. Braun, T. Bretz, A. Cañellas, E. Carmona , et al. (139 additional authors not shown)

Abstract: Galaxy clusters are being assembled today in the most energetic phase of hierarchical structure formation which manifests itself in powerful shocks that contribute to a substantial energy density of cosmic rays (CRs). Hence, clusters are expected to be luminous gamma-ray emitters since they also act as energy reservoirs for additional CR sources, such as active galactic nuclei and supernova-driven… ▽ More Galaxy clusters are being assembled today in the most energetic phase of hierarchical structure formation which manifests itself in powerful shocks that contribute to a substantial energy density of cosmic rays (CRs). Hence, clusters are expected to be luminous gamma-ray emitters since they also act as energy reservoirs for additional CR sources, such as active galactic nuclei and supernova-driven galactic winds. To detect the gamma-ray emission from CR interactions with the ambient cluster gas, we conducted the deepest to date observational campaign targeting a galaxy cluster at very high-energy gamma-rays and observed the Perseus cluster with the MAGIC Cherenkov telescopes for a total of ~85 hr of effective observing time. This campaign resulted in the detection of the central radio galaxy NGC 1275 at energies E > 100 GeV with a very steep energy spectrum. Here, we restrict our analysis to energies E > 630 GeV and detect no significant gamma-ray excess. This constrains the average CR-to-thermal pressure ratio to be <= 1-2%, depending on assumptions and the model for CR emission. Comparing these gamma-ray upper limits to predictions from cosmological cluster simulations that include CRs constrains the maximum CR acceleration efficiency at structure formation shocks to be < 50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Finally, we derive lower limits on the magnetic field distribution assuming that the Perseus radio mini-halo is generated by secondary electrons/positrons that are created in hadronic CR interactions: assuming a spectrum of E^-2.2 around TeV energies as implied by cluster simulations, we limit the central magnetic field to be > 4-9 microG, depending on the rate of decline of the magnetic field strength toward larger radii. △ Less

Submitted 4 May, 2012; v1 submitted 23 November, 2011; originally announced November 2011.

Comments: Corresponding authors are F. Zandanel, C. Pfrommer, P. Colin, A. Pinzke and S. Lombardi. Accepted for publication in A&A on 06/03/2012

Journal ref: A&A 541, A99, 2012

Observation of the Perseus cluster of galaxies with the MAGIC telescopes

Authors: Saverio Lombardi, Pierre Colin, Dorothee Hildebrand, Fabio Zandanel, Francisco Prada, Christoph Pfrommer, Anders Pinzke

Abstract: The MAGIC telescopes performed a deep observation of the central region of the Perseus galaxy cluster in stereoscopic mode between October 2009 and February 2011. The nearly 85 hr of collected data (after quality selection) represent the deepest observation of a cluster of galaxies at very high energies (VHE, E > 100 GeV) ever. The survey resulted in the detection of VHE gamma-ray emissions from i… ▽ More The MAGIC telescopes performed a deep observation of the central region of the Perseus galaxy cluster in stereoscopic mode between October 2009 and February 2011. The nearly 85 hr of collected data (after quality selection) represent the deepest observation of a cluster of galaxies at very high energies (VHE, E > 100 GeV) ever. The survey resulted in the detection of VHE gamma-ray emissions from its central galaxy NGC 1275 and from the radio galaxy IC 310. In addition, the deep survey also permits for the first time to constrain emission models predicting VHE gamma-rays from cosmic-ray acceleration in the cluster. In this contribution we report the latest MAGIC results concerning these topics. △ Less

Submitted 1 November, 2011; originally announced November 2011.

Comments: 2011 Fermi Symposium proceedings - eConf C110509

MAGIC detection of VHE Gamma-ray emission from NGC 1275 and IC 310

Authors: Dorothee Hildebrand, Saverio Lombardi, Pierre Colin, Julian Sitarek, Fabio Zandanel, Francisco Prada, Christoph Pfrommer, Anders Pinzke

Abstract: The MAGIC Cherenkov telescopes observed the Perseus cluster sky region in stereo mode for nearly 90 hr from October 2009 to February 2011. This campaign led to the discovery of very high energy Gamma-ray emission from the central radio galaxy NGC 1275 and the head-tail radio galaxy IC 310. Here we report the results on the most recent discovery of NGC 1275 which was detected at low energies in the… ▽ More The MAGIC Cherenkov telescopes observed the Perseus cluster sky region in stereo mode for nearly 90 hr from October 2009 to February 2011. This campaign led to the discovery of very high energy Gamma-ray emission from the central radio galaxy NGC 1275 and the head-tail radio galaxy IC 310. Here we report the results on the most recent discovery of NGC 1275 which was detected at low energies in the 2010/2011 data. We also present latest results on IC 310, which had been detected in the 2009/2010 data. △ Less

Submitted 27 October, 2011; v1 submitted 24 October, 2011; originally announced October 2011.

Comments: Proceedings of the 32nd International Cosmic Ray Conference (Beijing)

Observation of the Perseus galaxy cluster with the MAGIC telescopes

Authors: S. Lombardi, F. Zandanel, P. Colin, M. Doro, D. Hildebrand, F. Prada, C. Pfrommer, A. Pinzke

Abstract: The MAGIC ground-based Imaging Cherenkov experiment observed the Perseus galaxy cluster for a total of about 25 hr between November and December 2008 in single telescope mode and for nearly 90 hr between October 2009 and February 2011 in stereoscopic mode. This survey represents the deepest observation of a cluster of galaxies at very high energies ever. It resulted in the detection of the central… ▽ More The MAGIC ground-based Imaging Cherenkov experiment observed the Perseus galaxy cluster for a total of about 25 hr between November and December 2008 in single telescope mode and for nearly 90 hr between October 2009 and February 2011 in stereoscopic mode. This survey represents the deepest observation of a cluster of galaxies at very high energies ever. It resulted in the detection of the central radio galaxy NGC 1275 and the head-tail radio galaxy IC 310. It also permits for the first time to put constraints on emission models predicting gamma-rays from cosmic ray acceleration in the cluster and to investigate dark matter scenarios. Here, we will report the latest MAGIC results on these studies. △ Less

Submitted 28 September, 2011; originally announced September 2011.

Comments: 4 pages, 3 figures, contribution to the 32nd International Cosmic Ray Conference, Beijing, China, August 2011

Cosmic ray transport in galaxy clusters: implications for radio halos and gamma-rays

Authors: Christoph Pfrommer, Torsten Ensslin, Francesco Miniati, Kandaswamy Subramanian

Abstract: Observations of giant radio halos provide unambiguous evidence for the existence of cosmic ray (CR) electrons and magnetic fields in galaxy clusters. The physical mechanism generating radio halos is still heavily debated. We critically discuss the proposed models for the radio halo emission and highlight the weaknesses underlying each explanation. We present an idea how the interplay of CR propaga… ▽ More Observations of giant radio halos provide unambiguous evidence for the existence of cosmic ray (CR) electrons and magnetic fields in galaxy clusters. The physical mechanism generating radio halos is still heavily debated. We critically discuss the proposed models for the radio halo emission and highlight the weaknesses underlying each explanation. We present an idea how the interplay of CR propagation and turbulent advection selects a bimodal spatial CR distribution that is characteristic for the dynamical state of a cluster. As a result, strongly turbulent, merging clusters should have a more centrally concentrated CR energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected diffuse radio and gamma ray emission of clusters. Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR electrons. △ Less

Submitted 27 September, 2011; originally announced September 2011.

Comments: 7 pages, 5 figures, proceedings of an invited talk at "Non-Thermal Phenomena in Colliding Galaxy Clusters", held in Nice, 15-18 November 2010. Published in Mem. S.A.It

On the Cluster Physics of Sunyaev-Zel'dovich Surveys II: Deconstructing the Thermal SZ Power Spectrum

Authors: N. Battaglia, J. R. Bond, C. Pfrommer, J. L. Sievers

Abstract: Secondary anisotropies in the cosmic microwave background are a treasure-trove of cosmological information. Interpreting current experiments probing them are limited by theoretical uncertainties rather than by measurement errors. Here we focus on the secondary anisotropies resulting from the thermal Sunyaev-Zel'dovich (tSZ) effect; the amplitude of which depends critically on the average thermal p… ▽ More Secondary anisotropies in the cosmic microwave background are a treasure-trove of cosmological information. Interpreting current experiments probing them are limited by theoretical uncertainties rather than by measurement errors. Here we focus on the secondary anisotropies resulting from the thermal Sunyaev-Zel'dovich (tSZ) effect; the amplitude of which depends critically on the average thermal pressure profile of galaxy groups and clusters. To this end, we use a suite of SPH simulations that include radiative cooling, star formation, supernova feedback, and energetic feedback from active galactic nuclei (AGN). We examine in detail how the pressure profile depends on cluster radius, mass, and redshift and provide an empirical fitting function. We employ three different approaches for calculating the tSZ power spectrum: an analytical approach that uses our pressure profile fit, a semi-analytical method of pasting our pressure fit onto simulated clusters, and a direct numerical integration of our simulated volumes. We demonstrate that the detailed structure of the intracluster medium and cosmic web affect the tSZ power spectrum. In particular, the substructure and asphericity of clusters increase the tSZ power spectrum by 10-20% at ell ~2000-8000, with most of the additional power being contributed by substructures. The contributions to the power spectrum from radii larger than R_500 is ~ 20% at ell = 3000, thus clusters interiors (r < R_500) dominate the power spectrum amplitude at these angular scales. △ Less

Submitted 16 September, 2011; originally announced September 2011.

Comments: To be submitted to ApJ, 13 pages, 11 figures

On the Cluster Physics of Sunyaev-Zel'dovich Surveys I: The Influence of Feedback, Non-thermal Pressure and Cluster Shapes on Y-M Scaling Relations

Authors: N. Battaglia, J. R. Bond, C. Pfrommer, J. L. Sievers

Abstract: The utility of large Sunyaev Zel'dovich (SZ) surveys for determining cosmological parameters from cluster abundances is limited by the theoretical uncertainties in the integrated SZ-flux-to-mass relation, Y-M. We explore how non-thermal pressure and the anisotropic shape of the gas distribution of the intracluster medium (ICM) impacts Y-M scaling using a suite of SPH simulations of the cosmic web.… ▽ More The utility of large Sunyaev Zel'dovich (SZ) surveys for determining cosmological parameters from cluster abundances is limited by the theoretical uncertainties in the integrated SZ-flux-to-mass relation, Y-M. We explore how non-thermal pressure and the anisotropic shape of the gas distribution of the intracluster medium (ICM) impacts Y-M scaling using a suite of SPH simulations of the cosmic web. We contrast results for models with different treatments of entropy injection and transport, varying radiative cooling, star formation and accompanying supernova feedback, cosmic rays, and energetic feedback from active galactic nuclei (AGN). We find that the gas kinetic-to-thermal pressure ratio from internal bulk motions depends on the cluster mass, and increases in the outer-cluster due to enhanced substructure, as does the asphericity of the ICM gas. With only a ~5-10% correction to projected (observable) ellipticities, we can infer the 3D ellipticities. Our simulated Y-M-slope roughly follows the self-similar prediction, except for a steepening due to a deficit of gas in lower mass clusters at low redshift in our AGN-feedback simulations. AGN feedback enhances slightly the overall Y-M-scatter, from ~11% to ~13%, a reflection of accretion history variations due to cluster merging. If we split the cluster system into lower, middle and upper bands of both P_kin/P_th and long-to-short axis ratio, we find a ~10% effect on Y-M. Identifying observable second parameters related to internal bulk flows and anisotropy for cluster-selection to minimize Y-M scatter in a "fundamental plane" would allow tighter cosmological parameter constraints. △ Less

Submitted 16 September, 2011; originally announced September 2011.

Comments: To be submitted to ApJ, 22 pages, 22 figures

The Lyman-alpha forest in a blazar-heated Universe

Authors: Ewald Puchwein, Christoph Pfrommer, Volker Springel, Avery E. Broderick, Philip Chang

Abstract: It has been realised only recently that TeV emission from blazars can significantly heat the intergalactic medium (IGM) by pair-producing high-energy electrons and positrons, which in turn excite vigorous plasma instabilities, leading to a local dissipation of the pairs' kinetic energy. In this work, we use cosmological hydrodynamical simulations to model the impact of this blazar heating on the L… ▽ More It has been realised only recently that TeV emission from blazars can significantly heat the intergalactic medium (IGM) by pair-producing high-energy electrons and positrons, which in turn excite vigorous plasma instabilities, leading to a local dissipation of the pairs' kinetic energy. In this work, we use cosmological hydrodynamical simulations to model the impact of this blazar heating on the Lyman-alpha forest at redshifts z~2-3. We find that blazar heating produces an inverted temperature-density relation in the IGM and naturally resolves many of the problems present in previous simulations of the forest that included photoheating alone. In particular, our simulations with blazar heating simultaneously reproduce the observed effective optical depth and temperature as a function of redshift, the observed probability distribution functions of the transmitted flux, and the observed flux power spectra, over the full redshift range 2<z<3 analysed here. Additionally, by deblending the Lyman-alpha forest into a sum of thermally broadened individual lines, we find superb agreement with the observed lower cutoff of the line-width distribution and abundances of neutral hydrogen column densities. Using the most recent constraints on the cosmic ultraviolet (UV) background, this excellent agreement with observations does not require rescaling the amplitude of the UV background; a procedure that was routinely used in the past to match the observed level of transmitted flux. We also show that our blazar-heated model matches the data better than standard simulations even when such a rescaling is allowed. This concordance between Lyman-alpha data and simulations, which are based on the most recent cosmological parameters, suggests that the inclusion of blazar heating alleviates previous tensions on constraints for sigma_8 derived from Lyman-alpha measurements and other cosmological data. [abridged] △ Less

Submitted 16 February, 2012; v1 submitted 19 July, 2011; originally announced July 2011.

Comments: 18 pages, 13 figures, added an analysis in which all simulations are scaled to the observed effective optical depth, added some discussion about the local Lyman-alpha forest, accepted for publication in MNRAS

The Cosmological Impact of Luminous TeV Blazars III: Implications for Galaxy Clusters and the Formation of Dwarf Galaxies

Authors: Christoph Pfrommer, Philip Chang, Avery E. Broderick

Abstract: 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 suggest… ▽ More 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-Zel'dovich 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) △ Less

Submitted 9 May, 2012; v1 submitted 27 June, 2011; originally announced June 2011.

Comments: 22 pages, 7 figures, update that matches the final version in the ApJ

The Cosmological Impact of Luminous TeV Blazars II: Rewriting the Thermal History of the Intergalactic Medium

Authors: Philip Chang, Avery E. Broderick, Christoph Pfrommer

Abstract: The Universe is opaque to extragalactic very high-energy gamma rays (VHEGRs, E>100 GeV) because they annihilate and pair produce on the extragalactic background light. The resulting ultra-relativistic pairs are assumed to lose energy through inverse Compton scattering of CMB photons. In Broderick et al. (2011, Paper I of this three paper series), we argued that instead powerful plasma instabilitie… ▽ More The Universe is opaque to extragalactic very high-energy gamma rays (VHEGRs, E>100 GeV) because they annihilate and pair produce on the extragalactic background light. The resulting ultra-relativistic pairs are assumed to lose energy through inverse Compton scattering of CMB photons. In Broderick et al. (2011, Paper I of this three paper series), we argued that instead powerful plasma instabilities in the ultra-relativistic pair beam dissipate the kinetic energy of the TeV-generated pairs locally, heating the intergalactic medium (IGM). Here, we explore the effect of this heating upon the thermal history of the IGM. We collate the observed extragalactic VHEGR sources to determine a local VHEGR heating rate and correct for the pointed nature of VHEGR observations using Fermi observations of high and intermediate peaked BL Lacs. Because the local extragalactic VHEGR flux is dominated by TeV blazars, we tie the TeV blazar luminosity density to the quasar luminosity density, and produce a VHEGR heating rate as a function of redshift. This heating is relatively homogeneous for z<~4 with increasing spatial variation at higher redshift (order unity at z~6). This new heating process dominates photoheating at low redshift and the inclusion of TeV blazar heating qualitatively and quantitatively changes the structure and history of the IGM. TeV blazars produce a uniform volumetric heating rate that is sufficient to increase the temperature of the mean density IGM by nearly an order of magnitude, and at low densities by substantially more, naturally producing an inverted equation of state inferred by observations of the Ly-alpha forest, a feature that is difficult to reconcile with standard reionization models. Finally, we close with a discussion on the possibility of detecting this hot low-density IGM, but find that such measurements are currently not feasible. (abridged) △ Less

Submitted 9 May, 2012; v1 submitted 27 June, 2011; originally announced June 2011.

Comments: 23 pages, 14 figures, update that matches the final version in the ApJ

The Cosmological Impact of Luminous TeV Blazars I: Implications of Plasma Instabilities for the Intergalactic Magnetic Field and Extragalactic Gamma-Ray Background

Authors: Avery E. Broderick, Philip Chang, Christoph Pfrommer

Abstract: Inverse-Compton cascades initiated by energetic gamma rays (E>100 GeV) enhance the GeV emission from bright, extragalactic TeV sources. The absence of this emission from bright TeV blazars has been used to constrain the intergalactic magnetic field (IGMF), and the stringent limits placed upon the unresolved extragalactic gamma-ray background (EGRB) by Fermi has been used to argue against a large n… ▽ More Inverse-Compton cascades initiated by energetic gamma rays (E>100 GeV) enhance the GeV emission from bright, extragalactic TeV sources. The absence of this emission from bright TeV blazars has been used to constrain the intergalactic magnetic field (IGMF), and the stringent limits placed upon the unresolved extragalactic gamma-ray background (EGRB) by Fermi has been used to argue against a large number of such objects at high redshifts. However, these are predicated upon the assumption that inverse-Compton scattering is the primary energy-loss mechanism for the ultra-relativistic pairs produced by the annihilation of the energetic gamma rays on extragalactic background light photons. Here we show that for sufficiently bright TeV sources (isotropic-equivalent luminosities >10^{42} erg/s) plasma beam instabilities, specifically the "oblique" instability, present a plausible mechanism by which the energy of these pairs can be dissipated locally, heating the intergalactic medium. Since these instabilities typically grow on timescales short in comparison to the inverse-Compton cooling rate, they necessarily suppress the inverse-Compton cascades. As a consequence, this places a severe constraint upon efforts to limit the IGMF from the lack of a discernible GeV bump in TeV sources. Similarly, it considerably weakens the Fermi limits upon the evolution of blazar populations. Specifically, we construct a TeV-blazar luminosity function from those objects presently observed and find that it is very well described by the quasar luminosity function at z~0.1, shifted to lower luminosities and number densities, suggesting that both classes of sources are regulated by similar processes. Extending this relationship to higher redshifts, we show that the magnitude and shape of the EGRB above ~10 GeV is naturally reproduced with this particular example of a rapidly evolving TeV-blazar luminosity function. △ Less

Submitted 9 May, 2012; v1 submitted 27 June, 2011; originally announced June 2011.

Comments: 23 pages, 10 figures, update that matches the final version in the ApJ

A Comparison of Cosmological Codes: Properties of Thermal Gas and Shock Waves in Large Scale Structures

Authors: F. Vazza, K. Dolag, D. Ryu, G. Brunetti, C. Gheller, H. Kang, C. Pfrommer

Abstract: [...] We present results for the statistics of thermal gas and the shock wave properties for a large volume simulated with three different cosmological numerical codes: the Eulerian total variations diminishing code TVD, the Eulerian piecewise parabolic method-based code ENZO, and the Lagrangian smoothed-particle hydrodynamics code GADGET. Starting from a shared set of initial conditions, we prese… ▽ More [...] We present results for the statistics of thermal gas and the shock wave properties for a large volume simulated with three different cosmological numerical codes: the Eulerian total variations diminishing code TVD, the Eulerian piecewise parabolic method-based code ENZO, and the Lagrangian smoothed-particle hydrodynamics code GADGET. Starting from a shared set of initial conditions, we present convergence tests for a cosmological volume of side-length 100 Mpc/h, studying in detail the morphological and statistical properties of the thermal gas as a function of mass and spatial resolution in all codes. By applying shock finding methods to each code, we measure the statistics of shock waves and the related cosmic ray acceleration efficiencies, within the sample of simulations and for the results of the different approaches. We discuss the regimes of uncertainties and disagreement among codes, with a particular focus on the results at the scale of galaxy clusters. Even if the bulk of thermal and shock properties are reasonably in agreement among the three codes, yet some significant differences exist (especially between Eulerian methods and smoothed particle hydrodynamics). In particular, we report: a) differences of huge factors (10-100) in the values of average gas density, temperature, entropy, Mach number and shock thermal energy flux in the most rarefied regions of the simulations between grid and SPH methods; b) the hint of an entropy core inside clusters simulated in grid codes; c) significantly different phase diagrams of shocked cells in grid codes compared to SPH; d) sizable differences in the morphologies of accretion shocks between grid and SPH methods. △ Less

Submitted 3 August, 2011; v1 submitted 10 June, 2011; originally announced June 2011.

Comments: 30 pages, 24 figures. Version accepted by MNRAS. Only minor changes and updated references list compared to the first submitted version

Prospects of detecting gamma-ray emission from galaxy clusters: cosmic rays and dark matter annihilations

Authors: Anders Pinzke, Christoph Pfrommer, Lars Bergstrom

Abstract: We study the possibility for detecting gamma-ray emission from galaxy clusters. We consider 1) leptophilic models of dark matter (DM) annihilation that include a Sommerfeld enhancement (SFE), 2) different representative benchmark models of supersymmetric DM, and 3) cosmic ray (CR) induced pion decay. Among all clusters/groups of a flux-limited X-ray sample, we predict Virgo, Fornax and M49 to be t… ▽ More We study the possibility for detecting gamma-ray emission from galaxy clusters. We consider 1) leptophilic models of dark matter (DM) annihilation that include a Sommerfeld enhancement (SFE), 2) different representative benchmark models of supersymmetric DM, and 3) cosmic ray (CR) induced pion decay. Among all clusters/groups of a flux-limited X-ray sample, we predict Virgo, Fornax and M49 to be the brightest DM sources and find a particularly low CR-induced background for Fornax. For a minimum substructure mass given by the DM free-streaming scale, cluster halos maximize the substructure boost for which we find a factor above 1000. Since regions around the virial radius dominate the annihilation flux of substructures, the resulting surface brightness profiles are almost flat. This makes it very challenging to detect this flux with imaging atmospheric Cherenkov telescopes. Assuming cold dark matter with a substructure mass distribution down to an Earth mass and using extended Fermi upper limits, we rule out the leptophilic models in their present form in 28 clusters, and limit the boost from SFE in M49 and Fornax to be < 5. This corresponds to a limit on SFE in the Milky Way of < 3, which is too small to account for the increasing positron fraction with energy as seen by PAMELA and challenges the DM interpretation. Alternatively, if SFE is realized in Nature, this would imply a limiting substructure mass of M_lim > 10^4 M_sol - a problem for structure formation. Using individual cluster observations, it will be challenging for Fermi to constrain our selection of DM benchmark models without SFE. The Fermi upper limits are, however, closing in on our predictions for the CR flux using an analytic model based on cosmological hydrodynamical cluster simulations. We limit the CR-to-thermal pressure in nearby bright galaxy clusters of the Fermi sample to < 10% and in Norma and Coma to < 3%. △ Less

Submitted 18 October, 2011; v1 submitted 16 May, 2011; originally announced May 2011.

Comments: 43 pages, 23 figures, 10 tables. Accepted for publication in Phys. Rev. D: streamlined paper, added a paragraph about detectability to introduction, few references added, and few typos corrected

The Curious Case of Abell 2256

Authors: Tracy E. Clarke, Torsten Enßlin, Alexis Finoguenov, Huib Intema, Christoph Pfrommer, Reinout van Weeren, Huub Röttgering, Raymond Oonk

Abstract: Abell 2256 is a rich, nearby (z=0.0594) galaxy cluster that has significant evidence of merger activity. We present new radio and X-ray observations of this system. The low-frequency radio images trace the diffuse synchrotron emission of the Mpc-scale radio halo and relics as well as a number of recently discovered, more compact, steep spectrum sources. The spectral index across the relics steepen… ▽ More Abell 2256 is a rich, nearby (z=0.0594) galaxy cluster that has significant evidence of merger activity. We present new radio and X-ray observations of this system. The low-frequency radio images trace the diffuse synchrotron emission of the Mpc-scale radio halo and relics as well as a number of recently discovered, more compact, steep spectrum sources. The spectral index across the relics steepens from the north-west toward the south-east. Analysis of the spectral index gradients between low and and high-frequencies shows spectral differences away from the north-west relic edge such that the low-frequency index is significantly flatter than the high frequency spectral index near the cluster core. This trend would be consistent with an outgoing merger shock as the origin of the relic emission. New X-ray data from XMM-Newton reveal interesting structures in the intracluster medium pressure, entropy and temperature maps. The pressure maps show an overall low pressure core co-incident with the radio halo emission, while the temperature maps reveal multiple regions of cool emission within the central regions of Abell 2256. The two cold fronts in Abell 2256 both appear to have motion in similar directions. △ Less

Submitted 27 February, 2011; originally announced February 2011.

Comments: 4 pages, 3 figures, proceedings of contributed talk at "Non-Thermal Phenomena in Colliding Galaxy Clusters", held in Nice, 15-18 November 2010. To be published in Mem. S.A.It

Detection of very high energy gamma-ray emission from the Perseus cluster head-tail galaxy IC 310 by the MAGIC telescopes

Authors: MAGIC Collaboration, J. Aleksić, L. A. Antonelli, P. Antoranz, M. Backes, J. A. Barrio, D. Bastieri, J. Becerra González, W. Bednarek, A. Berdyugin, K. Berger, E. Bernardini, A. Biland, O. Blanch, R. K. Bock, A. Boller, G. Bonnoli, P. Bordas, D. Borla Tridon, V. Bosch-Ramon, D. Bose, I. Braun, T. Bretz, M. Camara, A. Cañellas , et al. (129 additional authors not shown)

Abstract: We report on the detection with the MAGIC telescopes of very high energy gamma-rays from IC 310, a head-tail radio galaxy in the Perseus galaxy cluster, observed during the interval November 2008 to February 2010. The Fermi satellite has also detected this galaxy. The source is detected by MAGIC at a high statistical significance of 7.6sigma in 20.6 hr of stereo data. The observed spectral energy… ▽ More We report on the detection with the MAGIC telescopes of very high energy gamma-rays from IC 310, a head-tail radio galaxy in the Perseus galaxy cluster, observed during the interval November 2008 to February 2010. The Fermi satellite has also detected this galaxy. The source is detected by MAGIC at a high statistical significance of 7.6sigma in 20.6 hr of stereo data. The observed spectral energy distribution is flat with a differential spectral index of -2.00 \pm 0.14. The mean flux above 300 GeV, between October 2009 and February 2010, (3.1 \pm 0.5)x10^{-12} cm^{-2} s^{-1}, corresponds to (2.5 \pm 0.4)% of Crab Nebula units. Only an upper limit, of 1.9% of Crab Nebula units above 300 GeV, was obtained with the 2008 data. This, together with strong hints (>3sigma) of flares in the middle of October and November 2009, implies that the emission is variable. The MAGIC results favour a scenario with the very high energy emission originating from the inner jet close to the central engine. More complicated models than a simple one-zone SSC scenario, e.g. multi-zone SSC, external Compton or hadronic, may be required to explain the very flat spectrum and its extension over more than three orders of magnitude in energy. △ Less

Submitted 11 September, 2010; originally announced September 2010.

Comments: 6 pages, 5 figures, submitted to ApJL

Cosmic ray transport in galaxy clusters: implications for radio halos, gamma-ray signatures, and cool core heating

Authors: Torsten A. Ensslin, Christoph Pfrommer, Francesco Miniati, Kandaswamy Subramanian

Abstract: We investigate the interplay of cosmic ray (CR) propagation and advection in galaxy clusters. Propagation in form of CR diffusion and streaming tends to drive the CR radial profiles towards being flat, with equal CR number density everywhere. Advection of CR by the turbulent gas motions tends to produce centrally enhanced profiles. We assume that the CR streaming velocity is of the order of the so… ▽ More We investigate the interplay of cosmic ray (CR) propagation and advection in galaxy clusters. Propagation in form of CR diffusion and streaming tends to drive the CR radial profiles towards being flat, with equal CR number density everywhere. Advection of CR by the turbulent gas motions tends to produce centrally enhanced profiles. We assume that the CR streaming velocity is of the order of the sound velocity. This is motivated by plasma physical arguments. The CR streaming is then usually larger than typical advection velocities and becomes comparable or lower than this only for periods with trans- and super-sonic cluster turbulence. As a consequence a bimodality of the CR spatial distribution results. Strongly turbulent, merging clusters should have a more centrally concentrated CR energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected diffuse radio and gamma-ray emission of clusters, since more centrally concentrated CR will find higher target densities for hadronic CR proton interactions, higher plasma wave energy densities for CR electron and proton re-acceleration, and stronger magnetic fields. Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR electrons. Energy dependence of the CR propagation should lead to spectral steepening of dying radio halos. Furthermore, we show that the interplay of CR diffusion with advection implies first order CR re-acceleration in the pressure-stratified atmospheres of galaxy clusters. Finally, we argue that CR streaming could be important in turbulent cool cores of galaxy clusters since it heats preferentially the central gas with highest cooling rate. △ Less

Submitted 9 December, 2010; v1 submitted 27 August, 2010; originally announced August 2010.

Comments: 21 pages, 8 figures, 1 table, revised manuscript

Dark Matter decay and annihilation in the Local Universe: CLUES from Fermi

Authors: A. J. Cuesta, T. E. Jeltema, F. Zandanel, S. Profumo, F. Prada, G. Yepes, A. Klypin, Y. Hoffman, S. Gottloeber, J. Primack, M. A. Sanchez-Conde, C. Pfrommer

Abstract: We present all-sky simulated Fermi maps of gamma-rays from dark matter decay and annihilation in the Local Universe. The dark matter distribution is obtained from a constrained cosmological simulation of the neighboring large-scale structure provided by the CLUES project. The dark matter fields of density and density squared are then taken as an input for the Fermi observation simulation tool to p… ▽ More We present all-sky simulated Fermi maps of gamma-rays from dark matter decay and annihilation in the Local Universe. The dark matter distribution is obtained from a constrained cosmological simulation of the neighboring large-scale structure provided by the CLUES project. The dark matter fields of density and density squared are then taken as an input for the Fermi observation simulation tool to predict the gamma-ray photon counts that Fermi would detect in 5 years of all-sky survey for given dark matter models. Signal-to-noise sky maps have also been obtained by adopting the current Galactic and isotropic diffuse background models released by the Fermi collaboration. We point out the possibility for Fermi to detect a dark matter gamma-ray signal in local extragalactic structures. In particular, we conclude here that Fermi observations of nearby clusters (e.g. Virgo and Coma) and filaments are expected to give stronger constraints on decaying dark matter compared to previous studies. As an example, we find a significant signal-to-noise ratio in dark matter models with a decay rate fitting the positron excess as measured by PAMELA. This is the first time that dark matter filaments are shown to be promising targets for indirect detection of dark matter. On the other hand, the prospects for detectability of annihilating dark matter in local extragalactic structures are less optimistic even with extreme cross-sections. We make the dark matter density and density squared maps available online at http://www.clues-project.org/articles/darkmattermaps.html △ Less

Submitted 16 January, 2012; v1 submitted 20 July, 2010; originally announced July 2010.

Comments: 7 pages, 2 figures, 1 table. Matches version published in ApJ Letters. High-resolution version of Figure 1 together with FITS files of the dark matter density and density squared maps at: http://www.clues-project.org/articles/darkmattermaps.html ; Version 3 includes 1-page ApJL Erratum: S/N values corrected, conclusions now point to a better detectability of Virgo and Coma in gamma-rays

Radio Galaxy NGC 1265 unveils the Accretion Shock onto the Perseus Galaxy Cluster

Authors: Christoph Pfrommer, Tom W. Jones

Abstract: We present a consistent 3D model for the head-tail radio galaxy NGC 1265 that explains the complex radio morphology and spectrum by a past passage of the galaxy and radio bubble through a shock wave. Using analytical solutions to the full Riemann problem and hydrodynamical simulations, we study how this passage transformed the plasma bubble into a toroidal vortex ring. Adiabatic compression of the… ▽ More We present a consistent 3D model for the head-tail radio galaxy NGC 1265 that explains the complex radio morphology and spectrum by a past passage of the galaxy and radio bubble through a shock wave. Using analytical solutions to the full Riemann problem and hydrodynamical simulations, we study how this passage transformed the plasma bubble into a toroidal vortex ring. Adiabatic compression of the aged electron population causes it to be energized and to emit low-surface brightness and steep-spectrum radio emission. The large infall velocity of NGC 1265 and the low Faraday rotation measure values and variance of the jet strongly argue that this transformation was due to the accretion shock onto Perseus situated roughly at R_200. Estimating the volume change of the radio bubble enables inferring a shock Mach number of M = 4.2_{-1.2}^{+0.8}, a density jump of 3.4_{-0.4}^{+0.2}, a temperature jump of 6.3_{-2.7}^{+2.5}, and a pressure jump of 21.5 +/- 10.5 while allowing for uncertainties in the equation of state of the radio plasma and volume of the torus. Extrapolating X-ray profiles, we obtain upper limits on the gas temperature and density in the infalling warm-hot intergalactic medium of kT < 0.4 keV and n < 5e-5 / cm^3. The orientation of the ellipsoidally shaped radio torus in combination with the direction of the galaxy's head and tail in the plane of the sky is impossible to reconcile with projection effects. Instead, this argues for post-shock shear flows that have been caused by curvature in the shock surface with a characteristic radius of 850 kpc. The energy density of the shear flow corresponds to a turbulent-to-thermal energy density of 14%. The shock-injected vorticity might be important in generating and amplifying magnetic fields in galaxy clusters. Future LOFAR observations of head-tail galaxies can be complementary probes of accretion shocks onto galaxy clusters. △ Less

Submitted 8 February, 2011; v1 submitted 20 April, 2010; originally announced April 2010.

Comments: 14 pages, 4 figures, ApJ, in print; v3: typos corrected to match the published version; v2: improved presentation, added 2D numerical simulations and exact solution to the 1D Riemann problem of a shock overrunning a spherical bubble that gets transformed into a vortex ring

Simulations of the Sunyaev-Zel'dovich Power Spectrum with AGN Feedback

Authors: N. Battaglia, J. R. Bond, C. Pfrommer, J. L. Sievers, D. Sijacki

Abstract: We explore how radiative cooling, supernova feedback, cosmic rays and a new model of the energetic feedback from active galactic nuclei (AGN) affect the thermal and kinetic Sunyaev-Zel'dovich (SZ) power spectra. To do this, we use a suite of hydrodynamical TreePM-SPH simulations of the cosmic web in large periodic boxes and tailored higher resolution simulations of individual galaxy clusters. Our… ▽ More We explore how radiative cooling, supernova feedback, cosmic rays and a new model of the energetic feedback from active galactic nuclei (AGN) affect the thermal and kinetic Sunyaev-Zel'dovich (SZ) power spectra. To do this, we use a suite of hydrodynamical TreePM-SPH simulations of the cosmic web in large periodic boxes and tailored higher resolution simulations of individual galaxy clusters. Our AGN feedback simulations match the recent universal pressure profile and cluster mass scaling relations of the REXCESS X-ray cluster sample better than previous analytical or numerical approaches. For multipoles $\ell\lesssim 2000$, our power spectra with and without enhanced feedback are similar, suggesting theoretical uncertainties over that range are relatively small, although current analytic and semi-analytic approaches overestimate this SZ power. We find the power at high 2000-10000 multipoles which ACT and SPT probe is sensitive to the feedback prescription, hence can constrain the theory of intracluster gas, in particular for the highly uncertain redshifts $>0.8$. The apparent tension between $σ_8$ from primary cosmic microwave background power and from analytic SZ spectra inferred using ACT and SPT data is lessened with our AGN feedback spectra. △ Less

Submitted 3 June, 2010; v1 submitted 22 March, 2010; originally announced March 2010.

Comments: 9 pages, 4 figures, 1 table. This modified version has been submitted to ApJ, expanding upon our original letter-size version, in response to comments we received. We have extended the discussion of our AGN energy input requirements; added kSZ power spectra; made MCMC comparisons with SPT data, shown in an extra figure, as well as the ACT data we originally used; added a few more references

Simulating the gamma-ray emission from galaxy clusters: a universal cosmic ray spectrum and spatial distribution

Authors: Anders Pinzke, Christoph Pfrommer

Abstract: Entering a new era of high-energy gamma-ray experiments, there is an exciting quest for the first detection of gamma-ray emission from clusters of galaxies. To complement these observational efforts, we use high-resolution simulations of a broad sample of galaxy clusters, and follow self-consistent cosmic ray (CR) physics using an improved spectral description. We study CR proton spectra as well a… ▽ More Entering a new era of high-energy gamma-ray experiments, there is an exciting quest for the first detection of gamma-ray emission from clusters of galaxies. To complement these observational efforts, we use high-resolution simulations of a broad sample of galaxy clusters, and follow self-consistent cosmic ray (CR) physics using an improved spectral description. We study CR proton spectra as well as the different contributions of the pion decay and inverse Compton emission to the total flux and present spectral index maps. We find a universal spectrum of the CR component in clusters with surprisingly little scatter across our cluster sample. The spatial CR distribution also shows approximate universality; it depends however on the cluster mass. This enables us to derive a semi-analytic model for both, the distribution of CRs as well as the pion-decay gamma-ray emission that results from hadronic CR interactions with ambient gas protons. In addition, we provide an analytic framework for the inverse Compton emission that is produced by shock-accelerated CR electrons and valid in the full gamma-ray energy range. Combining the complete sample of the brightest X-ray clusters observed by ROSAT with our gamma-ray scaling relations, we identify the brightest clusters for the gamma-ray space telescope Fermi and current imaging air Cherenkov telescopes (MAGIC, HESS, VERITAS). We reproduce the result in Pfrommer (2008), but provide somewhat more conservative predictions for the fluxes in the energy regimes of Fermi and imaging air Cherenkov telescopes when accounting for the bias of `artificial galaxies' in cosmological simulations. We find that it will be challenging to detect cluster gamma-ray emission with Fermi after the second year but this mission has the potential of constraining interesting values of the shock acceleration efficiency after several years of surveying. △ Less

Submitted 18 July, 2010; v1 submitted 28 January, 2010; originally announced January 2010.

Comments: 33 pages, 25 figures. Accepted for publication in MNRAS: Typos corrected and primary IC analysis now includes the Klein-Nishina effect and a simple-to-use semi-analytic formula

Detecting the orientation of magnetic fields in galaxy clusters

Authors: Christoph Pfrommer, L. Jonathan Dursi

Abstract: 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 cluster's thermal history; that some have not since cooled and collapsed is a mystery. In a seemingly unrelated pu… ▽ More 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 cluster's 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. △ Less

Submitted 28 September, 2011; v1 submitted 13 November, 2009; originally announced November 2009.

Comments: 48 pages, 21 figures, revised version to match published article in Nature Physics, high-resolution version available at http://www.cita.utoronto.ca/~pfrommer/Publications/pfrommer-dursi.pdf

Journal ref: Nature Phys.6:520, 2010

MAGIC Gamma-ray Telescope Observation of the Perseus Cluster of Galaxies: implications for cosmic rays, dark matter, and NGC1275

Authors: The MAGIC Collaboration, J. Aleksić, L. A. Antonelli, P. Antoranz, M. Backes, C. Baixeras, S. Balestra, J. A. Barrio, D. Bastieri, J. Becerra González, W. Bednarek, A. Berdyugin, K. Berger, E. Bernardini, A. Biland, R. K. Bock, G. Bonnoli, P. Bordas, D. Borla Tridon, V. Bosch-Ramon, D. Bose, I. Braun, T. Bretz, D. Britzger, M. Camara , et al. (130 additional authors not shown)

Abstract: The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, an… ▽ More The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and the central radio galaxy NGC1275. Results are compatible with cosmological cluster simulations for the cosmic-ray-induced gamma-ray emission, constraining the average cosmic ray-to-thermal pressure to <4% for the cluster core region (<8% for the entire cluster). Using simplified assumptions adopted in earlier work (a power-law spectrum with an index of -2.1, constant cosmic ray-to-thermal pressure for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow formation), we would limit the ratio of cosmic ray-to-thermal energy to E_CR/E_th<3%. The upper limit also translates into a level of gamma-ray emission from possible annihilations of the cluster dark matter (the dominant mass component) that is consistent with boost factors of ~10^4 for the typically expected dark matter annihilation-induced emission. Finally, the upper limits obtained for the gamma-ray emission of the central radio galaxy NGC1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to the extremely large Doppler factors required for the jet, a one-zone synchrotron self-Compton model is implausible in this case. We reproduce the observed spectral energy density by using the structured jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio galaxies. △ Less

Submitted 21 January, 2010; v1 submitted 17 September, 2009; originally announced September 2009.

Comments: 15 pages, 5 figures; accepted for publication in ApJ; updated to match the journal final version

Journal ref: Astrophys.J.710:634-647,2010

Gamma-rays from dark matter annihilations strongly constrain the substructure in halos

Authors: Anders Pinzke, Christoph Pfrommer, Lars Bergstrom

Abstract: 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 standard value to match the observations, we show that it immediately implies an observable level of gamma-ray emission for the Fermi telescope from… ▽ More 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 standard 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. △ Less

Submitted 20 November, 2009; v1 submitted 12 May, 2009; originally announced May 2009.

Comments: 5 pages, 2 figures (published version; minor corrections to figures and result, equation added)

Journal ref: Phys.Rev.Lett.103:181302,2009

Clusters and Large-Scale Structure: the Synchrotron Keys

Authors: L. Rudnick, P. Alexander, H. Andernach, N. Battaglia, S. Brown, Gf. Brunetti, J. Burns, T. Clarke, K. Dolag, D. Farnsworth, G. Giovannini, E. Hallman, M. Johnston-Hollitt, T. W. Jones, H. Kang, N. Kassim, A. Kravtsov, J. Lazio, C. Lonsdale, B. McNamara, S. Myers, F. Owen, C. Pfrommer, D. Ryu, C. Sarazin , et al. (3 additional authors not shown)

Abstract: For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studie… ▽ More For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studies at other wavebands. Over the next decade, they will fill essential gaps in both cluster astrophysics and the cosmological growth of structure in the universe, especially where the signatures of shocks and turbulence, or even the underlying thermal plasma itself, are otherwise undetectable. Simultaneously, synchrotron studies offer a unique tool for exploring the fundamental question of the origins of cosmic magnetic fields. This work will be based on the new generation of m/cm-wave radio telescopes now in construction, as well as major advances in the sophistication of 3-D MHD simulations. △ Less

Submitted 4 March, 2009; originally announced March 2009.

Comments: 8 pages, 5 figures. White Paper submitted to 2010 Decadal Survey

Galaxy Cluster Astrophysics and Cosmology: Questions and Opportunities for the Coming Decade

Authors: S. T. Myers, C. Pfrommer, J. Aguirre, J. R. Bond, J. O. Burns, T. Clarke, M. Devlin, A. Evrard, S. Golwala, S. Habib, K. Heitmann, W. L. Holzapfel, N. E. Kassim, A. Kravtsov, A. T. Lee, M. Markevich, D. Marrone, D. Nagai, L. Page, E. Pierpaoli, L. Rudnick, J. Sievers, G. Taylor, M. Voit

Abstract: We are learning much about how structure forms, in particular how clusters as nodes in the cosmic web evolve and accrete matter, and about the physical processes within these objects. In the next decade, the study of clusters will enable us to tackle important questions regarding the nature of Dark Matter and Dark Energy, how clusters co-evolve with super-massive black holes at their centers, an… ▽ More We are learning much about how structure forms, in particular how clusters as nodes in the cosmic web evolve and accrete matter, and about the physical processes within these objects. In the next decade, the study of clusters will enable us to tackle important questions regarding the nature of Dark Matter and Dark Energy, how clusters co-evolve with super-massive black holes at their centers, and to advance our knowledge about fundamental plasma astrophysics. This science white paper outlines the key questions and research opportunities in cluster astrophysics that are emerging in the coming decade and beyond, and serves as an overview to other cluster related white papers. △ Less

Submitted 3 March, 2009; originally announced March 2009.

Comments: 8 pages, 4 figures. Science White Paper submitted to CFP and GCT panels of Astro2010

Exploring the magnetized cosmic web through low frequency radio emission

Authors: N. Battaglia, C. Pfrommer, J. L. Sievers, J. R. Bond, T. A. Ensslin

Abstract: Recent improvements in the capabilities of low frequency radio telescopes provide a unique opportunity to study thermal and non-thermal properties of the cosmic web. We argue that the diffuse, polarized emission from giant radio relics traces structure formation shock waves and illuminates the large-scale magnetic field. To show this, we model the population of shock-accelerated relativistic ele… ▽ More Recent improvements in the capabilities of low frequency radio telescopes provide a unique opportunity to study thermal and non-thermal properties of the cosmic web. We argue that the diffuse, polarized emission from giant radio relics traces structure formation shock waves and illuminates the large-scale magnetic field. To show this, we model the population of shock-accelerated relativistic electrons in high-resolution cosmological simulations of galaxy clusters and calculate the resulting radio synchrotron emission. We find that individual shock waves correspond to localized peaks in the radio surface brightness map which enables us to measure Mach numbers for these shocks. We show that the luminosities and number counts of the relics strongly depend on the magnetic field properties, the cluster mass and dynamical state. By suitably combining different cluster data, including Faraday rotation measures, we are able to constrain some macroscopic parameters of the plasma at the structure formation shocks, such as models of turbulence. We also predict upper limits for the properties of the warm-hot intergalactic medium, such as its temperature and density. We predict that the current generation of radio telescopes (LOFAR, GMRT, MWA, LWA) have the potential to discover a substantially larger sample of radio relics, with multiple relics expected for each violently merging cluster. Future experiments (SKA) should enable us to further probe the macroscopic parameters of plasma physics in clusters. △ Less

Submitted 7 August, 2009; v1 submitted 19 June, 2008; originally announced June 2008.

Comments: 18 pages, 13 figures, accepted to MNRAS Oct 23 2008, minor changes, full resolution version available at http://www.astro.utoronto.ca/~battaglia/Publications/

Journal ref: Mon.Not.Roy.Astron.Soc.393:1073-1089,2009

Simulations of cosmic ray feedback by AGN in galaxy clusters

Authors: D. Sijacki, C. Pfrommer, V. Springel, T. A. Ensslin

Abstract: We investigate a numerical model for AGN feedback where for the first time a relativistic particle population in AGN-inflated bubbles is followed within a full cosmological context. In our high-resolution simulations of galaxy cluster formation, we assume that BH accretion is accompanied by energy feedback that occurs in two different modes, depending on the accretion rate itself. Unlike in prev… ▽ More We investigate a numerical model for AGN feedback where for the first time a relativistic particle population in AGN-inflated bubbles is followed within a full cosmological context. In our high-resolution simulations of galaxy cluster formation, we assume that BH accretion is accompanied by energy feedback that occurs in two different modes, depending on the accretion rate itself. Unlike in previous work, we inject a non-thermal particle population of relativistic protons into the AGN bubbles, instead of adopting a purely thermal heating. We then follow the subsequent evolution of the cosmic ray (CR) plasma inside the bubbles, considering both its hydrodynamical interactions and dissipation processes relevant for the CR population. Due to the different buoyancy of relativistic plasma and the comparatively long CR dissipation timescale we find substantial changes in the evolution of clusters as a result of CR feedback. In particular, the non-thermal population can provide significant pressure support in central cluster regions at low thermal temperatures, providing a natural explanation for the decreasing temperature profiles found in cool core clusters. At the same time, the morphologies of the bubbles and of the induced X-ray cavities show a striking similarity to observational findings. AGN feedback with CRs also proves efficient in regulating cluster cooling flows so that the total baryon fraction in stars becomes limited to realistic values of the order of 10%. We find that the partial CR support of the intracluster gas also affects the expected signal of the thermal Sunyaev-Zel'dovich effect, with typical modifications of the integrated Compton-y parameter within the virial radius of the order of 10%. [Abridged] △ Less

Submitted 28 May, 2008; v1 submitted 21 January, 2008; originally announced January 2008.

Comments: 15 pages, 7 figures, minor revisions, MNRAS accepted

Draping of Cluster Magnetic Fields over Bullets and Bubbles -- Morphology and Dynamic Effects

Authors: L. J. Dursi, C. Pfrommer

Abstract: High-resolution X-ray observations have revealed cavities and `cold fronts' with sharp edges in temperature, density, and metallicity within galaxy clusters. Their presence poses a puzzle since these features are not expected to be hydrodynamically stable, or to remain sharp in the presence of diffusion. However, a moving core or bubble in even a very weakly magnetized plasma necessarily sweeps… ▽ More High-resolution X-ray observations have revealed cavities and `cold fronts' with sharp edges in temperature, density, and metallicity within galaxy clusters. Their presence poses a puzzle since these features are not expected to be hydrodynamically stable, or to remain sharp in the presence of diffusion. However, a moving core or bubble in even a very weakly magnetized plasma necessarily sweeps up enough magnetic field to build up a dynamically important sheath around the object; the layer's strength is set by a competition between `plowing up' of field and field lines slipping around the core. We show that a two-dimensional approach to the problem is quite generally not possible. In three dimensions, we show with analytic arguments and in numerical experiments, that this magnetic layer modifies the dynamics of a plunging core, greatly modifies the effects of hydrodynamic instabilities on the core, modifies the geometry of stripped material, and even slows the fall of the core through magnetic tension. We derive an expression for the maximum magnetic field strength, the thickness of the layer, and the opening angle of the magnetic wake. The morphology of the magnetic draping layer implies the suppression of thermal conduction across the layer, thus conserving strong temperature gradients over the contact surface. The intermittent amplification of the magnetic field as well as the injection of MHD turbulence in the wake of the core is identified to be due to vorticity generation within the magnetic draping layer. These results have important consequences for understanding the physical properties and the complex gasdynamical processes of the intra-cluster medium, and apply quite generally to motions through other magnetized environments, e.g., the ISM. △ Less

Submitted 8 January, 2008; v1 submitted 1 November, 2007; originally announced November 2007.

Comments: For version of this paper with interactive 3D graphics and full-resolution figures, see http://www.cita.utoronto.ca/~ljdursi/draping/ . 19p, 26 figures, emulateapj format. Version accepted by ApJ - new references, improved figures

Simulating cosmic rays in clusters of galaxies - II. A unified scheme for radio halos and relics with predictions of the gamma-ray emission

Authors: Christoph Pfrommer, Torsten A. Ensslin, Volker Springel

Abstract: The thermal plasma of galaxy clusters lost most of its information on how structure formation proceeded as a result of dissipative processes. In contrast, non-equilibrium distributions of cosmic rays (CR) preserve the information about their injection and transport processes and provide thus a unique window of current and past structure formation processes. This information can be unveiled by ob… ▽ More The thermal plasma of galaxy clusters lost most of its information on how structure formation proceeded as a result of dissipative processes. In contrast, non-equilibrium distributions of cosmic rays (CR) preserve the information about their injection and transport processes and provide thus a unique window of current and past structure formation processes. This information can be unveiled by observations of non-thermal radiative processes, including radio synchrotron, hard X-ray, and gamma-ray emission. To explore this, we use high-resolution simulations of a sample of galaxy clusters spanning a mass range of about two orders of magnitudes, and follow self-consistent CR physics on top of the radiative hydrodynamics. We model CR electrons that are accelerated at cosmological structure formation shocks and those that are produced in hadronic interactions of CRs with ambient gas protons. We find that CR protons trace the time integrated non-equilibrium activities of clusters while shock-accelerated CR electrons probe current accretion and merging shock waves. The resulting inhomogeneous synchrotron emission matches the properties of observed radio relics. We propose a unified model for the generation of radio halos. Giant radio halos are dominated in the centre by secondary synchrotron emission with a transition to the synchrotron radiation emitted from shock-accelerated electrons in the cluster periphery. This model is able to explain the observed correlation of mergers with radio halos, the larger peripheral variation of the spectral index, and the large scatter in the scaling relation between cluster mass and synchrotron emission. Future low-frequency radio telescopes (LOFAR, GMRT, MWA, LWA) are expected to probe the accretion shocks of clusters. [abridged] △ Less

Submitted 6 February, 2008; v1 submitted 12 July, 2007; originally announced July 2007.

Comments: 32 pages, 19 figures, small changes to match the version to be published by MNRAS, full resolution version available at http://www.cita.utoronto.ca/~pfrommer/Publications/CRs_non-thermal.pdf

Simulating cosmic rays in clusters of galaxies - III. Non-thermal scaling relations and comparison to observations

Authors: Christoph Pfrommer

Abstract: Complementary views of galaxy clusters in the radio synchrotron, hard X-ray inverse Compton, and high-energy gamma-ray regimes are critical in calibrating them as high-precision cosmological probes. We present predictions for scaling relations between cluster mass and these non-thermal observables. To this end, we use high-resolution simulations of a sample of galaxy clusters spanning a mass ran… ▽ More Complementary views of galaxy clusters in the radio synchrotron, hard X-ray inverse Compton, and high-energy gamma-ray regimes are critical in calibrating them as high-precision cosmological probes. We present predictions for scaling relations between cluster mass and these non-thermal observables. To this end, we use high-resolution simulations of a sample of galaxy clusters spanning a mass range of almost two orders of magnitudes, and follow self-consistent cosmic ray physics on top of the radiative hydrodynamics. Calibrating the magnetic fields of our model with Faraday rotation measurements (RM), the synchrotron emission of our relativistic electron populations matches the radio synchrotron luminosities and morphologies of observed giant radio halos and mini-halos surprisingly well. Using the complete sample of the brightest X-ray clusters observed by ROSAT in combination with our gamma-ray scaling relation, we predict GLAST will detect about ten clusters allowing for Eddington bias due to the scatter in the scaling relation. The brightest gamma-ray clusters are Ophiuchus, Fornax, Coma, A3627, Perseus, and Centaurus. We provide an absolute lower flux limit for the gamma-ray emission of Coma in the hadronic model which can be made tighter for magnetic field values derived from RM values to match the GLAST sensitivity, providing thus a unique test for the possible hadronic origin of radio halos. Our predicted hard X-ray emission, due to inverse Compton emission of shock accelerated and hadronically produced relativistic electrons, falls short of the detections in Coma and Perseus by a factor of 50. This casts doubts on the inverse Compton interpretation and reinforces the known discrepancy of magnetic field estimates from Faraday RM values and those obtained by combining synchrotron and inverse Compton emission. [abridged] △ Less

Submitted 6 February, 2008; v1 submitted 11 July, 2007; originally announced July 2007.

Comments: 15 pages, 4 figures, to be published by MNRAS, added concept of minimum gamma-ray emission from clusters in the hadronic model and improved discussion on hard X-ray emission, full resolution version available at http://www.cita.utoronto.ca/~pfrommer/Publications/CRs_scaling.pdf