We present a novel and fast application to generate glass-like initial conditions for Lagrangian hydrodynamic schemes (e.g. Smoothed Particle Hydrodynamics (SPH)) following arbitrary density models based on weighted Voronoi tessellations and combine it with improved initial configurations and an additional particle reshuffling scheme. We show our application's ability to sample different kinds of density features and to converge properly towards the given model density as well as a glass-like particle configuration. We analyse convergence with iterations as well as with varying particle number. Additionally, we demonstrate the versatility of the implemented algorithms by providing an extensive test suite for standard (magneto-) hydrodynamic test cases as well as a few common astrophysical applications. We indicate the potential to bridge further between observational astronomy and simulations as well as applicability to other fields of science by advanced features such as describing a density model using gridded data for exampling from an image file instead of an analytic model.

The archetype FR-II galaxy Cygnus A lies in a moderately rich cluster about to undergo a major merger. We study the pre-merger Cygnus cluster environment using smoothed particle hydrodynamics simulations constrained by 2Ms of Chandra observations of the hot intracluster medium. The observations constrain the total gravitating mass and concentration parameter, and the simulations provide the quiescent and merger-enhanced temperature profiles of the pre- and post merger of the cluster excluding the central active galactic nucleus. We present the first detailed model of the sub cluster north west of Cygnus A, named CygNW. We find a lower baryon fraction and higher concentration parameter for CygA than expected from known scaling relations in the literature. The model suggests the Cygnus cluster hosts a pre-merger with a progenitor mass ratio of about 1.5:1 at the virial radius. We notice that the intra cluster medium is heated as a result of the merger, but we find no evidence for a (pre-)merger shock in the interstitial region between both cluster haloes. We attribute the merger-induced heating to compression of the cluster outskirts. The smooth model obtained from our simulations is subtracted from the observed cluster state and shows residual temperature structure that is neither hydrostatic nor merger-heated cluster gas. We speculate that this residual heating may result from previous AGN activity over the last ~100 Myr.

Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a fifth order weighted essentially non-oscillatory scheme for constrained-transport magnetohydrodynamics into the code WOMBAT. We establish the correctness of our implementation with an extensive number tests. We find that the fifth order scheme performs as accurately as a common second order scheme at half the resolution. We argue that for a given solution quality the new scheme is more computationally efficient than lower order schemes in three dimensions. We also establish the performance characteristics of the solver in the WOMBAT framework. Our implementation fully vectorizes using flattened arrays in thread-local memory. It performs at about 0.6 Million zones per second per node on Intel Broadwell. We present scaling tests of the code up to 98 thousand cores on the Cray XC40 machine "Hazel Hen", with a sustained performance of about 5 percent of peak at scale.

In galaxy clusters, modern radio interferometers observe non-thermal radio sources with unprecedented spatial and spectral resolution. For the first time, the new data allows to infer the structure of the intra-cluster magnetic fields on small scales via Faraday tomography. This leap forward demands new numerical models for the amplification of magnetic fields in cosmic structure formation - the cosmological magnetic dynamo. Here we present a novel numerical approach to astrophyiscal MHD simulations aimed to resolve this small-scale dynamo in future cosmological simulations. As a first step, we implement a fifth order WENO scheme in the new code WOMBAT. We show that this scheme doubles the effective resolution of the simulation and is thus less expensive than common second order schemes. WOMBAT uses a novel approach to parallelization and load balancing developed in collaboration with performance engineers at Cray Inc. This will allow us scale simulation to the exaflop regime and achieve kpc resolution in future cosmological simulations of galaxy clusters. Here we demonstrate the excellent scaling properties of the code and argue that resolved simulations of the cosmological small scale dynamo within the whole virial radius are possible in the next years.

