We investigate the impact of cosmic rays (CR) and different modes of CR transport on the properties of Milky Way-mass galaxies in cosmological magneto-hydrodynamical simulations in the context of the AURIGA project. We systematically study how advection, anisotropic diffusion and additional Alfv\'en-wave cooling affect the galactic disc and the circum-galactic medium (CGM). Global properties such as stellar mass and star formation rate vary little between simulations with and without various CR transport physics, whereas structural properties such as disc sizes, CGM densities or temperatures can be strongly affected. In our simulations, CRs affect the accretion of gas onto galaxies by modifying the CGM flow structure. This alters the angular momentum distribution which manifests itself as a difference in stellar and gaseous disc size. The strength of this effect depends on the CR transport model: CR advection results in the most compact discs while the Alfv\'en-wave model resembles more the AURIGA model. The advection and diffusion models exhibit large ($r\sim50$ kpc) CR pressure-dominated gas haloes causing a smoother and partly cooler CGM. The additional CR pressure smoothes small-scale density peaks and compensates for the missing thermal pressure support at lower CGM temperatures. In contrast, the Alfv\'en-wave model is only CR pressure dominated at the disc-halo interface and only in this model the gamma-ray emission from hadronic interactions agrees with observations. In contrast to previous findings, we conclude that details of CR transport are critical for accurately predicting the impact of CR feedback on galaxy formation.

In the outskirts of the intracluster medium (ICM) in galaxy clusters, the temperature decreases with radius. Due to the weakly collisional nature of the plasma, these regions are susceptible to the magneto-thermal instability (MTI), which can sustain turbulence and provide turbulent pressure support in the ICM. This instability arises due to heat conduction directed along the magnetic field, with a heat conductivity which is normally assumed to be given by the Spitzer value. Recent numerical studies of the ion mirror and the electron whistler instability using particle-in-cell codes have shown that microscale instabilities can lead to a reduced value for the heat conductivity in the ICM. This could in turn influence the efficiency with which the MTI drives turbulence. In this paper we investigate the influence of reduced heat transport on the nonlinear evolution of the MTI by employing a subgrid model that mimics the influence of the mirror instability on the heat conductivity. We use this subgrid model to assess the effect of microscales on the large scale dynamics of the ICM. We find that the nonlinear saturation of the MTI is surprisingly robust. Over a factor of $\sim 10^3$ the thermal-to-magnetic pressure ratio and collisionality we find at most modest changes to the saturation of the MTI with respect to reference simulations where heat transport is unsuppressed.

Supernova remnants (SNRs) are believed to be the source of Galactic cosmic rays (CRs). SNR shocks accelerate CR protons and electrons which reveal key insights into the non-thermal physics by means of their synchrotron and $\gamma$-ray emission. The remnant SN 1006 is an ideal particle acceleration laboratory because it is observed across all electromagnetic wavelengths from radio to $\gamma$-rays. We perform three-dimensional (3D) magnetohydrodynamics (MHD) simulations where we include CR protons and follow the CR electron spectrum. By matching the observed morphology and non-thermal spectrum of SN 1006 in radio, X-rays and $\gamma$-rays, we gain new insight into CR electron acceleration and magnetic field amplification. 1. We show that a mixed leptonic-hadronic model is responsible for the $\gamma$-ray radiation: while leptonic inverse-Compton emission and hadronic pion-decay emission contribute equally at GeV energies observed by Fermi, TeV energies observed by imaging air Cherenkov telescopes are hadronically dominated. 2. We show that quasi-parallel acceleration (i.e., when the shock propagates at a narrow angle to the upstream magnetic field) is preferred for CR electrons and that the electron acceleration efficiency of radio-emitting GeV electrons at quasi-perpendicular shocks is suppressed at least by a factor ten. This poses a challenge for the interpretation of plasma particle-in-cell simulations of particle acceleration at shocks. 3. To match the radial emission profiles and the $\gamma$-ray spectrum, we require a volume-filling, turbulently amplified magnetic field and that the Bell-amplified magnetic field is damped in the immediate post-shock region. Our work connects micro-scale plasma physics simulations to the scale of SNRs.

Galactic cosmic rays (CRs) are believed to be accelerated at supernova remnant (SNR) shocks. In the hadronic scenario the TeV gamma-ray emission from SNRs originates from decaying pions that are produced in collisions of the interstellar gas and CRs. Using CR-magnetohydrodynamic simulations, we show that magnetic obliquity-dependent shock acceleration is able to reproduce the observed TeV gamma-ray morphology of SNRs such as Vela Jr. and SN1006 solely by varying the magnetic morphology. This implies that gamma-ray bright regions result from quasi-parallel shocks (i.e., when the shock propagates at a narrow angle to the upstream magnetic field), which are known to efficiently accelerate CR protons, and that gamma-ray dark regions point to quasi-perpendicular shock configurations. Comparison of the simulated gamma-ray morphology to observations allows us to constrain the magnetic coherence scale $\lambda_B$ around Vela Jr. and SN1006 to $\lambda_B \simeq 13_{-4.3}^{+13}$ pc and $\lambda_B >200_{-10}^{+80}$ pc, respectively, where the ambient magnetic field of SN1006 is consistent with being largely homogeneous. We find consistent pure hadronic and mixed hadronic-leptonic models that both reproduce the multi-frequency spectra from the radio to TeV gamma rays and match the observed gamma-ray morphology. Finally, to capture the propagation of a SNR shock in a clumpy interstellar medium, we study the interaction of a shock with a dense cloud with numerical simulations and analytics. We construct an analytical gamma-ray model for a core collapse SNR propagating through a structured interstellar medium, and show that the gamma-ray luminosity is only biased by 30% for realistic parameters.

High redshift, massive halos are observed to have sustained, high star formation rates, which require that the amount of cold gas in the halo is continuously replenished. The cooling time scale for the hot virialized halo gas is too long to provide the source of cold gas. Supersonic, cold streams have been invoked as a mechanism for feeding massive halos at high redshift and deliver the cold gas required for continued star formation at the rates observed. This mechanism for replenishing the cold gas reservoir is motivated by some cosmological simulations. However, the cold streams are likely to be subject to the supersonic version of the Kelvin-Helmholtz instability (KHI), which eventually leads to stream disruption. Cosmological simulations have yet to obtain the spatial resolution required for understanding the detailed stability properties of cold streams. In this paper, we consider instead an idealized model of magnetized cold streams that we spatially resolve. Using linear theory we show how magnetic fields with dynamically important field strengths do not inhibit the KHI but rather enhance its growth rate. We perform nonlinear simulations of magnetized stream disruption and find that magnetic fields can nevertheless increase stream survival times by suppressing the mixing rate of cold gas with the circumgalactic medium. We find that magnetic fields can allow streams to survive $\sim 2-8$ times longer and, consequently, that streams $\sim 2-8$ times thinner can reach the central galaxy if the magnetic field strength is $\sim 0.3-0.8 \mu$G.

We propose two effective parameters that fully characterise galactic-scale structure formation at high redshifts ($z\gtrsim5$) for a variety of dark matter (DM) models that have a primordial cutoff in the matter power spectrum. Our description is within the recently proposed ETHOS framework and includes standard thermal Warm DM (WDM) and models with dark acoustic oscillations (DAOs). To define and explore this parameter space, we use high-redshift zoom-in simulations that cover a wide range of non-linear scales from those where DM should behave as CDM ($k\sim10\,h\,{\rm Mpc}^{-1}$), down to those characterised by the onset of galaxy formation ($k\sim500\,h\,{\rm Mpc}^{-1}$). We show that the two physically motivated parameters $h_{\rm peak}$ and $k_{\rm peak}$, the amplitude and scale of the first DAO peak, respectively, are sufficient to parametrize the linear matter power spectrum and classify the DM models as belonging to effective non-linear structure formation regions. These are defined by their relative departure from Cold DM ($k_{\rm peak}\rightarrow\infty$) and WDM ($h_{\rm peak}=0$) according to the non-linear matter power spectrum and halo mass function. We identify a region where the DAOs still leave a distinct signature from WDM down to $z=5$, while a large part of the DAO parameter space is shown to be degenerate with WDM. Our framework can then be used to seamlessly connect a broad class of particle DM models to their structure formation properties at high redshift without the need of additional $N$-body simulations.

We study the longitudinal stability of beam-plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully-relativistic way. We do not assume the existence of Langmuir waves, and we focus on the rate and the mechanism by which waves are excited in such systems from an initial perturbation. We derive the structure of the unstable modes and compute an analytical approximation for their growth rates. Our computation is limited to dilute and cold beams, and shows an excellent agreement with particle-in-cell simulations performed using the SHARP code. We show that, due to such an inhomogeneity, the virulent beam-plasma instabilities in the intergalactic medium are not suppressed but their counterparts in the solar wind can be suppressed as evidenced by propagating type-III solar radio bursts.

Pair creation on the cosmic infrared background and subsequent inverse-Compton scattering on the CMB potentially reprocesses the TeV emission of blazars into faint GeV halos with structures sensitive to intergalactic magnetic fields (IGMF). Previous work has shown that these halos are then highly-anisotropic and extended. If the coherence length of the IGMF is greater than the inverse-Compton cooling length of the pairs, then the orientation of the gamma-ray halo will be correlated with the direction of the magnetic field which is unknown and expected to change for each source. In order to constructively add each source we then use angular power spectra which are insensitive to the jet orientation. By looking at known GeV blazars detected by Fermi, we exclude the existence of an IGMF with coherence lengths $>100$ Mpc at greater than $3.9\sigma$ with current-day strengths in the range $10^{-16}$ to $10^{-15}$ G, and at 2$\sigma$ from $10^{-17}$ to $10^{-14}$ G. This provides a direct measurement of the non-existence of gamma-ray halos, providing an important check on previous results.

