The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is over-simplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds (PhEW), that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a sub-grid prescription to model galactic winds more robustly both physically and numerically.

Dwarf galaxies are dark matter-dominated and therefore promising targets for the search for weakly interacting massive particles (WIMPs), which are well-known candidates for dark matter. Annihilation of WIMPs produce ultra-relativistic cosmic-ray electrons and positrons that emit synchrotron radiation in the presence of magnetic fields. For typical magnetic field strengths (few $\mu $G) and $\mathcal O$(GeV--TeV) WIMP masses, this emission peaks at hundreds of MHz. Here, we use the non-detection of 150-MHz radio continuum emission from the dwarf spheroidal galaxy Canes Venatici I with the LOw-Frequency ARray (LOFAR) to derive constraints on the annihilation cross section of WIMPs into primary electron-positron and other fundamental particle-antiparticle pairs. Our main underlying assumption is that the transport of the cosmic rays can be described by the diffusion approximation, thus requiring a non-zero magnetic field strength with small-scale structure. In particular, by adopting magnetic field strengths of $\mathcal O(1\,\mu$G) and diffusion coefficients $\sim 10^{27}~\rm cm^2\,s^{-1}$, we obtain limits that are comparable with those set by \emph{Fermi} Large Area Telescope using gamma-ray observations of this particular galaxy. Assuming s-wave annihilation and WIMPs making up 100 per cent of the DM density, our benchmark limits exclude several thermal WIMP realisations in the $[2,20]$-GeV mass range. We caution, however, that our limits for the cross section are subject to enormous uncertainties which we also quantitatively assess. In particular, variations on the propagation parameters or on the DM halo can shift our limits up by several orders of magnitude (in the pessimistic scenario).

We present LOFAR observations at 120-168 MHz of 42 systems with possible X-ray cavities in their hot atmosphere, of which 17 are groups or ellipticals, 19 are nearby clusters (z<0.3), and six are higher-redshift clusters (z>0.3). The X-ray cavities, formed by the radio lobes of the central active galactic nucleus (AGN), are evidence of radio-mode AGN feedback. In the groups and ellipticals sample, more than half of the systems have X-ray cavities for which no associated lobe emission was detected. Conversely, we report the discovery of large radio lobes in NGC 6338, extending far beyond the cavities reported previously in the literature. In the case of the nearby clusters, our observations show that there is little low-frequency radio emission that extends beyond the cavities (e.g., MS 0735.6+7421 and A2052). For the first time, we report secure cavity-radio associations in 2A 0335+096, ZwCl 2701, and ZwCl 8276 that strengthens their interpretation as AGN-created cavities. However, in some known cavity systems (e.g., A1795 and ZwCl 3146) we report the lack of detectable low-frequency radio emission associated with the cavities. Our sample of higher-redshift systems is small, and unfortunately the present LOFAR observations are not able to resolve the lobes in many of them. Nevertheless, our sample represents one of the best available for investigating the connection between radio and jet power in radio-mode AGN feedback.

We present LOFAR observations of one of the most spectacular objects in the radio sky: Abell 2255. This is a nearby ($z = 0.0806$) merging galaxy cluster hosting one of the first radio halos ever detected in the intra-cluster medium (ICM). The deep LOFAR images at 144 MHz of the central $\sim10$ Mpc$^2$ region show a plethora of emission on different scales, from tens of kpc to above Mpc sizes. In this work, we focus on the innermost region of the cluster. Among the numerous interesting features observed, we discover remarkable bright and filamentary structures embedded in the radio halo. We incorporate archival WSRT 1.2 GHz data to study the spectral properties of the diffuse synchrotron emission and find a very complex spectral index distribution in the halo spanning a wide range of values. We combine the radio data with Chandra observations to investigate the connection between the thermal and non-thermal components by quantitatively comparing the radio and X-ray surface brightness and the spectral index of the radio emission with the thermodynamical quantities of the ICM. Despite the multitude of structures observed in the radio halo, we find that the X-ray and radio emission are overall well correlated. The fact that the steepest spectrum emission is located in the cluster center and traces regions with high entropy possibly suggests the presence of seed particles injected by radio galaxies that are spread in the ICM by turbulence generating the extended radio halo.

The large-scale magnetic field of the Milky Way is thought to be created by an alpha-Omega dynamo, which implies that it should have opposite handedness North and South of the Galactic midplane. Here we attempt to detect a variation in handedness using polarization data from the Wilkinson Microwave Anisotropy Probe. Previous analyzes of the parity-even and parity-odd parts of linear polarization of the global dust and synchrotron emission have focused on quadratic correlations in spectral space of, and between, these two components. Here, by contrast, we analyze the parity-odd polarization itself and show that it has, on average, opposite signs in Northern and Southern Galactic hemispheres. Comparison with a Galactic mean-field dynamo model shows broad qualitative agreement and reveals that the sign of the observed hemispheric dependence of the azimuthally averaged parity-odd polarization is not determined by the sign of alpha, but by the sense of differential rotation.

Measuring the properties of extragalactic magnetic fields through the effect of Faraday rotation provides a means to understand the origin and evolution of cosmic magnetism. Here we use data from the LOFAR Two-Metre Sky Survey (LoTSS) to calculate the Faraday rotation measure (RM) of close pairs of extragalactic radio sources. By considering the RM difference ($\Delta$RM) between physical pairs (e.g. double-lobed radio galaxies) and non-physical pairs (i.e. close projected sources on the sky), we statistically isolate the contribution of extragalactic magnetic fields to $\Delta$RM along the line of sight between non-physical pairs. From our analysis, we find no significant difference between the $\Delta$RM distributions of the physical and non-physical pairs, limiting the excess Faraday rotation contribution to $< 1.9$ rad/m$^2$ ($\sim$$95\%$ confidence). We use this limit with a simple model of an inhomogeneous universe to place an upper limit of 4 nG on the cosmological co-moving magnetic field strength on Mpc scales. We also compare the RM data with a more realistic suite of cosmological MHD simulations, that explore different magnetogenesis scenarios. Both magnetization of the large scale structure by astrophysical processes such as galactic and AGN outflows, and simple primordial scenarios with seed magnetic field strengths $< 0.5$ nG cannot be rejected by the current data; while stronger primordial fields or models with dynamo amplification in filaments are disfavoured.

The origin of radio relics is usually explained via diffusive shock acceleration (DSA) or re-acceleration of electrons at/from merger shocks in galaxy clusters. The case of acceleration is challenged by the low predicted efficiency of low-Mach number merger shocks, unable to explain the power observed in most radio relics. In this Letter we present the discovery of a new giant radio relic around the galaxy cluster Abell 2249 ($z=0.0838$) using LOFAR. It is special since it has the lowest surface brightness of all known radio relics. We study its radio and X-ray properties combinig LOFAR data with uGMRT, JVLA and XMM. This object has a total power of $L_{1.4\rm GHz}=4.1\pm 0.8 \times 10^{23}$ W Hz$^{-1}$ and integrated spectral index $\alpha = 1.15\pm 0.23$. We infer for this radio relic a lower bound on the magnetisation of $B\geq 0.4\, \mu$G, a shock Mach number of $\mathcal{M}\approx 3.79$, and a low acceleration efficiency consistent with DSA. This result suggests that a missing population of relics may become visible thanks to the unprecedented sensitivity of the new generation of radio telescopes.

Giant radio galaxies (GRGs) are physically large radio sources that extend well beyond their host galaxy environment. Their polarization properties are affected by the poorly constrained magnetic field that permeates the intergalactic medium on Mpc scales. A low frequency ($<$ 200 MHz) polarization study of this class of radio sources is now possible with LOFAR. Here we investigate the polarization properties and Faraday rotation measure (RM) of a catalog of GRGs detected in the LoTSS. This is the first low-frequency polarization study of a large sample of radio galaxies selected on their physical size. We explore the magneto-ionic properties of their under-dense environment and probe intergalactic magnetic fields using the Faraday rotation properties of their radio lobes. We use RM synthesis in the 120-168 MHz band to search for polarized emission and to derive the RM and fractional polarization of each detected source component. We study the depolarization between 1.4 GHz and 144 MHz using images from the NVSS. From a sample of 240 GRGs, we detected 37 sources in polarization, all with a total flux density above 56 mJy. The fractional polarization of the detected GRGs at 1.4 GHz and 144 MHz is consistent with a small amount of Faraday depolarization (a Faraday dispersion $<$ 0.3 rad m$^{-2}$). Our analysis shows that the lobes are expanding into a low-density ($<10^{-5}$ cm$^{-3}$) local environment permeated by weak magnetic fields ($<$0.1 $\mu$G) with fluctuations on scales of 3 to 25 kpc. The presence of foreground galaxy clusters appears to influence the polarization detection rate up to 2R$_{500}$. In general, this work demonstrates the ability of LOFAR to quantify the rarefied environments in which these GRGs exist and highlights them as an excellent statistical sample to use as high precision probes of magnetic fields in the intergalactic medium and the Milky Way.

The degree of turbulent pressure support by residual gas motions in galaxy clusters is not well known. Mass modelling of combined X-ray and Sunyaev Zel'dovich observations provides an estimate of turbulent pressure support in the outer regions of several galaxy clusters. Here, we test two different filtering techniques to disentangle bulk from turbulent motions in non-radiative high-resolution cosmological simulations of galaxy clusters using the cosmological hydro code ENZO. We find that the radial behavior of the ratio of non-thermal pressure to total gas pressure as a function of cluster-centric distance can be described by a simple polynomial function. The typical non-thermal pressure support in the centre of clusters is $\sim$5%, increasing to $\sim$15% in the outskirts, in line with the pressure excess found in recent X-ray observations. While the complex dynamics of the ICM makes it impossible to reconstruct a simple correlation between turbulent motions and hydrostatic bias, we find that a relation between them can be established using the median properties of a sample of objects. Moreover, we estimate the contribution of radial accelerations to the non-thermal pressure support and conclude that it decreases moving outwards from 40% (in the core) to 15% (in the cluster's outskirts). Adding this contribution to one provided by turbulence, we show that it might account for the entire observed hydrostatic bias in the innermost regions of the clusters, and for less than 80% of it at $r > 0.8 r_{200, m}$.

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

Although a merging galaxy cluster is a useful laboratory to study many important astrophysical processes and fundamental physics, only limited interpretations are possible without careful analysis of the merger scenario. However, the study is demanding because a thorough comparison of multi-wavelength observations with detailed numerical simulations is required. In this paper, we present such a study for the off-axis binary merger A115. The system possesses a number of remarkable observational features, but no convincing merger scenario, explaining the shape and location of the radio relic in harmony with the orientation of the cometary X-ray tails, has been presented. Our hydrodynamical simulation, with adaptive mesh refinement, suggests that the cometary X-ray tail of A115 might be a slingshot tail and can arise $\sim0.3$ Gyrs after the impact and before the two subclusters reach their apocenters. This scenario can predict the location and orientation of the giant radio relic, which is parallel to the northern X-ray tail. In addition, our study indicates that diffusive shock acceleration alone cannot generate the observed radio power unless aided by stronger magnetic fields and/or more significant presence of fossil electrons.

Radio relics are vast synchrotron sources that sit on the outskirts of merging galaxy clusters. In this work we model their formation using a Press-Schechter formalism to simulate merger rates, analytical models for the intracluster medium and the shock dynamics, as well as a simple model for the cosmic-ray electrons at the merger shocks. We show that the statistical properties of the population of radio relics are strongly dependent on key physical parameters, such as the acceleration efficiency, the magnetic field strength at the relic, the geometry of the relic and the duration of the electron acceleration at merger shocks. It turns out that the flux distribution as well as the power-mass relation can constrain key parameters of the intracluster medium. With the advent of new large-area radio surveys, statistical analyses of radio relics will complement what we have learned from observations of individual objects.

Massive merging galaxy clusters often host diffuse Mpc-scale radio synchrotron emission. This emission originates from relativistic electrons in the ionized intracluster medium (ICM). An important question is how these synchrotron emitting relativistic electrons are accelerated. Our aim is to search for diffuse emission in the Frontier Fields clusters Abell S1063 and Abell 370 and characterize its properties. While these clusters are very massive and well studied at some other wavelengths, no diffuse emission has been reported for these clusters so far. We obtained 325 MHz Giant Metrewave Radio Telescope (GMRT) and 1--4 GHz Jansky Very Large Array (VLA) observations of Abell S1063 and Abell 370. We complement these data with Chandra and XMM-Newton X-ray observations. In our sensitive images, we discover radio halos in both clusters. In Abell S1063, a giant radio halo is found with a size of $\sim 1.2$ Mpc. The integrated spectral index between 325 MHz and 1.5 GHz is $-0.94\pm0.08$ and it steepens to $-1.77 \pm 0.20$ between 1.5 and 3.0 GHz. This spectral steepening provides support for the turbulent re-acceleration model for radio halo formation. Abell 370 hosts a faint radio halo mostly centred on the southern part of this binary merging cluster, with a size of $\sim 500-700$ kpc. The spectral index between 325 MHz and 1.5 GHz is $-1.10\pm0.09$. Both radio halos follow the known scaling relation between the cluster mass proxy $Y_{500}$ and radio power, consistent with the idea that they are related to ongoing cluster merger events.

Finding and classifying astronomical sources is key in the scientific exploitation of radio surveys. Source-finding usually involves identifying the parts of an image belonging to an astronomical source, against some estimated background. This can be problematic in the radio regime, owing to the presence of correlated noise, which can interfere with the signal from the source. In the current work, we present ConvoSource, a novel method based on a deep learning technique, to identify the positions of radio sources, and compare the results to a Gaussian-fitting method. Since the deep learning approach allows the generation of more training images, it should perform well in the source-finding task. We test the source-finding methods on artificial data created for the data challenge of the Square Kilometer Array (SKA). We investigate sources that are divided into three classes: star forming galaxies (SFGs) and two classes of active galactic nuclei (AGN). The~artificial data are given at two different frequencies (560~MHz and 1400~MHz), three total integration times (8 h, 100 h, 1000 h), and three signal-to-noise ratios (SNRs) of 1, 2, and 5. At~lower SNRs, ConvoSource tends to outperform a Gaussian-fitting approach in the recovery of SFGs and all sources, although at the lowest SNR of one, the better performance is likely due to chance matches. The~Gaussian-fitting method performs better in the recovery of the AGN-type sources at lower SNRs. At~a higher SNR, ConvoSource performs better on average in the recovery of AGN sources, whereas the Gaussian-fitting method performs better in the recovery of SFGs and all sources. ConvoSource usually performs better at shorter total integration times and detects more true positives and misses fewer sources compared to the Gaussian-fitting method; however, it detects more false positives.

The new generation of high-resolution broad-band radio telescopes, like the Low Frequency Array (LOFAR), produces, depending on the level of compression, between 1 to 10 TB of data per hour after correlation. Such a large amount of scientific data demand powerful computing resources and efficient data handling strategies to be mastered. The LOFAR Two-metre Sky Survey (LoTSS) is a Key Science Project (KSP) of the LOFAR telescope. It aims to map the entire northern hemisphere at unprecedented sensitivity and resolution. The survey consist of 3 168 pointings, requiring about 30 PBytes of storage space. As a member of the German Long Wavelength Consortioum (GLOW) the Forschungszentrum J\"ulich (FSZ) stores in the Long Term Archive (LTA) about 50% of all LoTSS observations conducted to date. In collaboration with SURFsara in Amsterdam we developed service tools that enable the KSP to process LOFAR data stored in the J\"ulich LTA at the supercomputer JUWELS in an automated and robust fashion. Through our system more than 500 out of 800 existing LoTSS observations have already been processed with the prefactor pipeline. This pipeline calibrates the direction-independent instrumental and ionospheric effects and furthermore reduces the data size significantly. For continuum imaging, this processing pipeline is the standard pipeline that is executed before more advanced processing and image reconstruction methods are applied.

It is well established that particle acceleration by shocks and turbulence in the intra-cluster medium can produce cluster-scale synchrotron emitting sources. However, the detailed physics of these particle acceleration processes is still not well understood. One of the main open questions is the role of fossil relativistic electrons that have been deposited in the intra-cluster medium by radio galaxies. These synchrotron-emitting electrons are very difficult to study, as their radiative life time is only tens of Myrs at GHz frequencies, and are therefore a relatively unexplored population. Despite the typical steep radio spectrum due to synchrotron losses, these fossil electrons are barely visible even at radio frequencies well below a GHz. However, when a pocket of fossil radio plasma is compressed, it boosts the visibility at sub-GHz frequencies, creating so-called radio phoenices. This compression can be the result of bulk motion and shocks in the ICM due to merger activity. In this paper, we demonstrate the discovery potential of low frequency radio sky surveys to find and study revived fossil plasma sources in galaxy clusters. We used the 150~MHz TGSS and 1.4 GHz NVSS sky surveys to identify candidate radio phoenices. A subset of three candidates were studied in detail using deep multi-band radio observations (LOFAR and GMRT), X-ray (\textit{Chandra} or \textit{XMM-Newton}) and archival optical observations. Two of the three sources are new discoveries. Using these observations, we identified common observational properties (radio morphology, ultra-steep spectrum, X-ray luminosity, dynamical state) that will enable us to identify this class of sources more easily, and help to understand the physical origin of these sources.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

We present the first $N$-body simulations that adapt the equations of smoothed particle hydrodynamics to capture the effect of dark matter self-interactions which are too frequent to be resolved explicitly. The relevant energy transfer equations are derived, the appropriate thermal conductivity is determined and the effects of different smoothing kernels are studied. We apply our framework to simulate the formation of isothermal cores in isolated dark matter haloes and determine the core growth rate as a function of the self-scattering cross section. Our approach may be combined with explicit simulations of rare scatterings in order to simulate accurately the effects of arbitrary dark matter self-interactions in future cosmological simulations.

We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.14 \pm 0.01$. We compute the shock strength using density jumps to conclude that the Mach number of the merging subcluster is small ($M \leq 1.15_{-0.09}^{+0.14}$). We also present pseudo-density, projected temperature, pseudo-pressure, and pseudo-entropy maps. Based on the pseudo-entropy map we conclude that the cluster is undergoing a mild merger, consistent with the small Mach number. On the other hand, radio relics extend over Mpc scale in the A3411-12 system, which strongly suggests that a population of energetic electrons already existed over extended regions of the cluster.

With ENZO simulations run on the J\"ulich supercomputers, we have investigated the evolution of magnetic fields in the largest cosmic structures (namely galaxy clusters and filaments connecting them) with unprecedented dynamical range. These simulations revealed the full development of the small-scale dynamo in Eulerian cosmological magneto-hydrodynamical simulations. The turbulent motions developed during the formation of clusters are energetic enough to foster the growth of magnetic fields by several orders of magnitude, starting from weak magnetic fields up strengths of $\sim \rm \mu G$ as observed. Furthermore, shock waves are launched during cluster formation and they are able to accelerate cosmic-ray electrons, that emit in the radio wavelengths. Radio observations of this emission provide information on the local magnetic field strength. We have incorporated, for the first time, the cooling of cosmic-ray electrons when modelling this emission. In this contribution, we present our advances in modelling these physical processes. Here, we mostly focus on the most interesting object in our sample of galaxy clusters, which shows the complexity of magnetic fields and the potential of existing and future multi-wavelengths observations in the study of the weakly collisional plasma on $\sim$ Megaparsecs scales.

Radio emission in the form of giant radio relics is observed at the periphery of galaxy clusters. This non-thermal emission is an important tracer for cosmic-ray electrons and intracluster magnetic fields. One striking observational feature of these objects is their high degree of polarisation which provides information on the magnetic fields at the relics' positions. In this contribution, we test if state-of-the-art high resolution cosmological simulations are able to reproduce the polarisation features of radio relics. Therefore, we present a new analysis of high-resolution cosmological simulations to study the polarisation properties of radio relics in detail. In order to compare our results with current and future radio observations, we create mock radio observations of the diffuse polarised emission from a massive galaxy clusters using six different projections, for different observing frequencies and for different telescopes. Our simulations suggest that, due to the effect of Faraday rotation, it is extremely difficult to relate the morphology of the polarised emission for observing frequencies below $1.4 \ \mathrm{GHz}$ to the real magnetic field structure in relics. We can reproduce the observed degree of polarisation and also several small-scale structures observed in real radio relics, but further work would be needed to reproduce some large-scale spectacular features as observed in real radio relics, such as the "Sausage" and the "Toothbrush" relics.