We review the present theoretical and numerical understanding of magnetic field amplification in cosmic large-scale structure, on length scales of galaxy clusters and beyond. Structure formation drives compression and turbulence, which amplify tiny magnetic seed fields to the microGauss values that are observed in the intracluster medium. This process is intimately connected to the properties of turbulence and the microphysics of the intra-cluster medium. Additional roles are played by merger induced shocks that sweep through the intra-cluster medium and motions induced by sloshing cool cores. The accurate simulation of magnetic field amplification in clusters still poses a serious challenge for simulations of cosmological structure formation. We review the current literature on cosmological simulations that include magnetic fields and outline theoretical as well as numerical challenges.

Previous studies have shown that CIZA J2242.8+5301 (the 'Sausage' cluster, $z=0.192$) is a massive merging galaxy cluster that hosts a radio halo and multiple relics. In this paper we present deep, high fidelity, low-frequency images made with the LOw-Frequency Array (LOFAR) between 115.5 and 179 MHz. These images, with a noise of 140 mJy/beam and a resolution of $\theta_{\text{beam}}=7.3"\times5.3"$, are an order of magnitude more sensitive and five times higher resolution than previous low-frequency images of this cluster. We combined the LOFAR data with the existing GMRT (153, 323, 608 MHz) and WSRT (1.2, 1.4, 1.7, 2.3 GHz) data to study the spectral properties of the radio emission from the cluster. Assuming diffusive shock acceleration (DSA), we found Mach numbers of $\mathcal{M}_{n}=2.7{}_{-0.3}^{+0.6}$ and $\mathcal{M}_{s}=1.9_{-0.2}^{+0.3}$ for the northern and southern shocks. The derived Mach number for the northern shock requires an acceleration efficiency of several percent to accelerate electrons from the thermal pool, which is challenging for DSA. Using the radio data, we characterised the eastern relic as a shock wave propagating outwards with a Mach number of $\mathcal{M}_{e}=2.4_{-0.3}^{+0.5}$, which is in agreement with $\mathcal{M}_{e}^{X}=2.5{}_{-0.2}^{+0.6}$ that we derived from Suzaku data. The eastern shock is likely to be associated with the major cluster merger. The radio halo was measured with a flux of $346\pm64\,\text{mJy}$ at $145\,\text{MHz}$. Across the halo, we observed a spectral index that remains approximately constant ($\alpha^{\text{145 MHz-2.3 GHz}}_{\text{across \(\sim\)1 Mpc}^2}=-1.01\pm0.10$) after the steepening in the post-shock region of the northern relic. This suggests a generation of post-shock turbulence that re-energies aged electrons.

The giant radio relic in CIZA J2242.8+5301 is likely evidence of a Mpc sized shock in a massive merging galaxy cluster. However, the exact shock properties are still not clearly determined. In particular, the Mach number derived from the integrated radio spectrum exceeds the Mach number derived from the X-ray temperature jump by a factor of two. We present here a numerical study, aiming for a model that is consistent with the majority of observations of this galaxy cluster. We first show that in the northern shock upstream X-ray temperature and radio data are consistent with each other. We then derive progenitor masses for the system using standard density profiles, X-ray properties and the assumption of hydrostatic equilibrium. We find a class of models that is roughly consistent with weak lensing data, radio data and some of the X-ray data. Assuming a cool-core versus non-cool-core merger, we find a fiducial model with a total mass of $1.6 \times 10^{15}\,M_\odot$, a mass ratio of 1.76 and a Mach number that is consistent with estimates from the radio spectrum. We are not able to match X-ray derived Mach numbers, because even low mass models over-predict the X-ray derived shock speeds. We argue that deep X-ray observations of CIZA J2242.8+5301 will be able to test our model and potentially reconcile X-ray and radio derived Mach numbers in relics.