Recent observations with the MeerKAT radio telescope reveal a unique population of faint non-thermal filaments pervading the central molecular zone (CMZ). Some of those filaments are organized into groups of almost parallel filaments, seemingly sorted by their length, so that their morphology resembles a harp with radio emitting "strings". We argue that the synchrotron emitting GeV electrons of these radio harps have been consecutively injected by the same source (a massive star or pulsar) into spatially intermittent magnetic fiber bundles within a magnetic flux tube or via time-dependent injection events. After escaping from this source, the propagation of cosmic ray (CR) electrons inside a flux tube is governed by the theory of CR transport. We propose to use observations of radio harp filaments to gain insight into the specifics of CR propagation along magnetic fields of which there are two principle modes: CRs could either stream with self-excited magneto-hydrodynamical waves or diffuse along the magnetic field. To disentangle these possibilities, we conduct hydrodynamical simulations of either purely diffusing or streaming CR electrons and compare the resulting brightness distributions to the observed synchrotron profiles of the radio harps. We find compelling evidence that CR streaming is the dominant propagation mode for GeV CRs in one of the radio harps. Observations at higher angular resolution should detect more radio harps and may help to disentangle projection effects of the possibly three-dimensional flux-tube structure of the other radio harps.

The circumgalactic medium (CGM) is one of the frontiers of galaxy formation and intimately connected to the galaxy via accretion of gas on to the galaxy and gaseous outflows from the galaxy. Here we analyse the magnetic field in the CGM of the Milky Way-like galaxies simulated as part of the \textsc{Auriga} project that constitutes a set of high resolution cosmological magnetohydrodynamical zoom simulations. We show that at high redshift the CGM becomes magnetised via galactic outflows that transport magnetised gas from the disk into the halo. At this time the magnetisation of the CGM closely follows its metal enrichment. We then show that at low redshift an in-situ turbulent dynamo that operates on a timescale of Gigayears further amplifies the magnetic field in the CGM and sets its strength down $z=0$. The magnetic field strength reaches a typical value of $0.1\,\mu G$ at the virial radius at $z=0$ and becomes mostly uniform within the virial radius. Its Faraday rotation signal is in excellent agreement with recent observations. For most of its evolution the magnetic field in the CGM is an unordered small scale field. Only strong coherent outflows at low redshift are able to order the magnetic field in parts of the CGM that are directly displaced by these outflows.

Cosmic ray (CR) protons are an important component in many astrophysical systems. Processes like CR injection, cooling, adiabatic changes as well as active CR transport through the medium strongly modify the CR momentum distribution and have to be taken into account in hydrodynamical simulations. We present an efficient novel numerical scheme to accurately compute the evolution of the particle distribution function by solving the Fokker-Planck equation with a low number of spectral bins (10 - 20), which is required to include a full spectrum for every computational fluid element. The distribution function is represented by piecewise power laws and is not forced to be continuous, which enables an optimal representation of the spectrum. The Fokker-Planck equation is solved with a two-moment approach evolving the CR number and energy density. The low numerical diffusion of the scheme reduces the numerical errors by orders of magnitude in comparison to classical schemes with piecewise constant spectral representations. With this method not only the spectral evolution of CRs can be computed accurately in magnetohydrodynamic simulations but also their dynamical impact as well as CR ionisation. This allows for more accurate models for astrophysical plasmas, like the interstellar medium, and direct comparisons with observations.

We present a method for efficiently modelling Braginskii viscosity on an unstructured, moving mesh. Braginskii viscosity, i.e., anisotropic transport of momentum with respect to the direction of the magnetic field, is thought to be of prime importance for studies of the weakly collisional plasma that comprises the intracluster medium (ICM) of galaxy clusters. Here anisotropic transport of heat and momentum has been shown to have profound consequences for the stability properties of the ICM. Our new method for modelling Braginskii viscosity has been implemented in the moving mesh code Arepo. We present a number of examples that serve to test the implementation and illustrate the modified dynamics found when including Braginskii viscosity in simulations. These include (but are not limited to) damping of fast magneto-sonic waves, interruption of linearly polarized Alfv\'en waves by the firehose instability and the inhibition of the Kelvin-Helmholtz instability by Braginskii viscosity. An explicit update of Braginskii viscosity is associated with a severe time step constraint that scales with $(\Delta x)^2$ where $\Delta x$ is the grid size. In our implementation, this restrictive time step constraint is alleviated by employing 2nd order accurate Runge-Kutta-Legendre super-time-stepping. We envision including Braginskii viscosity in future large-scale simulations of Kelvin-Helmholtz unstable cold fronts in cluster mergers and AGN-generated bubbles in central cluster regions.

Cosmic ray (CR) electrons reveal key insights into the non-thermal physics of the interstellar medium, galaxies, galaxy clusters, and active galactic nuclei by means of their inverse Compton {\gamma}-ray emission and synchrotron emission in magnetic fields. While magnetohydrodynamical (MHD) simulations with CR protons capture their dynamical impact on these systems, only few computational studies include CR electron physics because of the short cooling time-scales and complex hysteresis effects, which require a numerically expensive, high-resolution spectral treatment. Since CR electrons produce important non-thermal observational signatures, such a spectral CR electron treatment is important to link MHD simulations to observations. We present an efficient post-processing code for Cosmic Ray Electron Spectra that are evolved in Time (CREST) on Lagrangian tracer particles. The CR electron spectra are very accurately evolved on comparably large MHD time steps owing to an innovative hybrid numerical-analytical scheme. CREST is coupled to the cosmological MHD code AREPO and treats all important aspects of spectral CR electron evolution such as adiabatic expansion and compression, Coulomb losses, radiative losses in form of inverse Compton, bremsstrahlung and synchrotron processes, diffusive shock acceleration and reacceleration, and Fermi-II reacceleration. After showing various code validations of idealized one-zone simulations, we study the coupling of CREST to MHD simulations. We demonstrate that the CR electron spectra are efficiently and accurately evolved in shock-tube and Sedov-Taylor blast wave simulations. This opens up the possibility to produce self-consistent synthetic observables of non-thermal emission processes in various astrophysical environments.

We perform a series of cosmological hydrodynamic simulations to investigate the effects of non-gravitational dark matter (DM) interactions on the intergalactic medium (IGM). In particular, we use the ETHOS framework (Cyr-Racine et al. 2016; Vogelsberger et al. 2016) to compare statistics of the Lyman-$\alpha$ forest in cold dark matter (CDM) with an alternative model in which the DM couples strongly with a relativistic species in the early universe. These models are characterised by a cutoff in the linear power spectrum, followed by a series of 'dark acoustic oscillations' (DAOs) on sub-dwarf scales. While the primordial cutoff delays the formation of the first galaxies, structure builds-up more rapidly in the interacting DM model compared to CDM. We show that although DAOs are quickly washed away in the non-linear clustering of DM at $z\lesssim10$, their signature can be imprinted prominently in the Lyman-$\alpha$ flux power spectrum at $z>5$. On scales larger than the cutoff ($k\sim0.08$ s/km for the specific model considered here), the relative difference to CDM is reminiscent of a warm dark matter (WDM) model with a similar initial cutoff; however, the redshift evolution on smaller scales is distinctly different. The appearance and disappearance of DAOs in the Lyman-$\alpha$ flux spectrum provides a powerful way to distinguish interacting DM models from WDM and, indeed, variations in the thermal history of the IGM that may also induce a small-scale cutoff.

The Kelvin-Helmholtz instability (KHI) is a standard test of hydrodynamic and magnetohydrodynamic (MHD) simulation codes and finds many applications in astrophysics. The classic linear theory considers a discontinuity in density and velocity at the interface of two fluids. However, for numerical simulations of the KHI such initial conditions do not yield converged results even at the linear stage of the instability. Instead, smooth profiles of velocity and density are required for convergence. This renders analytical theory to be only approximately valid and hinders quantitative comparisons between the classical theory and simulations. In this paper we derive a linear theory for the KHI with smooth profiles and illustrate code testing with the MHD code Athena. We provide the linear solution for the KHI with smooth initial conditions in three different limits: inviscid hydrodynamics, ideal MHD and Braginskii-MHD. These linear solutions are obtained numerically with the framework Psecas (Pseudo-Spectral Eigenvalue Calculator with an Automated Solver), which generates and solves numerical eigenvalue problems using an equation-parser and pseudo-spectral methods. The Athena simulations are carried out on a periodic, Cartesian domain which is useful for code testing purposes. Using Psecas and analytic theory, we outline the differences between this artificial numerical setup and the KHI on an infinite Cartesian domain and the KHI in cylindrical geometry. We discuss several astrophysical applications, such as cold flows in galaxy formation and cold fronts in galaxy cluster mergers. Psecas, and the linear solutions used for code testing, are publicly available and can be downloaded from the web.

Star formation in galaxies appears to be self-regulated by energetic feedback processes. Among the most promising agents of feedback are cosmic rays (CRs), the relativistic ion population of interstellar and intergalactic plasmas. In these environments, energetic CRs are virtually collisionless and interact via collective phenomena mediated by kinetic-scale plasma waves and large-scale magnetic fields. The enormous separation of kinetic and global astrophysical scales requires a hydrodynamic description. Here, we develop a new macroscopic theory for CR transport in the self-confinement picture, which includes CR diffusion and streaming. The interaction between CRs and electromagnetic fields of Alfv\'enic turbulence provides the main source of CR scattering, and causes CRs to stream along the magnetic field with the Alfv\'en velocity if resonant waves are sufficiently energetic. However, numerical simulations struggle to capture this effect with current transport formalisms and adopt regularization schemes to ensure numerical stability. We extent the theory by deriving an equation for the CR momentum density along the mean magnetic field and include a transport equation for the Alfv\'en-wave energy. We account for energy exchange of CRs and Alfv\'en waves via the gyroresonant instability and include other wave damping mechanisms. Using numerical simulations we demonstrate that our new theory enables stable, self-regulated CR transport. The theory is coupled to magneto-hydrodynamics, conserves the total energy and momentum, and correctly recovers previous macroscopic CR transport formalisms in the steady-state flux limit. Because it is free of tunable parameters, it holds the promise to provide predictable simulations of CR feedback in galaxy formation.

