Showing 51–100 of 168 results for author: Pfrommer, C

ETHOS -- An effective parametrization and classification for structure formation: the non-linear regime at $z\gtrsim5$

Authors: Sebastian Bohr, Jesús Zavala, Francis-Yan Cyr-Racine, Mark Vogelsberger, Torsten Bringmann, Christoph Pfrommer

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

Submitted 7 September, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: 19 pages, 17 figures; Updated version accepted for publication in MNRAS

Journal ref: MNRAS 498, 3403 (2020)

The growth of the longitudinal beam-plasma instability in the presence of an inhomogeneous background

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

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

Submitted 5 March, 2020; originally announced March 2020.

Comments: 23 pages, 5 figures, accepted for publication in the Journal of Plasma Physics; comments are welcome

Probing Cosmic Ray Transport with Radio Synchrotron Harps in the Galactic Center

Authors: Timon Thomas, Christoph Pfrommer, Torsten A. Enßlin

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

Submitted 18 December, 2019; originally announced December 2019.

Comments: 7 pages, 5 figures, submitted to ApJL. Comments welcome!

Magnetising the circumgalactic medium of disk galaxies

Authors: Ruediger Pakmor, Freeke van de Voort, Rebekka Bieri, Facundo A. Gomez, Robert J. J. Grand, Thomas Guillet, Federico Marinacci, Christoph Pfrommer, Christine M. Simpson, Volker Springel

Abstract: 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 s… ▽ More 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 before $z=1$ 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 saturates before $z=0$. The magnetic field strength reaches a typical value of $0.1\,μ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. △ Less

Submitted 7 September, 2020; v1 submitted 25 November, 2019; originally announced November 2019.

Comments: 14 pages, 15 figures, accepted by MNRAS, updated Fig. 15 in proofs with errorbars

Braginskii viscosity on an unstructured, moving mesh accelerated with super-time-stepping

Authors: Thomas Berlok, Ruediger Pakmor, Christoph Pfrommer

Abstract: 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… ▽ More 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én 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 $(Δx)^2$ where $Δ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. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: 21 pages, 11 figures, accepted for publication in MNRAS

The effects of cosmic rays on the formation of Milky Way-mass galaxies in a cosmological context

Authors: Tobias Buck, Christoph Pfrommer, Rüdiger Pakmor, Robert J. J. Grand, Volker Springel

Abstract: 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én-wave cooling affect the galactic disc and the circum-galactic medium (CGM). Global properties such as… ▽ More 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én-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én-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én-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. △ Less

Submitted 2 July, 2020; v1 submitted 31 October, 2019; originally announced November 2019.

Comments: accepted by MNRAS, 26 pages (5 appendix), 20 figures (5 appendix)

Spectrally resolved cosmic ray hydrodynamics -- I. Spectral scheme

Authors: Philipp Girichidis, Christoph Pfrommer, Michal Hanasz, Thorsten Naab

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

Submitted 19 November, 2019; v1 submitted 27 September, 2019; originally announced September 2019.

Comments: 16 pages, MNRAS accepted, corrected typo in eq. (63)

The impact of magnetic fields on cold streams feeding galaxies

Authors: Thomas Berlok, Christoph Pfrommer

Abstract: 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 r… ▽ More 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 μ$G. △ Less

Submitted 2 June, 2020; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 19 pages, 15 figures, accepted for publication in MNRAS, animations of Figures 4 and 5 at http://tberlok.dk/movies/coldstream.html

Evolution of cosmic ray electron spectra in magnetohydrodynamical simulations

Authors: Georg Winner, Christoph Pfrommer, Philipp Girichidis, Rüdiger Pakmor

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

Submitted 23 September, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: 19 pages, 14 figures, 11 pages, 9 figures, published in MNRAS

Journal ref: MNRAS, 488, 2235W (2019)

On the Kelvin-Helmholtz instability with smooth initial conditions -- Linear theory and simulations

Authors: Thomas Berlok, Christoph Pfrommer

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

Submitted 4 February, 2019; originally announced February 2019.

Comments: 18 pages, 10 figures, accepted for publication in MNRAS. Code available at http://github.com/tberlok/psecas

The Sunyaev-Zel'dovich effect of simulated jet-inflated bubbles in clusters

Authors: Kristian Ehlert, Christoph Pfrommer, Rainer Weinberger, Rüdiger Pakmor, Volker Springel

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

Submitted 17 December, 2018; originally announced December 2018.