We present a new code for astrophysical magneto-hydrodynamics specifically designed and optimized for high performance and scaling on modern and future supercomputers. We describe a novel hybrid OpenMP/MPI programming model that emerged from a collaboration between Cray, Inc. and the University of Minnesota. This design utilizes MPI-RMA optimized for thread scaling, which allows the code to run extremely efficiently at very high thread counts ideal for the latest generation of the multi-core and many-core architectures. Such performance characteristics are needed in the era of "exascale" computing. We describe and demonstrate our high-performance design in detail with the intent that it may be used as a model for other, future astrophysical codes intended for applications demanding exceptional performance.

Giant radio relics are the arc-shaped diffuse radio emission regions observed in the outskirts of some merging galaxy clusters. They are believed to trace shock-waves in the intra-cluster medium. Recent observations demonstrated that some prominent radio relics exhibit a steepening above 2 GHz in their radio spectrum. This challenges standard theoretical models because shock acceleration is expected to accelerate electrons to very high energies with a power-law distribution in momentum. In this work we attempt to reconcile these data with the shock-acceleration scenario. We propose that the spectral steepening may be caused by the highest energy electrons emitting preferentially in lower magnetic fields than the bulk of synchrotron bright electrons in relics. Here, we focus on a model with an increasing mag- netic field behind the shock front, which quickly saturates and then declines. We derive the time-evolution of cosmic-ray electron spectra in time variable magnetic fields and an expanding medium. We then apply the formalism on the large radio relic in the cluster CIZA J2242.8+5301 (the Sausage relic). We show that under favourable circumstances of magnetic field amplification downstream, our model can explain the observed radio spectrum, the brightness profile and the spectral index profile of the relic. A possible interpretation for the required amplification of the magnetic field downstream is a dynamo acting behind the shock with an injection scale of magnetic turbulence of about 10 kpc. Our models require injection efficiencies of CRe - which are in tension with simple diffusive shock acceleration from the thermal pool. We show that this problem can likely be alleviated considering pre-existing CRe.

We present an on-the-fly geometrical approach for shock detection and Mach number calculation in simulations employing smoothed particle hydrodynamics (SPH). We utilize pressure gradients to select shock candidates and define up- and downstream positions. We obtain hydrodynamical states in the up- and downstream regimes with a series of normal and inverted kernel weightings parallel and perpendicular to the shock normals. Our on-the-fly geometrical Mach detector incorporates well within the SPH formalism and has low computational cost. We implement our Mach detector into the simulation code GADGET and alongside many SPH improvements. We test our shock finder in a sequence of shock-tube tests with successively increasing Mach numbers exceeding by far the typical values inside galaxy clusters. For all shocks, we resolve the shocks well and the correct Mach numbers are assigned. An application to a strong magnetized shock-tube gives stable results in full magnetohydrodynamic set-ups. We simulate a merger of two idealized galaxy clusters and study the shock front. Shock structures within the merging clusters as well as the cluster shock are well-captured by our algorithm and assigned correct Mach numbers.

Radio relics are patches of diffuse synchrotron radio emission that trace shock waves. Relics are thought to form when intra-cluster medium electrons are accelerated by cluster merger induced shock waves through the diffusive shock acceleration mechanism. In this paper, we present observations spanning 150 MHz to 30 GHz of the `Sausage' and `Toothbrush' relics from the Giant Metrewave and Westerbork telescopes, the Karl G. Jansky Very Large Array, the Effelsberg telescope, the Arcminute Microkelvin Imager and Combined Array for Research in Millimeter-wave Astronomy. We detect both relics at 30 GHz, where the previous highest frequency detection was at 16 GHz. The integrated radio spectra of both sources clearly steepen above 2 GHz, at the >6$\sigma$ significance level, supports the spectral steepening previously found in the `Sausage' and the Abell 2256 relic. Our results challenge the widely adopted simple formation mechanism of radio relics and suggest more complicated models have to be developed that, for example, involve re-acceleration of aged seed electrons.

We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.

