Search

This study explores the dynamical impact of cosmic rays (CRs) in Milky Way-like galaxies using the Rhea simulation suite. Cosmic rays, with their substantial energy density, influence the interstellar medium (ISM) by supporting galactic winds, modulating star formation, and shaping ISM energetics. The simulations incorporate a multi-phase ISM, self-consistent CR transport in the advection-diffusion approximation, and interactions with magnetic fields to study their effect on galaxy evolution. Key findings reveal that CRs reduce star formation rates, and drive weak but sustained outflows with mass loading factors of $\sim0.2$, transporting a substantial fraction (20%-60%) of the injected CR energy. These CR-driven outflows are launched not just from the galactic center but across the entire disk, illustrating their pervasive dynamical influence. Galactic disks supported by CRs exhibit broader vertical structures compared to magnetic-field-dominated setups, though the scale heights are similar. CR feedback enhances magnetic flux transport to the circumgalactic medium (CGM), leading to a magnetically enriched CGM with field strengths of $\sim0.5\mu\mathrm{G}$ while reducing gas temperatures to $\lesssim10^5\,\mathrm{K}$. The CR energy is relatively smoothly distributed in the disk, with gradient lengths exceeding the typical size of molecular clouds, indicating that the CR behavior is not adiabatic., Comment: 23 pages, 17 figures, submitted to A&A. See additional Rhea companion paper G\"oller et al. today on astro-ph

The circumgalactic medium (CGM) plays a critical role in galaxy evolution, influencing gas flows, feedback processes, and galactic dynamics. Observations show a substantial cold gas reservoir in the CGM, but the mechanisms driving its formation and evolution remain unclear. Cosmic rays (CRs), as a source of non-thermal pressure, are increasingly recognized as key regulators of cold gas dynamics. This study explores how CRs affect cold clouds that condense from the hot CGM via thermal instability (TI). Using 3D CR-magnetohydrodynamic (CRMHD) simulations with AREPO, we assess the impact of various CR transport models on cold gas evolution. Under purely advective CR transport, CR pressure significantly suppresses the collapse of thermally unstable regions, altering the CGM's structure. In contrast, realistic CR transport models reveal that CRs escape collapsing regions via streaming and diffusion along magnetic fields, diminishing their influence on the thermal and dynamic structure of the cold CGM. The ratio of the CR escape time to the cloud collapse time emerges as a critical factor in determining the impact of CRs on TI. CRs remain confined within cold clouds when effective CR diffusion is slow which maximizes their pressure support and inhibits collapse. Fast effective CR diffusion, as realized in our 2-moment CRMHD model, facilitates rapid CR escape, reducing their stabilizing effect. This realistic CR transport model shows a wide dynamic range of the effective CR diffusion coefficient, ranging from $10^{29}$ to $10^{30}\,\mathrm{cm^{2}\,s^{-1}}$ for thermally- to CR-dominated atmospheres, respectively. In addition to these CR transport-related effects, we demonstrate that high numerical resolution is crucial to avoid spuriously large clouds formed in low-resolution simulations, which would result in overly long CR escape times and artificially amplified CR pressure support., Comment: 20 pages and 17 figures in main text; submitted to A&A

It is generally accepted that radio relics are the result of synchrotron emission from shock-accelerated electrons. Current models, however, are still unable to explain several aspects of their formation. In this paper, we focus on three outstanding problems: i) Mach number estimates derived from radio data do not agree with those derived from X-ray data, ii) cooling length arguments imply a magnetic field that is at least an order of magnitude larger than the surrounding intracluster medium (ICM), and iii) spectral index variations do not agree with standard cooling models. To solve these problems, we first identify typical shock conditions in cosmological simulations, using the results to inform significantly higher resolution shock-tube simulations. We apply the cosmic ray electron spectra code CREST and the emission code CRAYON+ to these, thereby generating mock observables ab-initio. We identify that upon running into an accretion shock, merger shocks generate a shock-compressed sheet, which, in turn, runs into upstream density fluctuations in pressure equilibrium. This mechanism directly gives rise to solutions to the three problems: it creates a distribution of Mach numbers at the shock-front, which flattens cosmic ray electron spectra, thereby biasing radio-derived Mach number estimates to higher values. We show that this effect is particularly strong in weaker shocks. Secondly, the density sheet becomes Rayleigh-Taylor unstable at the contact discontinuity, causing turbulence and additional compression downstream. This amplifies the magnetic field from ICM-like conditions up to microgauss levels. We show that synchrotron-based measurements are strongly biased by the tail of the distribution here too. Finally, the same instability also breaks the common assumption that matter is advected at the post-shock velocity downstream, thus invalidating laminar-flow based cooling models., Comment: 28 pages, 19 figures. Submitted to A&A. The abstract has been abridged to fit within the allowed character limit