Comments: 7 pages, 4 figures, submitted to ApJ letters. Comments welcome!

ETHOS - an Effective Theory of Structure Formation: detecting dark matter interactions through the Lyman-$α$ forest

Authors: Sownak Bose, Mark Vogelsberger, Jesús Zavala, Christoph Pfrommer, Francis-Yan Cyr-Racine, Sebastian Bohr, Torsten Bringmann

Abstract: 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-$α$ forest in cold dark matter (CDM) with an alternative model in which the DM couples strongly wi… ▽ More 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-$α$ 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-$α$ 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-$α$ 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. △ Less

Submitted 30 May, 2019; v1 submitted 26 November, 2018; originally announced November 2018.

Comments: 15 pages, 12 figures, 1 table. Updated version accepted for publication in MNRAS

Journal ref: MNRAS 487, 522 (2019)

Missing Gamma-ray Halos and the Need for New Physics in the Gamma-ray Sky

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

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

Submitted 26 October, 2018; v1 submitted 8 August, 2018; originally announced August 2018.

Comments: Accepted to ApJ, 20 pages, 15 figures

The physics of multiphase gas flows: fragmentation of a radiatively cooling gas cloud in a hot wind

Authors: Martin Sparre, Christoph Pfrommer, Mark Vogelsberger

Abstract: 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… ▽ More 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-$α$ haloes, high-velocity clouds, and broad-line regions of active galactic nuclei. △ Less

Submitted 7 November, 2018; v1 submitted 20 July, 2018; originally announced July 2018.

Comments: Accepted for publication in MNRAS

Faraday rotation maps of disk galaxies

Authors: Ruediger Pakmor, Thomas Guillet, Christoph Pfrommer, Facundo A. Gomez, Robert J. J. Grand, Federico Marinacci, Christine M. Simpson, Volker Springel

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

Submitted 19 September, 2018; v1 submitted 5 July, 2018; originally announced July 2018.

Comments: 10 pages, 10 figures, accepted by MNRAS

Simulations of the dynamics of magnetised jets and cosmic rays in galaxy clusters

Authors: Kristian Ehlert, Rainer Weinberger, Christoph Pfrommer, Rüdiger Pakmor, Volker Springel

Abstract: 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én waves that are generated by streaming cosmic rays (CRs). In order to study this idea, we use three-dimensional magneto-hydrody… ▽ More 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én 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. △ Less

Submitted 27 November, 2020; v1 submitted 14 June, 2018; originally announced June 2018.

Comments: 25 pages, 20 figures, accepted by MNRAS

Journal ref: Mon. Not. R. Astron. Soc., 481, 2878-2900 (2018)

Cosmic-ray hydrodynamics: Alfvén-wave regulated transport of cosmic rays

Authors: Timon Thomas, Christoph Pfrommer

Abstract: 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 fi… ▽ More 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énic turbulence provides the main source of CR scattering, and causes CRs to stream along the magnetic field with the Alfvén 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én-wave energy. We account for energy exchange of CRs and Alfvén 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. △ Less

Submitted 27 January, 2019; v1 submitted 28 May, 2018; originally announced May 2018.

Comments: 34 pages, 6 figures, minor revision to match the accepted version to be published in MNRAS

Enhancing AGN efficiency and cool-core formation with anisotropic thermal conduction

Authors: David J. Barnes, Rahul Kannan, Mark Vogelsberger, Christoph Pfrommer, Ewald Puchwein, Rainer Weinberger, Volker Springel, Rüdiger Pakmor, Dylan Nelson, Federico Marinacci, Annalisa Pillepich, Paul Torrey, Lars Hernquist

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

Submitted 14 May, 2018; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 5 pages, 4 figures, submitted to MNRAS, comments welcome

Constraining the coherence scale of the interstellar magnetic field using TeV gamma-ray observations of supernova remnants

Authors: Matteo Pais, Christoph Pfrommer, Kristian Ehlert, Maria Werhahn, Georg Winner

Abstract: 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-ra… ▽ More 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 $λ_B$ around Vela Jr. and SN1006 to $λ_B \simeq 13_{-4.3}^{+13}$ pc and $λ_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. △ Less

Submitted 9 June, 2020; v1 submitted 8 May, 2018; originally announced May 2018.

Comments: 16 pages, 9 figures

MSC Class: 85A04 (Primary); 85A30 (Secondary)

IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays

Authors: Francois Boulanger, Torsten Ensslin, Andrew Fletcher, Philipp Girichides, Stefan Hackstein, Marijke Haverkorn, Joerg R. Hoerandel, Tess R. Jaffe, Jens Jasche, Michael Kachelriess, Kumiko Kotera, Christoph Pfrommer, Jorg P. Rachen, Luiz F. S. Rodrigues, Beatriz Ruiz-Granados, Amit Seta, Anvar Shukurov, Gunter Sigl, Theo Steininger, Valentina Vacca, Ellert van der Velden, Arjen van Vliet, Jiaxin Wang

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

Submitted 6 September, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: 56 pages, 15 figures, minor updates to match published JCAP version

Journal ref: JCAP08(2018)049

The effect of cosmic-ray acceleration on supernova blast wave dynamics

Authors: M. Pais, C. Pfrommer, K. Ehlert, R. Pakmor

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

Submitted 30 April, 2018; originally announced May 2018.

Comments: 19 pages, 18 figures

Growth of beam-plasma instabilities in the presence of background inhomogeneity

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

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

Submitted 13 April, 2018; originally announced April 2018.

Comments: 10 pages, 6 figures, Accepted by ApJ, comments welcome

The dependence of cosmic ray driven galactic winds on halo mass

Authors: Svenja Jacob, Rüdiger Pakmor, Christine M. Simpson, Volker Springel, Christoph Pfrommer

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

Submitted 13 December, 2017; originally announced December 2017.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

ETHOS -- an effective theory of structure formation: predictions for the high-redshift Universe -- abundance of galaxies and reionization

Authors: Mark R. Lovell, Jesús Zavala, Mark Vogelsberger, Xuejian Shen, Francis-Yan Cyr-Racine, Christoph Pfrommer, Kris Sigurdson, Michael Boylan-Kolchin, Annalisa Pillepich

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

Submitted 9 January, 2019; v1 submitted 28 November, 2017; originally announced November 2017.

Comments: 14 pages, 8 figures, version accepted for MNRAS publication. Please contact lovell@hi.is with any comments

Simulating Gamma-ray Emission in Star-forming Galaxies

Authors: C. Pfrommer, R. Pakmor, C. M. Simpson, V. Springel

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

Submitted 15 September, 2017; originally announced September 2017.

Comments: Accepted by ApJL, 8 pages, 4 figures

MERGHERS: An SZ-selected cluster survey with MeerKAT

Authors: Kenda Knowles, Andrew Baker, Kaustuv Basu, Vijaysarathi Bharadwaj, Roger Deane, Mark Devlin, Simon Dicker, Francesco de Gasperin, Chiara Ferrari, Matt Hilton, John P. Hughes, Huib T. Intema, Sphesihle Makhathini, Kavilan Moodley, Nadeem Oozeer, Christoph Pfrommer, Jonathan Sievers, Sinenhlanhla P. Sikhosana, Oleg Smirnov, Martin W. Sommer, Sara Stanchfield, Kurt van der Heyden, Jonathan T. L. Zwart

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

Submitted 11 September, 2017; originally announced September 2017.

Comments: 12 pages, 3 figures, proceedings from MeerKAT Science: On the Pathway to the SKA, 25-27 May, 2016

Importance of resolving the spectral support of beam-plasma instabilities in simulations

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

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

Submitted 4 October, 2017; v1 submitted 31 March, 2017; originally announced April 2017.

Comments: 7 pages, 9 figures, accepted by ApJ

Simulating the interaction of jets with the intra-cluster medium

Authors: Rainer Weinberger, Kristian Ehlert, Christoph Pfrommer, Rüdiger Pakmor, Volker Springel

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

Submitted 27 March, 2017; originally announced March 2017.

Comments: 18 pages, 15 figures, submitted to MNRAS, comments welcome

SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code

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

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

Submitted 11 June, 2017; v1 submitted 15 February, 2017; originally announced February 2017.

Comments: 26 pages, 19 figures, discussion about performance is added, published in ApJ

Constraints on the Intergalactic Magnetic Field from Bow Ties in the Gamma-ray Sky

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

Abstract: 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 cool… ▽ More 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σ$ with current-day strengths in the range $10^{-16}$ to $10^{-15}$ G, and at 2$σ$ 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. △ Less

Submitted 13 February, 2020; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: 12 pages, 14 figures, 1 appendix. Accepted to ApJ

Bow Ties in the Sky II: Searching for Gamma-ray Halos in the Fermi Sky Using Anisotropy

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

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