Galaxy clusters are unique laboratories to investigate turbulent fluid motions and large scale magnetic fields. Synchrotron radio halos at the center of merging galaxy clusters provide the most spectacular and direct evidence of the presence of relativistic particles and magnetic fields associated with the intracluster medium. The study of polarized emission from radio halos is extremely important to constrain the properties of intracluster magnetic fields and the physics of the acceleration and transport of the relativistic particles. However, detecting this polarized signal is a very hard task with the current radio facilities.We use cosmological magneto-hydrodynamical simulations to predict the expected polarized surface brightness of radio halos at 1.4 GHz. We compare these expectations with the sensitivity and the resolution reachable with the SKA1. This allows us to evaluate the potential for studying intracluster magnetic fields in the surveys planned for SKA1.

Particle acceleration by turbulence plays a role in many astrophysical environments. The non- linear evolution of the underlying cosmic-ray spectrum is complex and can be described by a Fokker-Planck equation, which in general has to be solved numerically. We present here an implementation to compute the evolution of a cosmic-ray spectrum coupled to turbulence considering isotropic particle pitch-angle distributions and taking into account the relevant particle energy gains and losses. Our code can be used in run time and post-processing to very large astrophysical fluid simulations. We also propose a novel method to compress cosmic- ray spectra by a factor of ten, to ease the memory demand in very large simulations. We show a number of code tests, which firmly establish the correctness of the code. In this paper we focus on relativistic electrons, but our code and methods can be easily extended to the case of hadrons. We apply our pipeline to the relevant problem of particle acceleration in galaxy clusters. We define a sub-grid model for compressible MHD-turbulence in the intra- cluster-medium and calculate the corresponding reacceleration timescale from first principles. We then use a magneto-hydrodynamic simulation of an isolated cluster merger to follow the evolution of relativistic electron spectra and radio emission generated from the system over several Gyrs.

Simulations of isolated binary mergers of galaxy clusters are a useful tool to study the evolution of these objects. For exceptionally massive systems they even represent the only viable way of simulation, because these are rare in typical cosmological simulations. We present a new practical model for these simulations based on the Hernquist dark matter profile. The hydrostatic equation is solved for a beta-model with $\beta$ = 2/3 in this potential and approximate expressions for X-ray brightness and Compton-y parameter are derived. We show in detail how to setup such a system using SPH. The theoretical and several numerical models are compared to observed scaling relations of galaxy clusters and satisfactory agreement with the self-similar relations is found. The model is then applied to investigate the observed cluster ACT-CT J0102-4915 (El Gordo), a particularly massive merging high redshift cluster. We are able to reproduce the X-ray luminosity, SZ-effect and dark matter core distance as well as the rough shape of the observed cluster for reasonable model parameters. The lack of substruc- ture prevents us from obtaining the fluctuations observed in the wake of the system and we argue that the parent cluster of the system was highly disturbed even before the main merger observed today.

The spectrum of the radio halo in the Coma cluster is measured over almost two decades in frequency. The current radio data show a steepening of the spectrum at higher frequencies, which has implications for models of the radio halo origin. There is an on-going debate on the possibility that the observed steepening is not intrinsic to the emitted radiation, but is instead caused by the SZ effect. Recently, the Planck satellite measured the SZ signal and its spatial distribution in the Coma cluster allowing to test this hypothesis. Using the Planck results, we calculated the modification of the radio halo spectrum by the SZ effect in three different ways. With the first two methods we measured the SZ-decrement within the aperture radii used for flux measurements of the halo at the different frequencies. First we adopted the global compilation of data from Thierbach et al. and a reference aperture radius consistent with those used by the various authors. Second we used the available brightness profiles of the halo at different frequencies to derive the spectrum within two fixed apertures, and derived the SZ-decrement using these apertures. As a third method we used the quasi-linear correlation between the y and the radio-halo brightness at 330 MHz discovered by Planck to derive the modification of the radio spectrum by the SZ-decrement in a way that is almost independent of the adopted aperture radius. We found that the spectral modification induced by the SZ-decrement is 4-5 times smaller than that necessary to explain the observed steepening. Consequently a break or cut-off in the spectrum of the emitting electrons is necessary to explain current data. We also show that, if a steepening is absent from the emitted spectrum, future deep observations at 5 GHz with single dishes are expected to measure a halo flux in a 40 arcmin radius that would be 7-8 times higher than currently seen.