Our goal is to study the termination of an AGN jet in the young universe and to deduce physical parameters of the jet and the intergalactic medium. We use LOFAR to image the long-wavelength radio emission of the high-redshift blazar S5 0836+710 on arcsecond scales between 120 MHz and 160 MHz. The LOFAR image shows a compact unresolved core and a resolved emission region about 1.5 arcsec to the southwest of the radio core. This structure is in general agreement with previous higher-frequency radio observations with MERLIN and the VLA. The southern component shows a moderately steep spectrum with a spectral index of about $\gtrsim -1$ while the spectral index of the core is flat to slightly inverted. In addition, we detect for the first time a resolved steep-spectrum halo with a spectral index of about $-1$ surrounding the core. The arcsecond-scale radio structure of S5 0836+710 can be understood as an FR II-like radio galaxy observed at a small viewing angle. The southern component can be interpreted as the region of the approaching jet's terminal hotspot and the halo-like diffuse component near the core can be interpreted as the counter-hotspot region. From the differential Doppler boosting of both features, we can derive the hotspot advance speed to $(0.01-0.036)$ c. At a constant advance speed, the derived age of the source would exceed the total lifetime of such a powerful FR II-like radio galaxy substantially. Thus, the hotspot advance speed must have been higher in the past in agreement with a scenario in which the originally highly relativistic jet has lost collimation due to the growth of instabilities and has transformed into an only mildly relativistic flow. Our data suggest that the density of the intergalactic medium around this distant ($z=2.22$) AGN could be substantially higher than the values typically found in less distant FR II radio galaxies.

Radio relics are sites of electron (re)acceleration in merging galaxy clusters but the mechanism of acceleration and the topology of the magnetic field in and near relics are yet to be understood. We are carrying out an observational campaign on double relic galaxy clusters starting with RXC J1314.4-2515. With $Jansky Very Large Array$ multi-configuration observations in the frequency range 1-4 GHz, we perform both spectral and polarization analyses, using the Rotation Measure synthesis technique. We use archival $XMM-Newton$ observations to constrain the properties of the shocked region. We discover a possible connection between the activity of a radio galaxy and the emission of the eastern radio relic. In the northern elongated arc of the western radio relic, we detect polarized emission with an average polarization fraction of $31 \ \%$ at 3 GHz and we derive the Mach number of the underlying X-ray shock. Our observations reveal low levels of fractional polarization and Faraday-complex structures in the southern region of the relic, which point to the presence of thermal gas and filamentary magnetic field morphology inside the radio emitting volume. We measured largely different Rotation Measure dispersion from the two relics. Finally, we use cosmological magneto-hydrodynamical simulations to constrain the magnetic field, viewing angle, and to derive the acceleration efficiency of the shock. We find that the polarization properties of RXC J1314.4-2515 are consistent with a radio relic observed at $70^{\circ}$ with respect to the line of sight and that efficient re-acceleration of fossil electrons has taken place.

Diffuse radio emission associated with the intra-cluster medium (ICM) is observed in a number of merging galaxy clusters. It is currently believed that in mergers a fraction of the kinetic energy is channeled into non-thermal components, such as turbulence, cosmic rays and magnetic fields, that may lead to the formation of giant synchrotron sources in the ICM. Studying merging galaxy clusters in different evolutionary phases is fundamental to understanding the origin of radio emission in the ICM. We observed the nearby galaxy cluster pair RXC J1825.3+3026 ($z\sim0.065$) and CIZA J1824.1+3029 ($z\sim0.071$) at 120-168 MHz with the LOw Frequency ARray (LOFAR) and made use of a deep (240 ks) XMM-Newton dataset to study the non-thermal and thermal properties of the system. RXC J1825.3+3026 is in a complex dynamical state, with a primary on-going merger in the E-W direction and a secondary later stage merger with a group of galaxies in the SW, while CIZA J1824.1+3029 is dynamically relaxed. These two clusters are in a pre-merger phase. We report the discovery of a Mpc-scale radio halo with a low surface brightness extension in RXC J1825.3+3026 that follows the X-ray emission from the cluster center to the remnant of a galaxy group in the SW. This is among the least massive systems and the faintest giant radio halo known to date. Contrary to this, no diffuse radio emission is observed in CIZA J1824.1+3029 nor in the region between the pre-merger cluster pair. The power spectra of the X-ray surface brightness fluctuations of RXC J1825.3+3026 and CIZA J1824.1+3029 are in agreement with the findings for clusters exhibiting a radio halo and the ones where no radio emission has been detected, respectively. We provide quantitative support to the idea that cluster mergers play a crucial role in the generation of non-thermal components in the ICM.

We investigate the possibility of measuring intergalactic magnetic fields using the dispersion measures and rotation measures of fast radio bursts. With Bayesian methods, we produce probability density functions for values of these measures. We distinguish between contributions from the intergalactic medium, the host galaxy and the local environment of the progenitor. To this end, we use constrained, magnetohydrodynamic simulations of the local Universe to compute lines-of-sight integrals from the position of the Milky Way. In particular, we differentiate between predominantly astrophysical and primordial origins of magnetic fields in the intergalactic medium. We test different possible types of host galaxies and probe different distribution functions of fast radio burst progenitor locations inside the host galaxy. Under the assumption that fast radio bursts are produced by magnetars, we use analytic predictions to account for the contribution of the local environment. We find that less than 100 fast radio bursts from magnetars in stellar-wind environments hosted by starburst dwarf galaxies at redshift $z \gtrsim 0.5$ suffice to discriminate between predominantly primordial and astrophysical origins of intergalactic magnetic fields. However, this requires the contribution of the Milky Way to be removed with a precision of $\approx 1 \rm~rad~m^{-2}$. We show the potential existence of a subset of fast radio bursts whose rotation measure carry information on the strength of the intergalactic magnetic field and its origins.

In this Letter, we report the discovery of a radio halo in the high-redshift galaxy cluster PSZ2 G099.86+58.45 ($z=0.616$) with the LOw Frequency ARray (LOFAR) at 120-168 MHz. This is one of the most distant radio halos discovered so far. The diffuse emission extends over $\sim$ 1 Mpc and has a morphology similar to that of the X-ray emission as revealed by XMM-Newton data. The halo is very faint at higher frequencies and is barely detected by follow-up 1-2 GHz Karl G.~Jansky Very Large Array (JVLA) observations, which enable us to constrain the radio spectral index to be $\alpha\leq 1.5-1.6$, i.e.; with properties between canonical and ultra-steep spectrum radio halos. Radio halos are currently explained as synchrotron radiation from relativistic electrons that are re-accelerated in the intra-cluster medium (ICM) by turbulence driven by energetic mergers. We show that in such a framework radio halos are expected to be relatively common at $\sim150$ MHz ($\sim30-60\%$) in clusters with mass and redshift similar to PSZ2 G099.86+58.45; however, at least 2/3 of these radio halos should have steep spectrum and thus be very faint above $\sim 1$ GHz frequencies. Furthermore, since the luminosity of radio halos at high redshift depends strongly on the magnetic field strength in the hosting clusters, future LOFAR observations will also provide vital information on the origin and amplification of magnetic fields in galaxy clusters.

Cosmic rays play a pivotal role in launching galactic winds, particularly in quiescently star-forming galaxies where the hot gas alone is not sufficient to drive a wind. Except for the Milky Way, not much is known about the transport of cosmic rays in galaxies. In this Letter, we present low-frequency observations of the nearby edge-on spiral galaxy NGC 4565 using the LOw-Frequency ARray (LOFAR). With our deep 144-MHz observations, we obtain a clean estimate of the emission originating from old cosmic-ray electrons (CRe), which is almost free from contamination by thermal emission. We measured vertical profiles of the non-thermal radio continuum emission that we fitted with Gaussian and exponential functions. The different profile shapes correspond to 1D cosmic-ray transport models of pure diffusion and advection, respectively. We detect a warp in the radio continuum that is reminiscent of the previously known HI warp. Because the warp is not seen at GHz-frequencies in the radio continuum, its minimum age must be about 100 Myr. The warp also explains the slight flaring of the thick radio disc that can otherwise be well described by a Gaussian profile with an FWHM of 65 arcsec (3.7 kpc). The diffusive radio halo together with the extra-planar X-ray emission may be remnants of enhanced star-forming activity in the past where the galaxy had a galactic wind, as GHz-observations indicate only a weak outflow in the last 40 Myr. NGC 4565 could be in transition from an outflow- to an inflow-dominated phase.

Galaxy clusters are the most massive gravitationally bound structures in the Universe. They grow by accreting smaller structures in a merging process that produces shocks and turbulence in the intra-cluster gas. We observed a ridge of radio emission connecting the merging galaxy clusters Abell 0399 and Abell 0401 with the Low Frequency Array (LOFAR) at 140 MHz. This emission requires a population of relativistic electrons and a magnetic field located in a filament between the two galaxy clusters. We performed simulations to show that a volume-filling distribution of weak shocks may re-accelerate a pre-existing population of relativistic particles, producing emission at radio wavelengths that illuminates the magnetic ridge.

We present LOFAR data from 110--180~MHz of the merging galaxy cluster Abell 1682, alongside archival optical, radio and X-ray data. Our 6 arc-second resolution images at low frequencies reveal new structures associated with numerous radio galaxies in the cluster. At 20 arc-second resolution we see diffuse emission throughout the cluster over hundreds of kpc, indicating particle acceleration mechanisms are in play as a result of the cluster merger event and powerful active galactic nuclei. We show that a significant part of the cluster emission is from an old radio galaxy with very steep spectrum emission (having a spectral index of $\alpha < -2.5$). Furthermore we identify a new region of diffuse steep spectrum emission ($\alpha < -1.1$) as a candidate for a radio halo which is co-spatial with the centre of the cluster merger. We suggest its origin as a population of old and mildly relativistic electrons left over from radio galaxies throughout the cluster which have been re-accelerated to higher energies by shocks and turbulence induced by the cluster merger event. We also note the discovery of six new giant radio galaxies in the vicinity of Abell 1682.

The largest galaxy clusters are observed still to be forming through major cluster-cluster mergers, often showing observational signatures such as radio relics and giant radio haloes. Using LOFAR Two-meter Sky Survey data, we present new detections of both a radio halo (with a spectral index of $\alpha_{143}^{1400}=1.48^{+0.06}_{-0.23}$) and a likely radio relic in Abell 959, a massive cluster at a redshift of z=0.288. Using a sample of clusters with giant radio haloes from the literature (80 in total), we show that the radio halo in A959 lies reasonably well on the scaling relations between the thermal and non-thermal power of the system. Additionally, we find evidence that steep-spectrum haloes tend to reside in clusters with high X-ray luminosities relative to those expected from cluster LM scaling relations, indicating that such systems may preferentially lie at an earlier stage of the merger, consistent with the theory that some steep-spectrum haloes result from low-turbulence mergers. Lastly, we find that halo systems containing radio relics tend to lie at lower X-ray luminosities, relative to those expected from cluster LM scaling relations, for a given halo radio power than those without relics, suggesting that the presence of relics indicates a later stage of the merger, in line with simulations.

Detecting the thermal and non-thermal emission from the shocked cosmic gas surrounding large-scale structures represents a challenge for observations, as well as a unique window into the physics of the warm-hot intergalactic medium. In this work, we present synthetic radio and X-ray surveys of large cosmological simulations in order to assess the chances of jointly detecting the cosmic web in both frequency ranges. We then propose best observing strategies tailored for existing (LOFAR, MWA and XMM) or future instruments (SKA-LOW and SKA-MID, ATHENA and eROSITA). We find that the most promising targets are the extreme peripheries of galaxy clusters in an early merging stage, where the merger causes the fast compression of warm-hot gas onto the virial region. By taking advantage of a detection in the radio band, future deep X-ray observations will probe this gas in emission, and help us to study plasma conditions in the dynamic warm-hot intergalactic medium with unprecedented detail.

Next-generation radio surveys will yield an unprecedented amount of data, warranting analysis by use of machine learning techniques. Convolutional neural networks are the deep learning technique that has proven to be the most successful in classifying image data. Capsule networks are a more recently developed technique that use capsules comprised of groups of neurons, that describe properties of an image including the relative spatial locations of features. The current work explores the performance of different capsule network architectures against simpler convolutional neural network architectures, in reproducing the classifications into the classes of unresolved, FRI and FRII morphologies. We utilise images from a LOFAR survey which is the deepest, wide-area radio survey to date, revealing more complex radio-source structures compared to previous surveys, presenting further challenges for machine learning algorithms. The 4- and 8-layer convolutional networks attain an average precision of 93.3% and 94.3% respectively, compared to 89.7% obtained with the capsule network, when training on original and augmented images. Implementing transfer learning achieves a precision of 94.4%, that is within the confidence interval of the 8-layer convolutional network. The convolutional networks always outperform any variation of the capsule network, as they prove to be more robust to the presence of noise in images. The use of pooling appears to allow more freedom for the intra-class variability of radio galaxy morphologies, as well as reducing the impact of noise.

Recent observations suggest that dense gas clouds can survive even in hot galactic winds. Here we show that the inclusion of turbulent densities with different statistical properties has significant effects on the evolution of wind-swept clouds. We investigate how the initial standard deviation of the log-normal density field influences the dynamics of quasi-isothermal clouds embedded in supersonic winds. We compare uniform, fractal solenoidal, and fractal compressive cloud models in both 3D and 2D hydrodynamical simulations. We find that the processes of cloud disruption and dense gas entrainment are functions of the initial density distribution in the cloud. Fractal clouds accelerate, mix, and are disrupted earlier than uniform clouds. Within the fractal cloud sample, compressive clouds retain high-density nuclei, so they are more confined, less accelerated, and have lower velocity dispersions than their solenoidal counterparts. Compressive clouds are also less prone to Kelvin-Helmholtz and Rayleigh-Taylor instabilities, so they survive longer than solenoidal clouds. By comparing the cloud properties at the destruction time, we find that dense gas entrainment is more effective in uniform clouds than in either of the fractal clouds, and it is more effective in solenoidal than in compressive models. In contrast, mass loading into the wind is more efficient in compressive cloud models than in uniform or solenoidal models. Overall, wide density distributions lead to inefficient entrainment, but they facilitate mass loading and favour the survival of very dense gas in hot galactic winds.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We investigate the evolution of magnetic fields in galaxy clusters starting from constant primordial fields using highly resolved ($\approx \rm 4 ~kpc$) cosmological MHD simulations. The magnetic fields in our sample exhibit amplification via a small-scale dynamo and compression during structure formation. In particular, we study how the spectral properties of magnetic fields are affected by mergers, and we relate the measured magnetic energy spectra to the dynamical evolution of the intracluster medium. The magnetic energy grows by a factor of $\sim$ 40-50 in a time-span of $\sim 9$ Gyr and equipartition between kinetic and magnetic energy occurs on a range of scales ($< 160 \rm ~kpc$ at all epochs) depending on the turbulence state of the system. We also find that, in general, the outer scale of the magnetic field and the MHD scale are not simply correlated in time. The effect of major mergers is to shift the peak magnetic spectra to it smaller scales, whereas the magnetic amplification only starts after $\lesssim$ 1 Gyr. In contrast, continuous minor mergers promote the steady growth of the magnetic field. We discuss the implications of these findings in the interpretation of future radio observations of galaxy clusters.

We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We investigate the impact of different magneto-genesis scenarios, both, primordial and astrophysical, on the propagation of cosmic rays. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements. The distribution of nearby sources causes anisotropy at very high energies, independent of the magnetic field model. We compare our results to the dipole signal measured by the Pierre Auger Observatory. All our models could reproduce the observed dipole amplitude with a pure iron injection composition. This is due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of 4 - 8 EeV.