We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cell mass of $\sim10^5\mathrm{M}_{\odot}$ and minimum gas softening of $\sim180$~pc) within ETHOS to date -- plus a CDM counterpart -- to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM $\tau$ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep JWST surveys of strongly-lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.

Feedback by active galactic nuclei (AGNs) is essential for regulating the fast radiative cooling of low-entropy gas at the centers of galaxy clusters and for reducing star formation rates of central ellipticals. The details of self-regulation depend critically on the unknown contents of AGN-inflated bubbles. Observations of the Sunyaev-Zeldovich (SZ) signal of AGN bubbles provide us with the ability to directly measure the lobe electron pressure given a bubble morphology. Here we compute the SZ signal of jet-inflated bubbles in three-dimensional magnetohydrodynamical simulations of the galaxy cluster MS0735.6+7421 with the Arepo code, and compare our synthetic SZ results to inferences obtained with popular modelling approaches. We find that cutting out ellipsoidal bubbles from a double-beta pressure profile only matches the inner bubble edges in the simulations and fails to account for the emission of the shock-enhanced pressure cocoon outside the bubbles. This additional contribution significantly worsens the accuracy of the cut-out method for jets with small inclinations with respect to the line of sight. Also, the kinetic SZ effect of the bubbles, a previously neglected contribution, becomes relevant at these smaller inclinations due to entrainment and mixing of the intracluster medium with low-density jet material. Fortunately, the different signs of the kinetic SZ signal in opposite lobes allow modelling this effect. We present an approximate method to determine the jet inclination, which combines jet power and lifetime estimates, the stand-off distance between jet head and bow shock, and the kinetic SZ effect, thereby helping to correctly infer the bubble contents.

Galactic winds exhibit a multiphase structure that consists of hot-diffuse and cold-dense phases. Here we present high-resolution idealised simulations of the interaction of a hot supersonic wind with a cold cloud with the moving-mesh code arepo in setups with and without radiative cooling. We demonstrate that cooling causes clouds with sizes larger than the cooling length to fragment in two- and three-dimensional simulations (2D and 3D). We confirm earlier 2D simulations by McCourt et al. 2018 and highlight differences of the shattering processes of 3D clouds that are exposed to a hot wind. The fragmentation process is quantified with a friends-of-friends analysis of shattered cloudlets and density power spectra. Those show that radiative cooling causes the power spectral index to gradually increase when the initial cloud radius is larger than the cooling length and with increasing time until the cloud is fully dissolved in the hot wind. A resolution of around 1 pc is required to reveal the effect of cooling-induced fragmentation of a 100 pc outflowing cloud. Thus, state-of-the-art cosmological zoom simulations of the circumgalactic medium (CGM) fall short by orders of magnitudes from resolving this fragmentation process. This physics is, however, necessary to reliably model observed column densities and covering fractions of Lyman-$\alpha$ haloes, high-velocity clouds, and broad-line regions of active galactic nuclei.

An intergalactic magnetic field stronger than $3\times10^{-13}$~G would explain the lack of a bright, extended degree-scale, GeV-energy inverse Compton component in the gamma-ray spectra of TeV-blazars. A robustly predicted consequence of the presence of such a field is the existence of degree-scale GeV-energy gamma-ray halos -- gamma-ray bow ties -- about TeV-bright active galactic nuclei, corresponding to more than half of all radio galaxies. However, the emitting regions of these halos are confined to and aligned with the direction of the relativistic jets associated with gamma-ray sources. Based on the orientation of radio jets, we align and stack corresponding degree-scale gamma-ray images of isolated Fanaroff-Riley class I and II objects and exclude the existence of these halos at overwhelming confidence, limiting the intergalactic field strength to $<10^{-15}$~G for large-scale fields and progressively larger in the diffusive regime when the correlation length of the field becomes small in comparison to 1 Mpc. When combined with prior limits on the strength of the intergalactic magnetic field, this excludes a purely magnetic explanation for the absence of halos. Thus, it requires the existence of novel physical processes that preempt the creation of halos, e.g., the presence of beam-plasma instabilities in the intergalactic medium or a drastic cutoff of the very high energy spectrum of these sources.

Faraday rotation is one of the most widely used observables to infer the strength and configuration of the magnetic field in the ionised gas of the Milky Way and nearby spiral galaxies. Here we compute synthetic Faraday rotation maps at $z=0$ for a set of disk galaxies from the Auriga high-resolution cosmological simulations, for different observer positions within and outside the galaxy. We find that the strength of the Faraday rotation of our simulated galaxies for a hypothetic observer at the solar circle is broadly consistent with the Faraday rotation seen for the Milky Way. The same holds for an observer outside the galaxy and the observed signal of the nearby spiral galaxy M51. However, we also find that the structure and angular power spectra of the synthetic all-sky Faraday rotation maps vary strongly with azimuthal position along the solar circle. We argue that this variation is a result of the structure of the magnetic field of the galaxy that is dominated by an azimuthal magnetic field ordered scales of several kpc, but has radial and vertical magnetic field components that are only ordered on scales of 1-2 kpc. Because the magnetic field strength decreases exponentially with height above the disk, the Faraday rotation for an observer at the solar circle is dominated by the local environment. This represents a severe obstacle for attempts to reconstruct the global magnetic field of the Milky Way from Faraday rotation maps alone without including additional observables.

In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.)

Feedback processes by active galactic nuclei in the centres of galaxy clusters appear to prevent large-scale cooling flows and impede star formation. However, the detailed heating mechanism remains uncertain. One promising heating scenario invokes the dissipation of Alfv\'en waves that are generated by streaming cosmic rays (CRs). In order to study this idea, we use three-dimensional magneto-hydrodynamical simulations with the AREPO code that follow the evolution of jet-inflated bubbles that are filled with CRs in a turbulent cluster atmosphere. We find that a single injection event produces the CR distribution and heating rate required for a successful CR heating model. As a bubble rises buoyantly, cluster magnetic fields drape around the leading interface and are amplified to strengths that balance the ram pressure. Together with helical magnetic fields in the bubble, this initially confines the CRs and suppresses the formation of interface instabilities. But as the bubble continues to rise, bubble-scale eddies significantly amplify radial magnetic filaments in its wake and enable CR transport from the bubble to the cooling intracluster medium. By varying the jet parameters, we obtain a rich and diverse set of jet and bubble morphologies ranging from Fanaroff-Riley type I-like (FRI) to FRII-like jets. We identify jet energy as the leading order parameter (keeping the ambient density profiles fixed), whereas jet luminosity is primarily responsible for setting the Mach numbers of shocks around FRII-like sources. Our simulations also produce FRI-like jets that inflate bubbles without detectable shocks and show morphologies consistent with cluster observations.

Understanding how baryonic processes shape the intracluster medium (ICM) is of critical importance to the next generation of galaxy cluster surveys. However, most models of structure formation neglect potentially important physical processes, like anisotropic thermal conduction (ATC). In this letter, we explore the impact of ATC on the prevalence of cool-cores (CCs) using 12 pairs of magnetohydrodynamical galaxy cluster simulations, simulated using the IllustrisTNG model with and without ATC. Although the impact of ATC varies from cluster to cluster and with CC criterion, its inclusion produces a systematic shift to larger CC fractions at z = 0 for all CC criteria considered. Additionally, the inclusion of ATC yields a flatter CC fraction redshift evolution, easing the tension with the observed evolution. With ATC included, the energy required for the central black hole to achieve self-regulation is reduced and the gas fraction in the cluster core increases, resulting in larger CC fractions. ATC makes the ICM unstable to perturbations and the increased efficiency of AGN feedback suggests that its inclusion results in a greater level of mixing in the ICM. Therefore, ATC is potentially an important physical process in reproducing the thermal structure of the ICM.

Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh code {\sc arepo} for different magnetic field morphologies, ranging from homogeneous to turbulent configurations. A Sedov-Taylor explosion in a homogeneous magnetic field generates an oblate ellipsoidal shock surface due to the slower propagating blast wave in the direction of the magnetic field. This is because of the efficient cosmic ray (CR) production in the quasi-parallel polar cap regions, which softens the equation of state and increases the compressibility of the post-shock gas. We find that the solution remains self-similar because the ellipticity of the propagating blast wave stays constant in time. This enables us to derive an effective ratio of specific heats for a composite of thermal gas and CRs as a function of the maximum acceleration efficiency. We finally discuss the behavior of supernova remnants expanding into a turbulent magnetic field with varying coherence lengths. For a maximum CR acceleration efficiency of about 15 per cent at quasi-parallel shocks (as suggested by kinetic plasma simulations), we find an average efficiency of about 5 per cent, independent of the assumed magnetic coherence length.

We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be suppressed by inhomogeneities in the intergalactic medium as was recently claimed in the literature. Using high resolution Particle-In-Cell simulations with the SHARP code, we show that the growth of the instability is local, i.e., regions with almost homogeneous background density will support the growth of the Langmuir waves, even when they are separated by strongly inhomogeneous regions, resulting in an overall slower growth of the instability. We also show that if the background density is continuously varying, the growth rate of the instability is lower; though in all cases, the system remains within the linear regime longer and the instability is not extinguished. In all cases, the beam loses approximately the same fraction of its initial kinetic energy in comparison to the uniform case at non-linear saturation. Thus, inhomogeneities in the intergalactic medium are unlikely to suppress the growth of blazar-driven beam-plasma instabilities.