Submitted 8 February, 2017; originally announced February 2017.

Comments: Submitted to ApJ, 22 pages, 27 figures, 9 pages in appendix

Magnetic field formation in the Milky Way-like disk galaxies of the Auriga project

Authors: Ruediger Pakmor, Facundo A. Gomez, Robert J. J. Grand, Federico Marinacci, Christine M. Simpson, Volker Springel, David J. R. Campbell, Carlos S. Frenk, Thomas Guillet, Christoph Pfrommer, Simon D. M. White

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

Submitted 2 May, 2017; v1 submitted 24 January, 2017; originally announced January 2017.

Comments: 16 pages, 16 figures, accepted by MNRAS, section 4.1 (turbulent dynamo) extended substantially, now includes magnetic and kinetic power spectra

Increasing blackhole feedback induced quenching with anisotropic thermal conduction

Authors: Rahul Kannan, Mark Vogelsberger, Christoph Pfrommer, Rainer Weinberger, Volker Springel, Lars Hernquist, Ewald Puchwein, Ruediger Pakmor

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

Submitted 5 December, 2016; originally announced December 2016.

Comments: 5 pages, 5 figures, Submitted to ApJ Letters

Turbulence and Particle Acceleration in Giant Radio Haloes: the Origin of Seed Electrons

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

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

Submitted 22 November, 2016; originally announced November 2016.

Comments: 18 pages, 7 figures, 1 table. This article supersedes arXiv:1503.07870, and has a significantly enlarged analysis and presentation. To be published in MNRAS, ref. MN-15-2504-MJ.R3

The Linear Instability of Dilute Ultrarelativistic $e^{\pm}$ Pair Beams

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

Abstract: 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 inst… ▽ More 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{ü}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. △ Less

Submitted 6 October, 2016; originally announced October 2016.

Comments: 13 pages, 4 figure, accepted to ApJ

Cosmic ray heating in cool core clusters - II. Self-regulation cycle and non-thermal emission

Authors: Svenja Jacob, Christoph Pfrommer

Abstract: 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, rad… ▽ More 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{é}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. △ Less

Submitted 26 April, 2017; v1 submitted 20 September, 2016; originally announced September 2016.

Comments: 20 pages including 4 pages of appendix, 11 figures, 4 tables. Accepted for publication in MNRAS

Journal ref: MNRAS 467, 1478-1495 (2017)

Cosmic ray heating in cool core clusters - I. Diversity of steady state solutions

Authors: Svenja Jacob, Christoph Pfrommer

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

Submitted 26 April, 2017; v1 submitted 20 September, 2016; originally announced September 2016.

Comments: 31 pages including 14 pages of appendix, 21 figures (13 in appendix), 3 tables (1 in appendix). Accepted for publication in MNRAS

Journal ref: MNRAS 467, 1449-1477 (2017)

Bow Ties in the Sky I: The Angular Structure of Inverse Compton Gamma-ray Halos in the Fermi Sky

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

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

Submitted 1 September, 2016; originally announced September 2016.

Comments: Accepted by ApJ, 18 pages, 17 figures, 13 pages in appendix

Cosmic ray-driven galactic winds: streaming or diffusion?

Authors: Joshua Wiener, Christoph Pfrommer, S. Peng Oh

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

Submitted 18 January, 2017; v1 submitted 8 August, 2016; originally announced August 2016.

Comments: 18 pages, 11 figures

The role of cosmic ray pressure in accelerating galactic outflows

Authors: Christine M. Simpson, Ruediger Pakmor, Federico Marinacci, Christoph Pfrommer, Volker Springel, Simon C. O. Glover, Paul C. Clark, Rowan J. Smith

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

Submitted 16 August, 2016; v1 submitted 7 June, 2016; originally announced June 2016.

Comments: 8 pages, 4 figures, accepted for publication in ApJL; movie of simulated gas densities can be found here: http://www.h-its.org/tap-images/galactic-outflows/

Galactic winds driven by isotropic and anisotropic cosmic ray diffusion in disk galaxies

Authors: Ruediger Pakmor, Christoph Pfrommer, Christine M. Simpson, Volker Springel

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

Submitted 16 July, 2016; v1 submitted 2 May, 2016; originally announced May 2016.

Comments: 7 pages, 5 figures, accepted by ApJL

Semi-implicit anisotropic cosmic ray transport on an unstructured moving mesh