In a growing number of galaxy clusters diffuse extended radio sources have been found. These sources are not directly associated with individual cluster galaxies. The radio emission reveal the presence of cosmic rays and magnetic fields in the intracluster medium (ICM). We classify diffuse cluster radio sources into radio halos, cluster radio shocks (relics), and revived AGN fossil plasma sources. Radio halo sources can be further divided into giant halos, mini-halos, and possible `intermediate' sources. Halos are generally positioned at cluster center and their brightness approximately follows the distribution of the thermal ICM. Cluster radio shocks (relics) are polarized sources mostly found in the cluster's periphery. They trace merger induced shock waves. Revived fossil plasma sources are characterized by their radio steep-spectra and often irregular morphologies. In this review we give an overview of the properties of diffuse cluster radio sources, with an emphasis on recent observational results. We discuss the resulting implications for the underlying physical acceleration processes that operate in the ICM, the role of relativistic fossil plasma, and the properties of ICM shocks and magnetic fields. We also compile an updated list of diffuse cluster radio sources which will be available on-line http://galaxyclusters.com. We end this review with a discussion on the detection of diffuse radio emission from the cosmic web.

Faraday rotation of polarized emission from pulsars measured at radio frequencies provides a powerful tool to investigate the interstellar and interplanetary magnetic fields. However, besides being sensitive to the astrophysical media, pulsar observations in radio are affected by the highly time-variable ionosphere. In this article, the amount of ionospheric Faraday rotation has been computed by assuming a thin layer model. For this aim, ionospheric maps of the free electron density (based on Global Positioning System data) and semi-empirical geomagnetic models are needed. Through the data of five highly polarized pulsars observed with the individual German LOw-Frequency ARray stations, we investigate the performances of the ionospheric modelling. In addition, we estimate the parameters of the systematics and the correlated noise generated by the residual unmodelled ionospheric effects, and show the comparison of the different free-electron density maps. For the best ionospheric maps, we have found that the rotation measure corrections on one-year timescales after subtraction of diurnal periodicity are accurate to $\sim$ 0.06--0.07 rad m$^{-2}$.

The technique of Faraday tomography is a key tool for the study of magnetised plasmas in the new era of broadband radio polarisation observations. In particular, observations at metre-wavelengths provide significantly better Faraday depth accuracies compared to traditional cm-wavelength observations. However, the effect of Faraday depolarisation makes the polarised signal very challenging to detect at metre wavelengths (MHz frequencies). In this work, Faraday tomography is used to characterise the Faraday rotation properties of polarised sources found in data from the LOFAR Two-Metre Sky Survey (LoTSS). Of the 76 extragalactic polarised sources analysed here, we find that all host a radio-loud AGN. The majority of the sources ($\sim$64%) are large FRII radio galaxies with a median projected linear size of 710 kpc and median radio luminosity at 144 MHz of $4\times10^{26}$ W Hz$^{-1}$ (with $\sim$13% of all sources having a linear size $>1$ Mpc). In several cases, both hotspots are detected in polarisation at an angular resolution of $\sim$20". One such case allowed a study of intergalactic magnetic fields on scales of 3.4 Mpc. Other detected source types include an FRI radio galaxy and at least 8 blazars. Most sources display simple Faraday spectra, however, we highlight one blazar that displays a complex Faraday spectrum, with two close peaks in the Faraday dispersion function.

A number of radio observations have revealed the presence of large synchrotron-emitting sources associated with the intra-cluster medium. There is strong observational evidence that the emitting particles have been (re-)accelerated by shocks and turbulence generated during merger events. The particles that are accelerated are thought to have higher initial energies than those in the thermal pool but the origin of such mildly relativistic particles remains uncertain and needs to be further investigated. The galaxy cluster Abell 1914 is a massive galaxy cluster in which X-ray observations show clear evidence of merging activity. We carried out radio observations of this cluster with the LOw Frequency ARay (LOFAR) at 150 MHz and the Giant Metrewave Radio Telescope (GMRT) at 610 MHz. We also analysed Very Large Array (VLA) 1.4 GHz data, archival GMRT 325 MHz data, CFHT weak lensing data and Chandra observations. Our analysis shows that the ultra-steep spectrum source (4C38.39; $\alpha \lesssim -2$), previously thought to be part of a radio halo, is a distinct source with properties that are consistent with revived fossil plasma sources. Finally, we detect some diffuse emission to the west of the source 4C38.39 that could belong to a radio halo.

Massive, merging galaxy clusters often host giant, diffuse radio sources that arise from shocks and turbulence; hence, radio observations can be useful for determining the merger state of a cluster. In preparation for a larger study, we selected three clusters -- Abell 1319, Abell 1314, and RXC J1501.3+4220 (Z7215) -- making use of the new LOFAR Two-Metre Sky Survey (LoTSS) at 120-168 MHz, and together with archival data, show that these clusters appear to be in pre-merging, merging, and post-merging states, respectively. We argue that Abell 1319 is likely in its pre-merging phase, where three separate cluster components are about to merge. There are no radio halos nor radio relics detected in this system. Abell 1314 is a highly-disturbed, low-mass cluster which is likely in the process of merging. This low-mass system does not show a radio halo, however, we argue that the merger activates mechanisms that cause electron re-acceleration in the large 800 kpc radio tail associated with IC~711. In the cluster Z7215 we discover diffuse radio emission at the cluster center, and we classify this emission as a radio halo, although it is dimmer and smaller than expected by the radio halo power versus cluster mass correlation. We suggest that the disturbed cluster Z7215 is in its post-merging phase. Systematic studies of this kind over a larger sample of clusters observed with LoTSS will help constrain the time scales involved in turbulent re-acceleration and the subsequent energy losses of the underlying electrons.

Collisions of galaxy clusters generate shocks and turbulence in the intra-cluster medium (ICM). The presence of relativistic particles and magnetic fields is inferred through the detection of extended synchrotron radio sources such as haloes and relics and implies that merger shocks and turbulence are capable of (re-)accelerating particles to relativistic energies. However, the precise relationship between merger shocks, turbulence, and extended radio emission is still unclear. Studies of the most simple binary cluster mergers are important to help understand the particle acceleration in the ICM. Our main aim is to study the properties of the extended radio emission and particle acceleration mechanism(s) associated with the generation of relativistic particles in the ICM. We measure the low-frequency radio emission from the merging galaxy cluster Abell 2146 with LOFAR at 144 MHz. We characterize the spectral properties of the radio emission by combining these data with data from archival GMRT at 238 MHz and 612 MHz and VLA at 1.5 GHz. We observe extended radio emission at 144 MHz behind the NW and SE shocks. Across the NW extended source, the spectral index steepens from $-1.06\pm0.06$ to $-1.29\pm0.09$ in the direction of the cluster centre. This spectral behaviour suggests that a relic is associated with the NW upstream shock. The precise nature of the SE extended emission is unclear. It may be a radio halo bounded by a shock or a superposition of a relic and halo. At 144 MHz, we detect a faint emission that was not seen with high-frequency observations, implying a steep ($\alpha<-1.3$) spectrum nature of the bridge emission. Our results imply that the extended radio emission in Abell 2146 is probably associated with shocks and turbulence during cluster merger. The relativistic electrons in the NW and SE may originate from fossil plasma and thermal electrons, respectively.

Radio continuum (RC) emission in galaxies allows us to measure star formation rates (SFRs) unaffected by extinction due to dust, of which the low-frequency part is uncontaminated from thermal (free-free) emission. We calibrate the conversion from the spatially resolved 140 MHz RC emission to the SFR surface density ($\Sigma_{\rm SFR}$) at 1 kpc scale. We used recent observations of three galaxies (NGC 3184, 4736, and 5055) from the LOFAR Two-metre Sky Survey (LoTSS), and archival LOw-Frequency ARray (LOFAR) data of NGC 5194. Maps were created with the facet calibration technique and converted to radio $\Sigma_{\rm SFR}$ maps using the Condon relation. We compared these maps with hybrid $\Sigma_{\rm SFR}$ maps from a combination of GALEX far-ultraviolet and Spitzer 24 $\mu\rm m$ data using plots tracing the relation at $1.2\times 1.2$-kpc$^2$ resolution. The RC emission is smoothed with respect to the hybrid $\Sigma_{\rm SFR}$ owing to the transport of cosmic-ray electrons (CREs). This results in a sublinear relation $(\Sigma_{\rm SFR})_{\rm RC} \propto [(\Sigma_{\rm SFR})_{\rm hyb}]^{a}$, where $a=0.59\pm 0.13$ (140 MHz) and $a=0.75\pm 0.10$ (1365 MHz). Both relations have a scatter of $\sigma = 0.3~\rm dex$. If we restrict ourselves to areas of young CREs ($\alpha > -0.65$; $I_\nu \propto \nu^\alpha$), the relation becomes almost linear at both frequencies with $a\approx 0.9$ and a reduced scatter of $\sigma = 0.2~\rm dex$. We then simulate the effect of CRE transport by convolving the hybrid $\Sigma_{\rm SFR}$ maps with a Gaussian kernel until the RC-SFR relation is linearised; CRE transport lengths are $l=1$-5 kpc. Solving the CRE diffusion equation, we find diffusion coefficients of $D=(0.13$-$1.5) \times 10^{28} \rm cm^2\,s^{-1}$ at 1 GeV. A RC-SFR relation at $1.4$ GHz can be exploited to measure SFRs at redshift $z \approx 10$ using $140$ MHz observations.

Context: New generation low-frequency telescopes are exploring a new parameter space in terms of depth and resolution. The data taken with these interferometers, for example with the LOw Frequency ARray (LOFAR), are often calibrated in a low signal-to-noise ratio regime and the removal of critical systematic effects is challenging. The process requires an understanding of their origin and properties. Aim: In this paper we describe the major systematic effects inherent to next generation low-frequency telescopes, such as LOFAR. With this knowledge, we introduce a data processing pipeline that is able to isolate and correct these systematic effects. The pipeline will be used to calibrate calibrator observations as the first step of a full data reduction process. Methods: We processed two LOFAR observations of the calibrator 3C196: the first using the Low Band Antenna (LBA) system at 42-66 MHz and the second using the High Band Antenna (HBA) system at 115-189 MHz. Results: We were able to isolate and correct for the effects of clock drift, polarisation misalignment, ionospheric delay, Faraday rotation, ionospheric scintillation, beam shape, and bandpass. The designed calibration strategy produced the deepest image to date at 54 MHz. The image has been used to confirm that the spectral energy distribution of the average radio source population tends to flatten at low frequencies. Conclusions: We prove that LOFAR systematic effects can be described by a relatively small number of parameters. Furthermore, the identification of these parameters is fundamental to reducing the degrees of freedom when the calibration is carried out on fields that are not dominated by a strong calibrator.

A number of merging galaxy clusters show the presence of large-scale radio emission associated with the intra-cluster medium (ICM). These synchrotron sources are generally classified as radio haloes and radio relics. Whilst it is commonly accepted that mergers play a crucial role in the formation of radio haloes and relics, not all the merging clusters show the presence of giant diffuse radio sources and this provides important information concerning current models. The Abell 781 complex is a spectacular system composed of an apparent chain of clusters on the sky. Its main component is undergoing a merger and hosts peripheral emission that is classified as a candidate radio relic and a disputed radio halo. We used new LOw Frequency ARay (LOFAR) observations at 143 MHz and archival Giant Metrewave Radio Telescope (GMRT) observations at 325 and 610 MHz to study radio emission from non-thermal components in the ICM of Abell 781. Complementary information came from XMM-Newton data, which allowed us to investigate the connection with the thermal emission and its complex morphology. The origin of the peripheral emission is still uncertain. We speculate that it is related to the interaction between a head tail radio galaxy and shock. However, the current data allow us only to set an upper limit of $\mathcal{M} < 1.4$ on the Mach number of this putative shock. Instead, we successfully characterise the surface brightness and temperature jumps of a shock and two cold fronts in the main cluster component of Abell 781. Their positions suggest that the merger is involving three substructures. We do not find any evidence for a radio halo either at the centre of this system or in the other clusters of the chain. We place an upper limit to the diffuse radio emission in the main cluster of Abell 781 that is a factor of 2 below the current radio power-mass relation for giant radio haloes.

This paper presents a study of the local radio source population, by cross-comparing the data from the first data release (DR1) of the LOFAR Two-Metre Sky Survey (LoTSS) with the Sloan Digital Sky Survey (SDSS) DR7 main galaxy spectroscopic sample. The LoTSS DR1 provides deep data (median rms noise of 71 $\mathrm{\mu}$Jy at 150 MHz) over 424 square degrees of sky, which is sufficient to detect 10615 (32 per cent) of the SDSS galaxies over this sky area. An improved method to separate active galactic nuclei (AGN) accurately from sources with radio emission powered by star formation (SF) is developed and applied, leading to a sample of 2121 local ($z < 0.3$) radio AGN. The local 150 MHz luminosity function is derived for radio AGN and SF galaxies separately, and the good agreement with previous studies at 1.4 GHz suggests that the separation method presented is robust. The prevalence of radio AGN activity is confirmed to show a strong dependence on both stellar and black hole masses, remarkably reaching a fraction of 100 per cent of the most massive galaxies ($> 10^{11} \mathrm{M_{\odot}}$) displaying radio-AGN activity with $L_{\rm 150 MHz} \geq 10^{21}$W Hz$^{-1}$; thus, the most massive galaxies are always switched on at some level. The results allow the full Eddington-scaled accretion rate distribution (a proxy for the duty cycle) to be probed for massive galaxies. More than 50 per cent of the energy is released during the $\le 2$ per cent of the time spent at the highest accretion rates, $L_{\mathrm{mech}}/L_{\mathrm{Edd}} > 10^{-2.5}$. Stellar mass is shown to be a more important driver of radio-AGN activity than black hole mass, suggesting a possible connection between the fuelling gas and the surrounding halo. This result is in line with models in which these radio AGN are essential for maintaining the quenched state of galaxies at the centres of hot gas haloes.

We examine the possible acceleration mechanisms of the relativistic particles responsible for the extended radio emission in Abell 520. We used new LOFAR 145 MHz, archival GMRT 323 MHz and VLA 1.5 GHz data to study the morphological and spectral properties of extended cluster emission. The observational properties are discussed in the framework of particle acceleration models associated with cluster merger turbulence and shocks. In Abell 520, we confirm the presence of extended synchrotron radio emission that has been classified as a radio halo. The comparison between the radio and X-ray brightness suggests that the halo might originate in a cocoon rather than from the central X-ray bright regions of the cluster. The halo spectrum is roughly uniform on the scale of 66 kpc. There is a hint of spectral steepening from the SW edge towards the cluster centre. Assuming DSA, the radio data are suggestive of a shock of $\mathcal{M}_{SW}=2.6_{-0.2}^{+0.3}$ that is consistent with the X-ray derived estimates. This is in line with the scenario in which relativistic electrons in the SW radio edge gain their energies at the shock front via acceleration of either thermal or fossil electrons. We do not detect extended radio emission ahead of the SW shock that is predicted if the emission is the result of adiabatic compression. An X-ray surface brightness discontinuity is detected towards the NE region that may be a counter shock of $\mathcal{M}_{NE}^{X}=1.52\pm0.05$. This is lower than the value predicted from the radio emission ($\mathcal{M}_{NE}=2.1\pm0.2$). Our observations indicate that the SW radio emission in Abell 520 is likely effected by the prominent X-ray detected shock in which radio emitting particles are (re-)accelerated through the Fermi-I mechanism. The NE X-ray discontinuity that is approximately collocated with an edge in the radio emission hints at the presence of a counter shock.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

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.

We introduce a new sample of galaxy clusters obtained from a cosmological simulation covering an unprecedented dynamical range. All the clusters in our sample show a clear signature of small-scale dynamo amplification. We show that it is possible to use dynamo theory for studying the magnetic spectrum in the intracluster medium. We study if the intrinsic variations on the spectra depend on the dynamical history of each cluster or on some host cluster properties.

The joint analysis of the Dispersion and Faraday Rotation Measure from distant, polarised Fast Radio Bursts may be used to put constraints on the origin and distribution of extragalactic magnetic fields on cosmological scales. While the combination of Dispersion and Faraday Rotation Measure can in principle give the average magnetic fields along the line-of-sight, in practice this method must be used with care because it strongly depends on the assumed magnetisation model on large cosmological scales. Our simulations show that the observation of Rotation Measures with $\geq 1-10 ~\rm rad/m^2$ in $\sim 10^2$ Fast Radio Bursts will be able to discriminate between extreme scenarios for the origin of cosmic magnetic fields, independent of the exact distribution of sources with redshift. This represent a strong case for incoming (e.g. ALERT, CHIME) and future (e.g. with the Square Kilometer Array) radio polarisation surveys of the sky.

We have conducted a deep survey (with a central rms of $55\mu\textrm{Jy}$) with the LOw Frequency ARray (LOFAR) at 120-168 MHz of the Bo\"otes field, with an angular resolution of $3.98^{''}\times6.45^{''}$, and obtained a sample of 10091 radio sources ($5\sigma$ limit) over an area of $20\:\textrm{deg}^{2}$. The astrometry and flux scale accuracy of our source catalog is investigated. The resolution bias, incompleteness and other systematic effects that could affect our source counts are discussed and accounted for. The derived 150 MHz source counts present a flattening below sub-mJy flux densities, that is in agreement with previous results from high- and low- frequency surveys. This flattening has been argued to be due to an increasing contribution of star-forming galaxies and faint active galactic nuclei. Additionally, we use our observations to evaluate the contribution of cosmic variance to the scatter in source counts measurements. The latter is achieved by dividing our Bo\"otes mosaic into 10 non-overlapping circular sectors, each one with an approximate area of $2\:\textrm{deg}^{2}.$ The counts in each sector are computed in the same way as done for the entire mosaic. By comparing the induced scatter with that of counts obtained from depth observations scaled to 150MHz, we find that the $1\sigma$ scatter due to cosmic variance is larger than the Poissonian errors of the source counts, and it may explain the dispersion from previously reported depth source counts at flux densities $S<1\,\textrm{mJy}$. This work demonstrates the feasibility of achieving deep radio imaging at low-frequencies with LOFAR.

Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of sight velocity difference between likely merging subclusters. We find that the line-of-sight velocity components of the mergers are typically small compared with the maximum three dimensional relative velocity (usually $<1000$ km s$^{-1}$ and often consistent with zero). This suggests that the merger axes of these systems are generally in or near the plane of the sky matching findings in magneto-hydrodynamical simulations. In 28 of the 29 systems we identify substructures in the galaxy population aligned with the radio relic(s) and presumed associated merger induced shock. From this ensemble, we identify eight systems to include in a `gold' sample that is prime for further observation, modeling, and simulation study. Additional papers will present weak lensing mass maps and dynamical modeling for each merging system, ultimately leading to new insight into a wide range of astrophysical phenomena at some of the largest scales in the universe.

We present results from LOFAR and GMRT observations of the galaxy cluster MACS$\,$J0717.5$+$3745. The cluster is undergoing a violent merger involving at least four sub-clusters, and it is known to host a radio halo. LOFAR observations reveal new sources of radio emission in the Intra-Cluster Medium: (i) a radio bridge that connects the cluster to a head-tail radio galaxy located along a filament of galaxies falling into the main cluster, (ii) a 1.9 Mpc radio arc, that is located North West of the main mass component, (iii) radio emission along the X-ray bar, that traces the gas in the X-rays South West of the cluster centre. We use deep GMRT observations at 608 MHz to constrain the spectral indices of these new radio sources, and of the emission that was already studied in the literature at higher frequency. We find that the spectrum of the radio halo and of the relic at LOFAR frequency follows the same power law as observed at higher frequencies. The radio bridge, the radio arc, and the radio bar all have steep spectra, which can be used to constrain the particle acceleration mechanisms. We argue that the radio bridge could be caused by the re-acceleration of electrons by shock waves that are injected along the filament during the cluster mass assembly. Despite the sensitivity reached by our observations, the emission from the radio halo does not trace the emission of the gas revealed by X-ray observations. We argue that this could be due to the difference in the ratio of kinetic over thermal energy of the intra-cluster gas, suggested by X-ray observations.

We present LOFAR $120-168$ MHz images of the merging galaxy cluster Abell 1240 that hosts double radio relics. In combination with the GMRT $595-629$ MHz and VLA $2-4$ GHz data, we characterised the spectral and polarimetric properties of the radio emission. The spectral indices for the relics steepen from their outer edges towards the cluster centre and the electric field vectors are approximately perpendicular to the major axes of the relics. The results are consistent with the picture that these relics trace large-scale shocks propagating outwards during the merger. Assuming diffusive shock acceleration (DSA), we obtain shock Mach numbers of $\mathcal{M}=2.4$ and $2.3$ for the northern and southern shocks, respectively. For $\mathcal{M}\lesssim3$ shocks, a pre-existing population of mildly relativistic electrons is required to explain the brightness of the relics due to the high ($>10$ per cent) particle acceleration efficiency required. However, for $\mathcal{M}\gtrsim4$ shocks the required efficiency is $\gtrsim1\%$ and $\gtrsim0.5\%$, respectively, which is low enough for shock acceleration directly from the thermal pool. We used the fractional polarization to constrain the viewing angle to $\geqslant(53\pm3)^\circ$ and $\geqslant(39\pm5)^\circ$ for the northern and southern shocks, respectively. We found no evidence for diffuse emission in the cluster central region. If the halo spans the entire region between the relics ($\sim1.8\,\text{Mpc}$) our upper limit on the power is $P_\text{1.4 GHz}=(1.4\pm0.6)\times10^{23}\,\text{W}\text{Hz}^{-1}$ which is approximately equal to the anticipated flux from a cluster of this mass. However, if the halo is smaller than this, our constraints on the power imply that the halo is underluminous.

Radio halos and radio relics are diffuse synchrotron sources that extend over Mpc-scales and are found in a number of merger galaxy clusters. They are believed to form as a consequence of the energy that is dissipated by turbulence and shocks in the intra-cluster medium (ICM). However, the precise physical processes that generate these steep synchrotron spectrum sources are still poorly constrained. We present a new LOFAR observation of the double galaxy cluster Abell 1758. This system is composed of A1758N, a massive cluster hosting a known giant radio halo, and A1758S, which is a less massive cluster whose diffuse radio emission is confirmed here for the first time. Our observations have revealed a radio halo and a candidate radio relic in A1758S, and a suggestion of emission along the bridge connecting the two systems which deserves confirmation. We combined the LOFAR data with archival VLA and GMRT observations to constrain the spectral properties of the diffuse emission. We also analyzed a deep archival Chandra observation and used this to provide evidence that A1758N and A1758S are in a pre-merger phase. The ICM temperature across the bridge that connects the two systems shows a jump which might indicate the presence of a transversal shock generated in the initial stage of the merger.

Extended, steep-spectrum radio synchrotron sources are pre-dominantly found in massive galaxy clusters as opposed to groups. LOFAR Two-Metre Sky Survey images have revealed a diffuse, ultra-steep spectrum radio source in the low-mass cluster Abell 1931. The source has a fairly irregular morphology with a largest linear size of about 550 kpc. The source is only seen in LOFAR observations at 143 MHz and GMRT observations at 325 MHz. The spectral index of the total source between 143 MHz and 325 MHz is $\alpha_{143}^{325} = -2.86 \pm 0.36$. The source remains invisible in Very Large Array (1-2 GHz) observations as expected given the spectral index. Chandra X-ray observations of the cluster revealed a bolometric luminosity of $L_X = (1.65 \pm 0.39) \times 10^{43}$ erg s$^{-1}$ and a temperature of $2.92_{-0.87}^{+1.89}$ keV which implies a mass of around $\sim 10^{14} M_{\odot}$. We conclude that the source is a remnant radio galaxy that has shut off around 200 Myr ago. The brightest cluster galaxy, a radio-loud elliptical galaxy, could be the source for this extinct source. Unlike remnant sources studied in the literature, our source has a steep spectrum at low radio frequencies. Studying such remnant radio galaxies at low radio frequencies is important for understanding the scarcity of such sources and their role in feedback processes.