In principle, the angular anisotropy in the extragalactic gamma-ray background (EGRB) places severe constraints upon putative populations of unresolved gamma-ray point sources. Existing estimates of the EGRB anisotropy have been constructed by 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.

Galactic winds regulate star formation in disk galaxies and help to enrich the circum-galactic medium. They are therefore crucial for galaxy formation, but their driving mechanism is still poorly understood. Recent studies have demonstrated that cosmic rays (CRs) can drive outflows if active CR transport is taken into account. Using hydrodynamical simulations of isolated galaxies with virial masses between $10^{10}$ and $10^{13}\mathrm{~M_\odot}$, we study how the properties of CR-driven winds depend on halo mass. CRs are treated in a two-fluid approximation and their transport is modelled through isotropic or anisotropic diffusion. We find that CRs are only able to drive mass-loaded winds beyond the virial radius in haloes with masses below $10^{12}\mathrm{~M_\odot}$. For our lowest examined halo mass, the wind is roughly spherical and has velocities of $\sim20\mathrm{~km\;s^{-1}}$. With increasing halo mass, the wind becomes biconical and can reach ten times higher velocities. The mass loading factor drops rapidly with virial mass, a dependence that approximately follows a power-law with a slope between $-1$ and $-2$. This scaling is slightly steeper than observational inferences, and also steeper than commonly used prescriptions for wind feedback in cosmological simulations. The slope is quite robust to variations of the CR injection efficiency or the CR diffusion coefficient. In contrast to the mass loading, the energy loading shows no significant dependence on halo mass. While these scalings are close to successful heuristic models of wind feedback, the CR-driven winds in our present models are not yet powerful enough to fully account for the required feedback strength.

Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wave numbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion.

Star forming galaxies emit GeV- and TeV-gamma rays that are thought to originate from hadronic interactions of cosmic-ray (CR) nuclei with the interstellar medium. To understand the emission, we have used the moving mesh code Arepo to perform magneto-hydrodynamical galaxy formation simulations with self-consistent CR physics. Our galaxy models exhibit a first burst of star formation that injects CRs at supernovae. Once CRs have sufficiently accumulated in our Milky-Way like galaxy, their buoyancy force overcomes the magnetic tension of the toroidal disk field. As field lines open up, they enable anisotropically diffusing CRs to escape into the halo and to accelerate a bubble-like, CR-dominated outflow. However, these bubbles are invisible in our simulated gamma-ray maps of hadronic pion-decay and secondary inverse-Compton emission because of low gas density in the outflows. By adopting a phenomenological relation between star formation rate (SFR) and far-infrared emission and assuming that gamma rays mainly originate from decaying pions, our simulated galaxies can reproduce the observed tight relation between far-infrared and gamma-ray emission, independent of whether we account for anisotropic CR diffusion. This demonstrates that uncertainties in modeling active CR transport processes only play a minor role in predicting gamma-ray emission from galaxies. We find that in starbursts, most of the CR energy is "calorimetrically" lost to hadronic interactions. In contrast, the gamma-ray emission deviates from this calorimetric property at low SFRs due to adiabatic losses, which cannot be identified in traditional one-zone models.

We formulate an effective theory of structure formation (ETHOS) that enables cosmological structure formation to be computed in almost any microphysical model of dark matter physics. This framework maps the detailed microphysical theories of particle dark matter interactions into the physical effective parameters that shape the linear matter power spectrum and the self-interaction transfer cross section of non-relativistic dark matter. These are the input to structure formation simulations, which follow the evolution of the cosmological and galactic dark matter distributions. Models with similar effective parameters in ETHOS but with different dark particle physics would nevertheless result in similar dark matter distributions. We present a general method to map an ultraviolet complete or effective field theory of low energy dark matter physics into parameters that affect the linear matter power spectrum and carry out this mapping for several representative particle models. We further propose a simple but useful choice for characterizing the dark matter self-interaction transfer cross section that parametrizes self-scattering in structure formation simulations. Taken together, these effective parameters in ETHOS allow the classification of dark matter theories according to their structure formation properties rather than their intrinsic particle properties, paving the way for future simulations to span the space of viable dark matter physics relevant for structure formation.

The MeerKAT telescope will be one of the most sensitive radio arrays in the pre-SKA era. Here we discuss a low-frequency SZ-selected cluster survey with MeerKAT, the MeerKAT Extended Relics, Giant Halos, and Extragalactic Radio Sources (MERGHERS) survey. The primary goal of this survey is to detect faint signatures of diffuse cluster emission, specifically radio halos and relics. SZ-selected cluster samples offer a homogeneous, mass-limited set of targets out to higher redshift than X-ray samples. MeerKAT is sensitive enough to detect diffuse radio emission at the faint levels expected in low-mass and high-redshift clusters, thereby enabling radio halo and relic formation theories to be tested with a larger statistical sample over a significantly expanded phase space. Complementary multiwavelength follow-up observations will provide a more complete picture of any clusters found to host diffuse emission, thereby enhancing the scientific return of the MERGHERS survey.

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.

The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100\,\rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10\%$ of the turbulent energy density with a typical strength of $10-50\,\rm{\mu G}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.

The absence of large cooling flows in cool core clusters appears to require self-regulated energy feedback by active galactic nuclei (AGNs) but the exact heating mechanism has not yet been identified. Here, we analyse whether a combination of cosmic ray (CR) heating and thermal conduction can offset radiative cooling. To this end, we compile a large sample of 39 cool core clusters and determine steady state solutions of the hydrodynamic equations that are coupled to the CR energy equation. We find solutions that match the observed density and temperature profiles for all our clusters well. Radiative cooling is balanced by CR heating in the cluster centres and by thermal conduction on larger scales, thus demonstrating the relevance of both heating mechanisms. Our mass deposition rates vary by three orders of magnitude and are linearly correlated to the observed star formation rates. Clusters with large mass deposition rates show larger cooling radii and require a larger radial extent of the CR injection function. Interestingly, our sample shows a continuous sequence in cooling properties: clusters hosting radio mini halos are characterised by the largest cooling radii, star formation and mass deposition rates in our sample and thus signal the presence of a higher cooling activity. The steady state solutions support the structural differences between clusters hosting a radio mini halo and those that do not.