Authors: Ruediger Pakmor, Christoph Pfrommer, Christine M. Simpson, Rahul Kannan, Volker Springel

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

Submitted 16 July, 2016; v1 submitted 28 April, 2016; originally announced April 2016.

Comments: 16 pages, 16 figures, accepted by MNRAS

Shock finding on a moving-mesh: II. Hydrodynamic shocks in the Illustris universe

Authors: Kevin Schaal, Volker Springel, Rüdiger Pakmor, Christoph Pfrommer, Dylan Nelson, Mark Vogelsberger, Shy Genel, Annalisa Pillepich, Debora Sijacki, Lars Hernquist

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

Submitted 27 July, 2016; v1 submitted 25 April, 2016; originally announced April 2016.

Comments: 26 pages, 15 figures, 1 movie, published in mnras

Journal ref: MRNAS 461, 4441-4465 (2016)

Simulating cosmic ray physics on a moving mesh

Authors: C. Pfrommer, R. Pakmor, K. Schaal, C. M. Simpson, V. Springel

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

Submitted 24 November, 2016; v1 submitted 25 April, 2016; originally announced April 2016.

Comments: Accepted by MNRAS, 31 pages, 10 figures

Deep observation of the NGC 1275 region with MAGIC: search of diffuse gamma-ray emission from cosmic rays in the Perseus cluster

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, E. Carmona, A. Carosi, A. Chatterjee, R. Clavero , et al. (132 additional authors not shown)

Abstract: 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… ▽ More 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 $α$. For $α=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 $α=2.2$, $B_0\gtrsim5-8 μ$G, which is below the 25 $μ$G inferred from Faraday rotation measurements, whereas, for $α\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. △ Less

Submitted 9 February, 2016; originally announced February 2016.

Comments: 13 pages, 8 figures, accepted for publication in Section 3. Cosmology (including clusters of galaxies) of Astronomy and Astrophysics. Corresponding authors: P. Colin (colin@mppmu.mpg.de) and F. Zandanel (f.zandanel@uva.nl)

Journal ref: A&A 589, A33 (2016)

ETHOS - An Effective Theory of Structure Formation: Dark matter physics as a possible explanation of the small-scale CDM problems

Authors: Mark Vogelsberger, Jesus Zavala, Francis-Yan Cyr-Racine, Christoph Pfrommer, Torsten Bringmann, Kris Sigurdson

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

Submitted 4 June, 2016; v1 submitted 16 December, 2015; originally announced December 2015.

Comments: 17 pages, 13 figures, 4 tables. MNRAS accepted. This paper is part of a series of papers on constructing an effective theory of structure formation (ETHOS) that maps almost any microphysical model of dark matter physics to effective parameters for cosmological structure formation. Images and movies can be found at http://mvogelsb.scripts.mit.edu/ethos.php

Journal ref: MNRAS 460, 1399 (2016)

ETHOS - An Effective Theory of Structure Formation: From dark particle physics to the matter distribution of the Universe

Authors: Francis-Yan Cyr-Racine, Kris Sigurdson, Jesus Zavala, Torsten Bringmann, Mark Vogelsberger, Christoph Pfrommer

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

Submitted 11 September, 2017; v1 submitted 16 December, 2015; originally announced December 2015.

Comments: 16 pages + Appendix, 4 figures. Published in Phys. Rev. D. This paper is part of a series of papers on constructing an effective theory of structure formation (ETHOS) that maps almost any microphysical model of dark matter physics to effective parameters for cosmological structure formation. v3: Matches accepted version. v4: Updated definition of dark radiation perturbation variables

Journal ref: Phys. Rev. D 93, 123527 (2016)

Constraints on the cosmic ray cluster physics from a very deep observation of the Perseus cluster with MAGIC

Authors: Pierre Colin, Fabio Zandanel, Monica Vazquez Acosta, Joaquim Palacio, Christoph Pfrommer, Anders Pinzke

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

Submitted 20 August, 2015; originally announced August 2015.

Comments: Proceedings of The 34th International Cosmic Ray Conference (ICRC 2015), 8 pages, 8 figures

Turbulence and Particle Acceleration in Giant Radio Halos: the Origin of Seed Electrons

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

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

Submitted 26 March, 2015; originally announced March 2015.

Comments: 6 pages, 2 figures. Submitted

Patchy blazar heating: diversifying the thermal history of the intergalactic medium

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

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

Submitted 25 August, 2015; v1 submitted 27 February, 2015; originally announced February 2015.