Observations of Faraday rotation and synchrotron emission in galaxy clusters imply large-scale magnetic fields with $\mu\mathrm{G}$ strengths possibly extending back to $z=4$. Non-radiative cosmological simulations of galaxy clusters show a comparably slow magnetic field growth that only saturates at late times. We include galaxy formation physics and find a significantly accelerated magnetic field growth. We identify three crucial stages in the magnetic field evolution. 1) At high redshift, the central dominant galaxy serves as the prime agent that magnetizes not only its immediate vicinity but also most of the forming protocluster through a combination of a small-scale dynamo induced by gravitationally driven, compressive turbulence and stellar and active galactic nuclei feedback that distribute the magnetic field via outflows. 2) This process continues as other galaxies merge into the forming cluster in subsequent epochs, thereby transporting their previously amplified magnetic field to the intracluster medium (ICM) through ram pressure stripping and galactic winds. 3) At lower redshift, gas accretion and frequent cluster mergers trigger additional small-scale dynamo processes, thereby preventing the decay of the magnetic field and fostering the increase of the magnetic coherence scale. We show that the magnetic field observed today in the weakly collisional ICM is consistently amplified on collisional scales. Initially, this occurs in the collisional interstellar medium during protocluster assembly, and later in the ICM on the magnetic coherence scale, which always exceeds the particle mean free path, supporting the use of magneto-hydrodynamics for studying the cluster dynamo. We generate synthetic Faraday rotation measure observations of protoclusters, highlighting the potential for studying magnetic field growth during the onset of cluster formation at cosmic dawn., Comment: 23 pages, 15 figures, submitted to A&A

Magnetic fields and turbulence are fundamental to the evolution of galaxies, yet their precise measurement and analysis present significant challenges. The recently developed Velocity Gradient Technique (VGT), which capitalizes on the anisotropy inherent in magnetohydrodynamic (MHD) turbulence, represents a new method for mapping magnetic fields in galaxies using spectroscopic observations. Most validations of VGT thus far, however, have relied upon idealized MHD turbulence simulations, which lack the more complex dynamics found in galaxies and galaxy mergers. In this study, we scrutinize VGT using an AREPO-based cosmological galaxy merger simulation, testing its effectiveness across pre-merger, merging, and post-merger stages. We examine the underlying assumptions of VGT and probe the statistics of gas density, velocity, and magnetic fields over time. We find that velocity fluctuations are indeed anisotropic at each stage, being larger in the direction perpendicular to the local magnetic field, as required by VGT. We find, additionally, that galaxy mergers substantially intensify velocity and density fluctuations and amplify magnetic fields at all scales. The observed scaling behavior of the velocity fluctuations corresponds to $r^{1/2}$ up to 0.4~kpc, shifting to a steeper trend between 0.6 and 3~kpc, and to a shallower trend thereafter. The scaling of the magnetic field and density fluctuations at scales $\lesssim$ 1.0 kpc also predominantly aligns with $r^{1/2}$. Finally, we compare results from VGT to those derived from polarization-based magnetic field measurements, finding consistent and statistically significant global agreement in all cases. This opens the way to applying VGT to external galaxies., Comment: 19 pages, 10 figures, submitted to ApJ