Merging galaxy clusters produce low Mach number shocks in the intracluster medium. These shocks can accelerate electrons to relativistic energies that are detectable at radio frequencies. MACS J0744.9+3927 is a massive ($M_{500} = (11.8 \pm 2.8) \times 10^{14} M_{\odot}$), high-redshift ($z=0.6976$) cluster where a Bullet-type merger is presumed to have taken place. Sunyaev-Zel'dovich maps from MUSTANG indicate that a shock, with Mach number $\mathcal{M} = 1.0-2.9$ and an extension of $\sim 200$ kpc, sits near the centre of the cluster. The shock is also detected as a brightness and temperature discontinuity in X-ray observations. To search for diffuse radio emission associated with the merger, we have imaged the cluster with the LOw Frequency ARray (LOFAR) at 120-165 MHz. Our LOFAR radio images reveal previously undetected AGN emission, but do not show clear cluster-scale diffuse emission in the form of a radio relic nor a radio halo. The region of the shock is on the western edge of AGN lobe emission from the brightest cluster galaxy. Correlating the flux of known shock-induced radio relics versus their size, we find that the radio emission overlapping the shocked region in MACS J0744.9+3927 is likely of AGN origin. We argue against the presence of a relic caused by diffusive shock acceleration and suggest that the shock is too weak to accelerate electrons from the intracluster medium.

Low-mass galaxies are subject to strong galactic outflows, in which cosmic rays may play an important role, they can be best traced with low-frequency radio continuum observations, which are less affected by spectral ageing. We present a study of the nearby star burst dwarf irregular galaxy IC 10 using observations at 140 MHz with the LOw-Frequency ARray (LOFAR), at 1580 MHz with the Very Large Array (VLA) and at 6200 MHz with the VLA and the 100-m Effelsberg telescope. We find that IC 10 has a low-frequency radio halo, which manifests itself as a second component (thick disc) in the minor axis profiles of the non-thermal radio continuum emission at 140 and 1580 MHz. These profiles are then fitted with 1D cosmic-ray transport models for pure diffusion and advection. We find that a diffusion model fits best, with a diffusion coefficient of $D=(0.4$-$0.8) \times 10^{26}(E/{\rm GeV})^{0.5}~{\rm cm^2\,s^{-1}}$, which is at least an order of magnitude smaller than estimates both from anisotropic diffusion and the diffusion length. In contrast, advection models, which cannot be ruled out due to the mild inclination, while providing poorer fits, result in advection speeds close to the escape velocity of $\approx$$50~\rm km\,s^{-1}$, as expected for a cosmic-ray driven wind. Our favoured model with an accelerating wind provides a self-consistent solution, where the magnetic field is in energy equipartition with both the warm neutral and warm ionized medium with an important contribution from cosmic rays. Consequently, cosmic rays can play a vital role for the launching of galactic winds in the disc--halo interface.

Machine learning techniques have been increasingly useful in astronomical applications over the last few years, for example in the morphological classification of galaxies. Convolutional neural networks have proven to be highly effective in classifying objects in image data. The current work aims to establish when multiple components are present, in the astronomical context of synthesis imaging observations of radio sources. To this effect, we design a convolutional neural network to differentiate between different morphology classes using sources from the Radio Galaxy Zoo (RGZ) citizen science project. In this first step, we focus on exploring the factors that affect the performance of such neural networks, such as the amount of training data, number and nature of layers and the hyperparameters. We begin with a simple experiment in which we only differentiate between two extreme morphologies, using compact and multiple component extended sources. We found that a three convolutional layer architecture yielded very good results, achieving a classification accuracy of 97.4% on a test data set. The same architecture was then tested on a four-class problem where we let the network classify sources into compact and three classes of extended sources, achieving a test achieving a test accuracy of 93.5%. The best-performing convolutional neural network setup has been verified against RGZ Data Release 1 where a final test accuracy of 94.8% was obtained, using both original and augmented images. The use of sigma clipping does not offer a significant benefit overall, except in cases with a small number of training images.

We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We use constrained initial conditions with the cosmological magnetohydrodynamics code {\sc ENZO}. The resulting models of the distribution of magnetic fields in the local Universe are used in the \crpropa code to simulate the propagation of ultra-high energy cosmic rays. We investigate the impact of six different magneto-genesis scenarios, both primordial and astrophysical, on the propagation of cosmic rays over cosmological distances. Moreover, we study the influence of different source distributions around the Milky Way. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements of ultra-high energy cosmic rays. However, at high energies above the GZK-limit, there is anisotropy caused by the distribution of nearby sources, independent of the magnetic field model. This provides a chance to identify cosmic ray sources with future full-sky measurements and high number statistics at the highest energies. Finally, we compare our results to the dipole signal measured by the Pierre Auger Observatory. All our source models and magnetic field models could reproduce the observed dipole amplitude with a pure iron injection composition. Our results indicate that the dipole is observed due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of $4-8 \rm\ EeV$.

We report L-band VLA observations of 1RXS J0603.3+4214, a cluster that hosts a bright radio relic, known as the Toothbrush, and an elongated giant radio halo. These new observations allow us to study the surface brightness distribution down to one arcsec resolution with very high sensitivity. Our images provide an unprecedented detailed view of the Toothbrush, revealing enigmatic filamentary structures. To study the spectral index distribution, we complement our analysis with published LOFAR and GMRT observations. The bright `brush' of the Toothbrush shows a prominent narrow ridge to its north with a sharp outer edge. The spectral index at the ridge is in the range $-0.70\leq\alpha\leq-0.80$. We suggest that the ridge is caused by projection along the line of sight. With a simple toy model for the smallest region of the ridge, we conclude that the magnetic field is below $5\,\rm\mu G$ and varies significantly across the shock front. Our model indicates that the actual Mach number is higher than that obtained from the injection index and agrees well with the one derived from the overall spectrum, namely ${\cal M}=3.78^{+0.3}_{-0.2}$. The radio halo shows an average spectral index of $\alpha=-1.16\pm0.05$ and a slight gradient from north to south. The southernmost part of the halo is steeper and possibly related to a shock front. Excluding the southernmost part, the halo morphology agrees very well with the X-ray morphology. A power-law correlation is found between the radio and X-ray surface brightness

Faraday rotation and synchrotron emission from extragalactic radio sources give evidence for the presence of magnetic fields extending over ~Mpc scales. However, the origin of these fields remains elusive. With new high-resolution grid simulations we studied the growth of magnetic fields in a massive galaxy cluster that in several aspects is similar to the Coma cluster. We investigated models in which magnetic fields originate from primordial seed fields with comoving strengths of 0.1 nG at redshift z=30. The simulations show evidence of significant magnetic field amplification. At the best spatial resolution (3.95 kpc), we are able to resolve the scale where magnetic tension balances the bending of magnetic lines by turbulence. This allows us to observe the final growth stage of the small-scale dynamo. To our knowledge this is the first time that this is seen in cosmological simulations of the intracluster medium. Our mock observations of Faraday Rotation provide a good match to observations of the Coma cluster. However, the distribution of magnetic fields shows strong departures from a simple Maxwellian distribution, suggesting that the three-dimensional structure of magnetic fields in real clusters may be significantly different than what is usually assumed when inferring magnetic field values from rotation measure observations.

The origin of extragalactic magnetic fields is still poorly understood. Based on a dedicated suite of cosmological magneto-hydrodynamical simulations with the ENZO code we have performed a survey of different models that may have caused present-day magnetic fields in galaxies and galaxy clusters. The outcomes of these models differ in cluster outskirts, filaments, sheets and voids and we use these simulations to find observational signatures of magnetogenesis. With these simulations, we predict the signal of extragalactic magnetic fields in radio observations of synchrotron emission from the cosmic web, in Faraday Rotation, in the propagation of Ultra High Energy Cosmic Rays, in the polarized signal from Fast Radio Bursts at cosmological distance and in spectra of distant blazars. In general, primordial scenarios in which present-day magnetic fields originate from the amplification of weak (<nG) uniform seed fields result more homogeneous and relatively easier to observe magnetic fields than than astrophysical scenarios, in which present-day fields are the product of feedback processes triggered by stars and active galaxies. In the near future the best evidence for the origin of cosmic magnetic fields will most likely come from a combination of synchrotron emission and Faraday Rotation observed at the periphery of large-scale structures.

Feedback by radio-loud active galactic nuclei (AGN) in galaxy groups is not fully understood. Open questions include the duty cycle of the AGN, the spatial extent of the radio lobes, the effect they have on the intragroup medium, and the fate of the cosmic rays. We present the discovery of a 650 kpc-radio galaxy embedded in steep diffuse emission at $z = 0.18793 \pm 5 \times 10^{-5}$ located at the center of the galaxy group MaxBCG J199.31832+51.72503 using an observation from the LOFAR Two-meter Sky Survey (LoTSS) at the central frequency of 144 MHz. Subsequently, we performed a GMRT observation at the central frequency of 607 MHz to study the spectral properties of the source. The observations reveal a radio galaxy with a total radio power $P_{\rm tot, 1.4} \sim 2.1 \times 10^{24}$ W Hz$^{-1}$, exhibiting two asymmetrical jets and lobes. The derived spectral index map shows a steepening toward the inner regions and a steep-spectrum core region. We model the integrated radio spectrum, providing two possible interpretations: the radio source is evolved but still active or it is just at the end of its active phase. Finally, in the same field of view we have discovered Mpc-sized emission surrounding a close pair of AGN located at a redshift $z = 0.0587 \pm 2 \times 10^{-4}$ (SDSS J131544.56+521213.2 and SDSS J131543.99+521055.7) which could be a radio remnant source.

We present a deep, low-frequency radio continuum study of the nearby Fanaroff--Riley class I (FR I) radio galaxy 3C 31 using a combination of LOw Frequency ARray (LOFAR; 30--85 and 115--178 MHz), Very Large Array (VLA; 290--420 MHz), Westerbork Synthesis Radio Telescope (WSRT; 609 MHz) and Giant Metre Radio Telescope (GMRT; 615 MHz) observations. Our new LOFAR 145-MHz map shows that 3C 31 has a largest physical size of $1.1$ Mpc in projection, which means 3C 31 now falls in the class of giant radio galaxies. We model the radio continuum intensities with advective cosmic-ray transport, evolving the cosmic-ray electron population and magnetic field strength in the tails as functions of distance to the nucleus. We find that if there is no in-situ particle acceleration in the tails, then decelerating flows are required that depend on radius $r$ as $v\propto r^{\beta}$ ($\beta\approx -1$). This then compensates for the strong adiabatic losses due to the lateral expansion of the tails. We are able to find self-consistent solutions in agreement with the entrainment model of Croston & Hardcastle, where the magnetic field provides $\approx$$1/3$ of the pressure needed for equilibrium with the surrounding intra-cluster medium (ICM). We obtain an advective time-scale of $\approx$$190$ Myr, which, if equated to the source age, would require an average expansion Mach number ${\cal M} \approx 5$ over the source lifetime. Dynamical arguments suggest that instead, either the outer tail material does not represent the oldest jet plasma or else the particle ages are underestimated due to the effects of particle acceleration on large scales.

LOFAR observations at 144 MHz have revealed large-scale radio sources in the unrelaxed galaxy cluster Abell 1132. The cluster hosts diffuse radio emission on scales of $\sim$650 kpc near the cluster center and a head-tail (HT) radio galaxy, extending up to 1 Mpc, South of the cluster center. The central diffuse radio emission is not seen in NVSS, FIRST, WENSS, nor in C & D array VLA observations at 1.4 GHz, but is detected in our follow-up GMRT observations at 325 MHz. Using LOFAR and GMRT data, we determine the spectral index of the central diffuse emission to be $\alpha=-1.75\pm0.19$ ($S\propto\nu^{\alpha}$). We classify this emission as an ultra-steep spectrum radio halo and discuss the possible implications for the physical origin of radio halos. The HT radio galaxy shows narrow, collimated emission extending up to 1 Mpc and another 300 kpc of more diffuse, disturbed emission, giving a full projected linear size of 1.3 Mpc - classifying it as a giant radio galaxy (GRG) and making it the longest HT found to date. The head of the GRG coincides with an elliptical galaxy (SDSS J105851.01$+$564308.5) belonging to Abell 1132. In our LOFAR image, there appears to be a connection between the radio halo and the GRG. The turbulence that may have produced the halo may have also affected the tail of the GRG. In turn, the GRG may have provided seed electrons for the radio halo.

We present a study of AGN feedback at higher redshifts ($0.3<z<1.2$) using Sunyaev-Zel'dovich (SZ) selected samples of clusters from the South-Pole Telescope and Atacama Cosmology Telescope surveys. In contrast to studies of nearby systems, we do not find a separation between cooling flow clusters and non-cooling flow clusters based on the radio luminosity of the central radio source. This lack may be due to the increased incidence of galaxy-galaxy mergers at higher redshift that triggers AGN activity. In support of this scenario, we find evidence for evolution in the radio luminosity function of the central radio source: while the lower-luminosity sources do not evolve much, the higher-luminosity sources show a strong increase in the frequency of their occurrence at higher redshifts. We interpret this evolution as an increase in high-excitation radio galaxies (HERGs) in massive clusters at $z>0.6$, implying a transition from HERG-mode accretion to lower-power low-excitation radio galaxy (LERG)-mode accretion at intermediate redshifts. Additionally, we use local radio-to-jet power scaling relations to estimate feedback power and find that half of the cooling flow systems in our sample probably have enough heating to balance cooling. However, we postulate that the local relations are likely not well suited to predict feedback power in high-luminosity HERGs, as they are derived from samples composed mainly of lower-luminosity LERGs.

Cluster mergers leave distinct signatures in the ICM in the form of shocks and diffuse cluster radio sources that provide evidence for the acceleration of relativistic particles. However, the physics of particle acceleration in the ICM is still not fully understood. Here we present new 1-4 GHz Jansky Very Large Array (VLA) and archival Chandra observations of the HST Frontier Fields Cluster Abell 2744. In our new VLA images, we detect the previously known $\sim2.1$ Mpc radio halo and $\sim1.5$ Mpc radio relic. We carry out a radio spectral analysis from which we determine the relic's injection spectral index to be $\alpha_{\rm{inj}} = -1.12 \pm 0.19$. This corresponds to a shock Mach number of $\mathcal{M}$ = 2.05$^{+0.31}_{-0.19}$ under the assumption of diffusive shock acceleration. We also find evidence for spectral steepening in the post-shock region. We do not find evidence for a significant correlation between the radio halo's spectral index and ICM temperature. In addition, we observe three new polarized diffuse sources and determine two of these to be newly discovered giant radio relics. These two relics are located in the southeastern and northwestern outskirts of the cluster. The corresponding integrated spectral indices measure $-1.81 \pm 0.26$ and $-0.63 \pm 0.21$ for the SE and NW relics, respectively. From an X-ray surface brightness profile we also detect a possible density jump of $R=1.39^{+0.34}_{-0.22}$ co-located with the newly discovered SE relic. This density jump would correspond to a shock front Mach number of $\mathcal{M}=1.26^{+0.25}_{-0.15}$.

Turbulence generated by large-scale motions during structure formation affects the evolution of the thermal and non-thermal components of the intracluster medium. As enstrophy is a measure of the magnitude of vorticity, we study the generation and evolution of turbulence by analysing the Lagrangian history of enstrophy. For this purpose we combine cosmological simulations carried out with the ENZO-code with our Lagrangian post-processing tool CRaTer. This way we are able to quantify the individual source terms of enstrophy in the course of the accretion of groups onto galaxy clusters. Here we focus on the redshift range from $z=1$ to $z=0$. Finally, we measure the rate of dissipation of turbulence and estimate the resulting amplification of intracluster magnetic fields. We find that compressive and baroclinic motions are the main sources of enstrophy, while stretching motions and dissipation affect most of the ensuing enstrophy evolution. The rate of turbulent dissipation is able to sustain the amplification of intracluster magnetic fields to observed levels.

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.

Pairs of radio relics are believed to form during cluster mergers, and are best observed when the merger occurs in the plane of the sky. Mergers can also produce radio halos, through complex processes likely linked to turbulent re-acceleration of cosmic-ray electrons. However, only some clusters with double relics also show a radio halo. Here, we present a novel method to derive upper limits on the radio halo emission, and analyse archival X-ray Chandra data, as well as galaxy velocity dispersions and lensing data, in order to understand the key parameter that switches on radio halo emission. We place upper limits on the halo power below the $P_{\rm 1.4 \, GHz}\, M_{500}$ correlation for some clusters, confirming that clusters with double relics have different radio properties. Computing X-ray morphological indicators, we find that clusters with double relics are associated with the most disturbed clusters. We also investigate the role of different mass-ratios and time-since-merger. Data do not indicate that the merger mass ratio has an impact on the presence or absence of radio halos (the null hypothesis that the clusters belong to the same group cannot be rejected). However, the data suggests that the absence of radio halos could be associated with early and late mergers, but the sample is too small to perform a statistical test. Our study is limited by the small number of clusters with double relics. Future surveys with LOFAR, ASKAP, MeerKat and SKA will provide larger samples to better address this issue.

We present a study of the low-frequency radio properties of star forming (SF) galaxies and active galactic nuclei (AGN) up to redshift $z=2.5$. The new spectral window probed by the Low Frequency Array (LOFAR) allows us to reconstruct the radio continuum emission from 150 MHz to 1.4 GHz to an unprecedented depth for a radio-selected sample of $1542$ galaxies in $\sim 7~ \rm{deg}^2$ of the LOFAR Bo\"otes field. Using the extensive multi-wavelength dataset available in Bo\"otes and detailed modelling of the FIR to UV spectral energy distribution (SED), we are able to separate the star-formation (N=758) and the AGN (N=784) dominated populations. We study the shape of the radio SEDs and their evolution across cosmic time and find significant differences in the spectral curvature between the SF galaxy and AGN populations. While the radio spectra of SF galaxies exhibit a weak but statistically significant flattening, AGN SEDs show a clear trend to become steeper towards lower frequencies. No evolution of the spectral curvature as a function of redshift is found for SF galaxies or AGN. We investigate the redshift evolution of the infrared-radio correlation (IRC) for SF galaxies and find that the ratio of total infrared to 1.4 GHz radio luminosities decreases with increasing redshift: $ q_{\rm 1.4GHz} = (2.45 \pm 0.04) \times (1+z)^{-0.15 \pm 0.03} $. Similarly, $q_{\rm 150MHz}$ shows a redshift evolution following $ q_{\rm 150GHz} = (1.72 \pm 0.04) \times (1+z)^{-0.22 \pm 0.05}$. Calibration of the 150 MHz radio luminosity as a star formation rate tracer suggests that a single power-law extrapolation from $q_{\rm 1.4GHz}$ is not an accurate approximation at all redshifts.

In this paper, the second in a series investigating FR II radio galaxies at low frequencies, we use LOFAR and VLA observations between 117 and 456 MHz in addition to archival data to determine the dynamics and energetics of two radio galaxies, 3C452 and 3C223, through fitting of spectral ageing models on small spatial scales. We provide improved measurements for the physical extent of the two sources, including a previously unknown low surface brightness extension to the northern lobe of 3C223, and revised energetics based on these values. We find spectral ages of $77.05^{+9.22}_{-8.74}$ and $84.96^{+15.02}_{-13.83}$ Myr for 3C452 and 3C223 respectively suggesting a characteristic advance speed for the lobes of around one per cent the speed of light. For 3C452 we show that, even for a magnetic field strength not assumed to be in equipartition, a disparity of factor of approximately 2 exists between the spectral age and that determined from a dynamical standpoint. We confirm that the injection index of both sources (as derived from the lobe emission) remains steeper than classically assumed values even when considered on well resolved scales at low frequencies, but find an unexpected sharp discontinuity between the spectrum of the hotspots and the surrounding lobe emission. We suggest that this discrepancy is due to the absorption of hotspot emission and/or non-homogeneous and additional acceleration mechanisms and, as such, hotspots should not be used in the determination of the underlying initial electron energy distribution.

We have observed the steep spectrum radio source B2 0924+30 using the LOw Frequency ARray (LOFAR) telescope. Hosted by a z=0.026 elliptical galaxy, it has a relatively large angular size of 12' (corresponding to 360 kpc projected linear size) and a morphology reminiscent of a remnant Fanaroff-Riley type II (FRII) radio galaxy. We combine LOFAR and archival images to study the spectral properties at a spatial resolution of 1'. We derive low frequency spectral index maps and use synchrotron ageing models to infer ages for different regions of the source. Using radiative ageing model fitting we show that the AGN activity ceased around 50 Myr ago. We note that the outer regions of the lobes are younger than the inner regions which is interpreted as a sign that those regions are remnant hotspots. We demonstrate the usefulness of maps of AGN radio remnants taken at low frequencies and suggest caution over the interpretation of spectral ages derived from integrated flux density measurements versus age mapping. The spectral index properties as well as the derived ages of B2 0924+30 are consistent with it being an FRII AGN radio remnant. LOFAR data are proving to be instrumental in extending our studies to the lowest radio frequencies and enabling analysis of the oldest source regions.

To investigate the relationship between thermal and non-thermal components in merger galaxy clusters, we present deep JVLA and Chandra observations of the HST Frontier Fields cluster MACS J0717.5+3745. The Chandra image shows a complex merger event, with at least four components belonging to different merging subclusters. NW of the cluster, $\sim 0.7$ Mpc from the center, there is a ram-pressure-stripped core that appears to have traversed the densest parts of the cluster after entering the ICM from the direction of a galaxy filament to the SE. We detect a density discontinuity NNE of this core which we speculate is associated with a cold front. Our radio images reveal new details for the complex radio relic and radio halo in this cluster. In addition, we discover several new filamentary radio sources with sizes of 100-300 kpc. A few of these seem to be connected to the main radio relic, while others are either embedded within the radio halo or projected onto it. A narrow-angled-tailed (NAT) radio galaxy, a cluster member, is located at the center of the radio relic. The steep spectrum tails of this AGN leads into the large radio relic where the radio spectrum flattens again. This morphological connection between the NAT radio galaxy and relic provides evidence for re-acceleration (revival) of fossil electrons. The presence of hot $\gtrsim 20$ keV ICM gas detected by Chandra near the relic location provides additional support for this re-acceleration scenario.

On the largest scales, the Universe consists of voids and filaments making up the cosmic web. Galaxy clusters are located at the knots in this web, at the intersection of filaments. Clusters grow through accretion from these large-scale filaments and by mergers with other clusters and groups. In a growing number of galaxy clusters, elongated Mpc-size radio sources have been found, so-called radio relics. These relics are thought to trace relativistic electrons in the intracluster plasma accelerated by low-Mach number collisionless shocks generated by cluster-cluster merger events. A long-standing problem is how low-Mach number shocks can accelerate electrons so efficiently to explain the observed radio relics. Here we report on the discovery of a direct connection between a radio relic and a radio galaxy in the merging galaxy cluster Abell 3411-3412. This discovery indicates that fossil relativistic electrons from active galactic nuclei are re-accelerated at cluster shocks. It also implies that radio galaxies play an important role in governing the non-thermal component of the intracluster medium in merging clusters.

Observations of diffuse radio emission in galaxy clusters indicate that cosmic-ray electrons are accelerated on $\sim$Mpc scales. However, protons appear to be accelerated less efficiently since their associated hadronic $\gamma$-ray emission has not yet been detected. Inspired by recent particle-in-cell simulations, we study the cosmic-ray production and its signatures under the hypothesis that the efficiency of shock acceleration depends on the Mach number and on the shock obliquity. For this purpose, we combine ENZO cosmological magneto-hydrodynamical simulations with a Lagrangian tracer code to follow the properties of the cosmic rays. Our simulations suggest that the distribution of obliquities in galaxy clusters is random to first order. Quasi-perpendicular shocks are able to accelerate cosmic-ray electrons to the energies needed to produce observable radio emission. However, the $\gamma$-ray emission is lowered by a factor of a few, $\sim$3, if cosmic-ray protons are only accelerated by quasi-parallel shocks, reducing (yet not entirely solving) the tension with the non-detection of hadronic $\gamma$-ray emission by the \textit{Fermi}-satellite.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

1) Background: the budget of non-thermal energy in galaxy clusters is not well constrained, owing to the observational and theoretical difficulties in studying these diluted plasmas on large scales. 2) Method: we use recent cosmological simulations with complex physics in order to connect the emergence of non-thermal energy to the underlying evolution of gas and dark matter. 3) Results: the impact of non-thermal energy (e.g. cosmic rays, magnetic fields and turbulent motions) is found to increase in the outer region of galaxy clusters. Within numerical and theoretical uncertainties, turbulent motions dominate the budget of non-thermal energy in most of the cosmic volume. 4) Conclusion: assessing the distribution non-thermal energy in galaxy clusters is crucial to perform high-precision cosmology in the future. Constraining the level of non-thermal energy in cluster outskirts will improve our understanding of the acceleration of relativistic particles by cosmic shocks and of the origin of extragalactic magnetic fields.

Weak shocks in the intracluster medium may accelerate cosmic-ray protons and cosmic-ray electrons differently depending on the angle between the upstream magnetic field and the shock normal. In this work, we investigate how shock obliquity affects the production of cosmic rays in high-resolution simulations of galaxy clusters. For this purpose, we performed a magneto-hydrodynamical simulation of a galaxy cluster using the mesh refinement code \enzo. We use Lagrangian tracers to follow the properties of the thermal gas, the cosmic rays and the magnetic fields over time. We tested a number of different acceleration scenarios by varying the obliquity-dependent acceleration efficiencies of protons and electrons, and by examining the resulting hadronic $\gamma$-ray and radio emission. We find that the radio emission does not change significantly if only quasi-perpendicular shocks are able to accelerate cosmic-ray electrons. Our analysis suggests that radio emitting electrons found in relics have been typically shocked many times before $z=0$. On the other hand, the hadronic $\gamma$-ray emission from clusters is found to decrease significantly if only quasi-parallel shocks are allowed to accelerate cosmic-ray protons. This might reduce the tension with the low upper limits on $\gamma$-ray emission from clusters set by the \textit{Fermi}-satellite.

Turbulence is a key ingredient for the evolution of the intracluster medium, whose properties can be predicted with high resolution numerical simulations. We present initial results on the generation of solenoidal and compressive turbulence in the intracluster medium during the formation of a small-size cluster using highly resolved, non-radiative cosmological simulations, with a refined monitoring in time. In this first of a series of papers, we closely look at one simulated cluster whose formation was distinguished by a merger around $z \sim 0.3$. We separate laminar gas motions, turbulence and shocks with dedicated filtering strategies and distinguish the solenoidal and compressive components of the gas flows using Hodge-Helmholtz decomposition. Solenoidal turbulence dominates the dissipation of turbulent motions ($\sim 95\%$) in the central cluster volume at all epochs. The dissipation via compressive modes is found to be more important ($\sim 30 \%$ of the total) only at large radii ($\geq 0.5 ~r_{\rm vir}$) and close to merger events. We show that enstrophy (vorticity squared) is good proxy of solenoidal turbulence. All terms ruling the evolution of enstrophy (i.e. baroclinic, compressive, stretching and advective terms) are found to be significant, but in amounts that vary with time and location. Two important trends for the growth of enstrophy in our simulation are identified: first, enstrophy is continuously accreted into the cluster from the outside, and most of that accreted enstrophy is generated near the outer accretion shocks by baroclinic and compressive processes. Second, in the cluster interior vortex stretching is dominant, although the other terms also contribute substantially.

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

The Lockman Hole is a well-studied extragalactic field with extensive multi-band ancillary data covering a wide range in frequency, essential for characterising the physical and evolutionary properties of the various source populations detected in deep radio fields (mainly star-forming galaxies and AGNs). In this paper we present new 150-MHz observations carried out with the LOw Frequency ARray (LOFAR), allowing us to explore a new spectral window for the faint radio source population. This 150-MHz image covers an area of 34.7 square degrees with a resolution of 18.6$\times$14.7 arcsec and reaches an rms of 160 $\mu$Jy beam$^{-1}$ at the centre of the field. As expected for a low-frequency selected sample, the vast majority of sources exhibit steep spectra, with a median spectral index of $\alpha_{150}^{1400}=-0.78\pm0.015$. The median spectral index becomes slightly flatter (increasing from $\alpha_{150}^{1400}=-0.84$ to $\alpha_{150}^{1400}=-0.75$) with decreasing flux density down to $S_{150} \sim$10 mJy before flattening out and remaining constant below this flux level. For a bright subset of the 150-MHz selected sample we can trace the spectral properties down to lower frequencies using 60-MHz LOFAR observations, finding tentative evidence for sources to become flatter in spectrum between 60 and 150 MHz. Using the deep, multi-frequency data available in the Lockman Hole, we identify a sample of 100 Ultra-steep spectrum (USS) sources and 13 peaked spectrum sources. We estimate that up to 21 percent of these could have $z>4$ and are candidate high-$z$ radio galaxies, but further follow-up observations are required to confirm the physical nature of these objects.

We present LOFAR High-Band Array (HBA) observations of the Herschel-ATLAS North Galactic Pole survey area. The survey we have carried out, consisting of four pointings covering around 142 square degrees of sky in the frequency range 126--173 MHz, does not provide uniform noise coverage but otherwise is representative of the quality of data to be expected in the planned LOFAR wide-area surveys, and has been reduced using recently developed `facet calibration' methods at a resolution approaching the full resolution of the datasets ($\sim 10 \times 6$ arcsec) and an rms off-source noise that ranges from 100 $\mu$Jy beam$^{-1}$ in the centre of the best fields to around 2 mJy beam$^{-1}$ at the furthest extent of our imaging. We describe the imaging, cataloguing and source identification processes, and present some initial science results based on a 5-$\sigma$ source catalogue. These include (i) an initial look at the radio/far-infrared correlation at 150 MHz, showing that many Herschel sources are not yet detected by LOFAR; (ii) number counts at 150 MHz, including, for the first time, observational constraints on the numbers of star-forming galaxies; (iii) the 150-MHz luminosity functions for active and star-forming galaxies, which agree well with determinations at higher frequencies at low redshift, and show strong redshift evolution of the star-forming population; and (iv) some discussion of the implications of our observations for studies of radio galaxy life cycles.