Self-regulated feedback by active galactic nuclei (AGNs) appears to be critical in balancing radiative cooling of the low-entropy gas at the centres of galaxy clusters and in regulating star formation in central galaxies. In a companion paper, we found steady-state solutions of the hydrodynamic equations that are coupled to the CR energy equation for a large cluster sample. In those solutions, radiative cooling in the central region is balanced by streaming CRs through the generation and dissipation of resonantly generated Alfv{\'e}n waves and by thermal conduction at large radii. Here we demonstrate that the predicted non-thermal emission resulting from hadronic CR interactions in the intra-cluster medium exceeds observational radio (and gamma-ray) data in a subsample of clusters that host radio mini halos (RMHs). In contrast, the predicted non-thermal emission is well below observational data in cooling galaxy clusters without RMHs. These are characterised by exceptionally large AGN radio fluxes, indicating high CR yields and associated CR heating rates. We suggest a self-regulation cycle of AGN feedback in which non-RMH clusters are heated by streaming CRs homogeneously throughout the central cooling region. We predict {\em radio micro halos} surrounding the AGNs of these CR-heated clusters in which the primary emission may predominate the hadronically generated emission. Once the CR population has streamed sufficiently far and lost enough energy, the cooling rate increases, which explains the increased star formation rates in clusters hosting RMHs. Those could be powered hadronically by CRs that have previously heated the cluster core.

Jets from supermassive black holes in the centres of galaxy clusters are a potential candidate for moderating gas cooling and subsequent star formation through depositing energy in the intra-cluster gas. In this work, we simulate the jet-intra-cluster medium interaction using the moving-mesh magnetohydrodynamics code Arepo. Our model injects supersonic, low density, collimated and magnetised outflows in cluster centres, which are then stopped by the surrounding gas, thermalise and inflate low-density cavities filled with cosmic-rays. We perform high-resolution, non-radiative simulations of the lobe creation, expansion and disruption, and find that its dynamical evolution is in qualitative agreement with simulations of idealised low-density cavities that are dominated by a large-scale Rayleigh-Taylor instability. The buoyant rising of the lobe does not create energetically significant small-scale chaotic motion in a volume-filling fashion, but rather a systematic upward motion in the wake of the lobe and a corresponding back-flow perpendicular to it. We find that, overall, 50 per cent of the injected energy ends up in material which is not part of the lobe, and about 25 per cent remains in the inner 100 kpc. We conclude that jet-inflated, buoyantly rising cavities drive systematic gas motions which play an important role in heating the central regions, while mixing of lobe material is sub-dominant. Encouragingly, the main mechanisms responsible for this energy deposition can be modelled already at resolutions within reach in future, high-resolution cosmological simulations of galaxy clusters.

Many-degree-scale gamma-ray halos are expected to surround extragalactic high-energy gamma ray sources. These arise from the inverse Compton emission of an intergalactic population of relativistic electron/positron pairs generated by the annihilation of >100 GeV gamma rays on the extragalactic background light. These are typically anisotropic due to the jetted structure from which they originate or the presence of intergalactic magnetic fields. Here we propose a novel method for detecting these inverse-Compton gamma-ray halos based upon this anisotropic structure. Specifically, we show that by stacking suitably defined angular power spectra instead of images it is possible to robustly detect gamma-ray halos with existing Fermi Large Area Telescope (LAT) observations for a broad class of intergalactic magnetic fields. Importantly, these are largely insensitive to systematic uncertainties within the LAT instrumental response or associated with contaminating astronomical sources.

Cosmic rays (CRs) have recently re-emerged as attractive candidates for mediating feedback in galaxies because of their long cooling timescales. They can have energy densities comparable to the thermal gas, but do not suffer catastrophic cooling losses. Recent simulations have shown that the momentum and energy deposited by CRs moving with respect to the ambient medium can drive galactic winds. However, simulations are hampered by our ignorance of the details of CR transport. Two key limits previously considered model CR transport as a purely diffusive process (with a constant diffusion coefficient) and as an advective streaming process. With a series of GADGET simulations, we compare and contrast the results of these different assumptions. In idealised three-dimensional galaxy formation models, we show that these two cases result in significant differences for the galactic wind mass loss rates and star formation suppression in dwarf galaxies with halo masses $M\approx10^{10}\,\textrm{M}_\odot$: diffusive CR transport results in more than ten times larger mass loss rates compared to CR streaming models. We demonstrate that this is largely due to the excitation of Alfv\'en waves during the CR streaming process that drains energy from the CR population to the thermal gas, which is subsequently radiated away. By contrast, CR diffusion conserves the CR energy in the absence of adiabatic changes and if CRs are efficiently scattered by Alfv\'en waves that are propagating up the CR gradient. Moreover, because pressure gradients are preserved by CR streaming, but not diffusion, the two can have a significantly different dynamical evolution regardless of this energy exchange. In particular, the constant diffusion coefficients usually assumed can lead to unphysically high CR fluxes.

Feedback from central supermassive blackholes is often invoked to explain the low star formation rates in massive galaxies at the centers of galaxy clusters. However, the detailed physics of the coupling of the injected feedback energy with the intracluster medium is still unclear. Using high-resolution magnetohydrodynamic cosmological simulations of galaxy cluster formation, we investigate the role of anisotropic thermal conduction in shaping the thermodynamic structure of clusters, and, in particular, in modifying the impact of black hole feedback. Stratified anisotropically conducting plasmas are formally always unstable, and thus more prone to mixing, an expectation borne out by our results. The increased mixing efficiently isotropizes the injected feedback energy which in turn significantly improves the coupling between the feedback energy and the intracluster medium. This facilitates an earlier disruption of the cool core, reduces the star formation rate by more than an order of magnitude, and results in earlier quenching despite an overall lower amount of feedback energy injected into the cluster core. With conduction, the metallicity gradients and dispersions are lowered, aligning them better with observational constraints. These results highlight the important role of thermal conduction in establishing and maintaining quiescence of massive galaxies.

We discuss new methods to integrate the cosmic ray (CR) evolution equations coupled to magneto-hydrodynamics (MHD) on an unstructured moving mesh, as realised in the massively parallel AREPO code for cosmological simulations. We account for diffusive shock acceleration of CRs at resolved shocks and at supernova remnants in the interstellar medium (ISM), and follow the advective CR transport within the magnetised plasma, as well as anisotropic diffusive transport of CRs along the local magnetic field. CR losses are included in terms of Coulomb and hadronic interactions with the thermal plasma. We demonstrate the accuracy of our formalism for CR acceleration at shocks through simulations of plane-parallel shock tubes that are compared to newly derived exact solutions of the Riemann shock tube problem with CR acceleration. We find that the increased compressibility of the post-shock plasma due to the produced CRs decreases the shock speed. However, CR acceleration at spherically expanding blast waves does not significantly break the self-similarity of the Sedov-Taylor solution; the resulting modifications can be approximated by a suitably adjusted, but constant adiabatic index. In first applications of the new CR formalism to simulations of isolated galaxies and cosmic structure formation, we find that CRs add an important pressure component to the ISM that increases the vertical scale height of disk galaxies, and thus reduces the star formation rate. Strong external structure formation shocks inject CRs into the gas, but the relative pressure of this component decreases towards halo centres as adiabatic compression favours the thermal over the CR pressure.

About 1/3 of X-ray-luminous clusters show smooth, Mpc-scale radio emission, known as giant radio haloes. One promising model for radio haloes is Fermi-II acceleration of seed relativistic electrons by compressible turbulence. The origin of these seed electrons has never been fully explored. Here, we integrate the Fokker-Planck equation of the cosmic ray (CR) electron and proton distributions when post-processing cosmological simulations of cluster formation, and confront them with radio surface brightness and spectral data of Coma. For standard assumptions, structure formation shocks lead to a seed electron population which produces too centrally concentrated radio emission. Matching observations requires modifying properties of the CR population (rapid streaming; enhanced CR electron acceleration at shocks) or turbulence (increasing turbulent-to-thermal energy density with radius), but at the expense of fine-tuning. In a parameter study, we find that radio properties are exponentially sensitive to the amplitude of turbulence, which is inconsistent with small scatter in scaling relations. This sensitivity is removed if we relate the acceleration time to the turbulent dissipation time. In this case, turbulence above a threshold value provides a fixed amount of amplification; observations could thus potentially constrain the unknown CR seed population. To obtain sufficient acceleration, the turbulent magneto-hydrodynamics cascade has to terminate by transit time damping on CRs, i.e., thermal particles must be scattered by plasma instabilities. Understanding the small scatter in radio halo scaling relations may provide a rich source of insight on plasma processes in clusters.

The annihilation of TeV photons from extragalactic TeV sources and the extragalactic background light produces ultrarelativistic $e^{\pm}$ beams, which are subject to powerful plasma instabilities that sap their kinetic energy. Here we study the linear phase of the plasma instabilities that these pair beams drive. To this end, we calculate the linear growth rate of the beam plasma and oblique instability in the electrostatic approximation in both the reactive and kinetic regimes, assuming a Maxwell-J{\"u}ttner distribution for the pair beam. We reproduce the well-known reactive and kinetic growth rates for both the beam plasma and oblique mode. We demonstrate for the oblique instability that there is a broad spectrum of unstable modes that grow at the maximum rate for a wide range of beam temperatures and wave vector orientations relative to the beam. We also delineate the conditions for applicability for the reactive and kinetic regimes and find that the beam plasma mode transitions to the reactive regime at a lower Lorentz factor than the oblique mode due to a combination of their different scalings and the anisotropy of the velocity dispersions. Applying these results to the ultrarelativistic $e^{\pm}$ beams from TeV blazars, we confirm that these beams are unstable to both, the kinetic oblique mode and the reactive beam-plasma mode. These results are important in understanding how powerful plasma instabilities may sap the energy of the ultrarelativistic $e^{\pm}$ beams as they propagate through intergalactic space.

Extended inverse Compton halos are generally anticipated around extragalactic sources of gamma rays with energies above 100 GeV. These result from inverse Compton scattered cosmic microwave background photons by a population of high-energy electron/positron pairs produced by the annihilation of the high-energy gamma rays on the infrared background. Despite the observed attenuation of the high-energy gamma rays, the halo emission has yet to be directly detected. Here, we demonstrate that in most cases these halos are expected to be highly anisotropic, distributing the up-scattered gamma rays along axes defined either by the radio jets of the sources or oriented perpendicular to a global magnetic field. We present a pedagogical derivation of the angular structure in the inverse Compton halo and provide an analytic formalism that facilitates the generation of mock images. We discuss exploiting this fact for the purpose of detecting gamma-ray halos in a set of companion papers.

We study the formation of galactic outflows from supernova explosions (SNe) with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SNe placement and energy feedback, including cosmic rays (CR), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium.

Hydrodynamical shocks are a manifestation of the non-linearity of the Euler equations and play a fundamental role in cosmological gas dynamics. In this work, we identify and analyse shocks in the Illustris simulation, and contrast the results with those of non-radiative runs. We show that simulations with more comprehensive physical models of galaxy formation pose new challenges for shock finding algorithms due to radiative cooling and star-forming processes, prompting us to develop a number of methodology improvements. We find in Illustris a total shock surface area which is about 1.4 times larger at the present epoch compared to non-radiative runs, and an energy dissipation rate at shocks which is higher by a factor of around 7. Remarkably, shocks with Mach numbers above and below $\mathcal{M}\approx10$ contribute about equally to the total dissipation across cosmic time. This is in sharp contrast to non-radiative simulations, and we demonstrate that a large part of the difference arises due to strong black hole radio-mode feedback in Illustris. We also provide an overview of the large diversity of shock morphologies, which includes complex networks of halo-internal shocks, shocks on to cosmic sheets, feedback shocks due to black holes and galactic winds, as well as ubiquitous accretion shocks. In high redshift systems more massive than $10^{12}\,\mathrm{M}_\odot$ we discover the existence of a double accretion shock pattern in haloes. They are created when gas streams along filaments without being shocked at the outer accretion shock, but then forms a second, roughly spherical accretion shock further inside.

The physics of cosmic rays (CR) is a promising candidate for explaining the driving of galactic winds and outflows. Recent galaxy formation simulations have demonstrated the need for active CR transport either in the form of diffusion or streaming to successfully launch winds in galaxies. However, due to computational limitations, most previous simulations have modeled CR transport isotropically. Here, we discuss high resolution simulations of isolated disk galaxies in a $10^{11}\rm{M_\odot}$ halo with the moving mesh code {\sc Arepo} that include injection of CRs from supernovae, advective transport, CR cooling, and CR transport through isotropic or anisotropic diffusion. We show that either mode of diffusion leads to the formation of strong bipolar outflows. However, they develop significantly later in the simulation with anisotropic diffusion compared to the simulation with isotropic diffusion. Moreover, we find that isotropic diffusion allows most of the CRs to quickly diffuse out of the disk, while in the simulation with anisotropic diffusion, most CRs remain in the disk once the magnetic field becomes dominated by its azimuthal component, which occurs after $\sim 300\,{\rm Myrs}$. This has important consequences for the gas dynamics in the disk. In particular, we show that isotropic diffusion strongly suppresses the amplification of the magnetic field in the disk compared to anisotropic or no diffusion models. We therefore conclude that reliable simulations which include CR transport inevitably need to account for anisotropic diffusion.

In the interstellar medium of galaxies and the intracluster gas of galaxy clusters, the charged particles making up cosmic rays are moving almost exclusively along (but not across) magnetic field lines. The resulting anisotropic transport of cosmic rays in the form of diffusion or streaming not only affects the gas dynamics but also rearranges the magnetic fields themselves. The coupled dynamics of magnetic fields and cosmic rays can thus impact the formation and evolution of galaxies and the thermal evolution of galaxy clusters in critical ways. Numerically studying these effects requires solvers for anisotropic diffusion that are accurate, efficient, and robust, requirements that have proven difficult to satisfy in practice. Here, we present an anisotropic diffusion solver on an unstructured moving mesh that is conservative, does not violate the entropy condition, allows for semi-implicit time integration with individual timesteps, and only requires solving a single linear system of equations per timestep. We apply our new scheme to a large number of test problems and show that it works as well or better than previous implementations. Finally, we demonstrate for a numerically demanding simulation of the formation of an isolated disk galaxy that our local time-stepping scheme reproduces the results obtained with global time-stepping at a fraction of the computational cost.

We present the first simulations within an effective theory of structure formation (ETHOS), which includes the effect of interactions between dark matter and dark radiation on the linear initial power spectrum and dark matter self-interactions during non-linear structure formation. We simulate a Milky Way-like halo in four different dark matter models and the cold dark matter case. Our highest resolution simulation has a particle mass of $2.8\times 10^4\,{\rm M}_\odot$ and a softening length of $72.4\,{\rm pc}$. We demonstrate that all alternative models have only a negligible impact on large scale structure formation. On galactic scales, however, the models significantly affect the structure and abundance of subhaloes due to the combined effects of small scale primordial damping in the power spectrum and late time self-interactions. We derive an analytic mapping from the primordial damping scale in the power spectrum to the cutoff scale in the halo mass function and the kinetic decoupling temperature. We demonstrate that certain models within this extended effective framework that can alleviate the too-big-to-fail and missing satellite problems simultaneously, and possibly the core-cusp problem. The primordial power spectrum cutoff of our models naturally creates a diversity in the circular velocity profiles, which is larger than that found for cold dark matter simulations. We show that the parameter space of models can be constrained by contrasting model predictions to astrophysical observations. For example, some models may be challenged by the missing satellite problem if baryonic processes were to be included and even over-solve the too-big-to-fail problem; thus ruling them out.

Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse gamma-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E>100 GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6$^\circ$ from the centre, have been detected as point-like VHE gamma-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is well described by a power law with a photon index of $3.6\pm0.2_{stat}\pm0.2_{syst}$ between 90 GeV and 1.2 TeV. We do not detect any diffuse gamma-ray emission from the cluster and set stringent constraints on its CR population. In order to bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes $\alpha$. For $\alpha=2.2$, the CR-to-thermal pressure within the cluster virial radius is constrained to be below 1-2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to <20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, $B_0$, that depend on the CR distribution. For $\alpha=2.2$, $B_0\gtrsim5-8 \mu$G, which is below the 25 $\mu$G inferred from Faraday rotation measurements, whereas, for $\alpha\lesssim2.1$, the hadronic interpretation of the diffuse radio emission is in contrast with our gamma-ray flux upper limits independently of the magnetic field strength.

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.

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).

TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron-positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature-density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts on the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large scale simulations to determine the impact of inhomogeneous heating. We show that, because of blazar clustering, blazar heating is inhomogeneous for z>= 2. At high redshift, the temperature-density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature-density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large scale simulations and provides a cost-effective alternative to subgrid models.

Galaxy clusters are the largest and most massive gravitationally bound structures known in the Universe. Cosmic-Ray (CR) hadrons accelerated at structure formation shocks and injected by galaxies, are confined in galaxy clusters where they accumulate for cosmological times. The presence of diffuse synchrotron radio emission in several clusters proves the existence of high-energy electrons, and magnetic fields. However, a direct proof of CR proton acceleration is missing. The presence of CR protons can be probe through the diffuse gamma-ray emission induced by their hadronic interaction with the Intra-Cluster Medium (ICM). The Perseus cluster, a nearby cool-core cluster, has been identified to be among the best candidates to detect such emission. We present here the results of a very deep observation of the Perseus cluster with the MAGIC telescopes, accumulating about 250 hours of data from 2009 to 2014. No evidence of large-scale very-high-energy gamma-ray emission from CR-ICM interactions has been detected. The derived flux upper limits in the TeV regime allow us to put stringent constraints on the physics of cluster CRs, in particular on the CR-to-thermal pressure, the CR acceleration efficiency at formation shocks and the magnetic field of the central cluster region.

About 1/3 of X-ray-luminous clusters show smooth, unpolarized radio emission on ~Mpc scales, known as giant radio halos. One promising model for radio halos is Fermi-II acceleration of seed relativistic electrons by turbulence of the intracluster medium (ICM); Coulomb losses prohibit acceleration from the thermal pool. However, the origin of seed electrons has never been fully explored. Here, we integrate the Fokker-Planck equation of the cosmic ray (CR) electron and proton distributions in a cosmological simulations of cluster formation. For standard assumptions, structure formation shocks lead to a seed electron population which produces too centrally concentrated radio emission. Instead, we present three realistic scenarios that each can reproduce the spatially flat radio emission observed in the Coma cluster: (1) the ratio of injected turbulent energy density to thermal energy density increase significantly with radius, as seen in cosmological simulations. This generates a flat radio profile even if the seed population of CRs is steep with radius. (2) Self-confinement of energetic CR protons can be inefficient, and CR protons may stream at the Alfven speed to the cluster outskirts when the ICM is relatively quiescent. A spatially flat CR proton distribution develops and produces the required population of secondary seed electrons. (3) The CR proton to electron acceleration efficiency K_ep ~ 0.1 is assumed to be larger than in our Galaxy (K_ep ~ 0.01), due to the magnetic geometry at the shock. The resulting primary electron population dominates. Due to their weaker density dependence compared to secondary electrons, these primaries can also reproduce radio observations. These competing non-trivial solutions provide incisive probes of non thermal processes in the high-beta ICM.

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.

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.

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 $\gamma$ 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 $\gamma$-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\sigma$ 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 $\gamma$-ray flux from individual clusters in our sample.

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.

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]

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.

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 (\sigma_8 and \Omega_M) and scales as \sigma_8^7.4 and \Omega_M^1.9 around the multipole (ell) ~ 1000. We constrain \sigma_8 and \Omega_M from the cross-power spectrum to be \sigma_8 (\Omega_M/0.30)^0.26 = 0.8 +/- 0.02. Since this cross spectrum produces a tight constraint in the \sigma_8 and \Omega_M plane the errors on a \sigma_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 \sigma_8 and \Omega_M to constrain dark energy or the sum of the neutrino masses.

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.

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.

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.

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.

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.

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.

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.

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, <\sigma v>. (abr.)

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.

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)

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)

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.

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.

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 \Gamma=-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.

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]

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

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.

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.

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%.

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.

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.

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.

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.

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.