We review models for giant radio halos in clusters of galaxies, with a focus on numerical and theoretical work. After summarising the most important observations of these objects, we present an introduction to the theoretical aspects of hadronic models. We compare these models with observations using simulations and find severe problems for hadronic models. We give a short introduction to reacceleration models and show results from the first simulation of CRe reaccel- eration in cluster mergers. We find that in-line with previous theoretical work, reacceleration models are able to elegantly explain main observables of giant radio halos.

We present the first high resolution MHD simulation of cosmic-ray electron reacceleration by turbulence in cluster mergers. We use an idealised model for cluster mergers, combined with a numerical model for the injection, cooling and reacceleration of cosmic-ray electrons, to investigate the evolution of cluster scale radio emission in these objects. In line with theoretical expectations, we for the first time, show in a simulation that reacceleration of CRe has the potential to reproduce key observables of radio halos. In particular, we show that clusters evolve being radio loud or radio quiet, depending on their evolutionary stage during the merger. We thus recover the observed transient nature of radio halos. In the simulation the diffuse emission traces the complex interplay between spatial distribution of turbulence injected by the halo infall and the spatial distribution of the seed electrons to reaccelerate. During the formation and evolution of the halo the synchrotron emission spectra show the observed variety: from power-laws with spectral index of 1 to 1.3 to curved and ultra-steep spectra with index $> 1.5$.

We present high resolution simulations of a multiple merger of three disk galaxies including the evolution of magnetic fields performed with the N-body/SPH code Gadget. For the first time, we embed the galaxies in a magnetized, low-density medium, thus modeling an ambient IGM. The simulations include radiative cooling and a model for star formation and supernova feedback. The progenitor disks have initial magnetic seed fields in the range of 10e-9 to 10e-6 G and the IGM has initial fields of 10e-12 to 10e-9 G. The simulations are compared to a run excluding magnetic fields. We show that the propagation of interaction-driven shocks depends significantly on the initial magnetic field strength. The shocks propagate faster in simulations with stronger initial field, suggesting that the shocks are supported by magnetic pressure. The Mach numbers of the shocks range from approximately M=1.5 for the non-magnetized case up to M=6 for the highest initial magnetization, resulting in higher temperatures of the shock heated IGM gas. The magnetic field in the system saturates rapidly after the mergers at ~ 10e-6 G within the galaxies and ~ 10e-8 G in the IGM independent of the initial value. These field strengths agree with observed values and correspond to the equipartition value of the magnetic pressure with the turbulent pressure in the system. We also present synthetic radio and polarization maps for different phases of the evolution showing that shocks driven by the interaction produce a high amount of polarized emission. These idealized simulations indicate that magnetic fields play an important role for the hydrodynamics of the IGM during galactic interactions. We also show that even weak seed fields are efficiently strengthened during multiple galactic mergers. This interaction driven amplification might have been a key process for the magnetization of the Universe.

Using cosmological MHD simulations of the magnetic field in galaxy clusters and filaments we evaluate the possibility to infer the magnetic field strength in filaments by measuring cross-correlation functions between Faraday Rotation Measures (RM) and the galaxy density field. We also test the reliability of recent estimates considering the problem of data quality and Galactic foreground (GF) removal in current datasets. Besides the two self-consistent simulations of cosmological magnetic fields based on primordial seed fields and galactic outflows analyzed here, we also explore a larger range of models scaling up the resulting magnetic fields of one of the simulations. We find that, if an unnormalized estimator for the cross-correlation functions and a GF removal procedure is used, the detectability of the cosmological signal is only possible for future instruments (e.g. SKA and ASKAP). However, mapping of the observed RM signal to the underlying magnetization of the Universe (both in space and time) is an extremely challenging task which is limited by the ambiguities of our model parameters, as well as to the weak response of the RM signal in low density environments. Therefore, we conclude that current data cannot constrain the amplitude and distribution of magnetic fields within the large scale structure and a detailed theoretical understanding of the build up and distribution of magnetic fields within the Universe will be needed for the interpretation of future observations.