In weakly collisional, strongly magnetised plasmas such as the intracluster medium (ICM), hot accretion flows and the solar corona, the transport of heat and momentum occurs primarily along magnetic field lines. In this paper we present a new scheme for modelling anisotropic thermal conduction which we have implemented in the moving mesh code AREPO. Our implementation uses a semi-implicit time integration scheme which works accurately and efficiently with individual timestepping, making the scheme highly suitable for use in cosmological simulations. We apply the scheme to a number of test-problems including the diffusion of a hot patch of gas in a circular magnetic field, the progression of a point explosion in the presence of thermal conduction, and the evolution and saturation of buoyancy instabilities in anisotropically conducting plasmas. We use these idealised tests to demonstrate the accuracy and stability of the solver and highlight the ways in which anisotropic conduction can fundamentally change the behaviour of the system. Finally, we demonstrate the solver's capability when applied to highly non-linear problems with deep timestep hierarchies by performing high-resolution cosmological zoom-in simulations of a galaxy cluster with conduction. We show that anisotropic thermal conduction can have a significant impact on the temperature distribution of the ICM and that whistler suppression may be relevant on cluster scales. The new scheme is, therefore, well suited for future work which will explore the role of anisotropic thermal conduction in a range of astrophysical contexts including the ICM of clusters and the circumgalactic medium of galaxies., Comment: 20 pages, 14 figures, submitted to MNRAS, comments welcome

Feedback processes in galaxies dictate their structure and evolution. Baryons can be cycled through stars, which inject energy into the interstellar medium (ISM) in supernova explosions, fueling multiphase galactic winds. Cosmic rays (CRs) accelerated at supernova remnants are an important component of feedback. CRs can effectively contribute to wind driving; however, their impact heavily depends on the assumed CR transport model. We run high-resolution "tallbox" simulations of a patch of a galactic disk using the moving mesh magnetohydrodynamics code AREPO, including varied CR implementations and the CRISP non-equilibrium thermochemistry model. We characterize the impact of CR feedback on star formation and multiphase outflows. While CR-driven winds are able to supply energy to a global-scale wind, a purely thermal wind loses most of its energy by the time it reaches 3 kpc above the disk midplane. We further find that the adopted CR transport model significantly affects the steady-state of the wind. In the model with CR advection, streaming, diffusion, and nonlinear Landau damping, CRs provide very strong feedback. Additionally accounting for ion-neutral damping (IND) decouples CRs from the cold ISM, which reduces the impact of CRs on the star formation rate. Nevertheless, CRs in this most realistic model are able to accelerate warm gas and levitate cool gas in the wind but have little effect on cold gas and hot gas. This model displays moderate mass loading and significant CR energy loading, demonstrating that IND does not prevent CRs from providing effective feedback., Comment: 27 pages, 13 figures. Submitted to ApJ