[...] 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.

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.

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.

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.

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.

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.

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 $\sigma_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.

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.

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.

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.

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.

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.

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]

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]

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]

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.

Cosmic rays (CRs) play a decisive role within our own Galaxy. They provide partial pressure support against gravity, they trace past energetic events such as supernovae, and they reveal the underlying structure of the baryonic matter distribution through their interactions. To study the impact of CRs on galaxy and cosmic structure formation and evolution, we develop an approximative framework for treating dynamical and radiative effects of CRs in cosmological simulations. Our guiding principle is to try to find a balance between capturing as many physical properties of CR populations as possible while at the same time requiring as little extra computational resources as possible. We approximate the CR spectrum of each fluid element by a single power-law, with spatially and temporally varying normalisation, low-energy cut-off, and spectral index. Principles of conservation of particle number, energy, and pressure are then used to derive evolution equations for the basic variables describing the CR spectrum, both due to adiabatic and non-adiabatic processes. The processes considered include compression and rarefaction, CR injection via shocks in supernova remnants, injection in structure formation shock waves, in-situ re-acceleration of CRs, CR spatial diffusion, CR energy losses due to Coulomb interactions, ionisation losses, Bremsstrahlung losses, and, finally, hadronic interactions with the background gas, including the associated gamma-ray and radio emission due to subsequent pion decay. We show that the formalism reproduces CR energy densities, pressure, and cooling rates with an accuracy of ~10% in steady state conditions where CR injection balances cooling. Our framework is therefore well suited to be included into numerical simulation schemes of galaxy and structure formation. (abridged)

We performed high-resolution simulations of a sample of 14 galaxy clusters that span a mass range from 5 x 10^13 M_solar/h to 2 x 10^15 M_solar/h to study the effects of cosmic rays (CRs) on thermal cluster observables such as X-ray emission and the Sunyaev-Zel'dovich effect. We analyse the CR effects on the intra-cluster medium while simultaneously taking into account the cluster's dynamical state as well as the mass of the cluster. The modelling of the cosmic ray physics includes adiabatic CR transport processes, injection by supernovae and cosmological structure formation shocks, as well as CR thermalization by Coulomb interaction and catastrophic losses by hadronic interactions. While the relative pressure contained in CRs within the virial radius is of the order of 2 per cent in our non-radiative simulations, their contribution rises to 32 per cent in our simulations with dissipative gas physics including radiative cooling, star formation, and supernova feedback. Interestingly, in the radiative simulations the relative CR pressure reaches high values of the order of equipartition with the thermal gas in each cluster galaxy due to the fast thermal cooling of gas which diminishes the thermal pressure support relative to that in CRs. This also leads to a lower effective adiabatic index of the composite gas that increases the compressibility of the intra-cluster medium. This effect slightly increases the central density, thermal pressure and the gas fraction. While the X-ray luminosity in low mass cool core clusters is boosted by up to 40 per cent, the integrated Sunyaev-Zel'dovich effect appears to be remarkably robust and the total flux decrement only slightly reduced by typically 2 per cent. The resolved Sunyaev-Zel'dovich maps, however, show a larger variation with an increased central flux decrement. [abridged]

