ABSTRACT
Non-thermal emission from relativistic cosmic ray (CR) electrons gives insight into the strength and morphology of intra-cluster magnetic fields, as well as providing powerful tracers of ...structure formation shocks. Emission caused by CR protons on the other hand still challenges current observations and is therefore testing models of proton acceleration at intra-cluster shocks. Large-scale simulations including the effects of CRs have been difficult to achieve and have been mainly reduced to simulating an overall energy budget, or tracing CR populations in post-processing of simulation output and has often been done for either protons or electrons. We introduce crescendo: Cosmic Ray Evolution with SpeCtral Electrons aND prOtons, an efficient on-the-fly Fokker–Planck solver to evolve distributions of CR protons and electrons within every resolution element of our simulation. The solver accounts for CR (re-)acceleration at intra-cluster shocks, based on results of recent particle-in-cell simulations, adiabatic changes, and radiative losses of electrons. We show its performance in test cases as well as idealized galaxy cluster (GC) simulations. We apply the model to an idealized GC merger following best-fitting parameters for CIZA J2242.4 + 5301-1 and study CR injection, radio relic morphology, spectral steepening, and synchrotron emission.
Abstract
We investigate shock structures driven by merger events in high-resolution simulations that result in a galaxy with a virial mass
M
≈ 10
12
M
⊙
. We find that the sizes and morphologies of ...the internal shocks resemble remarkably well those of the newly detected class of odd radio circles (ORCs). This would highlight a so-far overlooked mechanism to form radio rings, shells, and even more complex structures around elliptical galaxies. Mach numbers of
= 2–3 for such internal shocks are in agreement with the spectral indices of the observed ORCs. We estimate that ∼5% of galaxies could undergo merger events, which occasionally lead to such prominent structures within the galactic halo during their lifetime, explaining the low number of observed ORCs. At the time when the shock structures are matching the physical sizes of the observed ORCs, the central galaxies are typically classified as early-type galaxies, with no ongoing star formation, in agreement with observational findings. Although the energy released by such mergers could potentially power the observed radio luminosity already in Milky Way–like halos, our predicted luminosity from a simple, direct shock acceleration model is much smaller than the observed one. Considering the estimated number of candidates from our cosmological simulations and the higher observed energies, we suggest that the proposed scenario is more likely for halo masses around 10
13
M
⊙
in agreement with the observed stellar masses of the galaxies at the center of ORCs. Such shocks might be detectable with next-generation X-ray instruments like the Line Emission Mapper (LEM).
Abstract
Radio relics are typically found to be arc-like regions of synchrotron emission in the outskirts of merging galaxy clusters, bowing out from the cluster center. In most cases they show ...synchrotron spectra that steepen toward the cluster center, indicating that they are caused by relativistic electrons being accelerated at outward traveling merger shocks. A number of radio relics break with this ideal picture and show morphologies that are bent the opposite way and show spectral index distributions that do not follow expectations from the ideal picture. We propose that these “wrong way” relics can form when an outward traveling shock wave is bent inward by an infalling galaxy cluster or group. We test this in an ultra-high-resolution zoom-in simulation of a massive galaxy cluster with an on-the-fly spectral cosmic-ray model. This allows us to study not only the synchrotron emission at colliding shocks, but also their synchrotron spectra to address the open question of relics with strongly varying spectral indices over the relic surface.
Abstract We present the first results of one extremely high-resolution, nonradiative magnetohydrodynamical cosmological zoom-in simulation of a massive cluster with a virial mass of M vir = 2.0 × 10 ...15 solar masses. We adopt a mass resolution of 4 × 10 5 M ⊙ with a maximum spatial resolution of around 250 pc in the central regions of the cluster. We follow the detailed amplification process in a resolved small-scale turbulent dynamo in the intracluster medium (ICM) with strong exponential growth until redshift 4, after which the field grows weakly in the adiabatic compression limit until redshift 2. The energy in the field is slightly reduced as the system approaches redshift zero in agreement with adiabatic decompression. The field structure is highly turbulent in the center and shows field reversals on a length scale of a few tens of kiloparsecs and an anticorrelation between the radial and angular field components in the central region that is ordered by small-scale turbulent dynamo action. The large-scale field on megaparsec scales is almost isotropic, indicating that the structure formation process in massive galaxy cluster formation suppresses any memory of both the initial field configuration and the amplified morphology via the turbulent dynamo. We demonstrate that extremely high-resolution simulations of the magnetized ICM are within reach that can simultaneously resolve the small-scale magnetic field structure, which is of major importance for the injection of and transport of cosmic rays in the ICM. This work is a major cornerstone for follow-up studies with an on-the-fly treatment of cosmic rays to model in detail electron-synchrotron and gamma-ray emissions.
Abstract
We present non-radiative, cosmological zoom-in simulations of galaxy-cluster formation with magnetic fields and (anisotropic) thermal conduction of one massive galaxy cluster with
M
vir
∼ 2 ...× 10
15
M
⊙
at
z
∼ 0. We run the cluster on three resolution levels (1×, 10×, 25×), starting with an effective mass resolution of 2 × 10
8
M
⊙
, subsequently increasing the particle number to reach 4 × 10
6
M
⊙
. The maximum spatial resolution obtained in the simulations is limited by the gravitational softening reaching
ϵ
= 1.0 kpc at the highest resolution level, allowing one to resolve the hierarchical assembly of the structures in fine detail. All simulations presented are carried out with the SPMHD code
gadget3
with an updated SPMHD prescription. The primary focus of this paper is to investigate magnetic field amplification in the intracluster medium. We show that the main amplification mechanism is the small-scale turbulent dynamo in the limit of reconnection diffusion. In our two highest resolution models we start to resolve the magnetic field amplification driven by the dynamo and we explicitly quantify this with the magnetic power spectra and the curvature of the magnetic field lines, consistent with dynamo theory. Furthermore, we investigate the ∇ ·
B
= 0 constraint within our simulations and show that we achieve comparable results to state-of-the-art AMR or moving-mesh techniques, used in codes such as
enzo
and
arepo
. Our results show for the first time in a cosmological simulation of a galaxy cluster that dynamo action can be resolved with modern numerical Lagrangian magnetohydrodynamic methods, a study that is currently missing in the literature.
Radio Relics are typically found to be arc-like regions of synchrotron emission in the outskirts of merging galaxy clusters, bowing out from the cluster center. In most cases they show synchrotron ...spectra that steepen towards the cluster center, indicating that they are caused by relativistic electrons being accelerated at outwards traveling merger shocks. A number of radio relics break with this ideal picture and show morphologies that are bent the opposite way and show spectral index distributions which do not follow expectations from the ideal picture. We propose that these `Wrong Way' Relics can form when an outwards travelling shock wave is bent inwards by an in-falling galaxy cluster or group. We test this in an ultra-high resolution zoom-in simulation of a massive galaxy cluster with an on-the-fly spectral Cosmic Ray model. This allows us to study not only the synchrotron emission at colliding shocks, but also their synchrotron spectra to adress the open question of relics with strongly varying spectral indices over the relic surface.
Non-thermal emission from relativistic electrons gives insight into the strength and morphology of intra-cluster magnetic fields, as well as providing powerful tracers of structure formation shocks. ...Emission caused by Cosmic Ray (CR) protons on the other hand still challenges current observations and is therefore testing models of proton acceleration at intra-cluster shocks. Large-scale simulations including the effects of CRs have been difficult to achieve and have been mainly reduced to simulating an overall energy budget, or tracing CR populations in post-processing of simulation output and has often been done for either protons or electrons. We use an efficient on-the-fly Fokker-Planck solver to evolve distributions of CR protons and electrons within every resolution element of our simulation. The solver accounts for CR acceleration at intra-cluster shocks, based on results of recent PIC simulations, re-acceleration due to shocks and MHD turbulence, adiabatic changes and radiative losses of electrons. We apply this model to zoom simulations of galaxy clusters, recently used to show the evolution of the small-scale turbulent dynamo on cluster scales. For these simulations we use a spectral resolution of 48 bins over 6 orders of magnitude in momentum for electrons and 12 bins over 6 orders of magnitude in momentum for protons. We present preliminary results about a possible formation mechanism for Wrong Way Radio Relics in our simulation.
Non-thermal emission from relativistic Cosmic Ray (CR) electrons gives insight into the strength and morphology of intra-cluster magnetic fields, as well as providing powerful tracers of structure ...formation shocks. Emission caused by CR protons on the other hand still challenges current observations and is therefore testing models of proton acceleration at intra-cluster shocks. Large-scale simulations including the effects of CRs have been difficult to achieve and have been mainly reduced to simulating an overall energy budget, or tracing CR populations in post-processing of simulation output and has often been done for either protons or electrons. We introduce CRESCENDO: Cosmic Ray Evolution with SpeCtral Electrons aND prOtons, an efficient on-the-fly Fokker-Planck solver to evolve distributions of CR protons and electrons within every resolution element of our simulation. The solver accounts for CR (re-)acceleration at intra-cluster shocks, based on results of recent PIC simulations, adiabatic changes and radiative losses of electrons. We show its performance in test cases as well as idealized galaxy cluster (GC) simulations. We apply the model to an idealized GC merger following best-fit parameters for CIZA J2242.4+5301-1 and study CR injection, radio relic morphology, spectral steepening and synchrotron emission.
Aims: Detecting diffuse synchrotron emission from the cosmic web is still a challenge for current radio telescopes. We aim to make predictions for the detectability of cosmic web filaments from ...simulations. Methods: We present the first cosmological MHD simulation of a 500 \(h^{-1} c\)Mpc volume with an on-the-fly spectral cosmic ray (CR) model. This allows us to follow the evolution of populations of CR electrons and protons within every resolution element of the simulation. We model CR injection at shocks, while accounting for adiabatic changes to the CR population and high energy loss processes of electrons. The synchrotron emission is then calculated from the aged electron population, using the simulated magnetic field, as well as different models for origin and amplification of magnetic fields. We use constrained initial conditions, which closely resemble the local Universe and compare the results of the cosmological volume to zoom-in simulation of the Coma cluster, to study the impact of resolution and turbulent re-acceleration of CRs on the results. Results: We find consistent injection of CRs at accretion shocks onto cosmic web filaments and galaxy clusters. This leads to diffuse emission from filaments of the order \(S_\nu \approx 0.1 \: \mu\)Jy beam\(^{-1}\) for a potential LOFAR observation at 144 MHz, when assuming the most optimistic magnetic field model and the inclusion of an on-the-fly treatment of re-acceleration of electrons by turbulence. The flux can be increased by up-to two orders of magnitude for different choices of CR injection parameters. This can bring the flux within a factor of 10 of the current limits for direct detection. We find a spectral index of the simulated synchrotron emission from filaments of {\alpha} {\approx} 1.0 - 1.5.