The formation of galaxies is significantly influenced by galactic winds, possibly driven by cosmic rays due to their long cooling times and better coupling to plasma compared to radiation. In this study, we compare the radio observations of the edge-on galaxy NGC 4217 from the CHANG-ES collaboration catalog with a mock observation of an isolated galaxy based on the arepo simulation that adopts the state-of-the-art two-moment cosmic ray transport treatment and multiphase interstellar medium model. We find significant agreement between the simulated and observed images and spectroscopic data for reasonable model parameters. Specifically, we find that (i) the shape of the intensity profiles depends weakly on the magnitude of the magnetic field, the distance of the simulated galaxy, and the normalization of the CR electron spectrum. The agreement between the mock and actual observations is degenerate with respect to these factors; (ii) the multi-wavelength spectrum above 0.1 GHz is in agreement with the radio observations and its slope is also only weakly sensitive to the magnetic field strength; (iii) the magnetic field direction exhibits X-shaped morphology, often seen in edge-on galaxies, which is consistent with the observations and indicates the presence of a galactic-scale outflow. Our results highlight the importance of incorporating advanced cosmic ray transport models in simulations and provide a deeper understanding of galactic wind dynamics and its impact on galaxy evolution., Comment: 21 pages. Accepted for publication in ApJ

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\'en 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., Comment: Comments are welcome

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., Comment: 16 pages, 9 figures, submitted to A&A, comments welcome

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., Comment: 9 pages, 5 figures, accepted for publication in A&A (Sect. 4. Extragalactic astronomy)

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-$\beta$ 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., Comment: 27 pages, 17 figures. Final version published in A&A

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., Comment: 14 pages, 7 figures, Accepted for publication in MNRAS

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 magnetothermal 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., Comment: 7 pages, 4 figures. Final version published in A&A

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., Comment: 19 pages, 12 figures, accepted for publication by MNRAS, comments welcome

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., Comment: invited A&ARv review; revised version; accepted for publication; 238 pages

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., Comment: 13 pages + references, 4 figures. To appear in JPP Letters

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\'{e}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., Comment: 28 pages, MNRAS, accepted

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., Comment: 22 pages, 12 figures. Published in MNRAS on 07 September 2023

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., Comment: 35 pages, 11 figures, submitted to JPP. Comments are welcome

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., Comment: 18 pages, 11 figures, submitted to MNRAS, comments are welcome!

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., Comment: 22 pages, 20 figures, accepted for publication in MNRAS, comments welcome

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., Comment: 18 pages, 14 figures, accepted for publication by MNRAS, comments welcome

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., Comment: 15 pages, 12 figures, 2 tables, submitted to ApJ

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., Comment: 13 pages, 12 figures

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<\rho\right>_V$. Assuming that the large-scale ISM density is statistically isotropic, we derive a relation between the three-dimensional clumping factor, $\mathcal{C}_\rho$, and the clumping factor of the $4\pi$ column density distribution on the cloud surface, $\mathcal{C}_\Sigma$, and find $\mathcal{C}_\Sigma=\mathcal{C}_\rho^{2/3}$. Applying our model to calculate the covering fraction, i.e., the $4\pi$ 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., Comment: 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

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., Comment: 27 pages, 21 Figures, submitted to MNRAS, comments welcome

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., Comment: 26 pages, 19 figures, published in MNRS

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\'en 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\'en-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\'en 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., Comment: 19 pages, 12 figures, submitted to MNRAS, comments welcome

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., Comment: 15 pages, 11 figures, 3 Tables. Comments are welcome. Updated based on referee's comments. Accepted to MNRAS on July 29, 2022

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., Comment: 12 pages, 8 figures, 1 table, Accepted to ApJ

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., Comment: 22 pages, 9 figures, revised version, Accepted in MNRAS

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 $\sigma$ away from $M_{\rm{dyn},500}$), showing that J0019 is a merging cluster and probably in a post-merging phase., Comment: 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"

Radio relics are arc-like synchrotron sources at the periphery of galaxy clusters, produced by cosmic-ray electrons in a $\mu$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\'en 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., Comment: 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"

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., Comment: accepted for publication in MNRAS

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., Comment: 15 pages, 13 figures, submitted to MNRAS

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., Comment: 21 pages, 13 figures, published in MNRAS, erratum added

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., Comment: 24 pages, 15 figures, submitted to MNRAS, comments are welcome

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), Comment: 37 pages, 26 figures, accepted for publication in MNRAS

The $\gamma$-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 $\gamma$-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, $\zeta_{\mathrm{SN}}=0.05$) can successfully reproduce the observed far infrared (FIR)-$\gamma$-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 $\gamma$-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 $\gamma$-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 $\gamma$-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., Comment: 20 pages, 11 figures, accepted to MNRAS

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., Comment: 23 pages, 17 figures, accepted to MNRAS

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\'en 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., Comment: 13 pages, 3 figure; comments are welcome!

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 $\alpha$-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 $[\alpha$/Fe] vs. [Fe/H] plane of rather low intrinsic scatter potentially in tension with the observed abundance scatter., Comment: 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

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., Comment: 19 pages, 9 figures, submitted to MNRS. Comments welcome!

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., Comment: 27 pages, 18 figures. Published in MNRAS on 19 May 2021

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\'en 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., Comment: 24 pages, 15 figures, submitted to MNRAS, comments are welcome!

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\'en 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., Comment: 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

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\alpha$-emitting tail of the jellyfish galaxy JO206. The tail has a large-scale magnetic field ($>4.1\,\mu$G), a steep radio spectral index ($\alpha \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., Comment: accepted for publication in Nature Astronomy, Sharedlt: https://rdcu.be/b9c0L (free read only)

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 $\delta\rho/\rho_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., Comment: 20 pages, 12 figures, accepted for publication on MNRAS

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., Comment: Accepted for publication in MNRAS