There is growing consensus that feedback from AGN is the main mechanism responsible for stopping cooling flows in clusters of galaxies. AGN are known to inflate buoyant bubbles that supply mechanical power to the intracluster gas (ICM). High Reynolds number hydrodynamical simulations show that such bubbles get entirely disrupted within 100 Myr, as they rise in cluster atmospheres, which is contrary to observations. This artificial mixing has consequences for models trying to quantify the amount of heating and star formation in cool core clusters of galaxies. It has been suggested that magnetic fields can stabilize bubbles against disruption. We perform MHD simulations of fossil bubbles in the presence of tangled magnetic fields using the high order PENCIL code. We focus on the physically-motivated case where thermal pressure dominates over magnetic pressure and consider randomly oriented fields with and without maximum helicity and a case where large scale external fields drape the bubble.We find that helicity has some stabilizing effect. However, unless the coherence length of magnetic fields exceeds the bubble size, the bubbles are quickly shredded. As observations of Hydra A suggest that lengthscale of magnetic fields may be smaller then typical bubble size, this may suggest that other mechanisms, such as viscosity, may be responsible for stabilizing the bubbles. However, since Faraday rotation observations of radio lobes do not constrain large scale ICM fields well if they are aligned with the bubble surface, the draping case may be a viable alternative solution to the problem. A generic feature found in our simulations is the formation of magnetic wakes where fields are ordered and amplified. We suggest that this effect could prevent evaporation by thermal conduction of cold Halpha filaments observed in the Perseus cluster.

Cosmological shock waves during structure formation not only play a decisive role for the thermalization of gas in virializing structures but also for the acceleration of relativistic cosmic rays (CRs) through diffusive shock acceleration. We discuss a novel numerical treatment of the physics of cosmic rays in combination with a formalism for identifying and measuring the shock strength on-the-fly during a smoothed particle hydrodynamics simulation. In our methodology, the non-thermal CR population is treated self-consistently in order to assess its dynamical impact on the thermal gas as well as other implications on cosmological observables. Using this formalism, we study the history of the thermalization process in high-resolution hydrodynamic simulations of the Lambda cold dark matter model. Collapsed cosmological structures are surrounded by shocks with high Mach numbers up to 1000, but they play only a minor role in the energy balance of thermalization. However, this finding has important consequences for our understanding of the spatial distribution of CRs in the large-scale structure. In high resolution simulations of galaxy clusters, we find a low contribution of the averaged CR pressure, due to the small acceleration efficiency of lower Mach numbers of flow shocks inside halos and the softer adiabatic index of CRs. However, within cool core regions, the CR pressure reaches equipartition with the thermal pressure leading there to a lower effective adiabatic index and thus to an enhanced compressibility of the central intracluster medium. This effect increases the central density and pressure of the cluster and thus the resulting X-ray emission and the central Sunyaev-Zel'dovich flux decrement. The integrated Sunyaev-Zel'dovich effect, however, is only slightly changed.

It is well known that cosmic rays (CRs) contribute significantly to the pressure of the interstellar medium in our own Galaxy, suggesting that they may play an important role in regulating star formation during the formation and evolution of galaxies. We will present a novel numerical treatment of the physics of CRs and its implementation in the parallel smoothed particle hydrodynamics (SPH) code GADGET-2. In our methodology, the non-thermal CR population is treated self-consistently in order to assess its dynamical impact on the thermal gas as well as other implications on cosmological observables. In simulations of galaxy formation, we find that CRs can significantly reduce the star formation efficiencies of small galaxies. This effect becomes progressively stronger towards low mass scales. In cosmological simulations of the formation of dwarf galaxies at high redshift, we find that the total mass-to-light ratio of small halos and the faint-end of the luminosity function are affected. In high resolution simulations of galaxy clusters, we find lower contributions of CR pressure, due to the smaller CR injection efficiencies at low Mach number flow shocks inside halos, and the softer adiabatic index of CRs, which disfavours them when a composite of thermal gas and CRs is adiabatically compressed. Within cool core regions, the CR pressure reaches equipartition with the thermal pressure leading to an enhanced compressibility of the central intra-cluster medium, an effect that increases the central density and pressure of the gas. While the X-ray luminosity in low mass cool core clusters is boosted, the integrated Sunyaev-Zel'dovich effect is only slightly changed. The resolved Sunyaev-Zel'dovich maps, however, show a larger variation with an increased central flux decrement.

We develop a formalism for the identification and accurate estimation of the strength of structure formation shocks during cosmological smoothed particle hydrodynamics simulations. Shocks not only play a decisive role for the thermalization of gas in virialising structures but also for the acceleration of relativistic cosmic rays (CRs) through diffusive shock acceleration. Our formalism is applicable both to ordinary non-relativistic thermal gas, and to plasmas composed of CRs and thermal gas. To this end, we derive an analytical solution to the one-dimensional Riemann shock tube problem for a composite plasma of CRs and thermal gas. We apply our methods to study the properties of structure formation shocks in high-resolution hydrodynamic simulations of the LCDM model. We find that most of the energy is dissipated in weak internal shocks with Mach numbers M~2 which are predominantly central flow shocks or merger shock waves traversing halo centres. Collapsed cosmological structures are surrounded by external shocks with much higher Mach numbers up to M~1000, but they play only a minor role in the energy balance of thermalization. We show that after the epoch of cosmic reionisation the Mach number distribution is significantly modified by an efficient suppression of strong external shock waves due to the associated increase of the sound speed of the diffuse gas. Invoking a model for CR acceleration in shock waves, we find that the average strength of shock waves responsible for CR energy injection is higher than that for shocks that dominate the thermalization of the gas. When combined with radiative dissipation and star formation, our formalism can also be used to study CR injection by supernova shocks, or to construct models for shock-induced star formation in the interstellar medium. (abridged)

It is well known that cosmic rays (CRs) contribute significantly to the pressure of the interstellar medium in our own Galaxy, suggesting that they may play an important role in regulating star formation during the formation and evolution of galaxies. We here discuss a novel numerical treatment of the physics of CRs and its implementation in the parallel smoothed particle hydrodynamics code GADGET-2. In our methodology, the non-thermal CR population of each gaseous fluid element is approximated by a simple power law spectrum in particle momentum, characterized by an amplitude, a cut-off, and a fixed slope. Adiabatic compression, and a number of physical source and sink terms are modelled which modify the CR pressure of each particle. The most important sources considered are injection by supernovae and diffusive shock acceleration, while the primary sinks are thermalization by Coulomb interactions, and catastrophic losses by hadronic interactions. We also include diffusion of CRs. Our scheme allows us to carry out the first cosmological structure formation simulations that self-consistently account for CR physics. In simulations of isolated galaxies, we find that CRs can significantly reduce the star formation efficiencies of small galaxies, with virial velocities below \~80 km/s, an effect that becomes progressively stronger towards low mass scales. In cosmological simulations at high redshift, the total mass-to-light ratio of small halos and the faint-end of the luminosity function are strongly affected. When CR acceleration in shocks is followed as well, up to ~40% of the energy dissipated at structure formation shocks can appear as CR pressure at z~3-6, but this fraction drops to ~10% at low redshifts when the shock distribution becomes increasingly dominated by lower Mach numbers. (abridged)

A brief overview about our knowledge on galaxy cluster magnetic fields is provided. Emphasize is given to the mutual dependence of our knowledge on relativistic particles in galaxy clusters and the magnetic field strength. Furthermore, we describe efforts to measure magnetic field strengths, characteristic length-scales, and power-spectra with reliable accuracy. An interpretation of these results in terms of non-helical dynamo theory is given. If this interpretation turns out to be correct, the understanding of cluster magnetic fields is directly connected to our understanding of intra-cluster turbulence.

There have been many speculations about the presence of cosmic ray protons (CRps) in galaxy clusters over the past two decades. However, no direct evidence such as the characteristic gamma-ray signature of decaying pions has been found so far. These pions would be a direct tracer of hadronic CRp interactions with the ambient thermal gas also yielding observable synchrotron and inverse Compton emission by additionally produced secondary electrons. The obvious question concerns the type of galaxy clusters most likely to yield a signal: Particularly suited sites should be cluster cooling cores due to their high gas and magnetic energy densities. We studied a nearby sample of clusters evincing cooling cores in order to place stringent limits on the cluster CRp population by using non-detections of EGRET. In this context, we examined the possibility of a hadronic origin of Coma-sized radio halos as well as radio mini-halos. Especially for mini-halos, strong clues are provided by the very plausible small amount of required CRp energy density and a matching radio profile. Introducing the hadronic minimum energy criterion, we show that the energetically favored CRp energy density is constrained to 2% +/- 1% of the thermal energy density in Perseus. We also studied the CRp population within the cooling core region of Virgo using the TeV gamma-ray detection of M 87 by HEGRA. Both the expected radial gamma-ray profile and the required amount of CRp support this hadronic scenario.