We report on high-resolution JVLA and Chandra observations of the HST Frontier Cluster MACS J0717.5+3745. MACS J0717.5+3745 offers the largest contiguous magnified area of any known cluster, making it a promising target to search for lensed radio and X-ray sources. With the high-resolution 1.0-6.5 GHz JVLA imaging in A and B configuration, we detect a total of 51 compact radio sources within the area covered by the HST imaging. Within this sample we find 7 lensed sources with amplification factors larger than $2$. None of these sources are identified as multiply-lensed. Based on the radio luminosities, the majority of these sources are likely star forming galaxies with star formation rates of 10-50 M$_\odot$ yr$^{-1}$ located at $1 \lesssim z \lesssim 2$. Two of the lensed radio sources are also detected in the Chandra image of the cluster. These two sources are likely AGN, given their $2-10$ keV X-ray luminosities of $\sim 10^{43-44}$ erg s$^{-1}$. From the derived radio luminosity function, we find evidence for an increase in the number density of radio sources at $0.6<z<2.0$, compared to a $z < 0.3$ sample. Our observations indicate that deep radio imaging of lensing clusters can be used to study star forming galaxies, with star formation rates as low as $\sim10$ M$_{\odot}$ yr$^{-1}$, at the peak of cosmic star formation history.

We present the first wide area (19 deg$^2$), deep ($\approx120-150$ {\mu}Jy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Bo\"otes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which includes full direction-dependent calibration, are described here. We present a radio source catalogue containing 6276 sources detected over an area of $19$\,deg$^2$, with a peak flux density threshold of $5\sigma$. As the first thorough test of the facet calibration strategy, introduced by van Weeren et al., we investigate the flux and positional accuracy of the catalogue. We present differential source counts that reach an order of magnitude deeper in flux density than previously achieved at these low frequencies, and show flattening at 150 MHz flux densities below 10 mJy associated with the rise of the low flux density star-forming galaxies and radio-quiet AGN.

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

With Low-Frequency Array (LOFAR) observations, we have discovered a diverse assembly of steep spectrum emission that is apparently associated with the intra cluster medium (ICM) of the merging galaxy cluster Abell 2034. Such a rich variety of complex emission associated with the ICM has been observed in few other clusters. This not only indicates that Abell 2034 is a more interesting and complex system than previously thought but it also demonstrates the importance of sensitive and high-resolution, low-frequency observations. These observations can reveal emission from relativistic particles which have been accelerated to sufficient energy to produce observable emission or have had their high energy maintained by mechanisms in the ICM. The most prominent feature in our maps is a bright bulb of emission connected to two steep spectrum filamentary structures, the longest of which extends perpendicular to the merger axis for 0.5Mpc across the south of the cluster. The origin of these objects is unclear, with no shock detected in the X-ray images and no obvious connection with cluster galaxies or AGNs. We also find that the X-ray bright region of the cluster coincides with a giant radio halo with an irregular morphology and a very steep spectrum. In addition, the cluster hosts up to three possible radio relics, which are misaligned with the cluster X-ray emission. Finally, we have identified multiple regions of emission with a very steep spectral index that seem to be associated with either tailed radio galaxies or a shock.

We explore the impact of electron thermal conduction on the evolution of radiatively-cooled cold clouds embedded in flows of hot and fast material, as occur in outflowing galaxies. Performing a parameter study of three-dimensional adaptive mesh refinement hydrodynamical simulations, we show that electron thermal conduction causes cold clouds to evaporate, but it can also extend their lifetimes by compressing them into dense filaments. We distinguish between low column-density clouds, which are disrupted on very short times, and high-column density clouds with much-longer disruption times that are set by a balance between impinging thermal energy and evaporation. We provide fits to the cloud lifetimes and velocities that can be used in galaxy-scale simulations of outflows, in which the evolution of individual clouds cannot be modeled with the required resolution. Moreover, we show that the clouds are only accelerated to a small fraction of the ambient velocity because compression by evaporation causes the clouds to present a small cross-section to the ambient flow. This means that either magnetic fields must suppress thermal conduction, or that the cold clouds observed in galaxy outflows are not formed of cold material carried out from the galaxy.

We present deep LOFAR observations between 120-181 MHz of the "Toothbrush" (RX J0603.3+4214), a cluster that contains one of the brightest radio relic sources known. Our LOFAR observations exploit a new and novel calibration scheme to probe 10 times deeper than any previous study in this relatively unexplored part of the spectrum. The LOFAR observations, when combined with VLA, GMRT, and Chandra X-ray data, provide new information about the nature of cluster merger shocks and their role in re-accelerating relativistic particles. We derive a spectral index of $\alpha = -0.8 \pm 0.1$ at the northern edge of the main radio relic, steepening towards the south to $\alpha \approx - 2$. The spectral index of the radio halo is remarkably uniform ($\alpha = -1.16$, with an intrinsic scatter of $\leq 0.04$). The observed radio relic spectral index gives a Mach number of $\mathcal{M} = 2.8^{+0.5}_{-0.3}$, assuming diffusive shock acceleration (DSA). However, the gas density jump at the northern edge of the large radio relic implies a much weaker shock ($\mathcal{M} \approx 1.2$, with an upper limit of $\mathcal{M} \approx 1.5$). The discrepancy between the Mach numbers calculated from the radio and X-rays can be explained if either (i) the relic traces a complex shock surface along the line of sight, or (ii) if the radio relic emission is produced by a re-accelerated population of fossil particles from a radio galaxy. Our results highlight the need for additional theoretical work and numerical simulations of particle acceleration and re-acceleration at cluster merger shocks.

LOFAR, the Low-Frequency Array, is a powerful new radio telescope operating between 10 and 240 MHz. LOFAR allows detailed sensitive high-resolution studies of the low-frequency radio sky. At the same time LOFAR also provides excellent short baseline coverage to map diffuse extended emission. However, producing high-quality deep images is challenging due to the presence of direction dependent calibration errors, caused by imperfect knowledge of the station beam shapes and the ionosphere. Furthermore, the large data volume and presence of station clock errors present additional difficulties. In this paper we present a new calibration scheme, which we name facet calibration, to obtain deep high-resolution LOFAR High Band Antenna images using the Dutch part of the array. This scheme solves and corrects the direction dependent errors in a number of facets that cover the observed field of view. Facet calibration provides close to thermal noise limited images for a typical 8 hr observing run at $\sim$ 5arcsec resolution, meeting the specifications of the LOFAR Tier-1 northern survey.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

Remnant radio galaxies represent the final dying phase of radio galaxy evolution, in which the jets are no longer active. Due to their rarity in flux limited samples and the difficulty of identification, this dying phase remains poorly understood and the luminosity evolution largely unconstrained. Here we present the discovery and detailed analysis of a large (700 kpc) remnant radio galaxy with a low surface brightness that has been identified in LOFAR images at 150 MHz. By combining LOFAR data with new follow-up Westerbork observations and archival data at higher frequencies, we investigated the source morphology and spectral properties from 116 to 4850 MHz. By modelling the radio spectrum we probed characteristic timescales of the radio activity. The source has a relatively smooth, diffuse, amorphous appearance together with a very weak central compact core which is associated with the host galaxy located at z=0.051. From our ageing and morphological analysis it is clear that the nuclear engine is currently switched off or, at most, active at a very low power state. The host galaxy is currently interacting with another galaxy located at a projected separation of 15 kpc and a radial velocity offset of 300 km/s. This interaction may have played a role in the triggering and/or shut down of the radio jets. The spectral shape of this remnant radio galaxy differs from the majority of the previously identified remnant sources, which show steep or curved spectra at low to intermediate frequencies. In light of this finding and in preparation for new-generation deep low-frequency surveys, we discuss the selection criteria to be used to select representative samples of these sources.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

Steep spectrum radio sources associated with active galactic nuclei (AGN) may contain remnants of past AGN activity episodes. Novel instruments like the LOw Frequency ARray (LOFAR) are enabling studies of these fascinating structures to be made at tens to hundreds of MHz with sufficient resolution to analyse their complex morphology. Our goal is to characterize the integrated and resolved spectral properties of VLSS J1431+1331 and estimate source ages based on synchrotron radio emission models, thus putting constraints on the AGN duty cycle. Using a broad spectral coverage, we have derived spectral and curvature maps, and used synchrotron ageing models to determine the time elapsed from the last time the source plasma was energized. We used LOFAR, Giant Metrewave Radio Telescope (GMRT) and Jansky Very Large Array (VLA) data. Based on our ageing analysis, we infer that the AGN that created this source currently has very low levels of activity or that it is switched off. The derived ages for the larger source component range from around 60 to 130 Myr, hinting that the AGN activity decreased or stopped around 60 Myr ago. Our analysis suggests that VLSS J1431.8+1331 is an intriguing, two-component source. The larger component seems to host a faint radio core, suggesting that the source may be an AGN radio relic. The spectral index we observe from the smaller component is distinctly flatter at lower frequencies than the spectral index of the larger component, suggesting the possibility that the smaller component may be a shocked plasma bubble. From the integrated source spectrum, we deduce that its shape and slope can be used as tracers of the activity history of this type of steep spectrum radio source.

The merger of two white dwarfs (WDs) has for many years not been considered as the favoured model for the progenitor system of type Ia supernovae (SNe Ia). But recent years have seen a change of opinion as a number of studies, both observational and theoretical, have concluded that they should contribute significantly to the observed type Ia supernova rate. In this paper, we study the ignition and propagation of detonation through post-merger remnants and we follow the resulting nucleosynthesis up to the point where a homologous expansion is reached. In our study we cover the entire range of WD masses and compositions. For the emergence of a detonation we study three different setups. The first two are guided by the merger remnants from our earlier simulations (Dan et al. 2014}, while for the third one the ignitions were set by placing hotspots with properties determined by spatially resolved calculations taken from the literature. There are some caveats to our approach which we investigate. We carefully compare the nucleosynthetic yields of successful explosions with SN Ia observations. Only three of our models are consistent with all the imposed constraints and potentially lead to a standard type Ia event. The first one, a $0.45\,M_\odot\ {\rm helium\ (He)} + 0.9 \,M_\odot$ carbon-oxygen (CO) WD system produces a sub-luminous, SN 1991bg-like event while the other two, a $0.45\,M_\odot\ {\rm He} + 1.1\,M_\odot$ oxygen-neon (ONe) WD system and a $1.05+1.05\,M_\odot$ system with two CO WDs, are good candidates for common SNe Ia.

Radio halos are synchrotron radio sources detected in some massive galaxy clusters. Their Mpc-size indicates that (re)acceleration processes are taking place in the host cluster. X-ray catalogues of galaxy clusters have been used in the past to search for radio halos and to understand their connection with cluster-cluster mergers and with the thermal component of the intra-cluster medium. More recently, the Sunyaev-Zel'dovich effect has been proven to be a better route to search for massive clusters in a wider redshift range. With the aim of discovering new radio halos and understanding their connection with cluster-cluster mergers, we have selected from the Planck Early source catalog the most massive clusters, and we have observed with the Giant Metrewave Radio Telescope at 323 MHz those objects for which deep observations were not available. We have discovered new peculiar radio emission in three of the observed clusters finding: (i) a radio halo in the cluster RXCJ0949.8+1708; (ii) extended emission in Abell 1443 that we classify as a radio halo plus a radio relic, with a bright filament embedded in the radio halo; (iii) low-power radio emission is found in CIZA J1938.3+5409 which is ten times below the radio - X-ray correlation, and represents the first direct detection of the radio emission in the "upper-limit" region of the radio - X-ray diagram. We discuss the properties of these new radio halos in the framework of theoretical models for the radio emission.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

Using a new parallel algorithm implemented within the VisIt framework, we analysed large cosmological grid simulations to study the properties of baryons in filaments. The procedure allows us to build large catalogues with up to $\sim 3 \cdot 10^4$ filaments per simulated volume and to investigate the properties of cosmic filaments for very large volumes at high resolution (up to $300^3 ~\rm Mpc^3$ simulated with $2048^3$ cells). We determined scaling relations for the mass, volume, length and temperature of filaments and compared them to those of galaxy clusters. The longest filaments have a total length of about $200 ~\rm Mpc$ with a mass of several $10^{15} M_{\odot}$. We also investigated the effects of different gas physics. Radiative cooling significantly modifies the thermal properties of the warm-hot-intergalactic medium of filaments, mainly by lowering their mean temperature via line cooling. On the other hand, powerful feedback from active galactic nuclei in surrounding halos can heat up the gas in filaments. The impact of shock-accelerated cosmic rays from diffusive shock acceleration on filaments is small and the ratio of between cosmic ray and gas pressure within filaments is of the order of $\sim 10-20$ percent.

Content: We present the results from $Suzaku$ observations of the merging cluster of galaxies CIZA J2242.8+5301 at $z$=0.192. Aims. To study the physics of gas heating and particle acceleration in cluster mergers, we investigated the X-ray emission from CIZA J2242.8+5301, which hosts two giant radio relics in the northern/southern part of the cluster. Methods. We analyzed data from three-pointed Suzaku observations of CIZA J2242.8+5301 to derive the temperature distribution in four different directions. Results: The Intra-Cluster Medium (ICM) temperature shows a remarkable drop from 8.5$_{-0.6}^{+0.8}$ keV to 2.7$_{-0.4}^{+0.7}$ keV across the northern radio relic. The temperature drop is consistent with a Mach number ${\cal M}_n=2.7^{+0.7}_{-0.4}$ and a shock velocity $v_{shock:n}=2300_{-400}^{+700}\rm\,km\,s^{-1}$. We also confirm the temperature drop across the southern radio relic. However, the ICM temperature beyond this relic is much higher than beyond the northern one, which gives a Mach number ${\cal M}_s=1.7^{+0.4}_{-0.3}$ and shock velocity $v_{shock:s}=2040_{-410}^{+550}\rm \,km\,s^{-1}$. These results agree with other systems showing a relationship between the radio relics and shock fronts which are induced by merging activity. We compare the X-ray derived Mach numbers with the radio derived Mach numbers from the radio spectral index under the assumption of diffusive shock acceleration in the linear test particle regime. For the northern radio relic, the Mach numbers derived from X-ray and radio observations agree with each other. Based on the shock velocities, we estimate that CIZA J2242.8+5301 is observed approximately 0.6 Gyr after core passage. The magnetic field pressure at the northern relic is estimated to be 9% of the thermal pressure.

The cosmic web contains a large fraction of the total gas mass in the universe but is difficult to detect at most wavelengths. Synchrotron emission from shock-accelerated electrons may offer the chance of imaging the cosmic web at radio wavelengths. In this work we use 3D cosmological ENZO-MHD simulations (combined with a post-processing renormalisation of the magnetic field to bracket for missing physical ingredients and resolution effects) to produce models of the radio emission from the cosmic web. In post-processing we study the capabilities of 13 large radio surveys to detect this emission. We find that surveys by LOFAR, SKA1-LOW and MWA have a chance of detecting the cosmic web, provided that the magnetisation level of the tenuous medium in filaments is of the order of 1% of the thermal gas energy.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

To better understand the nature of the multiphase material found in outflowing galaxies, we study the evolution of cold clouds embedded in flows of hot and fast material. Using a suite of adaptive-mesh refinement simulations that include radiative cooling, we investigate both cloud mass loss and cloud acceleration under the full range of conditions observed in galaxy outflows. The simulations are designed to track the cloud center of mass, enabling us to study the cloud evolution at long disruption times. For supersonic flows, a Mach cone forms around the cloud, which damps the Kelvin-Helmholtz instability but also establishes a streamwise pressure gradient that stretches the cloud apart. If time is expressed in units of the cloud crushing time, both the cloud lifetime and the cloud acceleration rate are independent of cloud radius, and we find simple scalings for these quantities as a function of the Mach number of the external medium. A resolution study suggests that our simulations have sufficient resolution to accurately describe the evolution of cold clouds in the absence of thermal conduction and magnetic fields, physical processes whose roles will be studied in forthcoming papers.

Context. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SKA

Cluster mergers may play a fundamental role in the formation and evolution of cluster galaxies. Stroe et al. (2014) revealed unexpected over-densities of candidate H$\alpha$ emitters near the ~1 Mpc-wide shock fronts of the massive (~2x10$^{15}$M$_{\odot}$) "Sausage" merging cluster, CIZA J2242.8+5301. We used Keck/DEIMOS and WHT/AF2 to confirm 83 H$\alpha$ emitters in and around the merging cluster. We find that cluster star-forming galaxies in the hottest X-ray gas and/or in the cluster sub-cores (away from the shock fronts) show high [SII]6716/[SII]6761 and high [SII]6716/H$\alpha$, implying very low electron densities (<30x lower than all other star-forming galaxies outside the cluster) and significant contribution from supernovae, respectively. All cluster star-forming galaxies near the cluster centre show evidence of significant outflows (blueshifted Na D~200-300km/s), likely driven by supernovae. Strong outflows are also found for the cluster H$\alpha$ AGN. H$\alpha$ star-forming galaxies in the merging cluster follow the z~0 mass-metallicity relation, showing systematically higher metallicity (~0.15-0.2 dex) than H$\alpha$ emitters outside the cluster (projected R>2.5 Mpc). This suggests that the shock front may have triggered remaining metal-rich gas which galaxies were able to retain into forming stars. Our observations show that the merger of impressively massive (~10$^{15}$M$_\odot$) clusters can provide the conditions for significant star-formation and AGN activity, but, as we witness strong feedback by star-forming galaxies and AGN (and given how massive the merging cluster is), such sources will likely quench in a few 100 Myrs.

We investigate the possibility for the SKA to detect and study the magnetic fields in galaxy clusters and in the less dense environments surrounding them using Faraday Rotation Measures. To this end, we produce 3-dimensional magnetic field models for galaxy clusters of different masses and in different stages of their evolution, and derive mock rotation measure observations of background radiogalaxies. According to our results, already in phase I, we will be able to infer the magnetic field properties in galaxy clusters as a function of the cluster mass, down to $10^{13}$ solar-masses. Moreover, using cosmological simulations to model the gas density, we have computed the expected rotation measure through shock-fronts that occur in the intra-cluster medium during cluster mergers. The enhancement in the rotation measure due to the density jump will permit to constraint the magnetic field strength and structure after the shock passage. SKA observations of polarised sources located behind galaxy clusters will answer several questions about the magnetic field strength and structure in galaxy clusters, and its evolution with cosmic time.

The detection of the diffuse gas component of the cosmic web remains a formidable challenge. In this work we study synchrotron emission from the cosmic web with simulated SKA1 observations, which can represent an fundamental probe of the warm-hot intergalactic medium. We investigate radio emission originated by relativistic electrons accelerated by shocks surrounding cosmic filaments, assuming diffusive shock acceleration and as a function of the (unknown) large-scale magnetic fields. The detection of the brightest parts of large ($>10 \rm Mpc$) filaments of the cosmic web should be within reach of the SKA1-LOW, if the magnetic field is at the level of a $\sim 10$ percent equipartition with the thermal gas, corresponding to $\sim 0.1 \mu G$ for the most massive filaments in simulations. In the course of a 2-years survey with SKA1-LOW, this will enable a first detection of the "tip of the iceberg" of the radio cosmic web, and allow for the use of the SKA as a powerful tool to study the origin of cosmic magnetism in large-scale structures. On the other hand, the SKA1-MID and SKA1-SUR seem less suited for this science case at low redshift ($z \leq 0.4$), owing to the missing short baselines and the consequent lack of signal from the large-scale brightness fluctuations associated with the filaments. In this case only very long exposures ($\sim 1000$ hr) may enable the detection of $\sim 1-2$ filament for field of view in the SKA1-SUR PAF Band1.

Extended steep-spectrum radio emission in a galaxy cluster is usually associated with a recent merger. However, given the complex scenario of galaxy cluster mergers, many of the discovered sources hardly fit into the strict boundaries of a precise taxonomy. This is especially true for radio phoenixes that do not have very well defined observational criteria. Radio phoenixes are aged radio galaxy lobes whose emission is reactivated by compression or other mechanisms. Here, we present the detection of a radio phoenix close to the moment of its formation. The source is located in Abell 1033, a peculiar galaxy cluster which underwent a recent merger. To support our claim, we present unpublished Westerbork Synthesis Radio Telescope and Chandra observations together with archival data from the Very Large Array and the Sloan Digital Sky Survey. We discover the presence of two sub-clusters displaced along the N-S direction. The two sub-clusters probably underwent a recent merger which is the cause of a moderately perturbed X-ray brightness distribution. A steep-spectrum extended radio source very close to an AGN is proposed to be a newly born radio phoenix: the AGN lobes have been displaced/compressed by shocks formed during the merger event. This scenario explains the source location, morphology, spectral index, and brightness. Finally, we show evidence of a density discontinuity close to the radio phoenix and discuss the consequences of its presence.

Giant Radio Halos (RH) are diffuse, Mpc-sized, synchrotron radio sources observed in a fraction of merging galaxy clusters. The current scenario for the origin of RHs assumes that turbulence generated during cluster mergers re-accelerates pre-existing fossil and/or secondary electrons in the intra-cluster-medium (ICM) to the energies necessary to produce the observed radio emission. Moreover, more relaxed clusters could host diffuse "off state" halos produced by secondary electrons. In this Chapter we use Monte Carlo simulations, that combine turbulent-acceleration physics and the generation of secondaries in the ICM, to calculate the occurrence of RHs in the Universe, their spectral properties and connection with properties of the hosting clusters. Predictions for SKA1 surveys are presented at low (100-300 MHz) and mid (1-2 GHz) frequencies assuming the expected sensitivities and spatial resolutions of SKA1. SKA1 will step into an unexplored territory allowing us to study the formation and evolution of RHs in a totally new range of cluster masses and redshift, allowing firm tests of the current theoretical hypothesis. In particular, the combination of SKA1-LOW and SUR will allow the discovery of ~1000 ultrasteep- spectrum halos and to detect for the very first time "off state" RHs. We expect that at least ~2500 giant RHs will be discovered by SKA1-LOW surveys up to z~0.6. Remarkably these surveys will be sensitive to RHs in a cluster mass range (down to ~10^14 solar masses) and redshifts (up to ~1) that are unexplored by current observations. SKA1 surveys will be highly competitive with present and future SZ-surveys in the detection of high-redshift massive objects.

CIZA J2242.8+5301 (`Sausage') and 1RXS J0603.3+4213 (`Toothbrush') are two low-redshift ($z\sim0.2$), massive ($\sim2\times10^{15}M_\odot$), post-core passage merging clusters, which host shock waves traced by diffuse radio emission. To study their star-formation properties, we uniformly survey the `Sausage' and `Toothbrush' clusters in broad and narrow band filters and select a sample of $201$ and $463$ line emitters, down to a rest-frame equivalent width ($13${\AA}). We robustly separate between H$\alpha$ and higher redshift emitters using a combination of optical multi-band (B, g, V, r, i, z) and spectroscopic data. We build H$\alpha$ luminosity functions for the entire cluster region, near the shock fronts, and away from the shock fronts and find striking differences between the two clusters. In the dynamically younger, $1$ Gyr old `Sausage' cluster we find numerous ($59$) H$\alpha$ emitters above a star-formation rate (SFR) of $0.17$ M$_{\sun}$ yr$^{-1}$ surprisingly located in close proximity to the shock fronts, embedded in very hot intra-cluster medium plasma. The SFR density for the cluster population is at least at the level of typical galaxies at $z\sim2$. Down to the same star-formation rate, the possibly dynamically more evolved `Toothbrush' cluster has only $9$ H$\alpha$ galaxies. The cluster H$\alpha$ galaxies fall on the SFR-stellar mass relation $z\sim0.2$ for the field. However, the `Sausage' cluster has an H$\alpha$ emitter density $>20$ times that of blank fields. If the shock passes through gas-rich cluster galaxies, the compressed gas could collapse into dense clouds and excite star-formation for a few $100$ Myr. This process ultimately leads to a rapid consumption of the molecular gas, accelerating the transformation of gas-rich field spirals into cluster S0s or ellipticals.

The nuclear starburst in the nearby galaxy M82 provides an excellent laboratory for understanding the physics of star formation. This galaxy has been extensively observed in the past, revealing tens of radio-bright compact objects embedded in a diffuse free-free absorbing medium. Our understanding of the structure and physics of this medium in M82 can be greatly improved by high-resolution images at low frequencies where the effects of free-free absorption are most prominent. The aims of this study are, firstly, to demonstrate imaging using international baselines of the Low Frequency Array (LOFAR), and secondly, to constrain low-frequency spectra of compact and diffuse emission in the central starburst region of M82 via high-resolution radio imaging at low frequencies. The international LOFAR telescope was used to observe M82 at 110-126MHz and 146-162MHz. Images were obtained using standard techniques from very long baseline interferometry. images were obtained at each frequency range: one only using international baselines, and one only using the longest Dutch (remote) baselines. The 154MHz image obtained using international baselines is a new imaging record in terms of combined image resolution (0.3$"$) and sensitivity ($\sigma$=0.15mJy/beam) at low frequencies ($<327$MHz). We detected 16 objects at 154MHz, six of these also at 118MHz. Four weaker but resolved features are also found: a linear (50pc) filament and three other resolved objects, of which two show a clear shell structure. We do not detect any emission from either supernova 2008iz or from the radio transient source 43.78+59.3. The images obtained using remote baselines show diffuse emission, associated with the outflow in M82, with reduced brightness in the region of the edge-on star-forming disk.

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

The galaxy cluster Abell 2256 hosts one of the most intriguing examples in the class of radio relics. It has been found that this radio relic has a rather flat integrated spectrum at low frequencies that would imply an injection spectral index for the electrons that is inconsistent with the flattest allowed by the test particle diffusive shock acceleration (DSA). We performed new high-frequency observations at 2273, 2640, and 4850 MHz. Combining these new observations with images available in the literature, we constrain the radio integrated spectrum of the radio relic in Abell 2256 over the widest sampled frequency range collected so far for this class of objects (63 -10450 MHz). Moreover, we used X-ray observations of the cluster to check the temperature structure in the regions around the radio relic. We find that the relic keeps an unusually flat behavior up to high frequencies. Although the relic integrated spectrum between 63 and 10450 MHz is not inconsistent with a single power law with $\alpha_{63}^{10450}= 0.92\pm 0.02$, we find hints of a steepening at frequencies > 1400 MHz. The two frequency ranges 63-1369 MHz and 1369-10450 MHz are, indeed, best represented by two different power laws, with $\alpha_{63}^{1369}= 0.85\pm 0.01$ and $\alpha_{1369}^{10450}= 1.00\pm 0.02$. This broken power law would require special conditions to be explained in terms of test-particle DSA, e.g., non-stationarity of the spectrum and/or non-stationarity of the shock. On the other hand, the single power law would make of this relic the one with the flattest integrated spectrum known so far, even flatter than what allowed in the test-particle approach to DSA. We find a rather low temperature ratio of $T_2/T_1 \sim 1.7$ across the G region of the radio relic and no temperature jump across the H region.

The amplification of primordial magnetic fields via a small-scale turbulent dynamo during structure formation might be able to explain the observed magnetic fields in galaxy clusters. The magnetisation of more tenuous large-scale structures such as cosmic filaments is more uncertain, as it is challenging for numerical simulations to achieve the required dynamical range. In this work, we present magneto-hydrodynamical cosmological simulations on large uniform grids to study the amplification of primordial seed fields in the intracluster medium (ICM) and in the warm-hot-intergalactic medium (WHIM). In the ICM, we confirm that turbulence caused by structure formation can produce a significant dynamo amplification, even if the amplification is smaller than what is reported in other papers. In the WHIM inside filaments, we do not observe significant dynamo amplification, even though we achieve Reynolds numbers of $R_{\rm e} \sim 200-300$. The maximal amplification for large filaments is of the order of $\sim 100$ for the magnetic energy, corresponding to a typical field of a few $\sim \rm nG$ starting from a primordial weak field of $10^{-10}$ G (comoving). In order to start a small-scale dynamo, we found that a minimum of $\sim 10^2$ resolution elements across the virial radius of galaxy clusters was necessary. In filaments we could not find a minimum resolution to set off a dynamo. This stems from the inefficiency of supersonic motions in the WHIM in triggering solenoidal modes and small-scale twisting of magnetic field structures. Magnetic fields this small will make it hard to detect filaments in radio observations.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

The grand-design spiral galaxy M51 was observed with the LOFAR High Frequency Antennas (HBA) and imaged in total intensity and polarisation. This observation covered the frequencies between 115 MHz and 175 MHz. We produced an image of total emission of M51 at the mean frequency of 151 MHz with 20 arcsec resolution and 0.3 mJy rms noise, which is the most sensitive image of a galaxy at frequencies below 300 MHz so far. The integrated spectrum of total radio emission is described well by a power law, while flat spectral indices in the central region indicate thermal absorption. We observe that the disk extends out to 16 kpc and see a break in the radial profile near the optical radius of the disk. Our main results, the scale lengths of the inner and outer disks at 151 MHz and 1.4 GHz, arm--interarm contrast, and the break scales of the radio--far-infrared correlations, can be explained consistently by CRE diffusion, leading to a longer propagation length of CRE of lower energy. The distribution of CRE sources drops sharply at about 10 kpc radius, where the star formation rate also decreases sharply. We find evidence that thermal absorption is primarily caused by HII regions. The non-detection of polarisation from M51 at 151 MHz is consistent with the estimates of Faraday depolarisation. Future searches for polarised emission in this frequency range should concentrate on regions with low star formation rates.

Giant radio halos are Mpc-size sources found in some merging galaxy clusters. The synchrotron emitting electrons are thought to be (re)accelerated by plasma turbulence induced by the merging of two massive clusters. Cool core galaxy clusters have a low temperature core, likely an indication that a major merger has not recently occurred. CL1821+643 is one of the strongest cool core clusters known so far. Surprisingly, we detect a giant radio halo with a largest linear size of $\sim$ 1.1 Mpc. We discuss the radio and X-ray properties of the cluster in the framework of the proposed models for giant radio halos. If a merger is causing the radio emission, despite the presence of a cool-core, we suggest that it should be off-axis, or in an early phase, or a minor one.

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<\nu<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.

Recent observations of ram pressure stripped spiral galaxies in clusters revealed details of the stripping process, i.e., the truncation of all interstellar medium (ISM) phases and of star formation (SF) in the disk, and multiphase star-forming tails. Some stripped galaxies, in particular in merging clusters, develop spectacular star-forming tails, giving them a jellyfish-like appearance. In merging clusters, merger shocks in the intra-cluster medium (ICM) are thought to have overrun these galaxies, enhancing the ambient ICM pressure and thus triggering SF, gas stripping and tail formation. We present idealised hydrodynamical simulations of this scenario, including standard descriptions for SF and stellar feedback. To aid the interpretation of recent and upcoming observations, we focus on particular structures and dynamics in SF patterns in the remaining gas disk and in the near tails, which are easiest to observe. The observed jellyfish morphology is qualitatively reproduced for, both, face-on and edge-on stripping. In edge-on stripping, the interplay between the ICM wind and the disk rotation leads to asymmetries along the ICM wind direction and perpendicular to it. The apparent tail is still part of a highly deformed gaseous and young stellar disk. In both geometries, SF takes place in knots throughout the tail, such that the stars in the tails show no ordered age gradients. Significant SF enhancement in the disk occurs only at radii where the gas will be stripped in due course.

The CIZA J2242.8+5301 cluster hosts an extraordinary, narrow, Mpc-wide, diffuse patch of radio emission tracing travelling shock waves, called a relic. We perform a multi-wavelength of observations of the cluster. We discover radio spectral steepening and increasing spectral curvature in the shock downstream area, indicative of emission produced by spectrally-aged electrons, following a JP ageing model. The curved integrated spectrum towards $16$ GHz challenges the favoured relic formation model. Further, we find a boosting the number density of H$\alpha$ emitters by a factor of $10$, possibly caused by interactions of the galaxies with the shock wave.

We report the discovery of extended radio emission in the Phoenix cluster (SPT-CL J2344-4243, z=0.596) with the GMRT at 610 MHz. The diffuse emission extends over a region of at least 400-500 kpc and surrounds the central radio source of the Brightest Cluster Galaxy, but does not appear to be directly associated with it. We classify the diffuse emission as a radio mini-halo, making it the currently most distant mini-halo known. Radio mini-halos have been explained by synchrotron emitting particles re-accelerated via turbulence, possibly induced by gas sloshing generated from a minor merger event. Chandra observations show a non-concentric X-ray surface brightness distribution, which is consistent with this sloshing interpretation. The mini-halo has a flux density of $17\pm5$ mJy, resulting in a 1.4 GHz radio power of ($10.4\pm3.5) \times 10^{24}$ W Hz$^{-1}$. The combined cluster emission, which includes the central compact radio source, is also detected in a shallow GMRT 156 MHz observation and together with the 610 MHz data we compute a spectral index of $-0.84\pm0.12$ for the overall cluster radio emission. Given that mini-halos typically have steeper radio spectra than cluster radio galaxies, this spectral index should be taken as an upper limit for the mini-halo.

CIZA J2242.8+5301, a merging galaxy cluster at z=0.19, hosts a double-relic system and a faint radio halo. Radio observations at frequencies ranging from a few MHz to several GHz have shown that the radio spectral index at the outer edge of the N relic corresponds to a shock of Mach number 4.6+/-1.1, under the assumptions of diffusive shock acceleration of thermal particles in the test particle regime. Here, we present results from new Chandra observations of the cluster. The Chandra surface brightness profile across the N relic only hints to a surface brightness discontinuity (<2-sigma detection). Nevertheless, our reanalysis of archival Suzaku data indicates a temperature discontinuity across the relic that is consistent with a Mach number of 2.5+/-0.5, in agreement with previously published results. This confirms that the Mach number at the shock traced by the N relic is much weaker than predicted from the radio. Puzzlingly, in the Chandra data we also identify additional inner small density discontinuities both on and off the merger axis. Temperature measurements on both sides of the discontinuities do not allow us to undoubtedly determine their nature, although a shock front interpretation seems more likely. We speculate that if the inner density discontinuities are indeed shock fronts, then they are the consequence of violent relaxation of the dark matter cores of the clusters involved in the merger.

Observations have revealed cold gas with large velocity dispersions (~300 km/s) within the hot outflows of ultra-luminous infrared galaxies (ULIRGs). This gas may trace its origin to the Rayleigh-Taylor (RT) fragmentation of a super-bubble or may arise on smaller scales. We model a ULIRG outflow at two scales to recreate this gas in three-dimensional hydrodynamic simulations using FLASH. Although resolution is limited, these models successfully produce cold gas in outflows with large velocity dispersions. Our small-scale models produce this cold gas through RT fragmentation of the super-bubble wall, but the large-scale models produce the cold gas after hot bubbles fragment the disc's gas into cold clouds which are then accelerated by thermal pressure, and supplemented by cooling within the outflow. We produce simple mock spectra to compare these simulations to observed absorption spectra and find line-widths of ~250 km/s, agreeing with the lower end of observations.

The galaxy pair NGC 5580 and NGC 5588 are part of a loose group of galaxies. They are surrounded by steep-spectrum, extended radio emission which was previously suggested to be a down-scaled example of Mpc-size radio haloes present in galaxies clusters. We present a multi-frequency study of the radio-emission aimed to clarify its nature. The source has been observed with the Giant Meterwave Radio Telescope at 235, 325 and 610 MHz and the images obtained were combined with archival data to cover the frequency range 150-1400 MHz. The new observations revealed the presence of a second, fainter lobe on the South-East of NGC 5580. The spectral index study of the source shows a flattening of the spectrum (which implies a younger particle population) close to the two galaxies. We argue that the extended radio emission is the remnant of a past activity cycle of the active galactic nucleus present in NGC 5580 and therefore a notable example of a dying radio galaxy located outside a dense environments.

Radio relics are diffuse radio sources observed in galaxy clusters, probably produced by shock acceleration during cluster-cluster mergers. Their large size, of the order of 1 Mpc, indicates that the emitting electrons need to be (re)accelerated locally. The usually invoked Diffusive Shock Acceleration models have been challenged by recent observations and theory. We report the discovery of complex radio emission in the Galaxy cluster PLCKG287.0 +32.9, which hosts two relics, a radio halo, and several radio filamentary emission. Optical observations suggest that the cluster is elongated, likely along an intergalactic filament, and displays a significant amount of substructure. The peculiar features of this radio relic are that (i) it appears to be connected to the lobes of a radio galaxy and (ii) the radio spectrum steepens on either side of the radio relic. We discuss the origins of these features in the context of particle re-acceleration.

Recent X-ray and Sunyaev-Zel'dovich (SZ) observations have dramatically increased the number of known distant galaxy clusters. In some merging, low-redshift (z<0.4) clusters, centrally-located, diffuse, extended radio emission (called radio halos) has been found. Using the Giant Metrewave Radio Telescope (GMRT), we report the detection of diffuse radio emission in the binary-merging cluster PLCK G147.3-16.6 located at z = 0.65. We classify the emission as a giant radio halo due to the large physical extent of about 0.9 Mpc and low-surface brightness. We measure an integrated flux density of $7.3\pm1.1$ mJy at 610 MHz for the radio halo, resulting in a 1.4 GHz radio power of $5.1 \times 10^{24}$ W/Hz. The radio halo power is consistent with that expected from the known correlation between X-ray luminosity or the cluster integrated SZ signal and radio power. Our observations also suggest that more of these distant radio halos could be discovered with the GMRT.

Non-thermal (relativistic) particles are injected into the cosmos by structure formation shock waves, active galactic nuclei and stellar explosions. We present a suite of unigrid cosmological simulations (up to $2048^3$) using a two-fluid model in the grid code ENZO. The simulations include the dynamical effects of cosmic-ray (CR) protons and cover a range of theoretically motivated acceleration efficiencies. For the bulk of the cosmic volume the modelling of CR processes is rather stable with respect to resolution, provided that a minimum (cell) resolution of $\approx 100 ~\rm kpc/h$ is employed. However, the results for the innermost cluster regions depend on the assumptions for the baryonic physics. Inside clusters, non-radiative runs at high resolution tend to produce an energy density of CRs that are below available upper limits from the FERMI satellite, while the radiative runs are found to produce a higher budget of CRs. We show that weak ($M \leq 3-5$) shocks and shock-reacceleration are crucial to set the level of CRs in the innermost region of clusters, while in the outer regions the level of CR energy is mainly set via direct injection by stronger shocks, and is less sensitive to cooling and feedback from active galactic nuclei and supernovae.

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

Mechanical feedback via Active Galactic Nuclei (AGN) jets in the centres of galaxy groups and clusters is a crucial ingredient in current models of galaxy formation and cluster evolution. Jet feedback is believed to regulate gas cooling and thus star formation in the most massive galaxies, but a robust physical understanding of this feedback mode is currently lacking. The large collecting area, excellent spectral resolution and high spatial resolution of Athena+ will provide the breakthrough diagnostic ability necessary to develop this understanding, via: (1) the first kinematic measurements on relevant spatial scales of the hot gas in galaxy, group and cluster haloes as it absorbs the impact of AGN jets, and (2) vastly improved ability to map thermodynamic conditions on scales well-matched to the jets, lobes and gas disturbances produced by them. Athena+ will therefore determine for the first time how jet energy is dissipated and distributed in group and cluster gas, and how a feedback loop operates in group/cluster cores to regulate gas cooling and AGN fuelling. Athena+ will also establish firmly the cumulative impact of powerful radio galaxies on the evolution of baryons from the epoch of group/cluster formation to the present day.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

The aim of this work is to probe the magnetic field properties in relics and infall regions of galaxy clusters using Faraday Rotation Measures. We present Very Large Array multi-frequency observations of seven sources in the region South-West of the Coma cluster, where the infalling group NGC4839 and the relic 1253+275 are located. The Faraday Rotation Measure maps for the observed sources are derived and analysed to study the magnetic field in the South-West region of Coma. We discuss how to interpret the data by comparing observed and mock rotation measures maps that are produced simulating different 3-dimensional magnetic field models. The magnetic field model that gives the best fit to the Coma central region underestimates the rotation measure in the South-West region by a factor ~6, and no significant jump in the rotation measure data is found at the position of the relic. We explore different possibilities to reconcile observed and mock rotation measure trends, and conclude that an amplification of the magnetic field along the South-West sector is the most plausible solution. Our data together with recent X-ray estimates of the gas density obtained with Suzaku suggest that a magnetic field amplification by a factor ~3 is required throughout the entire South-West region in order to reconcile real and mock rotation measures trends. The magnetic field in the relic region is inferred to be ~2 microG, consistent with Inverse Compton limits.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

Giant cluster radio relics are thought to form at shock fronts in the course of collisions between galaxy clusters. Via processes that are still poorly understood, these shocks accelerate or re-accelerate cosmic-ray electrons and might amplify magnetic fields. The best object to study this phenomenon is the galaxy cluster CIZA J2242.8+5301 as it shows the most undisturbed relic. By means of Giant Metrewave Radio Telescope (GMRT) and Westerbork Synthesis Radio Telescope (WSRT) data at seven frequencies spanning from 153 MHz to 2272 MHz, we study the synchrotron emission in this cluster. We aim at distinguishing between theoretical injection and acceleration models proposed for the formation of radio relics. We also study the head-tail radio sources to reveal the interplay between the merger and the cluster galaxies. We produced spectral index, curvature maps and radio colour-colour plots and compared our data with predictions from models. We present one of the deepest 153 MHz maps of a cluster ever produced, reaching a noise level of 1.5 mJy/beam. We derive integrated spectra for four relics in the cluster, discovering extremely steep spectrum diffuse emission concentrated in multiple patches. We find a possible radio phoenix embedded in the relic to the south of the cluster. The spectral index of the northern relic retains signs of steepening from the front towards the back of the shock also at the radio frequencies below 600 MHz. The spectral curvature in the same relic also increases in the downstream area. The data is consistent with the Komissarov-Gubanov injection models, meaning that the emission we observe is produced by a single burst of spectrally-aged accelerated radio electrons.

The cluster 1RXS J0603.3+4214 is a merging galaxy cluster that hosts three radio relics and a giant radio halo. The northern relic, the Toothbrush, is 1.9-Mpc long and has an unusual linear morphology. According to simple diffusive shock acceleration theory, its radio spectral index indicates a Mach number of 3.3-4.6. Here, we present results from a deep XMM-Newton observation of the cluster. We observe two distinct cluster cores that have survived the merger. The presence of three shocks at or near the locations of the radio relics is confirmed by density and temperature discontinuities. However, the observation poses several puzzles that challenge our understanding of radio relics: (i) at the Toothbrush, the shock Mach number is not larger than 2, in apparent conflict with the shock strength predicted from the radio spectrum; (ii) at the Toothbrush, the shock front is, in part, spatially offset from the radio emission; (iii) at the eastern relic, we detect a temperature jump corresponding to a Mach number of approximately 2.5, but there is no associated surface brightness discontinuity. We discuss possible explanations for these findings.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

The Coma cluster is one of the nearest galaxy clusters, and the first one in which a radio halo and a peripheral relic were discovered. While its halo and the central parts of the intracluster medium have been studied extensively, X-ray observations of the plasma near its relic have been scarce. Here, we present results from a re-analysis of a 22-ks archival XMM-Newton observation. Across the relic, we detect a shock of Mach number about 2. This excludes the previously suggested hypothesis that the relic was formed by turbulence. Furthermore, multiwavelenth observations and numerical models do not support the scenario in which the shock at the Coma relic is an outgoing cluster-merger shock. Instead, our results lend support to the idea that the relic coincides with an infall shock front formed just as the NGC 4839 group falls onto the cluster along a cosmic filament.

We studied the intracluster medium of the galaxy cluster CIZA J2242.8+5301 using deep XMM-Newton observations. The cluster hosts a remarkable 2-Mpc long, ~50-kpc wide radio relic that has been nicknamed the "Sausage". A smaller, more irregular counter-relic is also present, along with a faint giant radio halo. We analysed the distribution of the ICM physical properties, and searched for shocks by trying to identify density and temperature discontinuities. East of the southern relic, we find evidence of shock compression corresponding to a Mach number of 1.3, and speculate that the shock extends beyond the length of the radio structure. The ICM temperature increases at the northern relic. More puzzling, we find a "wall" of hot gas east of the cluster centre. A partial elliptical ring of hot plasma appears to be present around the merger. While radio observations and numerical simulations predict a simple merger geometry, the X-ray results point towards a more complex merger scenario.

Deep radio observations of galaxy clusters have revealed the existence of diffuse radio sources related to the presence of relativistic electrons and weak magnetic fields in the intracluster volume. The role played by this non-thermal intracluster component on the thermodynamical evolution of galaxy clusters is debated, with important implications for cosmological and astrophysical studies of the largest gravitationally bound structures of the Universe. The low surface brightness and steep spectra of diffuse cluster radio sources make them more easily detectable at low-frequencies. LOFAR is the first instrument able to detect diffuse radio emission in hundreds of massive galaxy clusters up to their formation epoch. We present the first observations of clusters imaged by LOFAR and the huge perspectives opened by this instrument for non-thermal cluster studies.

Multiwavelength studies of radio relics at merger shocks set powerful constraints on the relics origin and formation mechanism. However, for X-ray observations, a main difficulty is represented by the low X-ray surface brightness far out in the cluster outskirts, where relics are typically found. Here, we present XMM-Newton results from a 130-ks observation of CIZA J2242.8+5301, a cluster at z=0.19 that hosts a double radio relic. We focus on the well-defined northern relic. There is a difference of ~55% between the temperature we measure behind the relic, and the temperature measured with Suzaku. We analyse the reasons for this large discrepancy, and discuss the possibility of reliably measuring the temperature beyond the northern relic.

Merger shocks induce turbulence in the intra-cluster medium (ICM), and, under some circumstances, accelerate electrons to relativistic velocities to form so-called radio relics. Relics are mostly found at the periphery of galaxy clusters and appear to have magnetic fields at the microGauss level. Here we investigate the possible origins of these magnetic fields. Turbulence produced by the shock itself cannot explain the magnitude of these fields. However, we argue that if the turbulent pressure support in the ICM upstream of the merger shock is of the order of 10 to 30 percent of the total pressure on scales of a few times 100 kpc, then vorticity generated by compressive and baroclinic effects across the shock discontinuity can lead to a sufficient amplification of the magnetic field. Compressional amplification can explain the large polarisation of the radio emission more easily than dynamo turbulent amplification. Finally, clumping of the ICM is shown to have a negligible effect on magnetic field amplification.

We investigate the origin and nature of the multiple sloshing cold fronts in the core of Abell 496 by direct comparison between observations and dedicated hydrodynamical simulations. Our simulations model a minor merger with a 4{\times}10^13M{\circ} subcluster crossing A496 from the south-west to the north-north-east, passing the cluster core in the south-east at a pericentre distance 100 to a few 100 kpc about 0.6 to 0.8 Gyr ago. The gas sloshing triggered by the merger can reproduce almost all observed features, e.g. the characteristic spiral-like brightness residual distribution in the cluster centre and its asymmetry out to 500 kpc, also the positions of and contrasts across the cold fronts. If the subcluster passes close (100 kpc) to the cluster core, the resulting shear flows are strong enough to trigger Kelvin-Helmholtz instabilities that in projection resemble the peculiar kinks in the cold fronts of Abell 496. Finally, we show that sloshing does not lead to a significant modification of the global ICM profiles but a mild oscillation around the initial profiles.

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

Utilizing $\sim 50$ ks of Chandra X-ray Observatory imaging, we present an analysis of the intracluster medium (ICM) and cavity system in the galaxy cluster RBS 797. In addition to the two previously known cavities in the cluster core, the new and deeper X-ray image has revealed additional structure associated with the active galactic nucleus (AGN). The surface brightness decrements of the two cavities are unusually large, and are consistent with elongated cavities lying close to our line-of-sight. We estimate a total AGN outburst energy and mean jet power of $\approx 3 - 6 \times 10^{60}$ erg and $\approx 3 - 6 \times 10^{45}$ erg s$^{-1}$, respectively, depending on the assumed geometrical configuration of the cavities. Thus, RBS 797 is apparently among the the most powerful AGN outbursts known in a cluster. The average mass accretion rate needed to power the AGN by accretion alone is $\sim 1 M_{\odot}$ yr$^{-1}$. We show that accretion of cold gas onto the AGN at this level is plausible, but that Bondi accretion of the hot atmosphere is probably not. The BCG harbors an unresolved, non-thermal nuclear X-ray source with a bolometric luminosity of $\approx 2 \times 10^{44}$ erg s$^{-1}$. The nuclear emission is probably associated with a rapidly-accreting, radiatively inefficient accretion flow. We present tentative evidence that star formation in the BCG is being triggered by the radio jets and suggest that the cavities may be driving weak shocks ($M \sim 1.5$) into the ICM, similar to the process in the galaxy cluster MS 0735.6+7421.

We study the properties of chaotic motions in the intra cluster medium using a set of 20 galaxy clusters simulated with large dynamical range, using the Adaptive Mesh Refinement code ENZO (e.g. Norman et al.2007). The adopted setup allows us to study the spectral and spatial properties of turbulent motions in galaxy clusters with unprecedented detail, achieving an maximum available Reynolds number of the order of R=500-1000 for the largest eddies. The correlations between the energy of these motions in the Intra Cluster Medium and the dynamical state of the host systems are studied, and the statistical properties of turbulent motions and their evolution with time support that major merger events are responsible for the injection of the bulk of turbulent kinetic energy inside cluster. Turbulence is found to account for a 20-30 per cent of the thermal energy in merging clusters, while it accounts for a 5 per cent in relaxed clusters. A comparison of the energies of turbulence and motions in our simulated clusters with present upper-limits in real nearby clusters, recently derived with XMM-Newton (Sanders et al.2010), is provided. When the same spatial scales of turbulent motions are compared, the data from simulations result well within the range presently allowed by observations. Finally, we comment on the possibility that turbulence may accelerate relativistic particles leading to the formation of giant radio halos in turbulent (merging) clusters. Based on our simulations we confirm previous semi-analytical studies that suggest that the fraction of turbulent clusters is consistent with that of clusters hosting radio halos.

The temperature profile of hot gas in galaxies and galaxy clusters is largely determined by the depth of the total gravitational potential and thereby by the dark matter (DM) distribution. We use high-resolution hydrodynamical simulations of galaxy formation to derive a surprisingly simple relation between the gas temperature and DM properties. We show that this relation holds not just for galaxy clusters but also for equilibrated and relaxed galaxies at radii beyond the central stellar-dominated region of typically a few kpc. It is then clarified how a measurement of the temperature and density of the hot gas component can lead to an indirect measurement of the DM velocity anisotropy in galaxies. We also study the temperature relation for galaxy clusters in the presence of self-regulated, recurrent active galactic nuclei (AGN), and demonstrate that this temperature relation even holds outside the inner region of 30 kpc in clusters with an active AGN.

(abridged) We perform hydrodynamical simulations of minor-merger induced gas sloshing and the subsequent formation of cold fronts in the Virgo cluster. We show for the first time that sloshing reproduces all characteristics of the observed cold fronts quantitatively, and we suggest a third cold front at 20 kpc NW of the Virgo core. We identify several new features typical for sloshing cold fronts, most importantly a large-scale brightness asymmetry. We can trace these new features not only in Virgo, but also in other sloshing cold front clusters. By comparing synthetic and real observations, we estimate that the original minor merger event took place about 1.5 Gyr ago when a subcluster of 2-4 \times 10^13 M\odot passed the Virgo core at 100 to 400 kpc distance, where a smaller mass corresponds to a smaller pericentre distance, and vice versa. From the merger geometry, we derive the current location of the disturbing subcluster to be about 1-2 Mpc E of the Virgo core. A possible candidate is M60. Additionally, we quantify the metal redistribution by sloshing and discuss its importance. We verify that the subcluster required to produce the observed cold fronts could be completely ram pressure stripped before reaching the Virgo centre, and discuss the conditions required for this to be achieved. Finally, we demonstrate that the bow shock of a fast galaxy passing the Virgo cluster at ~ 400 kpc distance also causes sloshing and leads to very similar cold front structures. The responsible galaxy would be located about 2 Mpc north of the Virgo centre. A possible candidate is M85.

We present multi-wavelength observations of the centre of RXCJ1504.1-0248 - the galaxy cluster with the most luminous and relatively nearby cool core at z~0.2. Although there are several galaxies within 100 kpc of the cluster core, only the brightest cluster galaxy (BCG), which lies at the peak of the X-ray emission, has blue colours and strong line-emission. Approximately 80 Msun/yr of intracluster gas is cooling below X-ray emitting temperatures, similar to the observed UV star formation rate of ~140 Msun/yr. Most star formation occurs in the core of the BCG and in a 42 kpc long filament of blue continuum, line emission, and X-ray emission, that extends southwest of the galaxy. The surrounding filamentary nebula is the most luminous around any observed BCG. The number of ionizing stars in the BCG is barely sufficient to ionize and heat the nebula, and the line ratios indicate an additional heat source is needed. This heat source can contribute to the H\alpha-deduced star formation rates (SFRs) in BCGs and therefore the derived SFRs should only be considered upper limits. AGN feedback can slow down the cooling flow to the observed mass deposition rate if the black hole accretion rate is of the order of 0.5 Msun/yr at 10% energy output efficiency. The average turbulent velocity of the nebula is vturb ~325 km/s which, if shared by the hot gas, limits the ratio of turbulent to thermal energy of the intracluster medium to less than 6%.

We present three-dimensional, adaptive mesh simulations of dwarf galaxy out- flows driven by supersonic turbulence. Here we develop a subgrid model to track not only the thermal and bulk velocities of the gas, but also its turbulent velocities and length scales. This allows us to deposit energy from supernovae directly into supersonic turbulence, which acts on scales much larger than a particle mean free path, but much smaller than resolved large-scale flows. Unlike previous approaches, we are able to simulate a starbursting galaxy modeled after NGC 1569, with realistic radiative cooling throughout the simulation. Pockets of hot, diffuse gas around individual OB associations sweep up thick shells of material that persist for long times due to the cooling instability. The overlapping of high-pressure, rarefied regions leads to a collective central outflow that escapes the galaxy by eating away at the exterior gas through turbulent mixing, rather than gathering it into a thin, unstable shell. Supersonic, turbulent gas naturally avoids dense regions where turbulence decays quickly and cooling times are short, and this further enhances density contrasts throughout the galaxy- leading to a complex, chaotic distribution of bubbles, loops and filaments as observed in NGC 1569 and other outflowing starbursts.

Cool cores of galaxy clusters are thought to be heated by low-power active galactic nuclei (AGN), whose accretion is regulated by feedback. However, the interaction between the hot gas ejected by the AGN and the ambient intracluster medium is extremely difficult to simulate, as it involves a wide range of spatial scales and gas that is Rayleigh-Taylor (RT) unstable. Here we use a subgrid model for RT-driven turbulence to overcome these problems and present the first observationally-consistent hydrodynamical simulations of AGN self-regulation in galaxy clusters. For a wide range of parameter choices the cluster in our three-dimensional simulations regulates itself for at least several Gyrs years. Heating balances cooling through a string of outbreaks with a typical recurrence time of approximately 80 Myrs, a timescale that depends only on the global cluster properties.

The physics of X-ray cavities in galaxy clusters is constrained by their observed morphological evolution, which depends on such poorly-understood properties as the turbulent density field and magnetic fields. Here we combine numerical simulations that include subgrid turbulence and software that produces synthetic X-ray observations to examine the evolution of X-ray cavities in the the absence of magnetic fields. Our results reveal an anisotropic size evolution that is very different from simplified, analytical predictions. These differences highlight some of the key issues that must be accurately quantified when studying AGN-driven cavities, and help to explain why the inferred pV energy in these regions appears to be correlated with their distance from the cluster center. Interpreting X-ray observations will require detailed modeling of effects including mass-entrainment, distortion by drag forces, and projection. Current limitations do not allow a discrimination between purely hydrodynamic and magnetically-dominated models for X-ray cavities.

We present a description of the public code XIM, a virtual X-ray observatory. XIM can be used to convert hydrodynamic simulations of astrophysical objects, such as large scale structure, galaxy clusters, groups, galaxies, supernova remnants, and similar extended objects, into virtual X-ray observations for direct comparison with observations and for post-processing with standard X-ray analysis tools. By default, XIM simulates Chandra and the International X-ray Observatory (IXO), but can accommodate any user-specified telescope parameters and instrument responses. Examples of XIM applications include virtual Chandra imaging of simulated X-ray cavities from AGN feedback in galaxy clusters, kinematic mapping of cluster velocity fields (e.g., due to mergers or AGN feedback), as well as detailed spectral modeling of multi-phase, multi-temperature spectra from space plasmas.

A wide range of recent observations have shown that AGN-driven cavities may provide the energy source that balances the cooling observed in the centres of cool-core galaxy clusters. One tool for better understanding the physics of these cavities is their observed morphological evolution, which is dependent on such poorly-understood properties as the turbulent density field and the impact of magnetic fields. Here we combine numerical simulations that include subgrid turbulence and software that produces synthetic X-ray observations to examine the evolution of X-ray cavities in the absence of magnetic fields. Our results reveal an anisotropic size evolution of that is dramatically different from simplified, analytical predictions. These differences highlight some of the key issues that must be accurately quantified when studying AGN-driven cavities, and help to explain why the inferred pV energy in these regions appears to be correlated with their distance from the cluster center. Interpreting X- ray observations will require detailed modeling of effects including mass-entrainment, distortion by drag forces, and pro jection. Current limitations do not allow a discrimination between purely hydrodynamic and magnetically-dominated models for X-ray cavities.

We analyzed a deep XMM-Newton observation of the cluster of galaxies Hydra A, focusing on the large-scale shock discovered as a surface brightness discontinuity in Chandra images. The shock front can be seen both in the pressure map and in temperature profiles in several sectors. The Mach numbers determined from the temperature jumps are in good agreement with the Mach numbers derived from EPIC/pn surface brightness profiles and previously from Chandra data and are consistent with M~1.3. The estimated shock age in the different sectors using a spherically symmetric point explosion model ranges between 130 and 230 Myr and the outburst energy between 1.5 and 3e61 ergs. The shape of the shock seen in the pressure map can be approximated with an ellipse centered 70 kpc towards the NE from the cluster center. We aimed to develop a better model that can explain the offset between the shock center and the AGN and give a consistent result on the shock age and energy. To this end, we performed 3D hydrodynamical simulations in which the shock is produced by a symmetrical pair of AGN jets launched in a spherical galaxy cluster. As an explanation of the observed offset of the shock center, we consider large-scale bulk flows in the intracluster medium. The simulation successfully reproduces the size, ellipticity, and average Mach number of the observed shock front. The predicted age of the shock is 160 Myr and the total input energy 3e61 erg. Both values are within the range determined by the spherically symmetric model. Matching the observed 70 kpc offset of the shock ellipse from the cluster center requires large-scale coherent motions with a high velocity of 670 km/s. We discuss the feasibility of this scenario and offer alternative ways to produce the offset and to further improve the simulation.

We analyzed global properties, radial profiles and 2D maps of the metal abundances and temperature in the cool core cluster of galaxies Hydra A using a deep XMM-Newton exposure. The best fit among the available spectral models is provided by a Gaussian distribution of the emission measure (gdem). We can accurately determine abundances for 7 elements in the cluster core with EPIC and 3 elements with RGS. The gdem model gives lower Fe abundances than a single temperature model. The abundance profiles for Fe, Si, S, but also O are centrally peaked. Combining the Hydra A results with 5 other clusters for which detailed chemical abundance studies are available, we find a significant decrease of O with radius, while the increase in the O/Fe ratio with radius is small within 0.1 r_200. We compare the observed abundance ratios with the mixing of various supernova type Ia and core-collapse yield models in different relative amounts. Producing the estimated O, Si and S peaks in Hydra A requires either an amount of metals ejected by stellar winds 3-8 times higher than predicted by available models or a remaining peak in the enrichment by core-collapse supernovae from the protocluster phase. The temperature map shows cooler gas extending in arm-like structures towards the north and south. These structures appear to be richer in metals than the ambient medium and spatially correlated with the large-scale radio lobes. We estimate the mass of cool gas, which was probably uplifted by buoyant bubbles of relativistic plasma produced by the AGN, to 1.6-6.1x10^9 M_sun, and the energy associated with this uplift to 3.3-12.5x10^58 ergs. The best estimate of the mass of Fe uplifted together with the cool gas is 1.7x10^7 M_sun, 15% of the total mass of Fe in the central 0.5arcmin region.

We present a numerical study of the mass transfer dynamics prior to the gravitational wave-driven merger of a double white dwarf system. Recently, there has been some discussion about the dynamics of these last stages, different methods seemed to provide qualitatively different results. While earlier SPH simulations indicated a very quick disruption of the binary on roughly the orbital time scale, more recent grid-based calculations find long-lived mass transfer for many orbital periods. Here we demonstrate how sensitive the dynamics of this last stage is to the exact initial conditions. We show that, after a careful preparation of the initial conditions, the reportedly short-lived systems undergo mass transfer for many dozens of orbits. The reported numbers of orbits are resolution-biased and therefore represent only lower limits to what is realized in nature. Nevertheless, the study shows convincingly the convergence of different methods to very similar results.

Hot, underdense bubbles powered by active galactic nuclei (AGN) are likely to play a key role in halting catastrophic cooling in the centers of cool-core galaxy clusters. We present three-dimensional simulations that capture the evolution of such bubbles, using an adaptive-mesh hydrodynamic code, FLASH3, to which we have added a subgrid model of turbulence and mixing. While pure-hydro simulations indicate that AGN bubbles are disrupted into resolution-dependent pockets of underdense gas, proper modeling of subgrid turbulence indicates that this a poor approximation to a turbulent cascade that continues far beyond the resolution limit. Instead, Rayleigh-Taylor instabilities act to effectively mix the heated region with its surroundings, while at the same time preserving it as a coherent structure, consistent with observations. Thus bubbles are transformed into hot clouds of mixed material as they move outwards in the hydrostatic intracluster medium (ICM), much as large airbursts lead to a distinctive ``mushroom cloud'' structure as they rise in the hydrostatic atmosphere of Earth. Properly capturing the evolution of such clouds has important implications for many ICM properties. In particular, it significantly changes the impact of AGN on the distribution of entropy and metals in cool-core clusters such as Perseus.

The onset of the thermonuclear runaway in a Chandrasekhar-mass white dwarf, leading to the explosion as a type Ia supernova, is studied with hydrodynamical simulations. We investigate the evolution of temperature fluctuations (``bubbles'') in the WD's convective core by means of 2D numerical simulations. We show how the occurrence of the thermonuclear runaway depends on various bubble parameters. The relevance of the progenitor's composition for the ignition process is also discussed.

M87 is a key object whose study can reveal the complex phenomena in cooling cores. We use a deep XMM-Newton observation of M87 to produce detailed temperature, pressure and entropy maps in order to analyze the physical processes of cooling cores and of their heating mechanisms. We employed both broad-band fitting and full spectroscopical one-temperature model analysis to derive temperature and surface brightness maps, from which the pseudo-deprojected entropy and pressure were calculated. We discuss possible physical interpretations of small deviations from radial and elliptical symmetry in these maps. The most prominent features observed are the E and SW X-ray arms that coincide with powerful radio lobes, a weak shock at a radius of 3', an overall ellipticity in the pressure map and a NW/SE asymmetry in the entropy map which we associate with the motion of the galaxy towards the NW. For the first time we find evidence that cold, metal-rich gas is being transported out of the center, possibly through bubble-induced mixing. Several edges in the abundance map indicate an oscillation of the galaxy along the NW/SE direction. Furthermore, the radio lobes appear to rise along the short axis of the elliptical pressure distribution, following the steepest gradient of the gravitational potential, and seem to contain a nonthermal pressure component.

In the framework of the Chandrasekhar-mass deflagration model for Type Ia supernovae (SNe Ia), a persisting free parameter is the initial morphology of the flame front, which is linked to the ignition process in the progenitor white dwarf. Previous analytical models indicate that the thermal runaway is driven by temperature perturbations (''bubbles'') that develop in the white dwarf's convective core. In order to probe the conditions at ignition (diameters, temperatures and evolutionary timescales), we have performed hydrodynamical 2D simulations of buoyant bubbles in white dwarf interiors. Our results show that fragmentation occurring during the bubble rise affects the outcome of the bubble evolution. Possible implications for the ignition process of SNe Ia are discussed.

The physics of the intracluster medium, in particular the values for the thermal conductivity and the viscosity are largely unknown and subject to an ongoing debate. Here, we study the effect of viscosity on the thermal state of the intracluster medium using three-dimensional cosmological simulations of structure formation. It is shown that viscosity, provided it is not too far off from the unmagnetised Spitzer value, has a significant effect on cluster profiles. In particular, it aids in heating the cool cores of clusters. The central cooling time of the most massive clusters in our simulation is increased by more than an order of magnitude. In large clusters, viscous heating may help to establish an entropy floor and to prevent a cooling catastrophe.

The morphology of the X-ray and radio emitting features in the central $\sim$ 50 kpc region around the galaxy M87 strongly suggests that buoyant bubbles of cosmic rays (inflated by an earlier nuclear active phase of the galaxy) rise through the cooling gas at roughly half the sound speed. In the absence of strong surface tension, initially spherical bubbles will transform into tori as they rise through an external medium. Such structures can be identified in the radio images of the halo of M87. During their rise, bubbles will uplift relatively cool X-ray emitting gas from the central regions of the cooling flow to larger distances. This gas is colder than the ambient gas and has a higher volume emissivity. As a result, rising ``radio'' bubbles may be trailed by elongated X-ray features as indeed is observed in M87. We performed simple hydrodynamic simulations to qualitatively illustrate the evolution of buoyant bubbles in the M87 environment.

Radio galaxies are known to inflate lobes of hot relativistic plasmas into the intergalactic medium. Here we present hydrodynamical and magnetohydrodynamical simulations of these hot plasma bubbles in FRII objects. We focus on the later stages of their evolution after the jet has died down and after the bow shock that surrounded the lobes at earlier stages has vanished. We investigate the evolution of the plasma bubbles as they become subject to Rayleigh-Taylor instabilities. From our simulations we calculate the radio and X-ray emissivities of the bubbles and discuss their appearance in observations. Finally, we investigate the influence of large-scale magnetic fields on the evolution of the bubbles. The issues of re-acceleration and diffusion of relativistic particles are briefly discussed.

We present results of three-dimensional (3D) simulations of the magnetohydrodynamic Kelvin-Helmholtz instability in a stratified shear layer. The magnetic field is taken to be uniform and parallel to the shear flow. We describe the evolution of the fluid flow and the magnetic field for a range of initial conditions. In particular, we investigate how the mixing rate of the fluid depends on the Richardson number and the magnetic field strength. It was found that the magnetic field can enhance as well as suppress mixing. Moreover, we have performed two-dimensional (2D) simulations and discuss some interesting differences between the 2D and 3D results.

In this paper we present the results of numerical simulations of the Kelvin-Helmholtz instability in a stratified shear layer. This shear instability is believed to be responsible for extra mixing in differentially rotating stellar interiors and is the prime candidate to explain the abundance anomalies observed in many rotating stars. All mixing prescriptions currently in use are based on phenomenological and heuristic estimates whose validity is often unclear. Using three-dimensional numerical simulations, we study the mixing efficiency as a function of the Richardson number and compare our results with some semi-analytical formalisms of mixing.

Inhomogeneities in wave propagation conditions near and below the solar surface have been detected by means of time-distance helioseismology. Here we calculate the effect of temperature inhomogeneities on the travel times of sound waves. A temperature increase, e.g. in active regions, not only increases the sound speed but also lengthens the path along which the wave travels because the expansion of the heated layers shifts the upper turning of the waves upward. Using a ray tracing approximation we find that in many cases the net effect of a temperature enhancement is an {\it increase} of the travel times. We argue that the reduced travel times that are observed are caused by a combination of magnetic fields in the active region and {\it reduced} subsurface temperatures. Such a reduction may be related to the increased radiative energy loss from small magnetic flux tubes.

In the solar convection zone acoustic waves are scattered by turbulent sound speed fluctuations. In this paper the scattering of waves by convective cells is treated using Rytov's technique. Particular care is taken to include diffraction effects which are important especially for high-degree modes that are confined to the surface layers of the Sun. The scattering leads to damping of the waves and causes a phase shift. Damping manifests itself in the width of the spectral peak of p-mode eigenfrequencies. The contribution of scattering to the line widths is estimated and the sensitivity of the results on the assumed spectrum of the turbulence is studied. Finally the theoretical predictions are compared with recently measured line widths of high-degree modes.

In this paper it is investigated whether torsion, which arises naturally in most theories of quantum gravity, has observable implications for the Big Bang nucleosynthesis. Torsion can lead to spin flips amongst neutrinos thus turning them into sterile neutrinos. In the early Universe they can alter the helium abundance which is tightly constrained by observations. Here I calculate to what extent torsion of the string theory type leads to a disagreement with the Big Bang nucleosynthesis predictions.

There has been some controversy over the expression for the so-called `interaction energy' due to screening of charged particles in a plasma. Even in the relatively simple case of weak screening, first discussed in the context of astrophysical plasmas by Salpeter (1954), there is disagreement. In particular, Shaviv and Shaviv (1996) have claimed recently that by not considering explicitly in his calculation the complete screening cloud, Salpeter obtained a result for the interaction energy between two nuclei separated by a distance $r$ which in the limit $r \to 0$ is only 2/3 the correct value. It appears that this claim has arisen from a fundamental misconception concerning the dynamics of the interaction. We rectify this misconception, and show that Salpeter's formula is indeed correct.