Showing 1–50 of 168 results for author: Pfrommer, C

The theory of resonant cosmic ray-driven instabilities -- Growth and saturation of single modes

Authors: Rouven Lemmerz, Mohamad Shalaby, Christoph Pfrommer, Timon Thomas

Abstract: Cosmic ray (CR) feedback is critical for galaxy formation as CRs drive galactic winds, regularize star formation in galaxies, and escape from active galactic nuclei to heat the cooling cores of galaxy clusters. The feedback strength of CRs depends on their coupling to the background plasma and, as such, on the effective CR transport speed. Traditionally, this has been hypothesized to depend on the… ▽ More Cosmic ray (CR) feedback is critical for galaxy formation as CRs drive galactic winds, regularize star formation in galaxies, and escape from active galactic nuclei to heat the cooling cores of galaxy clusters. The feedback strength of CRs depends on their coupling to the background plasma and, as such, on the effective CR transport speed. Traditionally, this has been hypothesized to depend on the balance between wave growth of CR-driven instabilities and their damping. Here, we study the physics of CR-driven instabilities from first principles, starting from a gyrotropic distribution of CR ions that stream along a background magnetic field. We develop a theory of the underlying processes that organize the particles' orbits and in particular their gyrophases, which provides an intuitive physical picture of (i) wave growth as the CR gyrophases start to bunch up lopsidedly towards the local wave magnetic field, (ii) instability saturation as a result of CRs overtaking the wave and damping its amplitude without isotropizing CRs in the wave frame, and (iii) CR back-reaction onto the unstable plasma waves as the CR gyrophases follow a pendulum motion around the wave magnetic field. Using our novel fluid-particle-in-cell code fluid-SHARP, we validate our theory on the evolution and excitation of individual unstable modes, such as forward and backward propagating Alfvén and whistler waves. We show that these kinetic simulations support our theoretical considerations, thus potentially foreshadowing a fundamental revision of the theory of CR transport in galaxies and galaxy clusters. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: Comments are welcome

Why are thermally- and cosmic ray-driven galactic winds fundamentally different?

Authors: Timon Thomas, Christoph Pfrommer, Rüdiger Pakmor