The Chandra X-ray Observatory is finding a large number of cavities in the X-ray emitting intra-cluster medium which often coincide with the lobes of the central radio galaxy. We propose high-resolution Sunyaev-Zel'dovich (SZ) observations in order to infer the still unknown dynamically dominating component of the radio plasma bubbles. This work calculates the thermal and relativistic SZ emission of different compositions of these plasma bubbles while simultaneously allowing for the cluster's kinetic SZ effect. As examples, we present simulations of an Atacama Large Millimeter Array (ALMA) observation and of a Green Bank Telescope (GBT) observation of the cores of the Perseus cluster and Abell 2052. We predict a 5 sigma detection of the southern radio bubble of Perseus in a few hours with the GBT and ALMA while assuming a relativistic electron population within the bubble. In Abell 2052, a similar detection would require a few ten hours with either telescope, the longer exposures mainly being the result of the higher redshift and the lower central temperature of this cluster. Future high-sensitivity multi-frequency SZ observations will be able to infer the energy spectrum of the dynamically dominating electron population in order to measure its temperature or spectral characteristics. This knowledge can yield indirect indications for an underlying radio jet model.

We wish to constrain the cosmic-ray proton (CRp) population in galaxy clusters. By hadronic interactions with the thermal gas of the intra-cluster medium (ICM), the CRp produce gamma-rays for which we develop an analytic formalism to deduce their spectral distribution. Assuming the CRp-to-thermal energy density ratio X_CRp and the CRp spectral index to be spatially constant, we derive an analytic relation between the gamma-ray and bolometric X-ray fluxes, F_gamma and F_X. Based on our relation, we compile a sample of suitable clusters which are promising candidates for future detection of gamma-rays resulting from hadronic CRp interactions. Comparing to EGRET upper limits, we constrain the CRp population in the cooling flow clusters Perseus and Virgo to X_CRp < 20%. Assuming a plausible value for the CRp diffusion coefficient kappa, we find the central CRp injection luminosity of M 87 to be limited to 10^43 erg s^-1 kappa/(10^29 cm^2 s^-1). The synchrotron emission from secondary electrons generated in CRp hadronic interactions allows even tighter limits to be placed on the CRp population using radio observations. We obtain excellent agreement between the observed and theoretical radio brightness profiles for Perseus, but not for Coma without a radially increasing CRp-to-thermal energy density profile. Since the CRp and magnetic energy densities necessary to reproduce the observed radio flux are very plausible, we propose synchrotron emission from secondary electrons as an attractive explanation of the radio mini-halos found in cooling flow clusters. This model can be tested with future sensitive gamma-ray observations of the accompanying pi0-decays. We identify Perseus (A 426), Virgo, Ophiuchus, and Coma (A 1656) as the most promising candidate clusters for such observations.

We investigate statistical distributions of differences in gravitational-lensing deflections between two light rays, the so-called lensing excursion angles. A probability distribution function of the lensing excursion angles, which plays a key role in estimates of lensing effects on angular clustering of objects (such as galaxies, QSOs and also the CMB temperature map), is known to consist of two components; a Gaussian core and an exponential tail. We use numerical gravitational-lensing experiments in a Lambda-CDM cosmology for quantifying these two components. We especially focus on the physical processes responsible for generating those two components. We develop a simple empirical model for the exponential tail which allows us to explore its origin. We find that the tail is generated by the coherent lensing scatter by massive halos with M>10^14h^-1Msun at z<1 and that its exponential shape arises due to the exponential cut-off of the halo mass function at that mass range. On scales larger than 1 arc minute, the tail does not have a practical influence on the lensing effects on the angular clustering. Our model predicts that the coherent scatter may have non-negligible effects on angular clustering at sub-arcminute scales.

All-sky thermal and kinetic Sunyaev-Zel'dovich (SZ) maps are presented for assessing how well the PLANCK-mission can find and characterise clusters of galaxies, especially in the presence of primary anisotropies of the cosmic microwave background (CMB) and various galactic and ecliptic foregrounds. The maps have been constructed from numerical simulations of structure formation in a standard LCDM cosmology and contain all clusters out to redshifts of z = 1.46 with masses exceeding 5e13 M_solar/h. By construction, the maps properly account for the evolution of cosmic structure, the halo-halo correlation function, the evolving mass function, halo substructure and adiabatic gas physics. The velocities in the kinetic map correspond to the actual density environment at the cluster positions. We characterise the SZ-cluster sample by measuring the distribution of angular sizes, the integrated thermal and kinetic Comptonisations, the source counts in the three relevant PLANCK-channels, and give the angular power-spectra of the SZ-sky. While our results are broadly consistent with simple estimates based on scaling relations and spherically symmetric cluster models, some significant differences are seen which may affect the number of cluster detectable by PLANCK.

The PLANCK mission is the most sensitive all-sky CMB experiment currently planned. The High Frequency Instrument (HFI) will be especially suited for observing clusters of galaxies by their thermal Sunyaev-Zel'dovich (SZ) effect. In order to assess PLANCKs SZ-capabilities in the presence of spurious signals, a simulation is presented that combines maps of the thermal and kinetic SZ-effects with a realisation of the cosmic microwave background (CMB), in addition to Galactic foregrounds (synchrotron emission, free-free emission, thermal emission from dust, CO-line radiation) as well as the sub-millimetric emission from celestial bodies of our Solar system. Additionally, observational issues such as the finite angular resolution and spatially non-uniform instrumental noise of PLANCKs sky maps are taken into account, yielding a set of all-sky flux maps, the auto-correlation and cross-correlation properties of which are examined in detail. In the second part of the paper, filtering schemes based on scale-adaptive and matched filtering are extended to spherical data sets, that enable the amplification of the weak SZ-signal in the presence of all contaminations stated above. The theory of scale-adaptive and matched filtering in the framework of spherical maps is developed, the resulting filter kernel shapes are discussed and their functionality is verified.

We wish to estimate magnetic field strengths of radio emitting galaxy clusters by minimising the non-thermal energy density contained in cosmic ray electrons (CRe), protons (CRp), and magnetic fields. The classical minimum energy estimate can be constructed independently of the origin of the radio synchrotron emitting CRe yielding thus an absolute minimum of the non-thermal energy density. Provided the observed synchrotron emission is generated by a CRe population originating from hadronic interactions of CRp with the ambient thermal gas of the intra-cluster medium, the parameter space of the classical scenario can be tightened by means of the hadronic minimum energy criterion. For both approaches, we derive the theoretically expected tolerance regions for the inferred minimum energy densities. Application to the radio halo of the Coma cluster and the radio mini-halo of the Perseus cluster yields equipartition between cosmic rays and magnetic fields within the expected tolerance regions. In the hadronic scenario, the inferred central magnetic field strength ranges from 2.4 muG (Coma) to 8.8 muG (Perseus), while the optimal CRp energy density is constrained to 2% +/- 1% of the thermal energy density (Perseus). We discuss the possibility of a hadronic origin of the Coma radio halo while current observations favour such a scenario for the Perseus radio mini-halo. Combining future expected detections of radio synchrotron, hard X-ray inverse Compton, and hadronically induced gamma-ray emission should allow an estimate of volume averaged cluster magnetic fields and provide information about their dynamical state.

A method for estimating redshifts of galaxy clusters based solely on resolved Sunyaev-Zel'dovich (SZ) images is proposed. Given a high resolution SZ cluster image (with FWHM of approx. 1 arcmin), the method indirectly measures its structure related parameters (amplitude, size, etc.) by fitting a model function to the higher order wavelet momenents of the cluster's SZ morphology. The applicability and accuracy of the wavelet method is assessed by applying it to maps of a set of clusters extracted from hydrodynamical simulations of cosmic structure formation. The parameters, derived by a fit to the spectrum of wavelet moments as a function of scale, are found to show a dependence on redshift z that is of the type x(z) = x_1 exp(-z/x_2) + x_3, where the monotony of this functional behaviour and the non-degeneracy of those parameters allow inversion and estimation of the redshift z. The average attainable accuracy in the z-estimation relative to 1+z is approx. 4-5% out to z = 1.2, which is comparable to photometric redshifts. For single-frequency SZ interferometers, where the ambient fluctuating CMB is the main noise source, the accuracy of the method drops slightly to <Delta z/(1+z)> = 6-7%.

Cooling flows are regions where the importance of non-thermal intra-cluster medium components such as magnetic fields and cosmic rays may be strongest within a galaxy cluster. They are also regions where such components are best detectable due to the high gas density which influences Faraday rotation measurements of magnetic fields and secondary particle production in hadronic interactions of cosmic ray nuclei with the ambient thermal gas. New estimates of magnetic fields in cooling flow and non-cooling flow clusters are presented, which are based on a newly developed Fourier analysis of extended Faraday rotation maps. We further present new constraints on the cluster cosmic ray proton population using radio and gamma-ray observations measurements of cluster cooling flows, which are especially suited for this purpose due to their high gas and magnetic energy densities. We argue that radio synchrotron emission of cosmic ray electrons generated hadronically by cosmic ray protons is a very plausible explanation for the radio mini-halos observed in some cooling flows.

We examine the cosmic ray proton (CRp) population within the giant elliptical galaxy M 87 using the TeV gamma-ray detection of the HEGRA collaboration. In our scenario, the gamma-rays are produced by decaying pions which result from hadronic CRp interactions with thermal gas of the interstellar medium of M 87. By comparing the gamma-ray emission to upper limits of EGRET, we constrain the spectral index of the CRp population to alpha(GeV,TeV) < 2.275 within our scenario. Both the expected radial gamma-ray profile and the required amount of CRp support this hadronic scenario. The accompanying radio mini-halo of hadronically originating cosmic ray electrons is outshone by the synchrotron emission of the relativistic jet of M 87 by one order of magnitude. According to our predictions, the future GLAST mission should allow us to test this hadronic scenario.