We use results from a constrained, cosmological MHD simulation of the Local Universe to predict radio halos and their evolution for a volume limited set of galaxy clusters and compare to current observations. The simulated magnetic field inside the clusters is a result of turbulent amplification within them, with the magnetic seed originating from star-burst driven, galactic outflows. We evaluate three models, where we choose different normalizations for the Cosmic Ray proton population within clusters. Similar to our previous analysis of the Coma cluster (Donnert et al. 2010), the radial profile and the morphological properties of observed radio halos can not be reproduced, even with a radially increasing energy fraction within the cosmic ray proton population. Scaling relations between X-ray luminosity and radio power can be reproduced by all models, however all models fail in the prediction of clusters with no radio emission. Also the evolutionary tracks of our largest clusters in all models fail to reproduce the observed bi-modality in radio luminosity. This provides additional evidence that the framework of hadronic, secondary models is disfavored to reproduce the large scale diffuse radio emission of galaxy clusters. We also provide predictions for the unavoidable emission of $\gamma$-rays from the hadronic models for the full cluster set. None of such secondary models is yet excluded by the observed limits in $\gamma$-ray emission, emphasizing that large scale diffuse radio emission is a powerful tool to constrain the amount of cosmic ray protons in galaxy clusters.

We use the results from a constrained, cosmological MHD simulation of the Local Universe to predict the radio halo and the gamma-ray flux from the Coma cluster and compare it to current observations. The simulated magnetic field within the Coma cluster is the result of turbulent amplification of the magnetic field during build-up of the cluster. The magnetic seed field originates from star-burst driven, galactic outflows. The synchrotron emission is calculated assuming a hadronic model. We follow four approaches with different distributions for the cosmic-ray proton (CRp) population within galaxy clusters. The radial profile the radio halo can only be reproduced with a radially increasing energy fraction within the cosmic ray proton population, reaching $>$100% of the thermal energy content at $\approx$ 1Mpc, e.g. the edge of the radio emitting region. Additionally the spectral steepening of the observed radio halo in Coma cannot be reproduced, even when accounting for the negative flux from the thermal SZ effect at high frequencies. Therefore the hadronic models are disfavored from present analysis. The emission of $\gamma$-rays expected from our simulated coma is still below the current observational limits (by a factor of $\sim$6) but would be detectable in the near future.

We performed cosmological, magneto-hydrodynamical simulations to follow the evolution of magnetic fields in galaxy clusters, exploring the possibility that the origin of the magnetic seed fields are galactic outflows during the star-burst phase of galactic evolution. To do this we coupled a semi-analytical model for magnetized galactic winds as suggested by \citet{2006MNRAS.370..319B} to our cosmological simulation. We find that the strength and structure of magnetic fields observed in galaxy clusters are well reproduced for a wide range of model parameters for the magnetized, galactic winds and do only weakly depend on the exact magnetic structure within the assumed galactic outflows. Although the evolution of a primordial magnetic seed field shows no significant differences to that of galaxy clusters fields from previous studies, we find that the magnetic field pollution in the diffuse medium within filaments is below the level predicted by scenarios with pure primordial magnetic seed field. We therefore conclude that magnetized galactic outflows and their subsequent evolution within the intra-cluster medium can fully account for the observed magnetic fields in galaxy clusters. Our findings also suggest that measuring cosmological magnetic fields in low-density environments such as filaments is much more useful than observing cluster magnetic fields to infer their possible origin.