Abstract: Galactic outflows influence the evolution of galaxies not only by expelling gas from their disks but also by injecting energy into the circumgalactic medium (CGM). This alters or even prevents the inflow of fresh gas onto the disk and thus reduces the star formation rate. Supernovae (SNe) are the engines of galactic winds as they release thermal and kinetic energy into the interstellar medium (ISM… ▽ More Galactic outflows influence the evolution of galaxies not only by expelling gas from their disks but also by injecting energy into the circumgalactic medium (CGM). This alters or even prevents the inflow of fresh gas onto the disk and thus reduces the star formation rate. Supernovae (SNe) are the engines of galactic winds as they release thermal and kinetic energy into the interstellar medium (ISM). Cosmic rays (CRs) are accelerated at the shocks of SN remnants and only constitute a small fraction of the overall SN energy budget. However, their long live-times allow them to act far away from the original injection site and thereby to participate in the galactic wind launching process. Using high-resolution simulations of an isolated Milky Way-type galaxy with the moving-mesh code Arepo and the new multi-phase ISM model Crisp (Cosmic Rays and InterStellar Physics), we investigate how SNe and CRs launch galactic outflows and how the inclusion of CR-mediated feedback boosts the energy and mass entrained in the galactic wind. We find that the majority of thermal SN energy and momentum is used for stirring turbulence either directly or indirectly by causing fountain flows, thereby self-regulating the ISM and not for efficiently driving outflows to large heights. A simulation without CRs only launches a weak galactic outflow at uniformly high temperatures and low densities by means of the thermal pressure gradient. By contrast, most of the CR energy accelerated at SN remnants ($\sim80\%$) escapes the ISM and moves into the CGM. In the inner CGM, CRs dominate the overall pressure and are able to accelerate a large mass fraction in a galactic wind. This wind is turbulent and multi-phase with cold cloudlets embedded in dilute gas at intermediate temperatures ($\sim10^5$ K) and the CGM shows enhanced OVI and CIV absorption in comparison to a simulation without CRs. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 16 pages, 9 figures, submitted to A&A, comments welcome

Comparing the interstellar and circumgalactic origin of gas in the tails of jellyfish galaxies

Authors: Martin Sparre, Christoph Pfrommer, Ewald Puchwein

Abstract: Simulations and observations find long tails in jellyfish galaxies, which are commonly thought to originate from ram-pressure stripped gas of the interstellar medium (ISM) in the immediate galactic wake. While at larger distances from the galaxy, they have been claimed to form in situ owing to thermal instability and fast radiative cooling of mixed ISM and intracluster medium (ICM). In this paper,… ▽ More Simulations and observations find long tails in jellyfish galaxies, which are commonly thought to originate from ram-pressure stripped gas of the interstellar medium (ISM) in the immediate galactic wake. While at larger distances from the galaxy, they have been claimed to form in situ owing to thermal instability and fast radiative cooling of mixed ISM and intracluster medium (ICM). In this paper, we use magneto-hydrodynamical windtunnel simulations of a galaxy with the arepo code to study the origin of gas in the tails of jellyfish galaxies. To this end, we model the galaxy orbit in a cluster by accounting for a time-varying galaxy velocity, ICM density and turbulent magnetic field. Tracking gas flows between the ISM, the circumgalactic medium (CGM) and the ICM, we find contrary to popular opinion that the majority of gas in the tail originated in the CGM. Prior to the central passage of the jellyfish galaxy in the cluster, the CGM is directly transported to the clumpy jellyfish tail that has been shattered into small cloudlets. After the central cluster passage, gas in the tail originates both from the initial ISM and the CGM, but that from the latter was accreted to the galactic ISM before being ram-pressure stripped to form filamentary tentacles in the tail. Our simulation shows a declining gas metallicity in the tail as a function of downstream distance from the galaxy. We conclude that the CGM plays an important role in shaping the tails of jellyfish galaxies. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 9 pages, 5 figures, submitted to A&A (Main Journal)

Thermal conductivity with bells and whistlers: suppression of the magnetothermal instability in galaxy clusters

Authors: Lorenzo Maria Perrone, Thomas Berlok, Christoph Pfrommer

Abstract: In the hot intracluster medium (ICM) in galaxy clusters, plasma microinstabilities may play an important role in the transport of heat and momentum on the large scales. In this paper, we continue our investigation of the effect of whistler suppression of thermal conductivity on the magneto-thermal instability (MTI), which may be active in the periphery of galaxy clusters and contribute to the obse… ▽ More In the hot intracluster medium (ICM) in galaxy clusters, plasma microinstabilities may play an important role in the transport of heat and momentum on the large scales. In this paper, we continue our investigation of the effect of whistler suppression of thermal conductivity on the magneto-thermal instability (MTI), which may be active in the periphery of galaxy clusters and contribute to the observed turbulence. We use a closure for the heat flux inspired by kinetic simulations and show that MTI turbulence with whistler suppression exhibits a critical transition: for modest suppression of the conductivity, the MTI turbulent velocities decrease in agreement with previous MTI scaling laws. However, for suppression above a critical threshold, the MTI loses its ability to maintain equipartition-level magnetic fields through a small-scale dynamo, and the system enters a ``death-spiral''. We propose a model to explain this critical transition, and speculate that conditions in the hot ICM are favourable to the upkeep of the dynamo. Additionaly, with whistler suppression high-$β$ regions are brought out of thermal equilibrium while the efficiency of MTI turbulent driving is reduced. Finally, we show that external turbulence interferes with the MTI and leads to lower levels of turbulence. While individually both external turbulence and whistler suppression weaken the MTI, we find that they can exhibit a complex interplay when acting in conjunction, with external turbulence boosting the whistler-suppressed thermal conductivity and even reviving a ``dead'' MTI. Our study illustrates how extending magnetohydrodynamics with a simple prescription for microscale plasma physics can lead to the formation of a complicated dynamical system and demonstrates that further work is needed in order to bridge the gap between micro- and macro scales in galaxy clusters. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 26 pages, 17 figures. Submitted to A&A. Comments welcome!

The magnetised and thermally unstable tails of jellyfish galaxies

Authors: Martin Sparre, Christoph Pfrommer, Ewald Puchwein

Abstract: Jellyfish galaxies are promising laboratories for studying radiative cooling and magnetic fields in multiphase gas flows. Their long, dense tails are observed to be magnetised, and they extend up to 100 kpc into the intracluster medium (ICM), suggesting that their gas is thermally unstable so that the cold gas mass grows with time rather than being fully dissolved in the hot wind as a result of hy… ▽ More Jellyfish galaxies are promising laboratories for studying radiative cooling and magnetic fields in multiphase gas flows. Their long, dense tails are observed to be magnetised, and they extend up to 100 kpc into the intracluster medium (ICM), suggesting that their gas is thermally unstable so that the cold gas mass grows with time rather than being fully dissolved in the hot wind as a result of hydrodynamical interface instabilities. In this paper we use the AREPO code to perform magnetohydrodynamical windtunnel simulations of a jellyfish galaxy experiencing ram-pressure stripping by interacting with an ICM wind. The ICM density, temperature and velocity that the galaxy encounters are time-dependent and comparable to what a real jellyfish galaxy experiences while orbiting the ICM. In simulations with a turbulent magnetised wind we reproduce observations, which show that the magnetic field is aligned with the jellyfish tails. During the galaxy infall into the cluster with a near edge-on geometry, the gas flow in the tail is fountain-like, implying preferential stripping of gas where the rotational velocity vectors add up with the ram pressure while fall-back occurs in the opposite case. Hence, the tail velocity shows a memory of the rotation pattern of the disc. At the time of the nearest cluster passage, ram-pressure stripping is so strong that the fountain flow is destroyed and instead the tail is dominated by removal of gas. We show that gas in the tail is very fragmentative, which is a prediction of shattering due to radiative cooling. △ Less

Submitted 20 November, 2023; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 14 pages, 7 figures, Accepted for publication in MNRAS

Does the magneto-thermal instability survive whistler-suppression of thermal conductivity in galaxy clusters?

Authors: Lorenzo Maria Perrone, Thomas Berlok, Christoph Pfrommer

Abstract: The hot and dilute intracluster medium (ICM) plays a central role in many key processes that shape galaxy clusters. Nevertheless, the nature of plasma turbulence and particle transport in the ICM remain poorly understood and quantifying the effect of kinetic plasma instabilities on the macroscopic dynamics represents an outstanding problem. Here we focus on the impact of whistler-wave suppression… ▽ More The hot and dilute intracluster medium (ICM) plays a central role in many key processes that shape galaxy clusters. Nevertheless, the nature of plasma turbulence and particle transport in the ICM remain poorly understood and quantifying the effect of kinetic plasma instabilities on the macroscopic dynamics represents an outstanding problem. Here we focus on the impact of whistler-wave suppression of the heat flux on the magneto-thermal instability (MTI), which is expected to drive significant turbulent motions in the periphery of galaxy clusters. We perform small-scale Boussinesq simulations with a sub-grid closure for the thermal diffusivity in the regime of whistler-wave suppression. Our model is characterized by a single parameter that quantifies the collisionality of the ICM on the astrophysical scales of interest that we tune to explore a range appropriate for the periphery of galaxy clusters. We find that the MTI is qualitatively unchanged for weak whistler-suppression. Conversely, with strong suppression the magnetic dynamo is interrupted and MTI-turbulence dies out. In the astrophysically relevant limit, however, the MTI is likely to be supplemented by additional sources of turbulence. Investigating this scenario, we show that the inclusion of external forcing has a beneficial impact and revives even MTI simulations with strong whistler-suppression. As a result, the plasma remains buoyantly unstable, with important consequences for turbulent mixing in the ICM. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 7 pages, 4 figures, accepted for publication in A&A

Magnetic field amplification in cosmological zoom simulations from dwarf galaxies to galaxy groups

Authors: Ruediger Pakmor, Rebekka Bieri, Freeke van de Voort, Maria Werhahn, Azadeh Fattahi, Thomas Guillet, Christoph Pfrommer, Volker Springel, Rosie Y. Talbot

Abstract: Magnetic fields are ubiquitous in the Universe. Recently, cosmological simulations of galaxies have successfully begun to incorporate magnetic fields and their evolution in galaxies and their haloes. However, so far they have mostly focused on Milky Way-like galaxies. Here we analyse a sample of high resolution cosmological zoom simulations of disc galaxies in haloes with mass $M_\mathrm{200c}$ fr… ▽ More Magnetic fields are ubiquitous in the Universe. Recently, cosmological simulations of galaxies have successfully begun to incorporate magnetic fields and their evolution in galaxies and their haloes. However, so far they have mostly focused on Milky Way-like galaxies. Here we analyse a sample of high resolution cosmological zoom simulations of disc galaxies in haloes with mass $M_\mathrm{200c}$ from $10^{10}\,\mathrm{M}_\odot$ to $10^{13}\,\mathrm{M}_\odot$, simulated with the Auriga galaxy formation model. We show that with sufficient numerical resolution the magnetic field amplification and saturation is converged. The magnetic field strength reaches equipartition with turbulent energy density for galaxies in haloes with $M_\mathrm{200c}\gtrsim 10^{11.5}\,\mathrm{M_\odot}$. For galaxies in less massive haloes, the magnetic field strength saturates at a fraction of equipartition that decreases with decreasing halo mass. For our lowest mass haloes, the magnetic field saturates significantly below $10\%$ of equipartition. We quantify the resolution we need to obtain converged magnetic field strengths and discuss our resolution requirements also in the context of the IllustrisTNG cosmological box simulations. We show that, at $z=0$, rotation-dominated galaxies in our sample exhibit for the most part an ordered large scale magnetic field, with fewer field reversals in more massive galaxies. Finally, we compare the magnetic fields in our cosmological galaxies at $z=0$ with simulations of isolated galaxies in a collapsing halo setup. Our results pave the way for detailed studies of cosmic rays and other physical processes in similar cosmological galaxy simulations that crucially depend on the strength and structure of magnetic fields. △ Less

Submitted 9 January, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: 19 pages, 12 figures, accepted for publication by MNRAS, comments welcome

Cosmic ray feedback in galaxies and galaxy clusters -- A pedagogical introduction and a topical review of the acceleration, transport, observables, and dynamical impact of cosmic rays

Authors: Mateusz Ruszkowski, Christoph Pfrommer

Abstract: Understanding the physical mechanisms that control galaxy formation is a fundamental challenge in contemporary astrophysics. Recent advances in the field of astrophysical feedback strongly suggest that cosmic rays (CRs) may be crucially important for our understanding of cosmological galaxy formation and evolution. The appealing features of CRs are their relatively long cooling times and relativel… ▽ More Understanding the physical mechanisms that control galaxy formation is a fundamental challenge in contemporary astrophysics. Recent advances in the field of astrophysical feedback strongly suggest that cosmic rays (CRs) may be crucially important for our understanding of cosmological galaxy formation and evolution. The appealing features of CRs are their relatively long cooling times and relatively strong dynamical coupling to the gas. In galaxies, CRs can be close to equipartition with the thermal, magnetic, and turbulent energy density in the interstellar medium, and can be dynamically very important in driving large-scale galactic winds. Similarly, CRs may provide a significant contribution to the pressure in the circumgalactic medium. In galaxy clusters, CRs may play a key role in addressing the classic cooling flow problem by facilitating efficient heating of the intracluster medium and preventing excessive star formation. Overall, the underlying physics of CR interactions with plasmas exhibit broad parallels across the entire range of scales characteristic of the interstellar, circumgalactic, and intracluster media. Here we present a review of the state-of-the-art of this field and provide a pedagogical introduction to cosmic ray plasma physics, including the physics of wave-particle interactions, acceleration processes, CR spatial and spectral transport, and important cooling processes. The field is ripe for discovery and will remain the subject of intense theoretical, computational, and observational research over the next decade with profound implications for the interpretation of the observations of stellar and supermassive black hole feedback spanning the entire width of the electromagnetic spectrum and multi-messenger data. △ Less

Submitted 6 October, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: invited A&ARv review; revised version; accepted for publication; 238 pages

Deciphering the physical basis of the intermediate-scale instability

Authors: Mohamad Shalaby, Timon Thomas, Christoph Pfrommer, Rouven Lemmerz, Virginia Bresci

Abstract: We study the underlying physics of cosmic-ray (CR) driven instabilities that play a crucial role for CR transport across a wide range of scales, from interstellar to galaxy cluster environments. By examining the linear dispersion relation of CR-driven instabilities in a magnetised electron-ion background plasma, we establish that both, the intermediate and gyroscale instabilities have a resonant o… ▽ More We study the underlying physics of cosmic-ray (CR) driven instabilities that play a crucial role for CR transport across a wide range of scales, from interstellar to galaxy cluster environments. By examining the linear dispersion relation of CR-driven instabilities in a magnetised electron-ion background plasma, we establish that both, the intermediate and gyroscale instabilities have a resonant origin and show that these resonances can be understood via a simple graphical interpretation. These instabilities destabilise wave modes parallel to the large-scale background magnetic field at significantly distinct scales and with very different phase speeds. Furthermore, we show that approximating the electron-ion background plasma with either magnetohydrodynamics (MHD) or Hall-MHD fails to capture the fastest growing instability in the linear regime, namely the intermediate-scale instability. This finding highlights the importance of accurately characterising the background plasma for resolving the most unstable wave modes. Finally, we discuss the implications of the different phase speeds of unstable modes on particle-wave scattering. Further work is needed to investigate the relative importance of these two instabilities in the non-linear, saturated regime and to develop a physical understanding of the effective CR transport coefficients in large-scale CR hydrodynamics theories. △ Less

Submitted 9 November, 2023; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: 13 pages + references, 4 figures. To appear in JPP Letters

Spectrally resolved cosmic rays -- III. Dynamical impact and properties of the circumgalactic medium

Authors: Philipp Girichidis, Maria Werhahn, Christoph Pfrommer, Rüdiger Pakmor, Volker Springel

Abstract: Cosmic rays (CRs) are dynamically important for the formation and evolution of galaxies by regulating star formation and by powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed along the magnetic field. Moreover, both properties sensitively depend on the CR momentum, wh… ▽ More Cosmic rays (CRs) are dynamically important for the formation and evolution of galaxies by regulating star formation and by powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed along the magnetic field. Moreover, both properties sensitively depend on the CR momentum, which is largely unexplored in three-dimensional hydrodynamical simulations. We perform magneto-hydrodynamical simulations of entire galaxies with masses ranging from $10^{10}$ to $10^{12}\,\mathrm{M}_\odot$ and compare dynamically coupled CRs in the grey approximation with a spectrally resolved model that includes CR momenta from $0.1\,\mathrm{GeV}~c^{-1}$ to $100\,\mathrm{TeV}~c^{-1}$. We find that hadronic cooling of CRs dominates over Alfvén cooling, with the latter emulating CR losses as a result of streaming of CRs down their pressure gradient. While star formation rates and galaxy morphologies are only mildly affected by the spectral CR modelling, mass loading factors of galactic outflows can differ by up to a factor of four in dwarf galaxies. All simulated low-mass halos ($M=10^{10}$, $10^{11}$, and $3\times10^{11}\,\mathrm{M}_\odot$) drive strong outflows, where CR transport is temporally dominated by advection. In contrast, the Milky Way-mass galaxy with $M=10^{12}\,\mathrm{M}_\odot$ does not drive sustained outflows, so that CR transport is entirely dominated by diffusion. The effective energy weighted diffusion coefficients vary by two orders of magnitude from the canonical energy-weighted values of $\langle{D}\rangle_{e_\mathrm{cr}}\sim10^{28}\,\mathrm{cm^2\,s^{-1}}$ in the disc up to $3\times10^{29}\,\mathrm{cm^2\,s^{-1}}$ in the circumgalactic medium, where we observe substantial temperature and CR pressure differences between our grey and spectral CR models. △ Less

Submitted 23 December, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 28 pages, MNRAS, accepted

The impact of magnetic fields on cosmological galaxy mergers -- II. Modified angular momentum transport and feedback

Authors: Joseph Whittingham, Martin Sparre, Christoph Pfrommer, Rüdiger Pakmor

Abstract: The role of magnetic fields in galaxy evolution is still an unsolved question in astrophysics. We have previously shown that magnetic fields play a crucial role in major mergers between disc galaxies; in hydrodynamic simulations of such mergers, the Auriga model produces compact remnants with a distinctive bar and ring morphology. In contrast, in magnetohydrodynamic (MHD) simulations, remnants for… ▽ More The role of magnetic fields in galaxy evolution is still an unsolved question in astrophysics. We have previously shown that magnetic fields play a crucial role in major mergers between disc galaxies; in hydrodynamic simulations of such mergers, the Auriga model produces compact remnants with a distinctive bar and ring morphology. In contrast, in magnetohydrodynamic (MHD) simulations, remnants form radially-extended discs with prominent spiral arm structure. In this paper, we analyse a series of cosmological "zoom-in" simulations of major mergers and identify exactly \textit{how} magnetic fields are able to alter the outcome of the merger. We find that magnetic fields modify the transport of angular momentum, systematically hastening the merger progress. The impact of this altered transport depends on the orientation of the field, with a predominantly non-azimuthal (azimuthal) orientation increasing the central baryonic concentration (providing support against collapse). Both effects act to suppress an otherwise existent bar-instability, which in turn leads to a fundamentally different morphology and manifestation of feedback. We note, in particular, that stellar feedback is substantially less influential in MHD simulations, which allows for the later accretion of higher angular momentum gas and the subsequent rapid radial growth of the remnant disc. A corollary of the increased baryonic concentration in MHD simulations is that black holes are able to grow twice as large, although this turns out to have little impact on the remnant's development. Our results show that galaxy evolution cannot be modelled correctly without including magnetic fields. △ Less

Submitted 26 September, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 22 pages, 12 figures. Published in MNRAS on 07 September 2023

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 526, Issue 1, November 2023, Pages 224-245

Coupling multi-fluid dynamics equipped with Landau closures to the particle-in-cell method

Authors: Rouven Lemmerz, Mohamad Shalaby, Timon Thomas, Christoph Pfrommer

Abstract: The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run-times and often to setups in less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the backgr… ▽ More The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run-times and often to setups in less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conservative multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's Law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron, or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison to linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic timescales. △ Less

Submitted 23 January, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: 35 pages, 11 figures, submitted to JPP. Comments are welcome

Gamma-ray emission from spectrally resolved cosmic rays in galaxies

Authors: Maria Werhahn, Philipp Girichidis, Christoph Pfrommer, Joseph Whittingham

Abstract: Cosmic rays (CRs) are ubiquitous in the interstellar medium (ISM) of nearby galaxies, but many of their properties are not well-constrained. Gamma-ray observations provide a powerful tool in this respect, allowing us to constrain both the interaction of CR protons with the ISM and their transport properties. To help better understand the link between observational signatures and CR physics, we use… ▽ More Cosmic rays (CRs) are ubiquitous in the interstellar medium (ISM) of nearby galaxies, but many of their properties are not well-constrained. Gamma-ray observations provide a powerful tool in this respect, allowing us to constrain both the interaction of CR protons with the ISM and their transport properties. To help better understand the link between observational signatures and CR physics, we use a series of magneto-hydrodynamical (MHD) AREPO simulations of isolated galaxies performed using spectrally-resolved CR transport in every computational cell, with subsequent gamma-ray emission calculated using the CRAYON+ (Cosmic RAY emissiON) code. In each of our simulated halos, modelling the energy-dependent spatial diffusion of CRs leads to a more extended distribution of high-energy (~100 GeV) gamma rays compared to that predicted by a 'grey' steady-state model, which is especially visible in the corresponding emission maps and radial profiles. Despite this, the total gamma-ray spectra can often be well approximated by the steady-state model, although recovering the same spectral index typically requires a minor variation of the energy dependence of the diffusion coefficient. Our simulations reproduce the observed spectral indices and gamma-ray spectra of nearby star-forming galaxies and also match recent observations of the far infrared--gamma-ray relation. We find, however, that the spectrally resolved model yields marginally smaller luminosities for lower star formation rates compared to grey simulations of CRs. Our work highlights the importance of modelling spectrally resolved CR transport for an accurate prediction of spatially resolved high-energy gamma-ray emission, as will be probed by the upcoming Cherenkov Telescope Array observatory. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Comments: 18 pages, 11 figures, submitted to MNRAS, comments are welcome!

Active galactic nucleus jet feedback in hydrostatic halos

Authors: Rainer Weinberger, Kung-Yi Su, Kristian Ehlert, Christoph Pfrommer, Lars Hernquist, Greg L. Bryan, Volker Springel, Yuan Li, Blakesley Burkhart, Ena Choi, Claude-André Faucher-Giguère

Abstract: Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution funct… ▽ More Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a $50$ Myr timescale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in halos, without a need for jet precession or wide opening angles. Comparing the jet feedback with a `kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star forming phase in place, while the wind acts more locally on the star forming phase, which manifests itself in different time-variability properties. △ Less

Submitted 5 May, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 22 pages, 20 figures, accepted for publication in MNRAS, comments welcome

The MillenniumTNG Project: The hydrodynamical full physics simulation and a first look at its galaxy clusters

Authors: Ruediger Pakmor, Volker Springel, Jonathan P. Coles, Thomas Guillet, Christoph Pfrommer, Sownak Bose, Monica Barrera, Ana Maria Delgado, Fulvio Ferlito, Carlos Frenk, Boryana Hadzhiyska, César Hernández-Aguayo, Lars Hernquist, Rahul Kannan, Simon D. M. White

Abstract: Cosmological simulations are an important theoretical pillar for understanding nonlinear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large volume simulations of cosmic structure formation aiming to better understand physical processes o… ▽ More Cosmological simulations are an important theoretical pillar for understanding nonlinear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large volume simulations of cosmic structure formation aiming to better understand physical processes on large scales and to help interpreting upcoming large-scale galaxy surveys. We here focus on the full physics box MTNG740 that computes a volume of $(740\,\mathrm{Mpc})^3$ with a baryonic mass resolution of $3.1\times~10^7\,\mathrm{M_\odot}$ using \textsc{arepo} with $80.6$~billion cells and the IllustrisTNG galaxy formation model. We verify that the galaxy properties produced by MTNG740 are consistent with the TNG simulations, including more recent observations. We focus on galaxy clusters and analyse cluster scaling relations and radial profiles. We show that both are broadly consistent with various observational constraints. We demonstrate that the SZ-signal on a deep lightcone is consistent with Planck limits. Finally, we compare MTNG740 clusters with galaxy clusters found in Planck and the SDSS-8 RedMaPPer richness catalogue in observational space, finding very good agreement as well. However, {\it simultaneously} matching cluster masses, richness, and Compton-$y$ requires us to assume that the SZ mass estimates for Planck clusters are underestimated by $0.2$~dex on average. Thanks to its unprecedented volume for a high-resolution hydrodynamical calculation, the MTNG740 simulation offers rich possibilities to study baryons in galaxies, galaxy clusters, and in large scale structure, and in particular their impact on upcoming large cosmological surveys. △ Less

Submitted 6 December, 2022; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 18 pages, 14 figures, accepted for publication by MNRAS, comments welcome

LYRA III: The smallest Reionization survivors

Authors: Thales A. Gutcke, Christoph Pfrommer, Greg L. Bryan, Rüdiger Pakmor, Volker Springel, Thorsten Naab

Abstract: The dividing line between galaxies that are quenched by reionization ("relics") and galaxies that survive reionization (i.e. continue forming stars) is commonly discussed in terms of a halo mass threshold. We probe this threshold in a physically more complete and accurate way than has been possible to date, using five extremely high resolution ($M_\mathrm{target}=4M_\odot$) cosmological zoom-in si… ▽ More The dividing line between galaxies that are quenched by reionization ("relics") and galaxies that survive reionization (i.e. continue forming stars) is commonly discussed in terms of a halo mass threshold. We probe this threshold in a physically more complete and accurate way than has been possible to date, using five extremely high resolution ($M_\mathrm{target}=4M_\odot$) cosmological zoom-in simulations of dwarf galaxies within the halo mass range $1-4\times10^9M_\odot$. The employed LYRA simulation model features resolved interstellar medium physics and individual, resolved supernova explosions. In our results, we discover an interesting intermediate population of dwarf galaxies close to the threshold mass but which are neither full reionization relics nor full reionization survivors. These galaxies initially quench at the time of reionization but merely remain quiescent for ~500Myr. At $z\approx5$ they recommence star formation in a synchronous way, and remain star-forming until the present day. These results demonstrate that the halo mass at $z=0$ is not a good indicator of survival close to the threshold. While the star formation histories we find are diverse, we show that they are directly related to the ability of a given halo to retain and cool gas. Whereas the latter is most strongly dependent on the mass (or virial temperature) of the host halo at the time of reionization, it also depends on its growth history, the UV background (and its decrease at late times) and the amount of metals retained within the halo. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: 15 pages, 12 figures, 2 tables, submitted to ApJ

Comparing energy and entropy formulations for cosmic ray hydrodynamics

Authors: Matthias Weber, Timon Thomas, Christoph Pfrommer

Abstract: Cosmic rays (CRs) play an important role in many astrophysical systems. Acting on plasma scales to galactic environments, CRs are usually modeled as a fluid, using the CR energy density as the evolving quantity. This method comes with the flaw that the corresponding CR evolution equation is not in conservative form as it contains an adiabatic source term that couples CRs to the thermal gas. In the… ▽ More Cosmic rays (CRs) play an important role in many astrophysical systems. Acting on plasma scales to galactic environments, CRs are usually modeled as a fluid, using the CR energy density as the evolving quantity. This method comes with the flaw that the corresponding CR evolution equation is not in conservative form as it contains an adiabatic source term that couples CRs to the thermal gas. In the absence of non-adiabatic changes, instead evolving the CR entropy density is a physically equivalent option that avoids this potential numerical inconsistency. In this work, we study both approaches for evolving CRs in the context of magneto-hydrodynamic (MHD) simulations using the massively parallel moving-mesh code AREPO. We investigate the performance of both methods in a sequence of shock-tube tests with various resolutions and shock Mach numbers. We find that the entropy-conserving scheme performs best for the idealized case of purely adiabatic CRs across the shock while both approaches yield similar results at lower resolution. In this setup, both schemes operate well and almost independently of the shock Mach number. Taking active CR acceleration at the shock into account, the energy-based method proves to be numerically much more stable and significantly more accurate in determining the shock velocity, in particular at low resolution, which is more typical for astrophysical large-scale simulations. For a more realistic application, we simulate the formation of several isolated galaxies at different halo masses and find that both numerical methods yield almost identical results with differences far below common astrophysical uncertainties. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 13 pages, 12 figures

How Cosmic Rays Mediate the Evolution of the Interstellar Medium

Authors: Christine M. Simpson, Rüdiger Pakmor, Christoph Pfrommer, Simon C. O. Glover, Rowan Smith

Abstract: We explore the impact of diffusive cosmic rays (CRs) on the evolution of the interstellar medium (ISM) under varying assumptions of supernova explosion environment. In practice, we systematically vary the relative fractions of supernovae (SN) occurring in star-forming high-density gas and those occurring in random locations decoupled from star-forming gas to account for SN from run-away stars or e… ▽ More We explore the impact of diffusive cosmic rays (CRs) on the evolution of the interstellar medium (ISM) under varying assumptions of supernova explosion environment. In practice, we systematically vary the relative fractions of supernovae (SN) occurring in star-forming high-density gas and those occurring in random locations decoupled from star-forming gas to account for SN from run-away stars or explosions in regions that have been cleared by prior SN, stellar winds, or radiation. We explore various mixed models by adjusting these fractions relative to each other. We find that in the simple system of a periodic stratified gas layer the ISM structure will evolve to one of two solutions: a "peak driving" state where warm gas is volume filling or a "thermal runaway" state where hot gas is volume filling. CR pressure and transport are important factors that strongly influence the solution state the ISM reaches and have the ability to flip the ISM between solutions. Observable signatures such as gamma ray emission and HI gas are explored. We find that gamma ray luminosity from pion decay is largely consistent with observations for a range of model parameters. The thickness of the HI gas layer may be too compact, however, this may be due to a large cold neutral fraction of midplane gas. The volume fraction of hot gas evolves to stable states in both solutions, but neither settles to a Milky Way-like configuration, suggesting that additional physics which is omitted here (e.g. a cosmological circum-galactic medium, radiation transport, or spectrally resolved and spatially varying CR transport) may be required. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 27 pages, 21 Figures, submitted to MNRAS, comments welcome

Escaping the maze: a statistical sub-grid model for cloud-scale density structures in the interstellar medium

Authors: Tobias Buck, Christoph Pfrommer, Philipp Girichidis, Bogdan Corobean

Abstract: The interstellar medium (ISM) is a turbulent, highly structured multi-phase medium. State-of-the-art cosmological simulations of the formation of galactic discs usually lack the resolution to accurately resolve those multi-phase structures. However, small-scale density structures play an important role in the life cycle of the ISM, and determine the fraction of cold, dense gas, the amount of star… ▽ More The interstellar medium (ISM) is a turbulent, highly structured multi-phase medium. State-of-the-art cosmological simulations of the formation of galactic discs usually lack the resolution to accurately resolve those multi-phase structures. However, small-scale density structures play an important role in the life cycle of the ISM, and determine the fraction of cold, dense gas, the amount of star formation and the amount of radiation and momentum leakage from cloud-embedded sources. Here, we derive a $statistical\, model$ to calculate the unresolved small-scale ISM density structure from coarse-grained, volume-averaged quantities such as the $gas\, clumping\, factor$, $\mathcal{C}$, and mean density $\left<ρ\right>_V$. Assuming that the large-scale ISM density is statistically isotropic, we derive a relation between the three-dimensional clumping factor, $\mathcal{C}_ρ$, and the clumping factor of the $4π$ column density distribution on the cloud surface, $\mathcal{C}_Σ$, and find $\mathcal{C}_Σ=\mathcal{C}_ρ^{2/3}$. Applying our model to calculate the covering fraction, i.e., the $4π$ sky distribution of optically thick sight-lines around sources inside interstellar gas clouds, we demonstrate that small-scale density structures lead to significant differences at fixed physical ISM density. Our model predicts that gas clumping increases the covering fraction by up to 30 per cent at low ISM densities compared to a uniform medium. On the other hand, at larger ISM densities, gas clumping suppresses the covering fraction and leads to increased scatter such that covering fractions can span a range from 20 to 100 per cent at fixed ISM density. All data and example code is publicly available at GitHub. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 16 pages with 8 figures, 14 pages main text with 7 figures, 1 page references, 1 page appendix with 1 figure, accepted by MNRAS on April 1st 2022

Self-regulated AGN feedback of light jets in cool-core galaxy clusters

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

Abstract: Heating from active galactic nuclei (AGN) is thought to stabilize cool-core clusters, limiting star formation and cooling flows. We employ radiative magneto-hydrodynamic (MHD) simulations to model light AGN jet feedback with different accretion modes (Bondi-Hoyle-Lyttleton and cold accretion) in an idealised Perseus-like cluster. Independent of the probed accretion model, accretion efficiency, jet… ▽ More Heating from active galactic nuclei (AGN) is thought to stabilize cool-core clusters, limiting star formation and cooling flows. We employ radiative magneto-hydrodynamic (MHD) simulations to model light AGN jet feedback with different accretion modes (Bondi-Hoyle-Lyttleton and cold accretion) in an idealised Perseus-like cluster. Independent of the probed accretion model, accretion efficiency, jet density and resolution, the cluster self-regulates with central entropies and cooling times consistent with observed cool-core clusters in this non-cosmological setting. We find that increased jet efficiencies lead to more intermittent jet powers and enhanced star formation rates. Our fiducial low-density jets can easily be deflected by orbiting cold gaseous filaments, which redistributes angular momentum and leads to more extended cold gas distributions and isotropic bubble distributions. In comparison to our fiducial low momentum-density jets, high momentum-density jets heat less efficiently and enable the formation of a persistent cold-gas disc perpendicular to the jets that is centrally confined. Cavity luminosities measured from our simulations generally reflect the cooling luminosities of the intracluster medium (ICM) and correspond to averaged jet powers that are relatively insensitive to short periods of low-luminosity jet injection. Cold gas structures in our MHD simulations with low momentum-density jets generally show a variety of morphologies ranging from discy to very extended filamentary structures. In particular, magnetic fields are crucial to inhibit the formation of unrealistically massive cold gas discs by redistributing angular momentum between the hot and cold phases and by fostering the formation of elongated cold filaments that are supported by magnetic pressure. △ Less

Submitted 25 March, 2023; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 26 pages, 19 figures, published in MNRS

Journal ref: MNRAS, 518, 4622 (2023)

Cosmic ray-driven galactic winds: transport modes of cosmic rays and Alfvén-wave dark regions

Authors: Timon Thomas, Christoph Pfrommer, Rüdiger Pakmor

Abstract: Feedback mediated by cosmic rays (CRs) is an important process in galaxy formation. Because CRs are long-lived and because they are transported along magnetic field lines independently of any gas flow, they can efficiently distribute their feedback energy within the galaxy. We present an in-depth investigation of (i) how CRs launch galactic winds from a disc that is forming in a… ▽ More Feedback mediated by cosmic rays (CRs) is an important process in galaxy formation. Because CRs are long-lived and because they are transported along magnetic field lines independently of any gas flow, they can efficiently distribute their feedback energy within the galaxy. We present an in-depth investigation of (i) how CRs launch galactic winds from a disc that is forming in a $10^{11} \mathrm{M}_\odot$ halo and (ii) how CR transport affects the dynamics in a galactic outflow. To this end, we use the Arepo moving-mesh code and model CR transport with the two-moment description of CR hydrodynamics. This model includes the CR interaction with gyroresonant Alfvén waves that enables us to self-consistently calculate the CR diffusion coefficient and CR transport speeds based on coarse-grained models for plasma physical effects. This delivers insight into key questions such as whether the effective CR transport is streaming-like or diffusive-like, how the CR diffusion coefficient and transport speed change inside the circumgalactic medium (CGM), and to what degree the two-moment approximation is needed to faithfully capture these effects. We find that the CR-diffusion coefficient reaches a steady-state in most environments with the notable exception of our newly discovered Alfvén-wave dark regions where the toroidal wind magnetic field is nearly perpendicular to the CR pressure gradient so that CRs are unable to excite gyroresonant Alfvén waves. However, CR transport itself cannot reach a steady-state and is not well described by either the CR streaming paradigm, the CR diffusion paradigm or a combination of both. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 19 pages, 12 figures, submitted to MNRAS, comments welcome

The feasibility of constraining DM interactions with high-redshift observations by JWST

Authors: Ali Kurmus, Sownak Bose, Mark Lovell, Francis-Yan Cyr-Racine, Mark Vogelsberger, Christoph Pfrommer, Jesús Zavala

Abstract: Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true with the advent of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as part of the Effective Theory of Structure Formation (ETHOS) project to compare high redshift galaxies in Cold (CDM) and alternative DM m… ▽ More Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true with the advent of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as part of the Effective Theory of Structure Formation (ETHOS) project to compare high redshift galaxies in Cold (CDM) and alternative DM models which have varying relativistic coupling and self-interaction strengths. The interacting DM scenarios produce a cutoff in the linear power spectrum on small-scales, followed by a series of "dark acoustic oscillations". We find that DM interactions suppress the abundance of galaxies below $M_\star \sim 10^8\,M_\odot$ for the models considered. The cutoff in the power spectrum delays structure formation relative to CDM. Objects in ETHOS that end up at the same final masses as their CDM counterparts are characterised by a more vigorous phase of early star formation. While galaxies with $M_\star \lesssim 10^6\,M_\odot$ make up more than 60 per cent of star formation in CDM at $z\approx 10$, they contribute only about half the star formation density in ETHOS. These differences diminish with decreasing redshift. We find that the effects of DM self-interactions are negligible compared to effects of relativistic coupling (i.e. the effective initial conditions for galaxy formation) in all properties of the galaxy population we examine. Finally, we show that the clustering strength of galaxies at high redshifts depends sensitively on DM physics, although these differences are manifest on scales that may be too small to be measurable by JWST. △ Less

Submitted 29 July, 2022; v1 submitted 9 March, 2022; originally announced March 2022.

Comments: 15 pages, 11 figures, 3 Tables. Comments are welcome. Updated based on referee's comments. Accepted to MNRAS on July 29, 2022

The mechanism of efficient electron acceleration at parallel non-relativistic shocks

Authors: Mohamad Shalaby, Rouven Lemmerz, Timon Thomas, Christoph Pfrommer

Abstract: Thermal electrons cannot directly participate in the process of diffusive acceleration at electron-ion shocks because their Larmor radii are smaller than the shock transition width: this is the well-known electron injection problem of diffusive shock acceleration. Instead, an efficient pre-acceleration process must exist that scatters electrons off of electromagnetic fluctuations on scales much sh… ▽ More Thermal electrons cannot directly participate in the process of diffusive acceleration at electron-ion shocks because their Larmor radii are smaller than the shock transition width: this is the well-known electron injection problem of diffusive shock acceleration. Instead, an efficient pre-acceleration process must exist that scatters electrons off of electromagnetic fluctuations on scales much shorter than the ion gyro radius. The recently found intermediate-scale instability provides a natural way to produce such fluctuations in parallel shocks. The instability drives comoving (with the upstream plasma) ion-cyclotron waves at the shock front and only operates when the drift speed is smaller than half of the electron Alfven speed. Here, we perform particle-in-cell simulations with the SHARP code to study the impact of this instability on electron acceleration at parallel non-relativistic, electron-ion shocks. To this end, we compare a shock simulation in which the intermediate-scale instability is expected to grow to simulations where it is suppressed. In particular, the simulation with an Alfvenic Mach number large enough to quench the intermediate instability shows a great reduction (by two orders of magnitude) of the electron acceleration efficiency. Moreover, the simulation with a reduced ion-to-electron mass ratio (where the intermediate instability is also suppressed) not only artificially precludes electron acceleration but also results in erroneous electron and ion heating in the downstream and shock transition regions. This finding opens up a promising route for a plasma physical understanding of diffusive shock acceleration of electrons, which necessarily requires realistic mass ratios in simulations of collisionless electron-ion shocks. △ Less

Submitted 4 May, 2022; v1 submitted 10 February, 2022; originally announced February 2022.

Comments: 12 pages, 8 figures, 1 table, Accepted to ApJ

Constraining blazar heating with the 2<z<3 Lyman-$α$ forest

Authors: A. Lamberts, E. Puchwein, C. Pfrommer, P. Chang, M. Shalaby, A. Broderick, P. Tiede, G. Rudie

Abstract: The intergalactic medium (IGM) acts like a calorimeter recording energy injection by cosmic structure formation, shocks and photoheating from stars and active galactic nuclei. It was recently proposed that spatially inhomogeneous TeV-blazars could significantly heat up the underdense IGM, resulting in patches of both cold and warm IGM around $z\simeq2-3$. The goal of this study is to compare predi… ▽ More The intergalactic medium (IGM) acts like a calorimeter recording energy injection by cosmic structure formation, shocks and photoheating from stars and active galactic nuclei. It was recently proposed that spatially inhomogeneous TeV-blazars could significantly heat up the underdense IGM, resulting in patches of both cold and warm IGM around $z\simeq2-3$. The goal of this study is to compare predictions of different blazar heating models with recent observations of the IGM. We perform a set of cosmological simulations and carefully compute mock observables of the Lyman-$α ($Ly$α$) forest. We perform a detailed assessment of different systematic uncertainties which typically impact this type of observables and find that they are smaller than the differences between our models. We find that our inhomogeneous blazar heating model is in good agreement with the Ly$α$ line properties and the rescaled flux probability distribution function at high redshift ($2.5<z<3$) but that our blazar heating models are challenged by lower redshift data ($2<z<2.5$). Our results could be explained by HeII reionisation although state-of-the-art models fall short on providing enough heating to the low-density IGM, thus motivating further radiative transfer studies of inhomogeneous HeII reionisation. If blazars are indeed hosted by group-mass halos of $2\times10^{13}{M}_\odot$, a later onset of blazar heating in comparison to previous models would be favoured, which could bring our findings here in agreement with the evidence of blazar heating from local gamma-ray observations. △ Less

Submitted 2 February, 2022; v1 submitted 31 January, 2022; originally announced January 2022.

Comments: Resubmitted after minor revisions

Cold and Hot gas distribution around the Milky-Way-M31 system in the HESTIA simulations

Authors: Mitali Damle, Martin Sparre, Philipp Richter, Maan H. Hani, Sebastián E. Nuza, Christoph Pfrommer, Robert J. J. Grand, Yehuda Hoffman, Noam Libeskind, Jenny G. Sorce, Matthias Steinmetz, Elmo Tempel, Mark Vogelsberger, Peng Wang

Abstract: Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three hig… ▽ More Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three high-resolution {\sc arepo}-based simulations of the LG, run using the Auriga galaxy formation model. For this paper, we focus only on the $z = 0$ simulation datasets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (HI and SiIII) are more clumpy in comparison to warmer gas tracers (OVI, OVII and OVIII). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies. HESTIA under-produces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way. △ Less

Submitted 16 March, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 22 pages, 9 figures, revised version, Accepted in MNRAS

A Multiwavelength Dynamical State Analysis of ACT-CL J0019.6+0336

Authors: Denisha S. Pillay, David J. Turner, Matt Hilton, Kenda Knowles, Kabelo C. Kesebonye, Kavilan Moodley, Tony Mroczkowski, Nadeem Oozeer, Christoph Pfrommer, Sinenhlanhla P. Sikhosana, Edward J. Wollack

Abstract: In our study, we show a multiwavelength view of ACT-CL J0019.6+0336 (which hosts a radio halo), to investigate the cluster dynamics, morphology, and ICM. We use a combination of XMM-Newton images, Dark Energy Survey (DES) imaging and photometry, SDSS spectroscopic information, and 1.16 GHz MeerKAT data to study the cluster properties. Various X-ray and optical morphology parameters are calculated… ▽ More In our study, we show a multiwavelength view of ACT-CL J0019.6+0336 (which hosts a radio halo), to investigate the cluster dynamics, morphology, and ICM. We use a combination of XMM-Newton images, Dark Energy Survey (DES) imaging and photometry, SDSS spectroscopic information, and 1.16 GHz MeerKAT data to study the cluster properties. Various X-ray and optical morphology parameters are calculated to investigate the level of disturbance. We find disturbances in two X-ray parameters and the optical density map shows elongated and axisymmetric structures with the main cluster component southeast of the cluster centre and another component northwest of the cluster centre. We also find a BCG offset of $\sim$950 km/s from the mean velocity of the cluster, and a discrepancy between the SZ mass, X-ray mass, and dynamical mass ($M_{X,500}$ and $M_{SZ,500}$ lies >3 $σ$ away from $M_{\rm{dyn},500}$), showing that J0019 is a merging cluster and probably in a post-merging phase. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 15 pages, 5 figures, Accepted for publication in special issue of Galaxies, from the conference "A New Window on the Radio Emission from Galaxies, Galaxy Clusters and Cosmic Web: Current Status and Perspectives"

EuCAPT White Paper: Opportunities and Challenges for Theoretical Astroparticle Physics in the Next Decade

Authors: R. Alves Batista, M. A. Amin, G. Barenboim, N. Bartolo, D. Baumann, A. Bauswein, E. Bellini, D. Benisty, G. Bertone, P. Blasi, C. G. Böhmer, Ž. Bošnjak, T. Bringmann, C. Burrage, M. Bustamante, J. Calderón Bustillo, C. T. Byrnes, F. Calore, R. Catena, D. G. Cerdeño, S. S. Cerri, M. Chianese, K. Clough, A. Cole, P. Coloma , et al. (112 additional authors not shown)

Abstract: Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, und… ▽ More Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: White paper of the European Consortium for Astroparticle Theory (EuCAPT). 135 authors, 400 endorsers, 133 pages, 1382 references

Turning AGN bubbles into radio relics with sloshing: modeling CR transport with realistic physics

Authors: John ZuHone, Kristian Ehlert, Rainer Weinberger, Christoph Pfrommer

Abstract: Radio relics are arc-like synchrotron sources at the periphery of galaxy clusters, produced by cosmic-ray electrons in a $μ$G magnetic field which are believed to have been (re-)accelerated by merger shock fronts. However, not all relics appear at the same location as shocks as seen in the X-ray. In a previous work, we suggested that the shape of some relics may result from the pre-existing spatia… ▽ More Radio relics are arc-like synchrotron sources at the periphery of galaxy clusters, produced by cosmic-ray electrons in a $μ$G magnetic field which are believed to have been (re-)accelerated by merger shock fronts. However, not all relics appear at the same location as shocks as seen in the X-ray. In a previous work, we suggested that the shape of some relics may result from the pre-existing spatial distribution of cosmic-ray electrons, and tested this hypothesis using simulations by launching AGN jets into a cluster atmosphere with sloshing gas motions generated by a previous merger event. We showed that these motions could transport the cosmic ray-enriched material of the AGN bubbles to large radii and stretch it in a tangential direction, producing a filamentary shape resembling a radio relic. In this work, we improve our physical description for the cosmic rays by modeling them as a separate fluid which undergoes diffusion and Alfvén losses. We find that including this additional cosmic ray physics significantly diminishes the appearance of these filamentary features, showing that our original hypothesis is sensitive to the modeling of cosmic ray physics in the intracluster medium. △ Less

Submitted 29 October, 2021; v1 submitted 8 October, 2021; originally announced October 2021.

Comments: 9 pages, 4 figures, accepted by the journal Galaxies for the Special Issue "A New Window on the Radio Emission from Galaxies, Galaxy Clusters and Cosmic Web: Current Status and Perspectives"

Two striking head-tail galaxies in the galaxy cluster IIZW108: insights into transition to turbulence, magnetic fields and particle re-acceleration

Authors: Ancla Müller, Christoph Pfrommer, Alessandro Ignesti, Alessia Moretti, Ana Lourenco, Rosita Paladino, Yara Jaffe, Myriam Gitti, Tiziana Venturi, Marco Gullieuszik, Bianca Poggianti, Benedetta Vulcani, Andrea Biviano, Björn Adebahr, Ralf-Jürgen Dettmar

Abstract: We present deep JVLA observations at 1.4 GHz and 2.7 GHz (full polarization), as well as optical OmegaWINGS/WINGS and X-ray observations of two extended radio galaxies in the IIZW108 galaxy cluster at z = 0.04889. They show a bent tail morphology in agreement with a radio lobed galaxy falling into the cluster potential. Both galaxies are found to possess properties comparable with {narrow-angle} t… ▽ More We present deep JVLA observations at 1.4 GHz and 2.7 GHz (full polarization), as well as optical OmegaWINGS/WINGS and X-ray observations of two extended radio galaxies in the IIZW108 galaxy cluster at z = 0.04889. They show a bent tail morphology in agreement with a radio lobed galaxy falling into the cluster potential. Both galaxies are found to possess properties comparable with {narrow-angle} tail galaxies in the literature even though they are part of a low mass cluster. We find a spectral index steepening and an increase in fractional polarization through the galaxy jets and an ordered magnetic field component mostly aligned with the jet direction. This is likely caused by either shear due to the velocity difference of the intracluster medium and the jet fluid and/or magnetic draping of the intracluster medium across the galaxy jets. We find clear evidence that one source is showing two active galactic nuclei (AGN) outbursts from which we expect the AGN has never turned off completely. We show that pure standard electron cooling cannot explain the jet length. We demonstrate therefore that these galaxies can be used as a laboratory to study gentle re-acceleration of relativistic electrons in galaxy jets via transition from laminar to turbulent motion. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: accepted for publication in MNRAS

Gas flows in galaxy mergers: supersonic turbulence in bridges, accretion from the circumgalactic medium, and metallicity dilution

Authors: Martin Sparre, Joseph Whittingham, Mitali Damle, Maan H. Hani, Philipp Richter, Sara L. Ellison, Christoph Pfrommer, Mark Vogelsberger

Abstract: In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magneto-hydrodynamical (MHD) zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated… ▽ More In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magneto-hydrodynamical (MHD) zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated by turbulent pressure, with turbulent Mach numbers peaking at values of 1.6-3.3. This implies that bridges are dominated by supersonic turbulence, and are thus ideal candidates for studying the impact of extreme environments on star formation. We also find that gas accreted from the circumgalactic medium (CGM) during the merger significantly contributes (27-51 per cent) to the star formation rate (SFR) at the time of coalescence and drives the subsequent reignition of star formation in the merger remnant. Indeed, 19-53 per cent of the SFR at z=0 originates from gas belonging to the CGM prior the merger. Finally, we investigate the origin of the metallicity-diluted gas at the centre of merging galaxies. We show that this gas is rapidly accreted onto the galactic centre with a time-scale much shorter than that of normal star-forming galaxies. This explains why coalescing galaxies are not well-captured by the fundamental metallicity relation. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 15 pages, 13 figures, submitted to MNRAS

Spectrally resolved cosmic rays: II -- Momentum-dependent cosmic ray diffusion drives powerful galactic winds

Authors: Philipp Girichidis, Christoph Pfrommer, Rüdiger Pakmor, Volker Springel

Abstract: Recently, cosmic ray (CR) feedback has been identified as a critical process in galaxy formation but most previous simulations have integrated out the energy-dependence of the CR distribution, despite its large extent over more than twelve decades in particle energy. To improve upon this simplification, we present the implementation and first application of spectrally resolved CRs which are couple… ▽ More Recently, cosmic ray (CR) feedback has been identified as a critical process in galaxy formation but most previous simulations have integrated out the energy-dependence of the CR distribution, despite its large extent over more than twelve decades in particle energy. To improve upon this simplification, we present the implementation and first application of spectrally resolved CRs which are coupled to the magneto-hydrodynamics in simulations of galaxy formation. The spectral model for the CRs enables more accurate cooling of CRs and allows for an energy-dependent spatial diffusion, for which we introduce a new stable numerical algorithm that proves essential in highly dynamical systems. We perform galaxy formation simulations with this new model and compare the results to a grey CR approach with a simplified diffusive transport and effective cooling that assumes steady-state spectra. We find that the galaxies with spectrally resolved CRs differ in morphology, star formation rate, and strength and structure of the outflows. Interestingly, the first outflow front is driven by CRs with average momenta of $\sim200-600\,\mathrm{Gev}~c^{-1}$. The subsequent formation of outflows, which reach mass loading factors of order unity, are primarily launched by CRs of progressively smaller average momenta of $\sim8-15\,\mathrm{GeV}~c^{-1}$. The CR spectra in the galactic centre quickly approach a steady state, which does not significantly vary over time. In the outer disc and outflow regions, the spectral shape approaches steady state only after $\sim2\,\mathrm{Gyr}$ of evolution. Furthermore, the shapes of the approximate steady state spectra differ for individual regions of the galaxy, which highlights the importance of actively including the full CR spectrum. △ Less

Submitted 14 March, 2023; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: 21 pages, 13 figures, published in MNRAS, erratum added

Cosmic rays and non-thermal emission in simulated galaxies: III. probing cosmic ray calorimetry with radio spectra and the FIR-radio correlation

Authors: Maria Werhahn, Christoph Pfrommer, Philipp Girichidis

Abstract: An extinction-free estimator of the star-formation rate (SFR) of galaxies is critical for understanding the high-redshift universe. To this end, the nearly linear, tight correlation of far-infrared (FIR) and radio luminosity of star-forming galaxies is widely used. While the FIR is linked to massive star formation, which also generates shock-accelerated cosmic ray (CR) electrons and radio synchrot… ▽ More An extinction-free estimator of the star-formation rate (SFR) of galaxies is critical for understanding the high-redshift universe. To this end, the nearly linear, tight correlation of far-infrared (FIR) and radio luminosity of star-forming galaxies is widely used. While the FIR is linked to massive star formation, which also generates shock-accelerated cosmic ray (CR) electrons and radio synchrotron emission, a detailed understanding of the underlying physics is still lacking. Hence, we perform three-dimensional magneto-hydrodynamical (MHD) simulations of isolated galaxies over a broad range of halo masses and SFRs using the moving-mesh code AREPO, and evolve the CR proton energy density self-consistently. In post-processing, we calculate the steady-state spectra of primary, shock-accelerated and secondary CR electrons, which result from hadronic CR proton interactions with the interstellar medium. The resulting total radio luminosities correlate with the FIR luminosities as observed and are dominated by primary CR electrons if we account for anisotropic CR diffusion. The increasing contribution of secondary emission up to 30 per cent in starbursts is compensated by the larger bremsstrahlung and Coulomb losses. CR electrons are in the calorimetric limit and lose most of their energy through inverse Compton interactions with star-light and cosmic microwave background (CMB) photons while less energy is converted to synchrotron emission. This implies steep steady-state synchrotron spectra in starbursts. Interestingly, we find that thermal free-free emission hardens the total radio spectra at high radio frequencies and reconciles calorimetric theory with observations while free-free absorption explains the observed low-frequency flattening towards the central regions of starbursts. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: 24 pages, 15 figures, submitted to MNRAS, comments are welcome

Simulating radio synchrotron emission in star-forming galaxies: small-scale magnetic dynamo and the origin of the far infrared-radio correlation

Authors: Christoph Pfrommer, Maria Werhahn, Rüdiger Pakmor, Philipp Girichidis, Christine M. Simpson

Abstract: In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). We examine magneto-hydrodynamic galaxy simulations with cosmic rays (CRs), accounting for their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shoc… ▽ More In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). We examine magneto-hydrodynamic galaxy simulations with cosmic rays (CRs), accounting for their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shock, which injects turbulence and drives a small-scale magnetic dynamo. As the shock propagates outwards and the associated turbulence decays, the large velocity shear between the supersonically rotating cool disc with respect to the (partially) pressure-supported hot circumgalactic medium excites Kelvin-Helmholtz surface and body modes. Those inject turbulence and drive multiple small-scale dynamos, which exponentially amplify magnetic fields. They grow in scale to reach equipartition with thermal and CR energies in Milky Way-mass galaxies. In small galaxies, the magnetic energy saturates at the turbulent energy while it fails to reach equipartition with thermal and CR energies. We solve for steady-state spectra of CR protons, secondary electrons/positrons from hadronic CR-proton interactions with the interstellar medium, and primary shock-accelerated electrons at supernovae. The radio-synchrotron emission is dominated by primary electrons, irradiates the magnetised disc, bulge, and bubble-shaped magnetically-loaded outflows of our simulated Milky Way-mass galaxy. Our star-forming and star-bursting galaxies with saturated magnetic fields match the global FIR-radio correlation (FRC) across four orders of magnitude. Its intrinsic scatter arises due to (i) different magnetic saturation levels that result from different seed magnetic fields, (ii) different radio synchrotron luminosities for different specific SFRs at fixed SFR and (iii) a varying radio intensity with galactic inclination. (abridged) △ Less

Submitted 5 July, 2022; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 37 pages, 26 figures, accepted for publication in MNRAS

Cosmic rays and non-thermal emission in simulated galaxies. II. $γ$-ray maps, spectra and the far infrared-$γ$-ray relation

Authors: Maria Werhahn, Christoph Pfrommer, Philipp Girichidis, Georg Winner

Abstract: The $γ$-ray emission of star-forming (SF) galaxies is attributed to hadronic interactions of cosmic ray (CR) protons with the interstellar gas and contributions from CR electrons via bremsstrahlung and inverse Compton (IC) scattering. The relative importance of these processes in different galaxy types is still unclear. We model these processes in three-dimensional magneto-hydrodynamical (MHD) sim… ▽ More The $γ$-ray emission of star-forming (SF) galaxies is attributed to hadronic interactions of cosmic ray (CR) protons with the interstellar gas and contributions from CR electrons via bremsstrahlung and inverse Compton (IC) scattering. The relative importance of these processes in different galaxy types is still unclear. We model these processes in three-dimensional magneto-hydrodynamical (MHD) simulations of the formation of isolated galactic discs using the moving-mesh code AREPO, including dynamically coupled CR protons and adopting different CR transport models. We calculate steady-state CR spectra and also account for the emergence of secondary electrons and positrons. This allows us to produce detailed $γ$-ray maps, luminosities and spectra of our simulated galaxies at different evolutionary stages. Our simulations with anisotropic CR diffusion and a low CR injection efficiency at supernovae (SNe, $ζ_{\mathrm{SN}}=0.05$) can successfully reproduce the observed far infrared (FIR)-$γ$-ray relation. Starburst galaxies are close to the calorimetric limit, where CR protons lose most of their energy due to hadronic interactions and hence, their $γ$-ray emission is dominated by neutral pion decay. However, in low SF galaxies, the increasing diffusive losses soften the CR proton spectra due to energy-dependent diffusion, and likewise steepen the pionic $γ$-ray spectra. In turn, IC emission hardens the total spectra and can contribute up to $\sim40$ per cent of the total luminosity in low SF galaxies. Furthermore, in order to match the observed $γ$-ray spectra of starburst galaxies, we require a weaker energy dependence of the CR diffusion coefficient, $D\propto{E}^{0.3}$, in comparison to Milky Way-like galaxies. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Comments: 20 pages, 11 figures, accepted to MNRAS

Cosmic rays and non-thermal emission in simulated galaxies. I. Electron and proton spectra compared to Voyager-1 data

Authors: Maria Werhahn, Christoph Pfrommer, Philipp Girichidis, Ewald Puchwein, Rüdiger Pakmor

Abstract: Current-day cosmic ray (CR) propagation studies use static Milky-Way models and fit parametrized source distributions to data. Instead, we use three-dimensional magneto-hydrodynamical (MHD) simulations of isolated galaxies with the moving-mesh code AREPO that self-consistently accounts for hydrodynamic effects of CR protons. In post-processing, we calculate their steady-state spectra, taking into… ▽ More Current-day cosmic ray (CR) propagation studies use static Milky-Way models and fit parametrized source distributions to data. Instead, we use three-dimensional magneto-hydrodynamical (MHD) simulations of isolated galaxies with the moving-mesh code AREPO that self-consistently accounts for hydrodynamic effects of CR protons. In post-processing, we calculate their steady-state spectra, taking into account all relevant loss processes. We show that this steady-state assumption is well justified in the disc and generally for regions that emit non-thermal radio and gamma rays. Additionally, we model the spectra of primary electrons, accelerated by supernova remnants, and secondary electrons and positrons produced in hadronic CR proton interactions with the gas. We find that proton spectra above 10 GeV only weakly depend on galactic radius, while they acquire a radial dependence at lower energies due to Coulomb interactions. Radiative losses steepen the spectra of primary CR electrons in the central galactic regions while diffusive losses dominate in the outskirts. Secondary electrons exhibit a steeper spectrum than primaries because they originate from the transported steeper CR proton spectra. Consistent with Voyager-1 and AMS-02 data, our models (i) show a turn-over of proton spectra below GeV energies due to Coulomb interactions so that electrons start to dominate the total particle spectra and (ii) match the shape of the positron fraction up to 10 GeV. We conclude that our steady-state CR modeling in MHD-CR galaxy simulations is sufficiently realistic to capture the dominant transport effects shaping their spectra, arguing for a full MHD treatment to accurately model CR transport in the future. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: 23 pages, 17 figures, accepted to MNRAS

Comparing different closure relations for cosmic ray hydrodynamics

Authors: Timon Thomas, Christoph Pfrommer

Abstract: Cosmic ray (CR) hydrodynamics is a (re-)emerging field of high interest due to the importance of CRs for the dynamical evolution of the interstellar, the circumgalactic, and the intracluster medium. In these environments, CRs with GeV energies can influence large-scale dynamics by regulating star formation, driving galactic winds or by altering the pressure balance of galactic halos. Recent effort… ▽ More Cosmic ray (CR) hydrodynamics is a (re-)emerging field of high interest due to the importance of CRs for the dynamical evolution of the interstellar, the circumgalactic, and the intracluster medium. In these environments, CRs with GeV energies can influence large-scale dynamics by regulating star formation, driving galactic winds or by altering the pressure balance of galactic halos. Recent efforts have moved the focus of the community from a one-moment description of CR transport towards a two-moment model as this allows for a more accurate description of the microphysics impacting the CR population. Like all hydrodynamical theories, these two-moment methods require a closure relation for a consistent and closed set of evolution equations. The goal of this paper is to quantify the impact of different closure relations on the resulting solutions. To this end, we review the common P1 and M1 closure relations, derive a new four-moment H1 description for CR transport and describe how to incorporate CR scattering by Alfvén waves into these three hydrodynamical models. While there are significant differences in the transport properties of radiation in the P1 and M1 approximations in comparison to more accurate radiative transfer simulations using the discrete ordinates approximation, we only find small differences between the three hydrodynamical CR transport models in the free streaming limit when we neglect CR scattering. Most importantly, for realistic applications in the interstellar, circumgalactic or intracluster medium where CR scattering is frequent, these differences vanish and all presented hydrodynamical models produce the same results. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Comments: 13 pages, 3 figure; comments are welcome!

The challenge of simultaneously matching the observed diversity of chemical abundance patterns in cosmological hydrodynamical simulations

Authors: Tobias Buck, Jan Rybizki, Sven Buder, Aura Obreja, Andrea V. Macciò, Christoph Pfrommer, Matthias Steinmetz, Melissa Ness

Abstract: With the advent of large spectroscopic surveys the amount of high quality chemo-dynamical data in the Milky Way (MW) increased tremendously. Accurately and correctly capturing and explaining the detailed features in the high-quality observational data is notoriously difficult for state-of-the-art numerical models. In order to keep up with the quantity and quality of observational datasets, improve… ▽ More With the advent of large spectroscopic surveys the amount of high quality chemo-dynamical data in the Milky Way (MW) increased tremendously. Accurately and correctly capturing and explaining the detailed features in the high-quality observational data is notoriously difficult for state-of-the-art numerical models. In order to keep up with the quantity and quality of observational datasets, improved prescriptions for galactic chemical evolution need to be incorporated into the simulations. Here we present a new, flexible, time resolved chemical enrichment model for cosmological simulations. Our model allows to easily change a number of stellar physics parameters such as the shape of the initial mass function (IMF), stellar lifetimes, chemical yields or SN Ia delay times. We implement our model into the Gasoline2 code and perform a series of cosmological simulations varying a number of key parameters, foremost evaluating different stellar yield sets for massive stars from the literature. We find that total metallicity, total iron abundance and gas phase oxygen abundance are robust predictions from different yield sets and in agreement with observational relations. On the other hand, individual element abundances, especially $α$-elements show significant differences across different yield sets and none of our models can simultaneously match constraints on the dwarf and MW mass scale. This offers a unique way of observationally constraining model parameters. For MW mass galaxies we find for most yield tables tested in this work a bimodality in the $[α$/Fe] vs. [Fe/H] plane of rather low intrinsic scatter potentially in tension with the observed abundance scatter. △ Less

Submitted 5 March, 2021; originally announced March 2021.

Comments: main text 19 pages and 11 figures, 4 pages of appendix, 23 pages total, python code and data at https://github.com/TobiBu/chemical_enrichment.git. submitted to MNRAS, comments very welcome

MERGHERS Pilot: MeerKAT discovery of diffuse emission in nine massive Sunyaev-Zel'dovich-selected galaxy clusters from ACT

Authors: K. Knowles, D. S. Pillay, S. Amodeo, A. J. Baker, K. Basu, D. Crichton, F. de Gasperin, M. Devlin, C. Ferrari, M. Hilton, K. M. Huffenberger, J. P. Hughes, B. J. Koopman, K. Moodley, T. Mroczkowski, S. Naess, F. Nati, L. B. Newburgh, N. Oozeer, L. Page, B. Partridge, C. Pfrommer, M. Salatino, A. Schillaci, C. Sifón , et al. (4 additional authors not shown)

Abstract: The MeerKAT Exploration of Relics, Giant Halos, and Extragalactic Radio Sources (MERGHERS) survey is a planned project to study a large statistical sample of galaxy clusters with the MeerKAT observatory. Here we present the results of a 16--hour pilot project, observed in response to the 2019 MeerKAT Shared Risk proposal call, to test the feasibility of using MeerKAT for a large cluster study usin… ▽ More The MeerKAT Exploration of Relics, Giant Halos, and Extragalactic Radio Sources (MERGHERS) survey is a planned project to study a large statistical sample of galaxy clusters with the MeerKAT observatory. Here we present the results of a 16--hour pilot project, observed in response to the 2019 MeerKAT Shared Risk proposal call, to test the feasibility of using MeerKAT for a large cluster study using short (0.2--2.1\,hour) integration times. The pilot focuses on 1.28\,GHz observations of 13 massive, low-to-intermediate redshift ($0.22 < z < 0.65$) clusters from the Sunyaev-Zel'dovich-selected Atacama Cosmology Telescope (ACT) DR5 catalogue that show multiwavelength indications of dynamical disturbance. With a 70 per cent detection rate (9/13 clusters), this pilot study validates our proposed MERGHERS observing strategy and provides twelve detections of diffuse emission, eleven of them new, indicating the strength of MeerKAT for such types of studies. The detections (signal-to-noise ratio $\gtrsim6$) are summarised as follows: two systems host both relic(s) and a giant radio halo, five systems host radio halos, and two have candidate radio halos. Power values, $k$-corrected to 1.4 GHz assuming a fiducial spectral index of $α= -1.3 \pm 0.4$, are consistent with known radio halo and relic scaling relations. △ Less

Submitted 15 April, 2021; v1 submitted 30 December, 2020; originally announced December 2020.

Comments: 9 pages, 4 figures, 2 tables. Accepted to MNRAS

Non-Kolmogorov turbulence in multiphase intracluster medium driven by cold gas precipitation and AGN jets

Authors: C. Wang, M. Ruszkowski, C. Pfrommer, S. Peng Oh, H. -Y. K. Yang

Abstract: AGN feedback is responsible for maintaining plasma in global thermal balance in extended halos of elliptical galaxies and galaxy clusters. Local thermal instability in the hot gas leads to the formation of precipitating cold gas clouds that feed the central supermassive black holes, thus heating the hot gas and maintaining global thermal equilibrium. We perform three dimensional MHD simulations of… ▽ More AGN feedback is responsible for maintaining plasma in global thermal balance in extended halos of elliptical galaxies and galaxy clusters. Local thermal instability in the hot gas leads to the formation of precipitating cold gas clouds that feed the central supermassive black holes, thus heating the hot gas and maintaining global thermal equilibrium. We perform three dimensional MHD simulations of self-regulated AGN feedback in a Perseus-like galaxy cluster with the aim of understanding the impact of the feedback physics on the turbulence properties of the hot and cold phases of the ICM. We find that, in general, the cold phase velocity structure function (VSF) is steeper than the prediction from Kolmogorov's theory. We attribute the physical origin of the steeper slope of the cold phase VSF to the driving of turbulent motions primarily by the gravitational acceleration acting on the ballistic clouds. We demonstrate that, in the pure hydrodynamical case, the precipitating cold filaments may be the dominant agent driving turbulence in the hot ICM. The arguments in favor of this hypothesis are that: (i) the cold phase mass dominates over hot gas mass in the inner cool core; (ii) hot and cold gas velocities are spatially correlated; (iii) both the cold and hot phase velocity distributions are radially biased. We show that, in the MHD case, the turbulence in the ambient hot medium (excluding the jet cone regions) can also be driven by the AGN jets. The driving is then facilitated by enhanced coupling due to magnetic fields of the ambient gas and the AGN jets. In the MHD case, turbulence may thus be driven by a combination of AGN jet stirring and filament motions. We conclude that future observations, including those from high spatial and spectral resolution X-ray missions, may help to constrain self-regulated AGN feedback by quantifying the multi-temperature VSF in the ICM. △ Less

Submitted 20 December, 2020; originally announced December 2020.

Comments: 12pages, 15figures, submitted to MNRAS

Connecting turbulent velocities and magnetic fields in galaxy cluster simulations with active galactic nuclei jets

Authors: Kristian Ehlert, Rainer Weinberger, Christoph Pfrommer, Volker Springel

Abstract: The study of velocity fields of the hot gas in galaxy clusters can help to unravel details of microphysics on small-scales and to decipher the nature of feedback by active galactic nuclei (AGN). Likewise, magnetic fields as traced by Faraday rotation measurements (RMs) inform about their impact on gas dynamics as well as on cosmic ray production and transport. We investigate the inherent relations… ▽ More The study of velocity fields of the hot gas in galaxy clusters can help to unravel details of microphysics on small-scales and to decipher the nature of feedback by active galactic nuclei (AGN). Likewise, magnetic fields as traced by Faraday rotation measurements (RMs) inform about their impact on gas dynamics as well as on cosmic ray production and transport. We investigate the inherent relationship between large-scale gas kinematics and magnetic fields through non-radiative magnetohydrodynamical simulations of the creation, evolution and disruption of AGN jet-inflated lobes in an isolated Perseus-like galaxy cluster, with and without pre-existing turbulence. In particular, we connect cluster velocity measurements with mock RM maps to highlight their underlying physical connection, which opens up the possibility of comparing turbulence levels in two different observables. For single jet outbursts, we find only a local impact on the velocity field, i.e. the associated increase in velocity dispersion is not volume-filling. Furthermore, in a setup with pre-existing turbulence, this increase in velocity dispersion is largely hidden. We use mock X-ray observations to show that at arcmin resolution, the velocity dispersion is therefore dominated by existing large-scale turbulence and is only minimally altered by the presence of a jet. For the velocity structure of central gas uplifted by buoyantly rising lobes, we find fast, coherent outflows with low velocity dispersion. Our results highlight that projected velocity distributions show complex structures which pose challenges for the interpretation of observations. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 19 pages, 9 figures, submitted to MNRS. Comments welcome!

The impact of magnetic fields on cosmological galaxy mergers. I: Reshaping gas and stellar discs

Authors: Joseph Whittingham, Martin Sparre, Christoph Pfrommer, Rüdiger Pakmor

Abstract: Mergers play an important role in galaxy evolution. In particular, major mergers are able to have a transformative effect on galaxy morphology. In this paper, we investigate the role of magnetic fields in gas-rich major mergers. To this end, we run a series of high-resolution magnetohydrodynamic (MHD) zoom-in simulations with the moving-mesh code Arepo and compare the outcome with hydrodynamic sim… ▽ More Mergers play an important role in galaxy evolution. In particular, major mergers are able to have a transformative effect on galaxy morphology. In this paper, we investigate the role of magnetic fields in gas-rich major mergers. To this end, we run a series of high-resolution magnetohydrodynamic (MHD) zoom-in simulations with the moving-mesh code Arepo and compare the outcome with hydrodynamic simulations run from the same initial conditions. This is the first time that the effect of magnetic fields in major mergers has been investigated in a cosmologically-consistent manner. In contrast to previous non-cosmological simulations, we find that the inclusion of magnetic fields has a substantial impact on the production of the merger remnant. Whilst magnetic fields do not strongly affect global properties, such as the star formation history, they are able to significantly influence structural properties. Indeed, MHD simulations consistently form remnants with extended discs and well-developed spiral structure, whilst hydrodynamic simulations form more compact remnants that display distinctive ring morphology. We support this work with a resolution study and show that whilst global properties are broadly converged across resolution and physics models, morphological differences only develop given sufficient resolution. We argue that this is due to the more efficient excitement of a small-scale dynamo in higher resolution simulations, resulting in a more strongly amplified field that is better able to influence gas dynamics. △ Less

Submitted 21 March, 2022; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 27 pages, 18 figures. Published in MNRAS on 19 May 2021

A finite volume method for two-moment cosmic-ray hydrodynamics on a moving mesh

Authors: Timon Thomas, Christoph Pfrommer, Rüdiger Pakmor

Abstract: We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh Arepo code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretised using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of Arepo. T… ▽ More We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh Arepo code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretised using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of Arepo. The interaction of CRs and gyroresonant Alfvén waves is described by short-timescale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magneto-hydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetised discontinuity, the hydrodynamic response to a CR overpressure, CR acceleration of a warm cloud, and a CR blast wave, which demonstrate that the coupling between CR and magneto-hydrodynamics is robust and accurate. We demonstrate the numerical convergence of the presented scheme using new linear and non-linear analytic solutions. △ Less

Submitted 22 October, 2020; originally announced October 2020.

Comments: 24 pages, 15 figures, submitted to MNRAS, comments are welcome!

A new cosmic ray-driven instability

Authors: Mohamad Shalaby, Timon Thomas, Christoph Pfrommer

Abstract: Cosmic ray (CR)-driven instabilities play a decisive role during particle acceleration at shocks and CR propagation in galaxies and galaxy clusters. These instabilities amplify magnetic fields and modulate CR transport so that the intrinsically collisionless CR population is tightly coupled to the thermal plasma and provides dynamical feedback. Here, we show that CRs with a finite pitch angle driv… ▽ More Cosmic ray (CR)-driven instabilities play a decisive role during particle acceleration at shocks and CR propagation in galaxies and galaxy clusters. These instabilities amplify magnetic fields and modulate CR transport so that the intrinsically collisionless CR population is tightly coupled to the thermal plasma and provides dynamical feedback. Here, we show that CRs with a finite pitch angle drive electromagnetic waves (along the background magnetic field) unstable on intermediate scales between the gyro-radii of CR ions and electrons as long as CRs are drifting with a velocity less than half of the Alfvén speed of electrons. By solving the linear dispersion relation, we show that this new instability typically grows faster by more than an order of magnitude in comparison to the commonly discussed resonant instability at the ion gyroscale. We find the growth rate for this intermediate-scale instability and identify the growing modes as background ion-cyclotron modes in the frame that is comoving with the CRs. We confirm the theoretical growth rate with a particle-in-cell (PIC) simulation and study the non-linear saturation of this instability. We identify three important astro-physical applications of this intermediate-scale instability, which is expected to 1. modulate CR transport and strengthen CR feedback in galaxies and galaxy clusters, 2. enable electron injection into the diffusive shock acceleration process, and 3. decelerate CR escape from the sites of particle acceleration which would generate gamma-ray halos surrounding CR sources such as supernova remnants. △ Less

Submitted 3 December, 2020; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 16 pages, 13 figures, 2 tables. Accepted to ApJ. Comments are welcome. An animation for the growth of the new instability can be found at http://www.youtube.com/watch?v=yP2dBCPOQYc

Highly ordered magnetic fields in the tail of the jellyfish galaxy JO206

Authors: Ancla Müller, Bianca Poggianti, Christoph Pfrommer, Björn Adebahr, Paolo Serra, Alessandro Ignesti, Martin Sparre, Myriam Gitti, Ralf-Jürgen Dettmar, Benedetta Vulcani, Alessia Moretti

Abstract: Jellyfish galaxies have long tails of gas that is stripped from the disc by ram pressure due to the motion of galaxies in the intracluster medium in galaxy clusters. We present the first measurement of the magnetic field strength and orientation within the disc and the (90$\,$kpc-long) $\rm Hα$-emitting tail of the jellyfish galaxy JO206. The tail has a large-scale magnetic field ($>4.1\,μ$G), a s… ▽ More Jellyfish galaxies have long tails of gas that is stripped from the disc by ram pressure due to the motion of galaxies in the intracluster medium in galaxy clusters. We present the first measurement of the magnetic field strength and orientation within the disc and the (90$\,$kpc-long) $\rm Hα$-emitting tail of the jellyfish galaxy JO206. The tail has a large-scale magnetic field ($>4.1\,μ$G), a steep radio spectral index ($α\sim -2.0$), indicating an aging of the electrons propagating away from the star-forming regions, and extremely high fractional polarisation ($>50\,$%), indicating low turbulent motions. The magnetic field vectors are aligned with (parallel to) the direction of the ionised-gas tail and stripping direction. High-resolution simulations of a large, cold gas cloud that is exposed to a hot, magnetised turbulent wind show that the high fractional polarisation and the ordered magnetic field can be explained by accretion of draped magnetised plasma from the hot wind that condenses onto the external layers of the tail, where it is adiabatically compressed and sheared. The ordered magnetic field, preventing heat and momentum exchange, may be a key factor in allowing in-situ star formation in the tail. △ Less

Submitted 26 October, 2020; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: accepted for publication in Nature Astronomy, Sharedlt: https://rdcu.be/b9c0L (free read only)

Simulating TeV gamma-ray morphologies of shell-type supernova remnants

Authors: Matteo Pais, Christoph Pfrommer

Abstract: Supernova remnant (SNR) shocks provide favourable sites of cosmic ray (CR) proton acceleration if the local magnetic field direction is quasi-parallel to the shock normal. Using the moving-mesh magneto-hydrodynamical (MHD) code AREPO we present a suite of SNR simulations with CR acceleration in the Sedov-Taylor phase that combine different magnetic field topologies, density distributions with grad… ▽ More Supernova remnant (SNR) shocks provide favourable sites of cosmic ray (CR) proton acceleration if the local magnetic field direction is quasi-parallel to the shock normal. Using the moving-mesh magneto-hydrodynamical (MHD) code AREPO we present a suite of SNR simulations with CR acceleration in the Sedov-Taylor phase that combine different magnetic field topologies, density distributions with gradients and large-scale fluctuations, and -- for our core-collapse SNRs -- a multi-phase interstellar medium with dense clumps with a contrast of $10^4$. Assuming the hadronic gamma-ray emission model for the TeV gamma-ray emission, we find that large-amplitude density fluctuations of $δρ/ρ_0\gtrsim75$ per cent are required to strongly modulate the gamma-ray emissivity in a straw man's model in which the acceleration efficiency is independent of magnetic obliquity. However, this causes strong corrugations of the shock surface that are ruled out by gamma-ray observations. By contrast, magnetic obliquity-dependent acceleration can easily explain the observed variance in gamma-ray morphologies ranging from SN1006 (with a homogeneous magnetic field) to Vela Junior and RX J1713 (with a turbulent field) in a single model that derives from plasma particle-in-cell simulations. Our best-fit model for SN1006 has a large-scale density gradient of $ \nabla{n}\simeq 0.0034~\mathrm{cm}^{-3}~\mathrm{pc}^{-1} $ pointing from south-west to north-east and a magnetic inclination with the plane of the sky of $\lesssim10^\circ$. Our best-fit model for Vela Junior and RX J1713 adopts a combination of turbulent magnetic field and dense clumps to explain their TeV gamma-ray morphologies and moderate shock corrugations. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 20 pages, 12 figures, accepted for publication on MNRAS

Interaction of a cold cloud with a hot wind: the regimes of cloud growth and destruction and the impact of magnetic fields

Authors: Martin Sparre, Christoph Pfrommer, Kristian Ehlert

Abstract: Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix… ▽ More Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix and dissolve in the hot wind. On the contrary, stripped cold gas from a large cloud mixes with the hot wind to intermediate temperatures, and then becomes thermally unstable and cools, causing a net accretion of hot gas to the cold tail. Using the magneto-hydrodynamical code AREPO, we perform cloud crushing simulations and test analytical criteria for the transition between the growth and destruction regimes to clarify a current debate in the literature. We find that the hot-wind cooling time sets the transition radius and not the cooling time of the mixed phase. Magnetic fields modify the wind-cloud interaction. Draping of wind magnetic field enhances the field upstream of the cloud and fluid instabilities are suppressed by a turbulently magnetised wind beyond what is seen for a wind with a uniform magnetic field. We furthermore predict jellyfish galaxies to have ordered magnetic fields aligned with their tails. We finally discuss how the results of idealised simulations can be used to provide input to subgrid models in cosmological (magneto-)hydrodynamical simulations, which cannot resolve the detailed small-scale structure of cold gas clouds in the circum-galactic medium. △ Less

Submitted 12 October, 2020; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: Accepted for publication in MNRAS

The Hestia project: simulations of the Local Group

Authors: Noam I. Libeskind, Edoardo Carlesi, Rob J. J. Grand, Arman Khalatyan, Alexander Knebe, Ruediger Pakmor, Sergey Pilipenko, Marcel S. Pawlowski, Martin Sparre, Elmo Tempel, Peng Wang, Helene M. Courtois, Stefan Gottloeber, Yehuda Hoffman, Ivan Minchev, Christoph Pfrommer, Jenny G. Sorce, Volker Springel, Matthias Steinmetz, R. Brent Tully, Mark Vogelsberger, Gustavo Yepes

Abstract: We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter on… ▽ More We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter only simulations, and selected for higher resolution magneto hydrodynamic simulation using the Arepo code. Baryonic physics follows the Auriga model of galaxy formation. The simulations contain a high resolution region of 3-5 Mpc in radius from the Local Group midpoint embedded in the correct cosmographic landscape. Within this region a simulated Local Group consisting of a Milky Way and Andromeda like galaxy forms, whose description is in excellent agreement with observations. The simulated Local Group galaxies resemble the Milky Way and Andromeda in terms of their halo mass, mass ratio, stellar disc mass, morphology separation, relative velocity, rotation curves, bulge-disc morphology, satellite galaxy stellar mass function, satellite radial distribution and in some cases, the presence of a Magellanic cloud like object. Because these simulations properly model the Local Group in their cosmographic context, they provide a testing ground for questions where environment is thought to play an important role. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: 20 pages, 13 figures, submitted to MNRAS

Whistler-regulated MHD: Transport equations for electron thermal conduction in the high $β$ intracluster medium of galaxy clusters

Authors: J. F. Drake, C. Pfrommer, C. S. Reynolds, M. Ruszkowski, M. Swisdak, A. Einarsson, T. Thomas, A. B. Hassam, G. T. Roberg-Clark

Abstract: Transport equations for electron thermal energy in the high $β_e$ intracluster medium (ICM) are developed that include scattering from both classical collisions and self-generated whistler waves. The calculation employs an expansion of the kinetic electron equation along the ambient magnetic field in the limit of strong scattering and assumes whistler waves with low phase speeds… ▽ More Transport equations for electron thermal energy in the high $β_e$ intracluster medium (ICM) are developed that include scattering from both classical collisions and self-generated whistler waves. The calculation employs an expansion of the kinetic electron equation along the ambient magnetic field in the limit of strong scattering and assumes whistler waves with low phase speeds $V_w\sim{v}_{te}/β_e\ll{v}_{te}$ dominate the turbulent spectrum, with $v_{te}$ the electron thermal speed and $β_e\gg1$ the ratio of electron thermal to magnetic pressure. We find: (1) temperature-gradient-driven whistlers dominate classical scattering when $L_c>L/β_e$, with $L_c$ the classical electron mean-free-path and $L$ the electron temperature scale length, and (2) in the whistler dominated regime the electron thermal flux is controlled by both advection at $V_w$ and a comparable diffusive term. The findings suggest whistlers limit electron heat flux over large regions of the ICM, including locations unstable to isobaric condensation. Consequences include: (1) the Field length decreases, extending the domain of thermal instability to smaller length-scales, (2) the heat flux temperature dependence changes from $T_e^{7/2}/L$ to $V_wnT_e\sim{T}_e^{1/2}$, (3) the magneto-thermal and heat-flux driven buoyancy instabilities are impaired or completely inhibited, and (4) sound waves in the ICM propagate greater distances, as inferred from observations. This description of thermal transport can be used in macroscale ICM models. △ Less

Submitted 27 August, 2021; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: Accepted for publication in The Astrophysical Journal

Suppressed heat conductivity in the intracluster medium: implications for the magneto-thermal instability

Authors: Thomas Berlok, Eliot Quataert, Martin E. Pessah, Christoph Pfrommer

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

Submitted 18 March, 2021; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 20 pages, 18 figures, accepted for publication in MNRAS

Evolution and observational signatures of the cosmic ray electron spectrum in SN 1006

Authors: Georg Winner, Christoph Pfrommer, Philipp Girichidis, Maria Werhahn, Matteo Pais

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

Submitted 5 December, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 18 pages, 18 figures, published in MNRAS

Journal ref: MNRAS, 499, 2785 (2020)