In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.)

We report the detection of diffuse radio emission which might be connected to a large-scale filament of the cosmic web covering a 8deg x 8deg area in the sky, likely associated with a z~0.1 over-density traced by nine massive galaxy clusters. In this work, we present radio observations of this region taken with the Sardinia Radio Telescope. Two of the clusters in the field host a powerful radio halo sustained by violent ongoing mergers and provide direct proof of intra-cluster magnetic fields. In order to investigate the presence of large-scale diffuse radio synchrotron emission in and beyond the galaxy clusters in this complex system, we combined the data taken at 1.4 GHz obtained with the Sardinia Radio Telescope with higher resolution data taken with the NRAO VLA Sky Survey. We found 28 candidate new sources with a size larger and X-ray emission fainter than known diffuse large-scale synchrotron cluster sources for a given radio power. This new population is potentially the tip of the iceberg of a class of diffuse large-scale synchrotron sources associated with the filaments of the cosmic web. In addition, we found in the field a candidate new giant radio galaxy.

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.

Galaxy clusters are the most massive constituents of the large-scale structure of the Universe. While the hot thermal gas that pervades galaxy clusters is relatively well understood through observations with X-ray satellites, our understanding of the non-thermal part of the intra-cluster medium remains incomplete. With LOFAR and GMRT observations, we have identified a phenomenon that can be unveiled only at extremely low radio-frequencies and offers new insights into the non-thermal component. We propose that the interplay between radio-emitting plasma and the perturbed intra-cluster medium can gently re-energise relativistic particles initially injected by active galactic nuclei. Sources powered through this mechanism can maintain electrons at higher energies than radiative ageing would allow. If this mechanism is common for aged plasma, a population of mildly relativistic electrons can be accumulated inside galaxy clusters providing the seed population for merger-induced re-acceleration mechanisms on larger scales such as turbulence and shock waves.

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be $n_\mathrm{s} = 0.968 \pm 0.006$ and tightly constrain its scale dependence to $d n_s/d \ln k =-0.003 \pm 0.007$ when combined with the Planck lensing likelihood. When the high-$\ell$ polarization data is included, the results are consistent and uncertainties are reduced. The upper bound on the tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode polarization constraint $r< 0.12$ (95% CL) obtained from a joint BICEP2/Keck Array and Planck analysis. These results imply that $V(\phi) \propto \phi^2$ and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as $R^2$ inflation. Three independent methods reconstructing the primordial power spectrum are investigated. The Planck data are consistent with adiabatic primordial perturbations. We investigate inflationary models producing an anisotropic modulation of the primordial curvature power spectrum as well as generalized models of inflation not governed by a scalar field with a canonical kinetic term. The 2015 results are consistent with the 2013 analysis based on the nominal mission data.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We study the intra-cluster magnetic field in the poor galaxy cluster Abell 194 by complementing radio data, at different frequencies, with data in the optical and X-ray bands. We analyze new total intensity and polarization observations of Abell 194 obtained with the Sardinia Radio Telescope (SRT). We use the SRT data in combination with archival Very Large Array observations to derive both the spectral aging and Rotation Measure (RM) images of the radio galaxies 3C40A and 3C40B embedded in Abell 194. The optical analysis indicates that Abell 194 does not show a major and recent cluster merger, but rather agrees with a scenario of accretion of small groups. Under the minimum energy assumption, the lifetimes of synchrotron electrons in 3C40B measured from the spectral break are found to be 157 Myrs. The break frequency image and the electron density profile inferred from the X-ray emission are used in combination with the RM data to constrain the intra-cluster magnetic field power spectrum. By assuming a Kolmogorov power law power spectrum, we find that the RM data in Abell 194 are well described by a magnetic field with a maximum scale of fluctuations of Lambda_max=64 kpc and a central magnetic field strength of <B0>=1.5 microG. Further out, the field decreases with the radius following the gas density to the power of eta=1.1. Comparing Abell 194 with a small sample of galaxy clusters, there is a hint of a trend between central electron densities and magnetic field strengths.

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

In this paper, we investigate numerically the stochastic ABC model, a toy model in the theory of astrophysical kinematic dynamos, within the recently proposed supersymmetric theory of stochastics (STS). STS characterises stochastic differential equations (SDEs) by the spectrum of the stochastic evolution operator (SEO) on elements of the exterior algebra or differentials forms over the system's phase space, X. STS can thereby classify SDEs as chaotic or non-chaotic by identifying the phenomenon of stochastic chaos with the spontaneously broken topological supersymmetry that all SDEs possess. We demonstrate the following three properties of the SEO, deduced previously analytically and from physical arguments: the SEO spectra for zeroth and top degree forms never break topological supersymmetry, all SDEs possesses pseudo-time-reversal symmetry, and each de Rahm cohomology class provides one supersymmetric eigenstate. Our results also suggests that the SEO spectra for forms of complementary degrees, i.e., k and dim X -k, may be isospectral.

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 present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

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 a new algorithm to reconstruct the Galactic free electron density from pulsar dispersion measures. The algorithm performs a nonparametric tomography for a density field with an arbitrary amount of degrees of freedom. It is based on approximating the Galactic free electron density as the product of a profile function with a statistically isotropic and homogeneous log-normal field. Under this approximation the algorithm generates a map of the free electron density as well as an uncertainty estimate without the need of information about the power spectrum. The uncertainties of the pulsar distances are treated consistently by an iterative procedure. We test the algorithm using the NE2001 model with modified fluctuations as a Galaxy model, pulsar populations generated from the Lorimer population model, and mock observations emulating the upcoming Square Kilometer Array. We show the quality of the reconstruction for mock data sets containing between 1000 and 10000 pulsars with distance uncertainties up to 25%. Our results show, that with the SKA nonparametric tomography of the Galactic free electron density becomes feasible, but the quality of the reconstruction is very sensitive to the distance uncertainties.

Determining magnetic field properties in different environments of the cosmic large-scale structure as well as their evolution over redshift is a fundamental step toward uncovering the origin of cosmic magnetic fields. Radio observations permit the study of extragalactic magnetic fields via measurements of the Faraday depth of extragalactic radio sources. Our aim is to investigate how much different extragalactic environments contribute to the Faraday depth variance of these sources. We develop a Bayesian algorithm to distinguish statistically Faraday depth variance contributions intrinsic to the source from those due to the medium between the source and the observer. In our algorithm the Galactic foreground and the measurement noise are taken into account as the uncertainty correlations of the galactic model. Additionally, our algorithm allows for the investigation of possible redshift evolution of the extragalactic contribution. This work presents the derivation of the algorithm and tests performed on mock observations. With cosmic magnetism being one of the key science projects of the new generation of radio interferometers we have made predictions for the algorithm's performance on data from the next generation of radio interferometers. Applications to real data are left for future work.

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70GHz and 100GHz channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70/100/143GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8deg-15arcmin are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement between 70 and 100GHz power spectra averaged over 70<l<390 at the 0.8% level, and agreement between 143 and 100GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100/143/217GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described...(abridged)

The Planck High Frequency Instrument (HFI) has observed the full sky at six frequencies (100, 143, 217, 353, 545, and 857 GHz) in intensity and at four frequencies in linear polarization (100, 143, 217, and 353 GHz). In order to obtain sky maps, the time-ordered information (TOI) containing the detector and pointing samples must be processed and the angular response must be assessed. The full mission TOI is included in the Planck 2015 release. This paper describes the HFI TOI and beam processing for the 2015 release. HFI calibration and map-making are described in a companion paper. The main pipeline has been modified since the last release (2013 nominal mission in intensity only), by including a correction for the non-linearity of the warm readout and by improving the model of the bolometer time response. The beam processing is an essential tool that derives the angular response used in all the Planck science papers and we report an improvement in the effective beam window function uncertainty of more than a factor 10 relative to the 2013 release. Noise correlations introduced by pipeline filtering function are assessed using dedicated simulations. Angular cross-power spectra using datasets that are decorrelated in time are immune to the main systematic effects.

Shortened abstract: Observations of the nearby Chamaeleon clouds in gamma rays with the Fermi Large Area Telescope and in thermal dust emission with Planck and IRAS have been used with the HI and CO radio data to (i) map the gas column densities in the different phases and at the dark neutral medium (DNM) transition between the HI-bright and CO-bright media; (ii) constrain the CO-to-$H_2$ conversion factor, $X_{CO}$; (iii) probe the dust properties per gas nucleon in each gas phase and spatially across the clouds. We have separated clouds in velocity in HI and CO emission and modelled the 0.4-100 GeV intensity, the dust optical depth at 353 GHz, the thermal radiance of the large grains, and an estimate of the dust extinction empirically corrected for the starlight intensity, $A_{VQ}$. The gamma-ray emissivity spectra confirm that the GeV-TeV cosmic rays uniformly permeate all gas phases up to the CO cores. The dust and cosmic rays reveal large amounts of DNM gas, with comparable spatial distributions and twice as much mass as in the CO-bright clouds. We give constraints on the HI-DNM-CO transitions and CO-dark $H_2$ fractions for separate clouds. The corrected extinction provides the best fit to the total gas traced by the gamma rays, but we find evidence for a rise in $A_{VQ}/N_H$ and a steep rise in opacity, with increasing $N_H$ and $H_2$ fraction, and with decreasing dust temperature. We observe less variations for the specific power of the grains, except for a decline by half in the CO cores. This combined information suggests grain evolution. The gamma rays and dust radiance yield consistent $X_{CO}$ estimates near $0.7\times10^{20}$ cm$^{-2}$ (K km/s)$^{-1}$. The other dust tracers yield biased values because of the grain opacity rise in the CO clouds. These results also confirm a factor of 2 difference between $X_{CO}$ estimates at pc and kpc scales.

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to $n_s = 0.9603 \pm 0.0073$, ruling out exact scale invariance at over 5 $\sigma$. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining $d n_s/d ln k = -0.0134 \pm 0.0090$. Several analyses dropping the slow-roll approximation are carried out, including detailed model comparison and inflationary potential reconstruction. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit $\chi^2$ by ~ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with bounds on $f_\mathrm{NL}$ measured by Planck. The fractional primordial contribution of CDM isocurvature modes in the curvaton and axion scenarios has upper bounds of 0.25% or 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelation can improve $\chi^2$ by approximatively 4 as a result of slightly lowering the theoretical prediction for the $\ell<40$ multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

To better understand the origin and properties of cosmological magnetic fields, a detailed knowledge of magnetic fields in the large-scale structure of the Universe (galaxy clusters, filaments) is crucial. We propose a new statistical approach to study magnetic fields on large scales with the rotation measure grid data that will be obtained with the new generation of radio interferometers.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

Using \Planck\ data combined with the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we address the study of peculiar motions by searching for evidence of the kinetic Sunyaev-Zeldovich effect (kSZ). By implementing various filters designed to extract the kSZ generated at the positions of the clusters, we obtain consistent constraints on the radial peculiar velocity average, root mean square (rms), and local bulk flow amplitude at different depths. For the whole cluster sample of average redshift 0.18, the measured average radial peculiar velocity with respect to the cosmic microwave background (CMB) radiation at that redshift, i.e., the kSZ monopole, amounts to $72 \pm 60$ km s$^{-1}$. This constitutes less than 1% of the relative Hubble velocity of the cluster sample with respect to our local CMB frame. While the linear $\Lambda$CDM prediction for the typical cluster radial velocity rms at $z=0.15$ is close to 230km s$^{-1}$, the upper limit imposed by \Planck\ data on the cluster subsample corresponds to 800 km s$^{-1}$ at 95% confidence level, i.e., about three times higher. \Planck\ data also set strong constraints on the local bulk flow in volumes centred on the Local Group. There is no detection of bulk flow as measured in any comoving sphere extending to the maximum redshift covered by the cluster sample. A blind search for bulk flows in this sample has an upper limit of 254 km s$^{-1}$ (95% confidence level) dominated by CMB confusion and instrumental noise, indicating that the Universe is largely homogeneous on Gpc scales. In this context, in conjunction with supernova observations, \Planck\ is able to rule out a large class of inhomogeneous void models as alternatives to dark energy or modified gravity. The \Planck\ constraints on peculiar velocities and bulk flows are thus consistent with the $\Lambda$CDM scenario.

The experimental situation of Dark Matter Direct Detection has reached an exciting cross-roads, with potential hints of a discovery of Dark Matter (DM) from the CDMS, CoGeNT, CRESST-II and DAMA experiments in tension with null-results from xenon-based experiments such as XENON100 and LUX. Given the present controversial experimental status, it is important that the analytical method used to search for DM in Direct Detection experiments is both robust and flexible enough to deal with data for which the distinction between signal and background points is difficult, and hence where the choice between setting a limit or defining a discovery region is debatable. In this article we propose a novel (Bayesian) analytical method, which can be applied to all Direct Detection experiments and which extracts the maximum amount of information from the data. We apply our method to the XENON100 experiment data as a worked example, and show that firstly our exclusion limit at 90% confidence is in agreement with their own for the 225 Live Days data, but is several times stronger for the 100 Live Days data. Secondly we find that, due to the two points at low values of S1 and S2 in the 225 days data-set, our analysis points to either weak consistency with low-mass Dark Matter or the possible presence of an unknown background. Given the null-result from LUX, the latter scenario seems the more plausible.

We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged)

Taking advantage of the all-sky coverage and broad frequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R500. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3 x R500, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and X-ray measurements inside R500. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from XMM-Newton in the region in which the profiles overlap (i.e., [0.1-1]R500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is slightly flatter than most predictions from numerical simulations. Combining the SZ and X-ray observed profiles into a joint fit to a generalised pressure profile gives best-fit parameters [P0, c500, gamma, alpha, beta] = [6.41, 1.81, 0.31, 1.33, 4.13]. Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3 x R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

In 2011, the XENON100 experiment has set unprecedented constraints on dark matter-nucleon interactions, excluding dark matter candidates with masses down to 6 GeV if the corresponding cross section is larger than 10^{-39} cm^2. The dependence of the exclusion limit in terms of the scintillation efficiency (Leff) has been debated at length. To overcome possible criticisms XENON100 performed an analysis in which Leff was considered as a nuisance parameter and its uncertainties were profiled out by using a Gaussian likelihood in which the mean value corresponds to the best fit Leff value smoothly extrapolated to zero below 3 keVnr. Although such a method seems fairly robust, it does not account for more extreme types of extrapolation nor does it enable to anticipate on how much the exclusion limit would vary if new data were to support a flat behaviour for Leff below 3 keVnr, for example. Yet, such a question is crucial for light dark matter models which are close to the published XENON100 limit. To answer this issue, we use a maximum Likelihood ratio analysis, as done by the XENON100 collaboration, but do not consider Leff as a nuisance parameter. Instead, Leff is obtained directly from the fits to the data. This enables us to define frequentist confidence intervals by marginalising over Leff.

Galactic winds are observed in many spiral galaxies with sizes from dwarfs up to the Milky Way, and they sometimes carry a mass in excess of that of newly formed stars by up to a factor of ten. Multiple driving processes of such winds have been proposed, including thermal pressure due to supernova-heating, UV radiation pressure on dust grains, or cosmic ray (CR) pressure. We here study wind formation due to CR physics using a numerical model that accounts for CR acceleration by supernovae, CR thermalization, and advective CR transport. In addition, we introduce a novel implementation of CR streaming relative to the rest frame of the gas. We find that CR streaming drives powerful and sustained winds in galaxies with virial masses M_200 < 10^{11} Msun. In dwarf galaxies (M_200 ~ 10^9 Msun) the winds reach a mass loading factor of ~5, expel ~60 per cent of the initial baryonic mass contained inside the halo's virial radius and suppress the star formation rate by a factor of ~5. In dwarfs, the winds are spherically symmetric while in larger galaxies the outflows transition to bi-conical morphologies that are aligned with the disc's angular momentum axis. We show that damping of Alfven waves excited by streaming CRs provides a means of heating the outflows to temperatures that scale with the square of the escape speed. In larger haloes (M_200 > 10^{11} Msun), CR streaming is able to drive fountain flows that excite turbulence. For halo masses M_200 > 10^{10} Msun, we predict an observable level of H-alpha and X-ray emission from the heated halo gas. We conclude that CR-driven winds should be crucial in suppressing and regulating the first epoch of galaxy formation, expelling a large fraction of baryons, and - by extension - aid in shaping the faint end of the galaxy luminosity function. They should then also be responsible for much of the metal enrichment of the intergalactic medium.

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

(abridged) Planck has observed the entire sky from 30 GHz to 857GHz. The observed foreground emission contains contributions from different phases of the interstellar medium (ISM). We have separated the observed Galactic emission into the different gaseous components (atomic, molecular and ionised) in each of a number of Galactocentric rings. Templates are created for various Galactocentric radii using velocity information from atomic (neutral hydrogen) and molecular (12CO) observations. The ionised template is assumed to be traced by free-free emission as observed by WMAP, while 408 MHz emission is used to trace the synchrotron component. Gas emission not traced by the above templates, namely "ark gas", as evidenced using Planck data, is included as an additional template, the first time such a component has been used in this way. These templates are then correlated with each of the Planck frequency bands, as well as other ancillary data. The emission per column density of the gas templates allows us to create distinct spectral energy distributions (SEDs) per Galactocentric ring and in each of the gaseous tracers from 1.4 GHz to 25 THz (12\mu m). Apart from the thermal dust and free-free emission, we have probed the Galaxy for anomalous (e.g., spinning) dust as well as synchrotron emission. We show that anomalous dust emission is present in the atomic, molecular and dark gas phases throughout the Galactic disk. The derived dust propeties associated with the dark gas phase are derived but do not allow us to reveal the nature of this phase. For all environments, the anomalous emission is consistent with rotation from polycyclic aromatic hydrocarbons (PAHs) and, according to our simple model, accounts for $(25\pm5)%$ (statistical) of the total emission at 30 GHz.

We derive a method to reconstruct Gaussian signals from linear measurements with Gaussian noise. This new algorithm is intended for applications in astrophysics and other sciences. The starting point of our considerations is the principle of minimum Gibbs free energy which was previously used to derive a signal reconstruction algorithm handling uncertainties in the signal covariance. We extend this algorithm to simultaneously uncertain noise and signal covariances using the same principles in the derivation. The resulting equations are general enough to be applied in many different contexts. We demonstrate the performance of the algorithm by applying it to specific example situations and compare it to algorithms not allowing for uncertainties in the noise covariance. The results show that the method we suggest performs very well under a variety of circumstances and is indeed qualitatively superior to the other methods in cases where uncertainty in the noise covariance is present.

This paper presents the first results of comparison of Planck along with IRAS data with Green Bank Telescope 21-cm observations in 14 fields covering more than 800 square degrees at high Galactic latitude. Galactic dust emission for fields with average HI column density lower than 2 x 10^20 cm^-2 is well correlated with 21-cm emission. The residual emission in these fields, once the HI-correlated emission is removed, is consistent with the expected statistical properties of the cosmic infrared background fluctuations. Fields with larger column densities show significant excess dust emission compared to the HI column density. Regions of excess lie in organized structures that suggest the presence of hydrogen in molecular form, though they are not always correlated with CO emission. Dust emission from intermediate-velocity clouds is detected with high significance. Its spectral properties are consistent with, compared to the local ISM values, significantly hotter dust (T~20 K), lower sub-millimeter dust opacity, and a relative abundance of very small grains to large grains about four times higher. These results are compatible with expectations for clouds that are part of the Galactic fountain in which there is dust shattering and fragmentation. Correlated dust emission in HVCs is not detected; the average of the 99.9% confidence upper limits to the emissivity is 0.15 times the local ISM value at 857 and 3000 GHz, in accordance with gas phase evidence for lower metallicity and depletion in these clouds. Unexpected anti-correlated variations of the dust temperature and emission cross-section per H atom are identified in the local ISM and IVCs, a trend that continues into molecular environments. This suggests that dust growth through aggregation, seen in molecular clouds, is active much earlier in the cloud condensation and star formation processes.

We present the XMM-Newton follow-up for confirmation of Planck cluster candidates. Twenty-five candidates have been observed to date using snapshot (~10 ksec) exposures, ten as part of a pilot programme to sample a low range of signal-to-noise ratios (4<S/N<6), and a further 15 in a programme to observe a sample of S/N>5 candidates. The sensitivity and spatial resolution of XMM-Newton allows unambiguous discrimination between clusters and false candidates. The 4 false candidates have S/N <= 4.1. A total of 21 candidates are confirmed as extended X-ray sources. Seventeen are single clusters, the majority of which are found to have highly irregular and disturbed morphologies (about ~70%). The remaining four sources are multiple systems, including the unexpected discovery of a supercluster at z=0.45. For 20 sources we are able to derive a redshift estimate from the X-ray Fe K line (albeit of variable quality). The new clusters span the redshift range 0.09 <= z <= 0.54, with a median redshift of z~0.37. A first determination is made of their X-ray properties including the characteristic size, which is used to improve the estimate of the SZ Compton parameter, Y_SZ. The follow-up validation programme has helped to optimise the Planck candidate selection process. It has also provided a preview of the X-ray properties of these newly-discovered clusters, allowing comparison with their SZ properties, and to the X-ray and SZ properties of known clusters observed in the Planck survey. Our results suggest that Planck may have started to reveal a non-negligible population of massive dynamically perturbed objects that is under-represented in X-ray surveys. However, despite their particular properties, these new clusters appear to follow the Y_SZ-Y_X relation established for X-ray selected objects, where Y_X is the product of the gas mass and temperature.

Observations of giant radio halos provide unambiguous evidence for the existence of cosmic ray (CR) electrons and magnetic fields in galaxy clusters. The physical mechanism generating radio halos is still heavily debated. We critically discuss the proposed models for the radio halo emission and highlight the weaknesses underlying each explanation. We present an idea how the interplay of CR propagation and turbulent advection selects a bimodal spatial CR distribution that is characteristic for the dynamical state of a cluster. As a result, strongly turbulent, merging clusters should have a more centrally concentrated CR energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected diffuse radio and gamma ray emission of clusters. Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR electrons.

The optimal reconstruction of cosmic metric perturbations and other signals requires knowledge of their power spectra and other parameters. If these are not known a priori, they have to be measured simultaneously from the same data used for the signal reconstruction. We formulate the general problem of signal inference in the presence of unknown parameters within the framework of information field theory. We develop a generic parameter uncertainty renormalized estimation (PURE) technique and address the problem of reconstructing Gaussian signals with unknown power-spectrum with five different approaches: (i) separate maximum-a-posteriori power spectrum measurement and subsequent reconstruction, (ii) maximum-a-posteriori power reconstruction with marginalized power-spectrum, (iii) maximizing the joint posterior of signal and spectrum, (iv) guessing the spectrum from the variance in the Wiener filter map, and (v) renormalization flow analysis of the field theoretical problem providing the PURE filter. In all cases, the reconstruction can be described or approximated as Wiener filter operations with assumed signal spectra derived from the data according to the same recipe, but with differing coefficients. All of these filters, except the renormalized one, exhibit a perception threshold in case of a Jeffreys prior for the unknown spectrum. Data modes, with variance below this threshold do not affect the signal reconstruction at all. Filter (iv) seems to be similar to the so called Karhune-Loeve and Feldman-Kaiser-Peacock estimators for galaxy power spectra used in cosmology, which therefore should also exhibit a marginal perception threshold if correctly implemented. We present statistical performance tests and show that the PURE filter is superior to the others.

COrE (Cosmic Origins Explorer) is a fourth-generation full-sky, microwave-band satellite recently proposed to ESA within Cosmic Vision 2015-2025. COrE will provide maps of the microwave sky in polarization and temperature in 15 frequency bands, ranging from 45 GHz to 795 GHz, with an angular resolution ranging from 23 arcmin (45 GHz) and 1.3 arcmin (795 GHz) and sensitivities roughly 10 to 30 times better than PLANCK (depending on the frequency channel). The COrE mission will lead to breakthrough science in a wide range of areas, ranging from primordial cosmology to galactic and extragalactic science. COrE is designed to detect the primordial gravitational waves generated during the epoch of cosmic inflation at more than $3\sigma $ for $r=(T/S)>=10^{-3}$. It will also measure the CMB gravitational lensing deflection power spectrum to the cosmic variance limit on all linear scales, allowing us to probe absolute neutrino masses better than laboratory experiments and down to plausible values suggested by the neutrino oscillation data. COrE will also search for primordial non-Gaussianity with significant improvements over Planck in its ability to constrain the shape (and amplitude) of non-Gaussianity. In the areas of galactic and extragalactic science, in its highest frequency channels COrE will provide maps of the galactic polarized dust emission allowing us to map the galactic magnetic field in areas of diffuse emission not otherwise accessible to probe the initial conditions for star formation. COrE will also map the galactic synchrotron emission thirty times better than PLANCK. This White Paper reviews the COrE science program, our simulations on foreground subtraction, and the proposed instrumental configuration.

We develop a novel maximum-a-posteriori method to measure magnetic power spectra from Faraday rotation data and implement it in the REALMAF code. A sophisticated model for the magnetic autocorrelation in real space permits us to alleviate previously required simplifying assumptions in the processing. We also introduce a way to treat the divergence relation of the magnetic field with a multiplicative factor in Fourier space, with which we can model the magnetic autocorrelation as a spherically symmetric function. Applied to the dataset of Hydra A north, we find a power law power spectrum on spatial scales between 0.3 kpc and 8 kpc, with no visible turnover at large scales within this range and a spectral index consistent with a Kolmogorov-like power law regime. The magnetic field strength profile seems to follow the electron density profile with an index alpha=1. A variation of alpha from 0.5 to 1.5 would lead to a spectral index between 1.55 and 2.05. The extrapolated magnetic field strength in the cluster centre highly depends on the assumed projection angle of the jet. For an angle of 45 degree we derive extrapolated 36 muG in the centre and directly probed 16 muG at 50 kpc radius.

Deep X-ray data from CHANDRA and XMM-Newton, and GMRT radio data are presented for ZwCl 2341.1+0000, an extremely unusual and complex merging cluster of galaxies at the intersection of optical filaments. We propose that energetics of multiple mergers and accretion flows has resulted in wide-spread shocks, acceleration of cosmic ray particles and amplification of weak magnetic fields. This results in Mpc-scale peripheral radio relics and halo like non-thermal emission observed near the merging center.

We investigate the interplay of cosmic ray (CR) propagation and advection in galaxy clusters. Propagation in form of CR diffusion and streaming tends to drive the CR radial profiles towards being flat, with equal CR number density everywhere. Advection of CR by the turbulent gas motions tends to produce centrally enhanced profiles. We assume that the CR streaming velocity is of the order of the sound velocity. This is motivated by plasma physical arguments. The CR streaming is then usually larger than typical advection velocities and becomes comparable or lower than this only for periods with trans- and super-sonic cluster turbulence. As a consequence a bimodality of the CR spatial distribution results. Strongly turbulent, merging clusters should have a more centrally concentrated CR energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected diffuse radio and gamma-ray emission of clusters, since more centrally concentrated CR will find higher target densities for hadronic CR proton interactions, higher plasma wave energy densities for CR electron and proton re-acceleration, and stronger magnetic fields. Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR electrons. Energy dependence of the CR propagation should lead to spectral steepening of dying radio halos. Furthermore, we show that the interplay of CR diffusion with advection implies first order CR re-acceleration in the pressure-stratified atmospheres of galaxy clusters. Finally, we argue that CR streaming could be important in turbulent cool cores of galaxy clusters since it heats preferentially the central gas with highest cooling rate.

We update our earlier calculations of gamma ray and radio observational constraints on annihilations of dark matter particles lighter than 10 GeV. We predict the synchrotron spectrum as well as the morphology of the radio emission associated with light decaying and annihilating dark matter candidates in both the Coma cluster and the Galactic Centre. Our new results basically confirm our previous findings: synchrotron emission in the very inner part of the Milky Way constrains or even excludes dark matter candidates if the magnetic field is larger than 50 micro Gauss. In fact, our results suggest that light annihilating candidates must have a S-wave suppressed pair annihilation cross section into electrons (or the branching ratio into electron positron must be small). If dark matter is decaying, it must have a life time that is larger than t = 3. 10^{25} s. Therefore, radio emission should always be considered when one proposes a "light" dark matter candidate.

Non-linear and non-Gaussian signal inference problems are difficult to tackle. Renormalization techniques permit us to construct good estimators for the posterior signal mean within information field theory (IFT), but the approximations and assumptions made are not very obvious. Here we introduce the simple concept of minimal Gibbs free energy to IFT, and show that previous renormalization results emerge naturally. They can be understood as being the Gaussian approximation to the full posterior probability, which has maximal cross information with it. We derive optimized estimators for three applications, to illustrate the usage of the framework: (i) reconstruction of a log-normal signal from Poissonian data with background counts and point spread function, as it is needed for gamma ray astronomy and for cosmography using photometric galaxy redshifts, (ii) inference of a Gaussian signal with unknown spectrum and (iii) inference of a Poissonian log-normal signal with unknown spectrum, the combination of (i) and (ii). Finally we explain how Gaussian knowledge states constructed by the minimal Gibbs free energy principle at different temperatures can be combined into a more accurate surrogate of the non-Gaussian posterior.

We present status and results of AstroGrid-D, a joint effort of astrophysicists and computer scientists to employ grid technology for scientific applications. AstroGrid-D provides access to a network of distributed machines with a set of commands as well as software interfaces. It allows simple use of computer and storage facilities and to schedule or monitor compute tasks and data management. It is based on the Globus Toolkit middleware (GT4). Chapter 1 describes the context which led to the demand for advanced software solutions in Astrophysics, and we state the goals of the project. We then present characteristic astrophysical applications that have been implemented on AstroGrid-D in chapter 2. We describe simulations of different complexity, compute-intensive calculations running on multiple sites, and advanced applications for specific scientific purposes, such as a connection to robotic telescopes. We can show from these examples how grid execution improves e.g. the scientific workflow. Chapter 3 explains the software tools and services that we adapted or newly developed. Section 3.1 is focused on the administrative aspects of the infrastructure, to manage users and monitor activity. Section 3.2 characterises the central components of our architecture: The AstroGrid-D information service to collect and store metadata, a file management system, the data management system, and a job manager for automatic submission of compute tasks. We summarise the successfully established infrastructure in chapter 4, concluding with our future plans to establish AstroGrid-D as a platform of modern e-Astronomy.

Reconstructing the matter density field from galaxy counts is a problem frequently addressed in current literature. Two main sources of error are shot noise from galaxy counts and insufficient knowledge of the correct galaxy position caused by peculiar velocities and redshift measurement uncertainty. Here we address the reconstruction problem of a Poissonian sampled log-normal density field with velocity distortions in a Bayesian way via a maximum a posteriory method. We test our algorithm on a 1D toy case and find significant improvement compared to simple data inversion. In particular, we address the following problems: photometric redshifts, mapping of extended sources in coded mask systems, real space reconstruction from redshift space galaxy distribution and combined analysis of data with different point spread functions.

We describe a technique for probing the statistical properties of cosmic magnetic fields based on radio polarimetry data. Second-order magnetic field statistics like the power spectrum cannot always distinguish between magnetic fields with essentially different spatial structure. Synchrotron polarimetry naturally allows certain 4th-order magnetic field statistics to be inferred from observational data, which lifts this degeneracy and can thereby help us gain a better picture of the structure of the cosmic fields and test theoretical scenarios describing magnetic turbulence. In this work we show that a 4th-order correlator of physical interest, the tension-force spectrum, can be recovered from the polarized synchrotron emission data. We develop an estimator for this quantity based on polarized-emission observations in the Faraday-rotation-free frequency regime. We consider two cases: a statistically isotropic field distribution, and a statistically isotropic field superimposed on a weak mean field. In both cases the tension force power spectrum is measurable; in the latter case, the magnetic power spectrum may also be obtainable. The method is exact in the idealized case of a homogeneous relativistic-electron distribution that has a power-law energy spectrum with a spectral index p=3, and assumes statistical isotropy of the turbulent field. We carry out tests of our method using synthetic data generated from numerically simulated magnetic fields. We show that the method is valid, that it is not prohibitively sensitive to the value of the electron spectral index, and that the observed tension-force spectrum allows one to distinguish between, e.g., a randomly tangled magnetic field (a default assumption in many studies) and a field organized in folded flux sheets or filaments.

We present a novel method to significantly speed up cosmological parameter sampling. The method relies on constructing an interpolation of the CMB-log-likelihood based on sparse grids, which is used as a shortcut for the likelihood-evaluation. We obtain excellent results over a large region in parameter space, comprising about 25 log-likelihoods around the peak, and we reproduce the one-dimensional projections of the likelihood almost perfectly. In speed and accuracy, our technique is competitive to existing approaches to accelerate parameter estimation based on polynomial interpolation or neural networks, while having some advantages over them. In our method, there is no danger of creating unphysical wiggles as it can be the case for polynomial fits of a high degree. Furthermore, we do not require a long training time as for neural networks, but the construction of the interpolation is determined by the time it takes to evaluate the likelihood at the sampling points, which can be parallelised to an arbitrary degree. Our approach is completely general, and it can adaptively exploit the properties of the underlying function. We can thus apply it to any problem where an accurate interpolation of a function is needed.

We search for an unusual alignment of the preferred axes of the quadrupole and octopole, the so-called axis of evil, in the CMB temperature and polarization data from WMAP. We use the part of the polarization map which is uncorrelated with the temperature map as a statistically independent probe of the axis of evil, which helps to assess whether the latter has a cosmological origin or if is a mere chance fluctuation in the temperature. Note, though, that for certain models creating a preferred axis in the temperature map, we would not expect to see the axis in the uncorrelated polarization map. We find that the axis of the quadrupole of the uncorrelated polarization map roughly aligns with the axis of evil within our measurement precision, whereas the axis of the octopole does not. However, with our measurement uncertainty, the probability of such a scenario to happen by chance in an isotropic universe is of the order of 50 per cent. We also find that the so-called cold spot present in the CMB temperature map is even colder in the part of the temperature map which is uncorrelated with the polarization, although there is still a large uncertainty in the latter. Therefore, our analysis of the axis of evil and a future analysis of the cold spot in the uncorrelated temperature data will strongly benefit from the polarization data expected from the Planck satellite.

In this work we present the first non-linear, non-Gaussian full Bayesian large scale structure analysis of the cosmic density field conducted so far. The density inference is based on the Sloan Digital Sky Survey data release 7, which covers the northern galactic cap. We employ a novel Bayesian sampling algorithm, which enables us to explore the extremely high dimensional non-Gaussian, non-linear log-normal Poissonian posterior of the three dimensional density field conditional on the data. These techniques are efficiently implemented in the HADES computer algorithm and permit the precise recovery of poorly sampled objects and non-linear density fields. The non-linear density inference is performed on a 750 Mpc cube with roughly 3 Mpc grid-resolution, while accounting for systematic effects, introduced by survey geometry and selection function of the SDSS, and the correct treatment of a Poissonian shot noise contribution. Our high resolution results represent remarkably well the cosmic web structure of the cosmic density field. Filaments, voids and clusters are clearly visible. Further, we also conduct a dynamical web classification, and estimated the web type posterior distribution conditional on the SDSS data.

We describe an exact, flexible, and computationally efficient algorithm for a joint estimation of the large-scale structure and its power-spectrum, building on a Gibbs sampling framework and present its implementation ARES (Algorithm for REconstruction and Sampling). ARES is designed to reconstruct the 3D power-spectrum together with the underlying dark matter density field in a Bayesian framework, under the reasonable assumption that the long wavelength Fourier components are Gaussian distributed. As a result ARES does not only provide a single estimate but samples from the joint posterior of the power-spectrum and density field conditional on a set of observations. This enables us to calculate any desired statistical summary, in particular we are able to provide joint uncertainty estimates. We apply our method to mock catalogs, with highly structured observational masks and selection functions, in order to demonstrate its ability to reconstruct the power-spectrum from real data sets, while fully accounting for any mask induced mode coupling.

We address the inverse problem of cosmic large-scale structure reconstruction from a Bayesian perspective. For a linear data model, a number of known and novel reconstruction schemes, which differ in terms of the underlying signal prior, data likelihood, and numerical inverse extra-regularization schemes are derived and classified. The Bayesian methodology presented in this paper tries to unify and extend the following methods: Wiener-filtering, Tikhonov regularization, Ridge regression, Maximum Entropy, and inverse regularization techniques. The inverse techniques considered here are the asymptotic regularization, the Jacobi, Steepest Descent, Newton-Raphson, Landweber-Fridman, and both linear and non-linear Krylov methods based on Fletcher-Reeves, Polak-Ribiere, and Hestenes-Stiefel Conjugate Gradients. The structures of the up-to-date highest-performing algorithms are presented, based on an operator scheme, which permits one to exploit the power of fast Fourier transforms. Using such an implementation of the generalized Wiener-filter in the novel ARGO-software package, the different numerical schemes are benchmarked with 1-, 2-, and 3-dimensional problems including structured white and Poissonian noise, data windowing and blurring effects. A novel numerical Krylov scheme is shown to be superior in terms of performance and fidelity. These fast inverse methods ultimately will enable the application of sampling techniques to explore complex joint posterior distributions. We outline how the space of the dark-matter density field, the peculiar velocity field, and the power spectrum can jointly be investigated by a Gibbs-sampling process. Such a method can be applied for the redshift distortions correction of the observed galaxies and for time-reversal reconstructions of the initial density field.

We present the largest Wiener reconstruction of the cosmic density field made to date. The reconstruction is based on the Sloan Digital Sky Survey data release 6 covering the northern Galactic cap. We use a novel supersampling algorithm to suppress aliasing effects and a Krylov-space inversion method to enable high performance with high resolution. These techniques are implemented in the ARGO computer code. We reconstruct the field over a 500 Mpc cube with Mpc grid-resolution while accounting both for the angular and radial selection functions of the SDSS, and the shot noise giving an effective resolution of the order of ~10 Mpc. In addition, we correct for the redshift distortions in the linear and nonlinear regimes in an approximate way. We show that the commonly used method of inverse weighting the galaxies by the corresponding selection function heads to excess noise in regions where the density of the observed galaxies is small. It is more accurate and conservative to adopt a Bayesian framework in which we model the galaxy selection/detection process to be Poisson-binomial. This results in heavier smoothing in regions of reduced sampling density. Our results show a complex cosmic web structure with huge void regions indicating that the recovered matter distribution is highly non-Gaussian. Filamentary structures are clearly visible on scales up to ~20 Mpc. We also calculate the statistical distribution of density after smoothing the reconstruction with Gaussian kernels of different radii r_S and find good agreement with a log-normal distribution for ~10 Mpc < r_S < ~30 Mpc.

Context: Hierarchal models of large scale structure (LSS) formation predict that galaxy clusters grow via gravitational infall and mergers of (smaller) mass concentrations, such as clusters and galaxy groups. Diffuse radio emission, in the form of radio halos and relics, is found in clusters undergoing a merger, indicating that shocks or turbulence associated with the merger are capable of accelerating electrons to highly relativistic energies. Here we report on radio observations of ZwCl 2341.1+0000, a complex merging structure of galaxies located at z=0.27, using Giant Metrewave Radio Telescope (GMRT) observations. Aims: The main aim of the observations is to study the nature of the diffuse radio emission in the galaxy cluster ZwCl 2341.1+0000. Methods: We have carried out GMRT 610, 241, and 157 MHz continuum observations of ZwCl 2341.1+0000. The radio observations are combined with X-ray and optical data of the cluster. Results: The GMRT observations show the presence of a double peripheral radio relic in the cluster ZwCl 2341.1+0000. The spectral index is -0.49 \pm 0.18 for the northern relic and -0.76 \pm 0.17 for the southern relic respectively. We have derived values of 0.48-0.93 microGauss for the equipartition magnetic field strength. The relics are probably associated with an outwards traveling merger shock waves.

We develop information field theory (IFT) as a means of Bayesian inference on spatially distributed signals, the information fields. A didactical approach is attempted. Starting from general considerations on the nature of measurements, signals, noise, and their relation to a physical reality, we derive the information Hamiltonian, the source field, propagator, and interaction terms. Free IFT reproduces the well known Wiener-filter theory. Interacting IFT can be diagrammatically expanded, for which we provide the Feynman rules in position-, Fourier-, and spherical harmonics space, and the Boltzmann-Shannon information measure. The theory should be applicable in many fields. However, here, two cosmological signal recovery problems are discussed in their IFT-formulation. 1) Reconstruction of the cosmic large-scale structure matter distribution from discrete galaxy counts in incomplete galaxy surveys within a simple model of galaxy formation. We show that a Gaussian signal, which should resemble the initial density perturbations of the Universe, observed with a strongly non-linear, incomplete and Poissonian-noise affected response, as the processes of structure and galaxy formation and observations provide, can be reconstructed thanks to the virtue of a response-renormalization flow equation. 2) We design a filter to detect local non-linearities in the cosmic microwave background, which are predicted from some Early-Universe inflationary scenarios, and expected due to measurement imperfections. This filter is the optimal Bayes' estimator up to linear order in the non-linearity parameter and can be used even to construct sky maps of non-linearities in the data.

Recent improvements in the capabilities of low frequency radio telescopes provide a unique opportunity to study thermal and non-thermal properties of the cosmic web. We argue that the diffuse, polarized emission from giant radio relics traces structure formation shock waves and illuminates the large-scale magnetic field. To show this, we model the population of shock-accelerated relativistic electrons in high-resolution cosmological simulations of galaxy clusters and calculate the resulting radio synchrotron emission. We find that individual shock waves correspond to localized peaks in the radio surface brightness map which enables us to measure Mach numbers for these shocks. We show that the luminosities and number counts of the relics strongly depend on the magnetic field properties, the cluster mass and dynamical state. By suitably combining different cluster data, including Faraday rotation measures, we are able to constrain some macroscopic parameters of the plasma at the structure formation shocks, such as models of turbulence. We also predict upper limits for the properties of the warm-hot intergalactic medium, such as its temperature and density. We predict that the current generation of radio telescopes (LOFAR, GMRT, MWA, LWA) have the potential to discover a substantially larger sample of radio relics, with multiple relics expected for each violently merging cluster. Future experiments (SKA) should enable us to further probe the macroscopic parameters of plasma physics in clusters.

We introduce a method to quantify the quality-of-fit between data and observables depending on the large scale Galactic magnetic field. We combine WMAP5 polarized synchrotron data and Rotation Measures of extragalactic sources in a joint analysis to obtain best fit parameters and confidence levels for GMF models common in the literature. None of the existing models provide a good fit in both the disk and halo regions, and in many instances best-fit parameters are quite different than the original values. We note that probing a very large parameter space is necessary to avoid false likelihood maxima. The thermal and relativistic electron densities are critical for determining the GMF from the observables but they are not well constrained. We show that some characteristics of the electron densities can already be constrained using our method and with future data it may be possible to carry out a self-consistent analysis in which models of the GMF and electron densities are simultaneously optimized.

Secondary anisotropies of the cosmic microwave background (CMB) can be detected by using the cross-correlation between the large-scale structure (LSS) and the CMB temperature fluctuations. In such studies, chance correlations of primordial CMB fluctuations with the LSS are the main source of uncertainty. We present a method for reducing this noise by exploiting information contained in the polarisation of CMB photons. The method is described in general terms and then applied to our recently proposed optimal method for measuring the integrated Sachs-Wolfe (ISW) effect. We obtain an expected signal-to-noise ratio of up to 8.5. This corresponds to an enhancement of the signal-to-noise by 23 per cent as compared to the standard method for ISW detection, and by 16 per cent w.r.t. our recently proposed method, both for the best-case scenario of having perfect (noiseless) CMB and LSS data.

We address the problem of discretizing continuous cosmological signals such as a galaxy distribution for further processing with Fast Fourier techniques. Discretizing, in particular representing continuous signals by discrete sets of sample points, introduces an enormous loss of information, which has to be understood in detail if one wants to make inference from the discretely sampled signal towards actual natural physical quantities. We therefore review the mathematics of discretizing signals and the application of Fast Fourier Transforms to demonstrate how the interpretation of the processed data can be affected by these procedures. It is also a well known fact that any practical sampling method introduces sampling artifacts and false information in the form of aliasing. These sampling artifacts, especially aliasing, make further processing of the sampled signal difficult. For this reason we introduce a fast and efficient supersampling method, frequently applied in 3D computer graphics, to cosmological applications such as matter power spectrum estimation. This method consists of two filtering steps which allow for a much better approximation of the ideal sampling procedure, while at the same time being computationally very efficient.Thus, it provides discretely sampled signals which are greately cleaned from aliasing contributions.

We analyse the local variance effect in the standard method for detecting the integrated Sachs-Wolfe effect (ISW) via cross-correlating the cosmic microwave background (CMB) with the large-scale structure (LSS). Local variance is defined as the systematic noise in the ISW detection that originates in the realisation of the matter distribution in the observed Universe. We show that the local variance contributes about 11 per cent to the total variance in the standard method, if a perfect and complete LSS survey up to z ~ 2 is assumed. Due to local variance, the estimated detection significance and cosmological parameter constraints in the standard method are biased. In this work, we present an optimal method of how to reduce the local variance effect in the ISW detection by working conditional on the LSS-data. The variance of the optimal method, and hence the signal-to-noise ratio, depends on the actual realisation of the matter distribution in the observed Universe. We show that for an ideal galaxy survey, the average signal-to-noise ratio is enhanced by about 7 per cent in the optimal method, as compared to the standard method. Furthermore, in the optimal method there is no need to estimate the covariance matrix by Monte Carlo simulations as in the standard method, which saves time and increases the accuracy. Finally, we derive the correct joint likelihood function for cosmological parameters given CMB- and LSS-data within the linear LSS formation regime, which includes a small coupling of the two datasets due to the ISW effect.

We present a publicly available code called Hammurabi for generating mock polarized observations of Galactic synchrotron emission for telescopes like LOFAR, SKA, Planck and WMAP, based on model inputs for the Galactic magnetic field (GMF), the cosmic-ray density distribution and the thermal electron density. We also present mock UHECR deflection measure (UDM) maps based on model inputs for the GMF. In future, when UHECR sources are identified, this will allow us to define UDM as a GMF probe in a similar way as polarized radio sources permit us to define rotation measures. To demonstrate the code's abilities mock observations are compared to real data as a means to constrain the input parameters of our simulations with a focus on large-scale magnetic field properties. As expected, attempts at trying to model the synchrotron, UHECR deflection and RM input parameters, show that any additional observational data set greatly increases the constraints on the models. The hammurabi code addresses this by allowing to perform simulations of several different data sets simultaneously, providing the means for a more reliable constraint of the magnetized inter-stellar-medium.

We investigate a numerical model for AGN feedback where for the first time a relativistic particle population in AGN-inflated bubbles is followed within a full cosmological context. In our high-resolution simulations of galaxy cluster formation, we assume that BH accretion is accompanied by energy feedback that occurs in two different modes, depending on the accretion rate itself. Unlike in previous work, we inject a non-thermal particle population of relativistic protons into the AGN bubbles, instead of adopting a purely thermal heating. We then follow the subsequent evolution of the cosmic ray (CR) plasma inside the bubbles, considering both its hydrodynamical interactions and dissipation processes relevant for the CR population. Due to the different buoyancy of relativistic plasma and the comparatively long CR dissipation timescale we find substantial changes in the evolution of clusters as a result of CR feedback. In particular, the non-thermal population can provide significant pressure support in central cluster regions at low thermal temperatures, providing a natural explanation for the decreasing temperature profiles found in cool core clusters. At the same time, the morphologies of the bubbles and of the induced X-ray cavities show a striking similarity to observational findings. AGN feedback with CRs also proves efficient in regulating cluster cooling flows so that the total baryon fraction in stars becomes limited to realistic values of the order of 10%. We find that the partial CR support of the intracluster gas also affects the expected signal of the thermal Sunyaev-Zel'dovich effect, with typical modifications of the integrated Compton-y parameter within the virial radius of the order of 10%. [Abridged]

The thermal plasma of galaxy clusters lost most of its information on how structure formation proceeded as a result of dissipative processes. In contrast, non-equilibrium distributions of cosmic rays (CR) preserve the information about their injection and transport processes and provide thus a unique window of current and past structure formation processes. This information can be unveiled by observations of non-thermal radiative processes, including radio synchrotron, hard X-ray, and gamma-ray emission. To explore this, we use high-resolution simulations of a sample of galaxy clusters spanning a mass range of about two orders of magnitudes, and follow self-consistent CR physics on top of the radiative hydrodynamics. We model CR electrons that are accelerated at cosmological structure formation shocks and those that are produced in hadronic interactions of CRs with ambient gas protons. We find that CR protons trace the time integrated non-equilibrium activities of clusters while shock-accelerated CR electrons probe current accretion and merging shock waves. The resulting inhomogeneous synchrotron emission matches the properties of observed radio relics. We propose a unified model for the generation of radio halos. Giant radio halos are dominated in the centre by secondary synchrotron emission with a transition to the synchrotron radiation emitted from shock-accelerated electrons in the cluster periphery. This model is able to explain the observed correlation of mergers with radio halos, the larger peripheral variation of the spectral index, and the large scatter in the scaling relation between cluster mass and synchrotron emission. Future low-frequency radio telescopes (LOFAR, GMRT, MWA, LWA) are expected to probe the accretion shocks of clusters. [abridged]

We report on initial results from a project to constrain the large-scale and turbulent magnetic fields of the Milky Way galaxy, which eventually will incorporate all of the relevant observational data. In this paper we fit popular large scale magnetic field models to WMAP3 polarization maps. We find that the polarization data can constrain certain model parameters but does not uniquely determine the best-fit parameters. We also find that the polarization data alone cannot distinguish between model symmetries, e.g., the existence of field reversals. We show how future UHECR data can break this degeneracy.

Observations of the Cosmic Microwave Background (CMB) have revealed an unexpected quadrupole-octopole alignment along a preferred axis pointing toward the Virgo cluster. We here investigate whether this feature can be explained in the framework of the concordance model by secondary anisotropies produced by the non-linear evolution of the gravitational potential, the so-called Rees-Sciama (RS) effect. We focus on the effect caused by the local superclusters, which we calculate using a constrained high-resolution hydrodynamical simulation, based on the IRAS 1.2-Jy all-sky galaxy redshift survey, which reproduces the main structures of our Universe out to a distance of 110 Mpc from our Galaxy. The resulting RS effect peaks at low multipoles and has a minimum/maximum amplitude of -6.6\mu K 1.9\mu K. Even though its quadrupole is well aligned with the one measured for the CMB, its amplitude is not sufficient to explain the observed magnitude of the quadrupole/octopole alignment. In addition, we analyze the WMAP-3 data with a linear matched filter in an attempt to determine an upper limit for the RS signal amplitude on large scales. We found that it is possible to infer a weak upper limit of 30\mu K for its maximum amplitude.

Cosmic rays (CRs) play a decisive role within our own Galaxy. They provide partial pressure support against gravity, they trace past energetic events such as supernovae, and they reveal the underlying structure of the baryonic matter distribution through their interactions. To study the impact of CRs on galaxy and cosmic structure formation and evolution, we develop an approximative framework for treating dynamical and radiative effects of CRs in cosmological simulations. Our guiding principle is to try to find a balance between capturing as many physical properties of CR populations as possible while at the same time requiring as little extra computational resources as possible. We approximate the CR spectrum of each fluid element by a single power-law, with spatially and temporally varying normalisation, low-energy cut-off, and spectral index. Principles of conservation of particle number, energy, and pressure are then used to derive evolution equations for the basic variables describing the CR spectrum, both due to adiabatic and non-adiabatic processes. The processes considered include compression and rarefaction, CR injection via shocks in supernova remnants, injection in structure formation shock waves, in-situ re-acceleration of CRs, CR spatial diffusion, CR energy losses due to Coulomb interactions, ionisation losses, Bremsstrahlung losses, and, finally, hadronic interactions with the background gas, including the associated gamma-ray and radio emission due to subsequent pion decay. We show that the formalism reproduces CR energy densities, pressure, and cooling rates with an accuracy of ~10% in steady state conditions where CR injection balances cooling. Our framework is therefore well suited to be included into numerical simulation schemes of galaxy and structure formation. (abridged)

We consider the role of the galactic kinetic Sunyaev Zeldovich (SZ) effect as a CMB foreground. While the galactic thermal Sunyaev Zeldovich effect has previously been studied and discarded as a potential CMB foreground, we find that the kinetic SZ effect is dominant in the galactic case. We analyse the detectability of the kinetic SZ effect by means of an optimally matched filter technique applied to a simulation of an ideal observation. We obtain no detection, getting a S/N ratio of 0.1, thereby demonstrating that the kinetic SZ effect can also safely be ignored as a CMB foreground. However we provide maps of the expected signal for inclusion in future high precision data processing. Furthermore, we rule out the significant contamination of the polarised CMB signal by second scattering of galactic kinetic Sunyaev-Zeldovich photons, since we show that the scattering of the CMB quadrupole photons by galactic electrons is a stronger effect than the Sunyaev Zeldovich second scattering, and has already been shown to produce no significant polarised contamination. We confirm the latter assessment also by means of an optimally matched filter.

One of the most outstanding problems in the gravitational collapse scenario of early structure formation is the cooling of primordial gas to allow for small mass objects to form. As the neutral primordial gas is a poor radiator at temperatures (T\le10^4\unit{K}), molecular hydrogen is needed for further cooling down to temperatures (T\sim100\unit{K}). The formation of molecular hydrogen is catalyzed by the presence of free electrons, which could be provided by the ionization due to an early population of cosmic rays. In order to investigate this possibility we developed a code to study the effects of ionizing cosmic rays on the thermal and chemical evolution of primordial gas. We found that cosmic rays can provide enough free electrons needed for the formation of molecular hydrogen, and therefore can increase the cooling ability of such primordial gas under following conditions: A dissociating photon flux with (F<10^{-18}\unit{erg cm^{-2} Hz^{-1} s^{-1}}), initial temperature of the gas (\sim10^{3}\unit{K}), total gas number densities (n\ge1\unit{cm^{-3}}), Cosmic ray sources with (\dot{\epsilon}_{CR}>10^{-33}\unit{erg cm^{-3} s^{-1}}).

Most cool core clusters of galaxies possess active galactic nuclei (AGN) in their centers. These AGN inflate buoyant bubbles containing non-thermal radio emitting particles. If such bubbles efficiently confine cosmic rays (CR) then this could explain ``radio ghosts'' seen far from cluster centers. We simulate the diffusion of cosmic rays from buoyant bubbles inflated by AGN. Our simulations include the effects of the anisotropic particle diffusion introduced by magnetic fields. Our models are consistent with the X-ray morphology of AGN bubbles, with disruption being suppressed by the magnetic draping effect. We conclude that for such magnetic field topologies, a substantial fraction of cosmic rays can be confined inside the bubbles on buoyant rise timescales even when the parallel diffusivity coefficient is very large. For isotropic diffusion at a comparable level, cosmic rays would leak out of the bubbles too rapidly to be consistent with radio observations. Thus, the long confinement times associated with the magnetic suppression of CR diffusion can explain the presence of radio ghosts. We show that the partial escape of cosmic rays is mostly confined to the wake of the rising bubbles, and speculate that this effect could: (1) account for the excitation of the H$\alpha$ filaments trailing behind the bubbles in the Perseus cluster, (2) inject entropy into the metal enriched material being lifted by the bubbles and, thus, help to displace it permanently from the cluster center and (3) produce observable $\gamma$-rays via the interaction of the diffusing cosmic rays with the thermal intracluster medium (ICM).

We performed high-resolution simulations of a sample of 14 galaxy clusters that span a mass range from 5 x 10^13 M_solar/h to 2 x 10^15 M_solar/h to study the effects of cosmic rays (CRs) on thermal cluster observables such as X-ray emission and the Sunyaev-Zel'dovich effect. We analyse the CR effects on the intra-cluster medium while simultaneously taking into account the cluster's dynamical state as well as the mass of the cluster. The modelling of the cosmic ray physics includes adiabatic CR transport processes, injection by supernovae and cosmological structure formation shocks, as well as CR thermalization by Coulomb interaction and catastrophic losses by hadronic interactions. While the relative pressure contained in CRs within the virial radius is of the order of 2 per cent in our non-radiative simulations, their contribution rises to 32 per cent in our simulations with dissipative gas physics including radiative cooling, star formation, and supernova feedback. Interestingly, in the radiative simulations the relative CR pressure reaches high values of the order of equipartition with the thermal gas in each cluster galaxy due to the fast thermal cooling of gas which diminishes the thermal pressure support relative to that in CRs. This also leads to a lower effective adiabatic index of the composite gas that increases the compressibility of the intra-cluster medium. This effect slightly increases the central density, thermal pressure and the gas fraction. While the X-ray luminosity in low mass cool core clusters is boosted by up to 40 per cent, the integrated Sunyaev-Zel'dovich effect appears to be remarkably robust and the total flux decrement only slightly reduced by typically 2 per cent. The resolved Sunyaev-Zel'dovich maps, however, show a larger variation with an increased central flux decrement. [abridged]

There is growing consensus that feedback from AGN is the main mechanism responsible for stopping cooling flows in clusters of galaxies. AGN are known to inflate buoyant bubbles that supply mechanical power to the intracluster gas (ICM). High Reynolds number hydrodynamical simulations show that such bubbles get entirely disrupted within 100 Myr, as they rise in cluster atmospheres, which is contrary to observations. This artificial mixing has consequences for models trying to quantify the amount of heating and star formation in cool core clusters of galaxies. It has been suggested that magnetic fields can stabilize bubbles against disruption. We perform MHD simulations of fossil bubbles in the presence of tangled magnetic fields using the high order PENCIL code. We focus on the physically-motivated case where thermal pressure dominates over magnetic pressure and consider randomly oriented fields with and without maximum helicity and a case where large scale external fields drape the bubble.We find that helicity has some stabilizing effect. However, unless the coherence length of magnetic fields exceeds the bubble size, the bubbles are quickly shredded. As observations of Hydra A suggest that lengthscale of magnetic fields may be smaller then typical bubble size, this may suggest that other mechanisms, such as viscosity, may be responsible for stabilizing the bubbles. However, since Faraday rotation observations of radio lobes do not constrain large scale ICM fields well if they are aligned with the bubble surface, the draping case may be a viable alternative solution to the problem. A generic feature found in our simulations is the formation of magnetic wakes where fields are ordered and amplified. We suggest that this effect could prevent evaporation by thermal conduction of cold Halpha filaments observed in the Perseus cluster.

Cosmological shock waves during structure formation not only play a decisive role for the thermalization of gas in virializing structures but also for the acceleration of relativistic cosmic rays (CRs) through diffusive shock acceleration. We discuss a novel numerical treatment of the physics of cosmic rays in combination with a formalism for identifying and measuring the shock strength on-the-fly during a smoothed particle hydrodynamics simulation. In our methodology, the non-thermal CR population is treated self-consistently in order to assess its dynamical impact on the thermal gas as well as other implications on cosmological observables. Using this formalism, we study the history of the thermalization process in high-resolution hydrodynamic simulations of the Lambda cold dark matter model. Collapsed cosmological structures are surrounded by shocks with high Mach numbers up to 1000, but they play only a minor role in the energy balance of thermalization. However, this finding has important consequences for our understanding of the spatial distribution of CRs in the large-scale structure. In high resolution simulations of galaxy clusters, we find a low contribution of the averaged CR pressure, due to the small acceleration efficiency of lower Mach numbers of flow shocks inside halos and the softer adiabatic index of CRs. However, within cool core regions, the CR pressure reaches equipartition with the thermal pressure leading there to a lower effective adiabatic index and thus to an enhanced compressibility of the central intracluster medium. This effect increases the central density and pressure of the cluster and thus the resulting X-ray emission and the central Sunyaev-Zel'dovich flux decrement. The integrated Sunyaev-Zel'dovich effect, however, is only slightly changed.

It is well known that cosmic rays (CRs) contribute significantly to the pressure of the interstellar medium in our own Galaxy, suggesting that they may play an important role in regulating star formation during the formation and evolution of galaxies. We will present a novel numerical treatment of the physics of CRs and its implementation in the parallel smoothed particle hydrodynamics (SPH) code GADGET-2. In our methodology, the non-thermal CR population is treated self-consistently in order to assess its dynamical impact on the thermal gas as well as other implications on cosmological observables. In simulations of galaxy formation, we find that CRs can significantly reduce the star formation efficiencies of small galaxies. This effect becomes progressively stronger towards low mass scales. In cosmological simulations of the formation of dwarf galaxies at high redshift, we find that the total mass-to-light ratio of small halos and the faint-end of the luminosity function are affected. In high resolution simulations of galaxy clusters, we find lower contributions of CR pressure, due to the smaller CR injection efficiencies at low Mach number flow shocks inside halos, and the softer adiabatic index of CRs, which disfavours them when a composite of thermal gas and CRs is adiabatically compressed. Within cool core regions, the CR pressure reaches equipartition with the thermal pressure leading to an enhanced compressibility of the central intra-cluster medium, an effect that increases the central density and pressure of the gas. While the X-ray luminosity in low mass cool core clusters is boosted, the integrated Sunyaev-Zel'dovich effect is only slightly changed. The resolved Sunyaev-Zel'dovich maps, however, show a larger variation with an increased central flux decrement.

The origin of intergalactic magnetic fields is still a mystery and several scenarios have been proposed so far: among them, primordial phase transitions, structure formation shocks and galactic outflows. In this work we investigate how efficiently galactic winds can provide an intense and widespread "seed" magnetisation. This may be used to explain the magnetic fields observed today in clusters of galaxies and in the intergalactic medium (IGM). We use semi-analytic simulations of magnetised galactic winds coupled to high resolution N-body simulations of structure formation to estimate lower and upper limits for the fraction of the IGM which can be magnetised up to a specified level. We find that galactic winds are able to seed a substantial fraction of the cosmic volume with magnetic fields. Most regions affected by winds have magnetic fields in the range -12 < Log B < -8 G, while higher seed fields can be obtained only rarely and in close proximity to wind-blowing galaxies. These seed fields are sufficiently intense for a moderately efficient turbulent dynamo to amplify them to the observed values. The volume filling factor of the magnetised regions strongly depends on the efficiency of winds to load mass from the ambient medium. However, winds never completely fill the whole Universe and pristine gas can be found in cosmic voids and regions unaffected by feedback even at z=0. This means that, in principle, there might be the possibility to probe the existence of primordial magnetic fields in such regions.

We argue that the recently reported Kolmogorov-like magnetic turbulence spectrum in the cool core of the Hydra A galaxy cluster can be understood by kinetic energy injection by active galaxies that drives a turbulent non-helical magnetic dynamo into its saturated state. Although dramatic differences exist between small-scale dynamo scenarios, their saturated state is expected to be similar, as we show for three scenarios: the flux rope dynamo, the fluctuation dynamo, and the explosive dynamo. Based on those scenarios, we develop an analytical model of the hydrodynamic and magnetic turbulence in cool cores. The model implies magnetic field strengths that fit well with Faraday rotation measurements and minimum energy estimates for the sample of cool core clusters having such data available. Predictions for magnetic fields in clusters for which the appropriate observational information is still missing, and for yet unobserved quantities like the hydrodynamical turbulence velocity and characteristic length-scale are provided. The underlying dynamo models suggest magnetic intermittency and possibly a large-scale hydrodynamic viscosity. We conclude that the success of the model to explain the field strength in cool core clusters indicates that in general cluster magnetic fields directly reflect hydrodynamical turbulence, also in clusters without cool cores.

We develop a formalism for the identification and accurate estimation of the strength of structure formation shocks during cosmological smoothed particle hydrodynamics simulations. Shocks not only play a decisive role for the thermalization of gas in virialising structures but also for the acceleration of relativistic cosmic rays (CRs) through diffusive shock acceleration. Our formalism is applicable both to ordinary non-relativistic thermal gas, and to plasmas composed of CRs and thermal gas. To this end, we derive an analytical solution to the one-dimensional Riemann shock tube problem for a composite plasma of CRs and thermal gas. We apply our methods to study the properties of structure formation shocks in high-resolution hydrodynamic simulations of the LCDM model. We find that most of the energy is dissipated in weak internal shocks with Mach numbers M~2 which are predominantly central flow shocks or merger shock waves traversing halo centres. Collapsed cosmological structures are surrounded by external shocks with much higher Mach numbers up to M~1000, but they play only a minor role in the energy balance of thermalization. We show that after the epoch of cosmic reionisation the Mach number distribution is significantly modified by an efficient suppression of strong external shock waves due to the associated increase of the sound speed of the diffuse gas. Invoking a model for CR acceleration in shock waves, we find that the average strength of shock waves responsible for CR energy injection is higher than that for shocks that dominate the thermalization of the gas. When combined with radiative dissipation and star formation, our formalism can also be used to study CR injection by supernova shocks, or to construct models for shock-induced star formation in the interstellar medium. (abridged)

It is well known that cosmic rays (CRs) contribute significantly to the pressure of the interstellar medium in our own Galaxy, suggesting that they may play an important role in regulating star formation during the formation and evolution of galaxies. We here discuss a novel numerical treatment of the physics of CRs and its implementation in the parallel smoothed particle hydrodynamics code GADGET-2. In our methodology, the non-thermal CR population of each gaseous fluid element is approximated by a simple power law spectrum in particle momentum, characterized by an amplitude, a cut-off, and a fixed slope. Adiabatic compression, and a number of physical source and sink terms are modelled which modify the CR pressure of each particle. The most important sources considered are injection by supernovae and diffusive shock acceleration, while the primary sinks are thermalization by Coulomb interactions, and catastrophic losses by hadronic interactions. We also include diffusion of CRs. Our scheme allows us to carry out the first cosmological structure formation simulations that self-consistently account for CR physics. In simulations of isolated galaxies, we find that CRs can significantly reduce the star formation efficiencies of small galaxies, with virial velocities below \~80 km/s, an effect that becomes progressively stronger towards low mass scales. In cosmological simulations at high redshift, the total mass-to-light ratio of small halos and the faint-end of the luminosity function are strongly affected. When CR acceleration in shocks is followed as well, up to ~40% of the energy dissipated at structure formation shocks can appear as CR pressure at z~3-6, but this fraction drops to ~10% at low redshifts when the shock distribution becomes increasingly dominated by lower Mach numbers. (abridged)

We present deep VLA observations of the merging galaxy cluster Abell 2256. This cluster is known to possess diffuse steep spectrum radio relic emission in the peripheral regions. Our new observations provide the first detailed image of the central diffuse radio halo emission in this cluster. The radio halo extends over more than 800 kpc in the cluster core, while the relic emission covers a region of ~1125 x 520 kpc. A spectral index map of the radio relic shows a spectral steepening from the northwest toward the southeast edge of the emission, with an average spectral index between 1369 MHz and 1703 MHz of -1.2 across the relic. Polarization maps reveal high fractional polarization of up to 45% in the relic region with an average polarization of 20% across the relic region. The observed Faraday rotation measure is consistent with the Galactic estimate and the dispersion in the rotation measure is small, suggesting that there is very little contribution to the rotation measure of the relic from the intracluster medium. We use these Faraday properties of the relic to argue that it is located on the front side of the cluster.

We discuss how to use the Rees-Sciama (RS) effect associated with merging clusters of galaxies to measure their kinematic properties. In a previous work (Rubino-Martin et al. 2004), the morphology and symmetries of the effect were examined by means of a simplified model. Here, we use realistic N-body simulations to better describe the effect, and to confirm that the signal has a characteristic quadrupole structure. From the amplitude of the signal obtained, we conclude that it is necessary to combine several cluster mergers in order to achieve a detection. Using the extended Press-Schechter formalism, we characterized the expected distribution of the parameters describing the mergers, and we used these results to generate realistic mock catalogues of cluster mergers. To optimize the extraction of the RS signal, we developed an extension of the spatial filtering method described in Haehnelt & Tegmark (1996). This extended filter has a general definition, so it can be applied in many other fields, such as gravitational lensing of the CMB or lensing of background galaxies. It has been applied to our mock catalogues, and we show that with the announced sensitivities of future experiments like the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT) or the Atacama Large Millimeter Array (ALMA), a detection of the signal will be possible if we consider of the order of 1,000 cluster mergers.

The Planck mission will permit measurements of the polarization of the cosmic microwave background and of polarized foregrounds such as our own Galaxy with an unprecedented combination of accuracy and completeness. This will provide information on cosmological and galactic magnetic fields. The latter can be studied in detail via nearly Faraday-rotation free synchrotron and polarized dust emission. Methods are discussed to extract physically relevant information on the magnetic turbulence from Planck data and other measurements.

We describe an assembly of numerical tools to model the output data of the Planck satellite. These start with the generation of a CMB sky in a chosen cosmology, add in various foreground sources, convolve the sky signal with arbitrary, even non-symmetric and polarised beam patterns, derive the time ordered data streams measured by the detectors depending on the chosen satellite-scanning strategy, and include noise signals for the individual detectors and electronic systems. The simulation products are needed to develop, verify, optimise, and characterise the accuracy and performance of all data processing and scientific analysis steps of the Planck mission, including data handling, data integrity checking, calibration, map making, physical component separation, and power spectrum estimation. In addition, the simulations allow detailed studies of the impact of many stochastic and systematic effects on the scientific results. The efficient implementation of the simulation allows the build-up of extended statistics of signal variances and co-variances. Although being developed specifically for the Planck mission, it is expected that the employed framework as well as most of the simulation tools will be of use for other experiments and CMB-related science in general.

We proposed a new method -- "Pacman" -- to calculate Faraday rotation measure (RM) maps from multi-frequency polarisation angle data (Dolag et al.) in order to avoid the so-called n-pi ambiguity. Here, we apply our "Pacman" algorithm to two polarisation data sets of extended radio sources in the Abell 2255 and the Hydra A cluster, and compare the RM maps obtained using "Pacman" to RM maps obtained employing already existing methods. Thereby, we provide a new high quality RM map of the Hydra north lobe which is in a good agreement with the existing one but find significant differences in the case of the south lobe of Hydra A. We demonstrate the reliability and the robustness of "Pacman". In order to study the influence of map making artefacts, which are imprinted by wrong solutions to the n-pi ambiguities, and of the error treatment of the data, we calculated and compared magnetic field power spectra from various RM maps. The power spectra were derived using the method recently proposed by Ensslin & Vogt (2003). We demonstrate the sensitivity of statistical analysis to artefacts and noise in the RM maps and thus, we demonstrate the importance of an unambiguous determination of RM maps and an understanding of the nature of the noise in the data. We introduce and perform statistical tests to estimate the quality of the derived RM maps, which demonstrate the quality improvements due to "Pacman".

We propose a new method to calculate Faraday rotation measure maps from multi-frequency polarisation angle data. In order to solve the so called n-pi ambiguity problem which arises from the observationally ambiguity of the polarisation angle which is only determined up to additions of +-n*pi, where n is an integer, we suggest using a global scheme. Instead of solving the n-pi ambiguity for each data point independently, our algorithm, which we chose to call "Pacman" (Polarisation Angle Correcting rotation Measure ANalysis), solves the n-pi ambiguity for a high signal-to-noise region ``democratically'' and uses this information to assist computations in adjacent low signal-to-noise areas.

A brief overview about our knowledge on galaxy cluster magnetic fields is provided. Emphasize is given to the mutual dependence of our knowledge on relativistic particles in galaxy clusters and the magnetic field strength. Furthermore, we describe efforts to measure magnetic field strengths, characteristic length-scales, and power-spectra with reliable accuracy. An interpretation of these results in terms of non-helical dynamo theory is given. If this interpretation turns out to be correct, the understanding of cluster magnetic fields is directly connected to our understanding of intra-cluster turbulence.

A critical discussion of our knowledge about extragalactic cosmic rays and magnetic fields is attempted. What do we know for sure? What are our prejudices? How do we confront our models with the observations? How can we assess the uncertainties in our modeling and in our observations? Unfortunately, perfect answers to these questions can not be given. Instead, I describe efforts I am involved in to gain reliable information about relativistic particles and magnetic fields in extragalactic space.

We present a Bayesian maximum likelihood analysis of Faraday rotation measure (RM) maps of extended radio sources to determine magnetic field power spectra in clusters of galaxies. Using this approach, it is possible to determine the uncertainties in the measurements. We apply this approach to the RM map of Hydra A and derive the power spectrum of the cluster magnetic field. For Hydra A, we measure a spectral index of -5/3 over at least one order of magnitude implying Kolmogorov type turbulence. We find a dominant scale of about 3 kpc on which the magnetic power is concentrated, since the magnetic autocorrelation length is lambda_B = 3 +/- 0.5 kpc. Furthermore, we investigate the influences of the assumption about the sampling volume (described by a window function) on the magnetic power spectrum. The central magnetic field strength was determined to be about 7 +/- 2 muG for the most likely geometries.

There have been many speculations about the presence of cosmic ray protons (CRps) in galaxy clusters over the past two decades. However, no direct evidence such as the characteristic gamma-ray signature of decaying pions has been found so far. These pions would be a direct tracer of hadronic CRp interactions with the ambient thermal gas also yielding observable synchrotron and inverse Compton emission by additionally produced secondary electrons. The obvious question concerns the type of galaxy clusters most likely to yield a signal: Particularly suited sites should be cluster cooling cores due to their high gas and magnetic energy densities. We studied a nearby sample of clusters evincing cooling cores in order to place stringent limits on the cluster CRp population by using non-detections of EGRET. In this context, we examined the possibility of a hadronic origin of Coma-sized radio halos as well as radio mini-halos. Especially for mini-halos, strong clues are provided by the very plausible small amount of required CRp energy density and a matching radio profile. Introducing the hadronic minimum energy criterion, we show that the energetically favored CRp energy density is constrained to 2% +/- 1% of the thermal energy density in Perseus. We also studied the CRp population within the cooling core region of Virgo using the TeV gamma-ray detection of M 87 by HEGRA. Both the expected radial gamma-ray profile and the required amount of CRp support this hadronic scenario.

The Chandra X-ray Observatory is finding a large number of cavities in the X-ray emitting intra-cluster medium which often coincide with the lobes of the central radio galaxy. We propose high-resolution Sunyaev-Zel'dovich (SZ) observations in order to infer the still unknown dynamically dominating component of the radio plasma bubbles. This work calculates the thermal and relativistic SZ emission of different compositions of these plasma bubbles while simultaneously allowing for the cluster's kinetic SZ effect. As examples, we present simulations of an Atacama Large Millimeter Array (ALMA) observation and of a Green Bank Telescope (GBT) observation of the cores of the Perseus cluster and Abell 2052. We predict a 5 sigma detection of the southern radio bubble of Perseus in a few hours with the GBT and ALMA while assuming a relativistic electron population within the bubble. In Abell 2052, a similar detection would require a few ten hours with either telescope, the longer exposures mainly being the result of the higher redshift and the lower central temperature of this cluster. Future high-sensitivity multi-frequency SZ observations will be able to infer the energy spectrum of the dynamically dominating electron population in order to measure its temperature or spectral characteristics. This knowledge can yield indirect indications for an underlying radio jet model.

An important aspect of the radio emission from galaxy clusters is represented by the diffuse radio sources associated with the intracluster medium: radio halos, relics and mini-halos. The radio halos and relics are indicators of cluster mergers, whereas mini-halos are detected at the center of cooling core clusters. SKA will dramatically improve the knowledge of these sources, thanks to the detection of new objects, and to detailed studies of their spectra and polarized emission. SKA will also provide the opportunity to investigate the presence of halos produced by radiation scattered by a powerful radio galaxy at the cluster centers.

We discuss the influence of large scale cosmic magnetic fields on the propagation of hadronic cosmic rays above 10^19 eV based on large scale structure simulations. Our simulations suggest that rather substantial deflection up to several tens of degrees at 10^20 eV are possible for nucleon primaries. Further, spectra and composition of cosmic rays from individual sources can depend on magnetic fields surrounding these sources in intrinsically unpredictable ways. This is true even if deflection from such individual sources is small. We conclude that the influence of large scale cosmic magnetic fields on ultra-high energy cosmic ray propagation is currently hard to quantify. We discuss possible reasons for discrepant results of simulations by Dolag et al. which predict deflections of at most a few degrees for nucleons. We finally point out that even in these latter simulations a possible heavy component would in general suffer substantial deflection.

We study signatures of a structured universe in the multi-pole moments, auto-correlation function, and cluster statistics of ultra-high energy cosmic rays above 10^19 eV. We compare scenarios where the sources are distributed homogeneously or according to the baryon density distribution obtained from a cosmological large scale structure simulation. The influence of extragalactic magnetic fields is studied by comparing the case of negligible fields with fields expected to be produced along large scale shocks with a maximal strength consistent with observations. We confirm that strongly magnetized observers would predict considerable anisotropy on large scales, which is already in conflict with current data. In the best fit scenario only the sources are strongly magnetized, although deflection can still be considerable, of order 20 degrees up to 10^20 eV, and a pronounced GZK cutoff is predicted. We then discuss signatures for future large scale full-sky detectors such as the Pierre Auger and EUSO projects. Auto-correlations are sensitive to the source density only if magnetic fields do not significantly affect propagation. In contrast, for a weakly magnetized observer, degree scale auto-correlations below a certain level indicate magnetized discrete sources. It may be difficult even for next generation experiments to distinguish between structured and unstructured source distributions.

We wish to constrain the cosmic-ray proton (CRp) population in galaxy clusters. By hadronic interactions with the thermal gas of the intra-cluster medium (ICM), the CRp produce gamma-rays for which we develop an analytic formalism to deduce their spectral distribution. Assuming the CRp-to-thermal energy density ratio X_CRp and the CRp spectral index to be spatially constant, we derive an analytic relation between the gamma-ray and bolometric X-ray fluxes, F_gamma and F_X. Based on our relation, we compile a sample of suitable clusters which are promising candidates for future detection of gamma-rays resulting from hadronic CRp interactions. Comparing to EGRET upper limits, we constrain the CRp population in the cooling flow clusters Perseus and Virgo to X_CRp < 20%. Assuming a plausible value for the CRp diffusion coefficient kappa, we find the central CRp injection luminosity of M 87 to be limited to 10^43 erg s^-1 kappa/(10^29 cm^2 s^-1). The synchrotron emission from secondary electrons generated in CRp hadronic interactions allows even tighter limits to be placed on the CRp population using radio observations. We obtain excellent agreement between the observed and theoretical radio brightness profiles for Perseus, but not for Coma without a radially increasing CRp-to-thermal energy density profile. Since the CRp and magnetic energy densities necessary to reproduce the observed radio flux are very plausible, we propose synchrotron emission from secondary electrons as an attractive explanation of the radio mini-halos found in cooling flow clusters. This model can be tested with future sensitive gamma-ray observations of the accompanying pi0-decays. We identify Perseus (A 426), Virgo, Ophiuchus, and Coma (A 1656) as the most promising candidate clusters for such observations.

We present X-ray spectroscopic and morphological analyses using Chandra ACIS and ROSAT observations of the giant elliptical galaxy NGC 4472 in the Virgo cluster. We discuss previously unobserved X-ray structures within the extended galactic corona. In the inner 2' of the galaxy, we find X-ray holes or cavities with radii of ~2 kpc, corresponding to the position of radio lobes. These holes were produced during a period of nuclear activity that began 1.2 x 10^7 years ago and may be ongoing. We also find an asymmetrical edge in the galaxy X-ray emission 3' (14 kpc) northeast of the core and an ~8' tail (36 kpc) extending southwest of the galaxy. These two features probably result from the interaction of NGC 4472 gas with the Virgo gas, which produces compression in the direction of NGC 4472's infall and an extended tail from ram pressure stripping. Assuming the tail is in pressure equilibrium with the surrounding gas, we compute its angle to our line of sight and estimate that its true extent exceeds 100 kpc. Finally, in addition to emission from the nucleus (first detected by Soldatenkov, Vikhlinin & Pavlinsky), we detect two small extended sources within 10'' of the nucleus of the galaxy, both of which have luminosities of ~7 x 10^38 erg/s.

We wish to estimate magnetic field strengths of radio emitting galaxy clusters by minimising the non-thermal energy density contained in cosmic ray electrons (CRe), protons (CRp), and magnetic fields. The classical minimum energy estimate can be constructed independently of the origin of the radio synchrotron emitting CRe yielding thus an absolute minimum of the non-thermal energy density. Provided the observed synchrotron emission is generated by a CRe population originating from hadronic interactions of CRp with the ambient thermal gas of the intra-cluster medium, the parameter space of the classical scenario can be tightened by means of the hadronic minimum energy criterion. For both approaches, we derive the theoretically expected tolerance regions for the inferred minimum energy densities. Application to the radio halo of the Coma cluster and the radio mini-halo of the Perseus cluster yields equipartition between cosmic rays and magnetic fields within the expected tolerance regions. In the hadronic scenario, the inferred central magnetic field strength ranges from 2.4 muG (Coma) to 8.8 muG (Perseus), while the optimal CRp energy density is constrained to 2% +/- 1% of the thermal energy density (Perseus). We discuss the possibility of a hadronic origin of the Coma radio halo while current observations favour such a scenario for the Perseus radio mini-halo. Combining future expected detections of radio synchrotron, hard X-ray inverse Compton, and hadronically induced gamma-ray emission should allow an estimate of volume averaged cluster magnetic fields and provide information about their dynamical state.

The Rees-Sciama effect produced in mergers of galaxy clusters is discussed, and an analytical approximation to compute this effect from numerical simulations is given. Using this approximation and a novel toy model describing the physics of the merger, we characterize the spatial properties and symmetries of the Rees-Sciama signal. Based on these properties, we propose a method to extract the physical parameters of the merger, which relies on the computation of the quadrupole moment of the observed brightness distribution on the sky. The relationships between the quadrupole coefficients and the physical parameters of the merger (physical separation, projection angle on the sky and angular momentum) are discussed. Finally, we propose a method to co-add coherently the RS signals from a sample of cluster mergers, in order to achieve an statistical detection of the effect for those cases where individual signals are masked by the kinetic SZ effect, the primordial CMB components, and by observational noise.

Populations of relativistic electrons in clusters of galaxies should imprint a characteristic spectral distortion signature on the Cosmic Microwave Background photons - the relativistic Sunyaev-Zeldovich (SZ) effect. We investigate how sensitive PLANCK and successor experiments will be in detecting or constraining the total optical depth tau of relativistic electrons. We find an expected sensitivity of PLANCK of sigma_tau ~ 2 10^-5 for an individual cluster, which is too insensitive to detect even optimistic scenarios, which predict tau ~ 10^-6 to 10^-7. However, by stacking the PLANCK-signal of all SZ detectable clusters the relativistic SZ effect can be statistically probed to a sensitivity level of sigma_tau ~ 10^-7. PLANCK successor experiments will be able to probe even conservative scenarios of relativistic electron populations in clusters of galaxies.

At the request of the conference attendees, we have compiled a classification of extended radio sources in clusters. These range from scales of tens of parsecs to over a megaparsec in scale, and include both sources associated with AGN and sources thought to derive from the electron population in the ionized ICM. We pay special attention to distinguishing between the types of AGN in the cores of cooling flow clusters and between the multiple classes of objects referred to over the years as ``radio relics.'' We suggest new names based on physical arguments for some of these classes of objects where their commonly used names are inappropriate or confusing.

Cooling flows are regions where the importance of non-thermal intra-cluster medium components such as magnetic fields and cosmic rays may be strongest within a galaxy cluster. They are also regions where such components are best detectable due to the high gas density which influences Faraday rotation measurements of magnetic fields and secondary particle production in hadronic interactions of cosmic ray nuclei with the ambient thermal gas. New estimates of magnetic fields in cooling flow and non-cooling flow clusters are presented, which are based on a newly developed Fourier analysis of extended Faraday rotation maps. We further present new constraints on the cluster cosmic ray proton population using radio and gamma-ray observations measurements of cluster cooling flows, which are especially suited for this purpose due to their high gas and magnetic energy densities. We argue that radio synchrotron emission of cosmic ray electrons generated hadronically by cosmic ray protons is a very plausible explanation for the radio mini-halos observed in some cooling flows.

We investigate the impact of a structured universe in the multi-pole moments, auto-correlation function, and cluster statistics of cosmic rays above 10^19 eV. We compare structured and uniform source distributions with and without magnetic fields obtained from a cosmological simulation. We find that current data marginally favor structured source distributions and magnetic fields reaching a few micro Gauss in galaxy clusters but below 0.1 micro Gauss in our local extragalactic neighborhood. A pronounced GZK cutoff is also predicted in this scenario. Future experiments will make the degree scale auto-correlation function a sensitive probe of micro Gauss fields surrounding the sources.

In this paper we compute the luminosity function of radio relics. In our calculation we include only those relics that are produced by the compression of former radio cocoons. This compression is provided by shocks that are generated in the process of structure formation. Starting from an analytical model for the luminosity evolution of ageing radio cocoons, the luminosity function of radio galaxies and the statistics of shocks as inferred from cosmological simulations, we are able to make the first estimates of the brightness distribution of radio relics. The computed luminosity function is consistent with current observations and predicts that more than $10^3$ radio relics should be discovered with the upcoming generation of low-frequency radio telescopes. Moreover, we predict that radio relics are predominantly found in low-pressure regions outside the cores of clusters.

We apply a novel technique of Faraday Rotation measure (RM) map analysis to three galaxy clusters, Abell 400, Abell 2634 and Hydra A, in order to estimate cluster magnetic field strengths, length scales and power spectra. This analysis is based on the assumption that the magnetic fields are statistically isotropically distributed across the Faraday screen. We investigate the difficulties involved in the application of the analysis to observational data. We derive magnetic power spectra for the three clusters and discuss influences on their shapes. We show that magnetic fluctuations are probed on length scales ranging over at least one order of magnitude. Using this range for the determination of central cluster magnetic field strength yields 3 muG in Abell 2634, 6 muG in Abell 400 and 12 muG in Hydra A. The magnetic field autocorrelation length was determined to be 4.9 kpc for Abell 2634, 3.6 kpc for Abell 400 and 0.9 kpc for Hydra A. We show that the RM autocorrelation length is larger than the magnetic field autocorrelation length. We investigate in a response analysis if it is possible to determine spectral slopes of the power spectra. We find that integrated numbers can be determined from this analysis but differential parameters such as spectral slopes have to be treated differently. Our response analysis results in spectral slopes of the power spectra of spectral indices alpha = 1.6 to 2.0 suggesting that Kolmogorov spectra are possible but flatter spectral slopes than alpha = 1.3 can be excluded.

We examine the cosmic ray proton (CRp) population within the giant elliptical galaxy M 87 using the TeV gamma-ray detection of the HEGRA collaboration. In our scenario, the gamma-rays are produced by decaying pions which result from hadronic CRp interactions with thermal gas of the interstellar medium of M 87. By comparing the gamma-ray emission to upper limits of EGRET, we constrain the spectral index of the CRp population to alpha(GeV,TeV) < 2.275 within our scenario. Both the expected radial gamma-ray profile and the required amount of CRp support this hadronic scenario. The accompanying radio mini-halo of hadronically originating cosmic ray electrons is outshone by the synchrotron emission of the relativistic jet of M 87 by one order of magnitude. According to our predictions, the future GLAST mission should allow us to test this hadronic scenario.

We investigate the origin of the high Faraday rotation measures (RMs) found for polarized radio galaxies in clusters. The two most likely origins are, magnetic fields local to the source, or magnetic fields located in the foreground intra-cluster medium (ICM). The latter is identified as the null hypothesis. Rudnick & Blundell (2003) have recently suggested that the presence of magnetic fields local to the source may be revealed in correlations of the position angles (PAs) of the source intrinsic linear polarization and the RMs. We investigate the claim of Rudnick & Blundell to have found a relationship between the intrinsic PA0 of the radio source PKS 1246-410 and its RM, by testing the clustering strength of the PA0-RM scatter plot. We show that the claimed relationship is an artifact of an improperly performed null-experiment. We describe a gradient alignment statistic aimed at finding correlations between PA0 and RM maps. This statistic does not require any null-experiment since it gives a unique (zero) result in the case of uncorrelated maps. We apply it to a number of extended radio sources in galaxy clusters (PKS 1246-410, Cygnus A, Hydra A, and 3C465). In no case is a significant large-scale alignment of PA0 and RM maps detected. We find significant small-scale co-alignment in all cases, but we are able to fully identify this with map making artifacts through a suitable statistical test. We conclude that there is presently no existing evidence for Faraday rotation local to radio lobes. Given the existing independent pieces of evidence, we favor the null hypothesis that the observed Faraday screens are produced by intracluster magnetic fields.

We simulate propagation of cosmic ray nucleons above 10^{19} eV in scenarios where both the source distribution and magnetic fields within about 50 Mpc from us are obtained from an unconstrained large scale structure simulation. We find that consistency of predicted sky distributions with current data above 4 x 10^{19} eV requires magnetic fields of ~0.1 microGauss in our immediate environment, and a nearby source density of ~10^{-4}-10^{-3} Mpc^{-3}. Radio galaxies could provide the required sources, but only if both high and low-luminosity radio galaxies are very efficient cosmic ray accelerators. Moreover, at ~10^{19} eV an additional isotropic flux component, presumably of cosmological origin, should dominate over the local flux component by about a factor three in order to explain the observed isotropy. This argues against the scenario in which local astrophysical sources of cosmic rays above ~10^{19} eV reside in strongly magnetized (B~0.1 microGauss) and structured intergalactic medium. Finally we discuss how future large scale full-sky detectors such as the Pierre Auger project will allow to put much more stringent constraints on source and magnetic field distributions.

We estimate the gamma ray fluxes from the residual annihilations of Dark Matter particles having a mass mdm \in [MeV, O(GeV)] and compare them to observations. We find that particles lighter than O(100 MeV) are excluded unless their cross section is S-wave suppressed.

Many compact radio sources like quasars, blazars, radio galaxies, and micro-quasars emit circular polarisation (CP) with surprising temporal persistent handedness. We propose that the CP is caused by Faraday conversion of linear polarisation synchrotron light which propagates along a line-of-sight through helical magnetic fields. Jet outflows from radio galaxies should have the required magnetic helicity in the emission region due to the magnetic torque of the accretion disc. Also advection dominated accretion flow (ADAF) should contain magnetic fields with the same helicity. However, a jet region seems to be the more plausible origin of CP. The proposed scenario requires Faraday rotation (FR) to be insignificant in the emission region. The proposed mechanism works in electron-positron e+/e- as well as electron-proton e/p plasma. In the latter case, the emission region should consist of individual flux tubes with independent polarities in order to suppress too strong FR - as it was already proposed for FR based CP generation models. The predominant CP is expected to mostly counter-rotate (rotation is measured here in sky-projection) with respect to the central engine in all cases (jet or ADAF, e+/e- or e/p plasma) and therefore allows to measure the sense of rotation of quasar engines. The engine of SgrA* is expected - in this scenario - to rotate clockwise and therefore counter-Galactic, as do the young hot stars in its vicinity, which are thought to feed SgrA* by their winds. Generally, sources with Stokes-V<0 (V>0) are expected to rotate clockwise

The autocorrelation function and similarly the Fourier-power spectrum of a rotation measure (RM) map of an extended background radio source can be used to measure components of the magnetic autocorrelation and power-spectrum tensor within a foreground Faraday screen. It is possible to reconstruct the full non-helical part of this tensor in the case of an isotropic magnetic field distribution statistics. The helical part is only accessible with additional information; e.g. the knowledge that the fields are force-free. The magnetic field strength, energy spectrum and autocorrelation length l_B can be obtained from the non-helical part alone. We demonstrate that l_B can differ substantially from l_RM, the observationally easily accessible autocorrelation length of an RM map. In typical astrophysical situation l_RM > l_B. Any RM study, which does not take this distinction into account, likely underestimates the magnetic field strength. For power-law magnetic power spectra, and for patchy magnetic field configurations the central RM autocorrelation function is shown to have characteristic asymptotic shapes. Ways to constrain the volume filling factor of a patchy field distribution are discussed. We discuss strategies to analyse observational data, taking into account - with the help of a window function - the limited extent of the polarised radio source, the spatial distribution of the electron density and average magnetic energy density in the screen, and allowing for noise reducing data weighting. We briefly discuss the effects of possible observational artefacts, and strategies to avoid them.

Many radio sources like quasars, blazars, radio galaxies, and micro-quasars exhibit circular polarisation (CP) with surprising temporal persistent handedness. As a possible explanation we propose that the CP is due to Faraday conversion (FC) of linear polarisation (LP) synchrotron light which propagates along a line-of-sight (LOS) through twisted magnetic fields. The rotational nature of accretion flows onto black holes naturally generates the required magnetic twist in the emission region, independent of whether it is a jet or an ADAF. The expected twist in both types of flows is of the order of what is required for optimal CP generation. This scenario requires that Faraday rotation (FR) is insignificant in the emission region. Although this is an assumption, it relaxes constraints on the plasma parameters, that were given in scenarios which rely on FR, since there the strength of FR can not be too far from the optimum. The proposed mechanism works in electron-positron (e+-) as well as electron-proton (e/p) plasma. In the latter case, the emission region should consist of individual flux tubes with independent polarities in order to suppress too strong FR. The predominant CP is expected to mostly counter-rotate (rotation is measured here in sky-projection) with respect to the central engine in all cases (jet or ADAF, e+- or e/p plasma). If the proposed mechanism is indeed operating, it will allow to measure the sense of rotation of quasar engines. The engine of SgrA* is then expected to rotate clockwise and therefore counter-Galactic, as do the young hot stars in its vicinity, which are thought to feed SgrA* by their winds. Similarly, we expect the microquasars SS 443 and GRS 1915+105 to rotate clockwise. Generally, sources with Stokes-V<0 (V>0) are expected to rotate clockwise (counter-clockwise) in this scenario.

(Abridged) A model for the escape of CR particles from radio galaxy cocoons is presented here. It is assumed that the radio cocoon is poorly magnetically connected to the environment. An extreme case of this kind is an insulating boundary layer of magnetic fields, which can efficiently suppress particle escape. More likely, magnetic field lines are less organised and allow the transport of CR particles from the source interior to the surface region. For such a scenario two transport regimes are analysed: diffusion of particles along inter-phase magnetic flux tubes (leaving the cocoon) and cross field transport of particles in flux tubes touching the cocoon surface. The cross field diffusion is likely the dominate escape path, unless a significant fraction of the surface is magnetically connected to the environment. Major cluster merger should strongly enhance the particle escape by two complementary mechanisms. i) The merger shock waves shred radio cocoons into filamentary structures, allowing the CRs to easily reach the radio cocoon boundary due to the changed morphology. ii) Also efficient particle losses can be expected for radio cocoons not compressed in shock waves. There, for a short period after the sudden injection of large scale turbulence, the (anomalous) cross field diffusion can be enhanced by several orders of magnitude. This lasts until the turbulent energy cascade has reached the microscopic scales, which determine the value of the microscopic diffusion coefficients.

We discuss the consequences of the spectral and morphological modification of galaxy cluster radio halos due to the Sunyaev-Zeldovich (SZ) effect for the interpretation of existing and upcoming high frequency radio observations. Likely these modifications have affected the interpretation of the existing Coma cluster radio data. The radio halo emission visible at low (< 5 GHz) frequencies is at higher (> 10 GHz) frequencies completely over-compensated by the thermal SZ decrement. Thus, the total radio emission of a galaxy cluster goes through zero (in comparison to the constant cosmic microwave background (CMB) emission) at a frequency of several GHz. Since the radio halo brightness has a narrow radial profile compared to the SZ decrement, a central emission region is surrounded by a decrement within the intermediate frequency range of several GHz. The size of this emission regions shrinks with increasing frequency until the decrement dominates everywhere in the cluster.

A significant fraction of galaxy clusters exhibits cluster wide radio halos. We give a simple prediction of the local and higher redshift radio halo luminosity function (RHLF) on the basis of (i) an observed and a theoretical X-ray cluster luminosity function (XCLF) (ii) the observed radio--X-ray luminosity correlation (RXLC) of galaxy clusters with radio halos (iii) an assumed fraction of 1/3 galaxy clusters to have radio halos as supported by observations. We then find 300-700 radio halos with S_1.4GHz > 1 mJy, and 10^5 - 10^6 radio halos with S_1.4GHz > 1 muJy should be visible on the sky. 14% of the S_1.4GHz > 1 mJy and 56% of the S_1.4GHz > 1 muJy halos are located at z>0.3. Subsequently, we give more realistic predictions taking into account (iv) a refined estimate of the radio halo fraction as a function of redshift and cluster mass, and (v) a decrease in intrinsic radio halo luminosity with redshift due to increased inverse Compton electron energy losses on the Cosmic Microwave Background (CMB). We find that this reduces the radio halo counts from the simple prediction by only 30 % totally, but the high redshift (z>0.3) counts are more strongly reduced by 50-70%. These calculations show that the new generation of sensitive radio telescopes like LOFAR, ATA, EVLA, SKA and the already operating GMRT should be able to detect large numbers of radio halos and will provide unique information for studies of galaxy cluster merger rates and associated non-thermal processes.

Radio plasma injected by active radio galaxies into clusters of galaxies quickly becomes invisible due to radiative losses of the relativistic electrons. In this talk, the fate of radio plasma and its role for the galaxy cluster is discussed: buoyancy removes it from the central regions and allows to transfer its energy into the ambient gas. The remaining low energy electron populations are still able to emit a low luminosity glow of observable radiation via synchrotron-self Comptonized emission. Shock waves in the ambient gas can re-ignite the radio emission.

We report the discovery of large-scale diffuse radio emission from what appears to be a large-scale filamentary network of galaxies in the region of cluster ZwCl 2341.1+0000, and stretching over an area of at least $6 h^{-1}_{50}$ Mpc in diameter. Multicolour CCD observations yield photometric redshifts indicating that a significant fraction of the optical galaxies in this region is at a redshift of z=0.3. This is supported by spectroscopic measurements of 4 galaxies in the SDSS survey at a mean z=0.27. We present VLA images at 20 cm (NVSS) and 90 cm wavelengths, showing the detailed radio structure of the filaments. Comparison with the VLA high resolution FIRST radio survey shows that the diffuse emission is not due to known individual point sources. The diffuse radio-emission has a spectral index $\alpha \lesssim -0.5$, and is most likely synchrotron emission from relativistic charged particles in an inter-galactic magnetic field. Furthermore, this optical/radio structure is detected in X-rays by the ROSAT all-sky survey. It has a 0.1--2.4 keV luminosity of about $10^{44}$ erg s$^{-1}$ and shows an extended highly non-relaxed morphology. These observations suggest that ZwCl 2341.1+0000 is possibly a proto-cluster of galaxies in which we are witnessing the process of structure formation. We show (both analytically and by numerical simulations) that the energetics of accretion shocks generated in forming large-scale structures are sufficient to produce enough high energy cosmic-ray (CR) electrons required to explain the observed radio emission, provided a magnetic field of strength of about 1 micro Gauss is present there.Thus it is the first evidence of cosmic-ray particle acceleration and magnetic fields occuring on a super-cluster scale. (Abridged)

A brief overview on the theory and observations of relativistic particle populations in clusters of galaxies is given. The following topics are addressed: (i) the diffuse relativistic electron population within the intra-cluster medium (ICM) as seen in the cluster wide radio halos and possibly also seen in the high energy X-ray and extreme ultraviolet excess emissions of some clusters, (ii) the observed confined relativistic electrons within fresh and old radio plasma and their connection to cluster radio relics at cluster merger shock waves, (iii) the relativistic proton population within the ICM, and its observable consequences (if it exists), and (iv) the confined relativistic proton population (if it exists) within radio plasma. The importance of upcoming, sensitive gamma-ray telescopes for this research area is highlighted.

We present a new wide-field map of the radio galaxy 3C 129 and its companion galaxy 3C 129.1 at lambda=90 cm. We see a distinct steep-spectrum feature near the head of 3C 129, extending in a direction perpendicular to the radio tails. We propose that this Crosspiece might consist of fossil radio plasma, which has been re-energized by the compression of the bow shock of the supersonically moving galaxy 3C 129. One possible origin of the fossil radio plasma could be the tail of a nearby head-tail radio galaxy. We discuss the implications of, and give testable predictions for this scenario.

The Chandra X-ray Observatory is finding a surprisingly large number of cavities in the X-ray emitting intracluster medium, produced by the release of radio plasma from active galactic nuclei. In this Letter, we present simple analytic formulae for the evolution of the X-ray deficit and for the radio spectrum of a buoyantly rising bubble. The aim of this work is to provide a theoretical framework for the planning and the analysis of X-ray and radio observations of galaxy clusters. We show that the cluster volume tested for the presence of cavities by X-ray observations is a strongly rising function of the sensitivity.

Powerful quasars and radio galaxies are injecting large amounts of energy in the form of radio plasma into the inter-galactic medium (IGM). Once this nonthermal component of the IGM has radiatively cooled the remaining radio emission is difficult to detect. Two scenarios in which the fossil radio plasma can be detected and thus be used to probe the IGM are discussed: a) re-illumination of the radio emission due to the compression in large-scale shock waves, and b) inverse Compton scattered radiation of the cosmic microwave background (CMB), cosmic radio background (CRB), and the internal very low frequency synchrotron emission of the still relativistic low energy electron population. We present 3-D magneto-hydrodynamical simulations of scenario a) and compare them to existing observations. Finally, we discuss the feasibility of the detection of process b) with upcoming instruments.

Most of the Universe's populations of low energy cosmic ray electrons in the energy range of 1-100 MeV still manage to elude from detection by our instruments, since their synchrotron emission is at too low frequencies. We investigate a mechanism which can lead to observable emission of such electron populations: synchrotron-self Comptonization (SSC). The inverse Compton (IC) scattering can shift an otherwise unobservable low-frequency 10 kHz-10 MHz photons into observable radio, infrared (IR) or even more energetic wave-bands. The resulting emission should be polarized. We also consider IC scattering of the cosmic microwave background (CMB) and the cosmic radio background (CRB). Electron spectral aging due to synchrotron, IC and adiabatic losses or gains influences the resulting spectrum. The predicted radiation spectrum is a strong function of the history of the source, of the low energy spectrum of relativistic electrons, and of redshift. It has typically two maxima, and a decrement in between at CMB frequencies. Detection will give a sensitive probe of the environment of radio galaxies. We demonstrate that the fossil remnants of powerful radio galaxies are promising detection candidates, especially when they are embedded in a dense intra-cluster medium (ICM). GHz peaked sources (GPS) have very low SSC luminosities, which may, however, extend into the X-ray or even the gamma ray regime. Clusters of galaxies with relativistic electron populations may be another detectable sources. Fossil radio plasma released by our own Galaxy could be revealed by its large angular scale SSC flux. We discuss the expected detectability of these sources with new and upcoming instruments as LOFAR, GMRT, EVLA, ATA, ALMA, PLANCK, and HERSCHEL.

(abridged) We present detailed 3-dimensional magneto-hydrodynamical simulations of the passage of a radio plasma cocoon filled with turbulent magnetic fields through a shock wave. Taking into account synchrotron, inverse Compton and adiabatic energy losses and gains we evolved the relativistic electron population to produce synthetic polarisation radio maps. On contact with the shock wave the radio cocoons are first compressed and finally torn into filamentary structures, as is observed in several cluster radio relics. In the synthetic radio maps the electric polarisation vectors are mostly perpendicular to the filamentary radio structures. If the magnetic field inside the cocoon is not too strong, the initially spherical radio cocoon is transformed into a torus after the passage of the shock wave. Very recent, high-resolution radio maps of cluster radio relics seem to exhibit such toroidal geometries in some cases. This supports the hypothesis that cluster radio relics are fossil radio cocoons that have been revived by a shock wave. For a late-stage relic the ratio of its global diameter to the filament diameter should correlate with the shock strength. Finally, we argue that the total radio polarisation of radio relic should be well correlated with the 3-dimensional orientation of the shock wave that produced the relic.

We present three-dimensional hydrodynamical simulations of buoyant gas in a typical cluster environment. The hot matter was injected continuously into a small region off-set from the cluster centre. In agreement with previous analytic estimates we found that the bubbles evolve very differently depending on their luminosity. Using tracer particles we computed radio maps of the bubbles based on different assumptions about the magnetic field. In the radio band the bubbles closely resemble FRI sources. For the bubbles to be detectable for long enough to account for FRI sources, we found that reacceleration has to take place. The bubbles are generally difficult to detect, both, in the radio and in the X-ray band. Thus it is possible to hide a significant amount of energy in the form of bubbles in clusters. Finally, we compute the efficiency of the bubbles to stir the ICM and find that recurrent low-power sources may be fairly effective in mixing the inner cluster region.

We present deep VLA observations of the galaxy clusters Abell 2256 and Abell 754, both of which appear to be in the violent stage of major cluster merger events. The complex nature of Abell 2256 is revealed through radio images which, in addition to the head-tail galaxies, show two extended, irregular, and sharp-edged regions of diffuse radio emission at the cluster periphery (so called radio relics), and a large-scale diffuse radio halo located in the central regions of the cluster. Polarimetry of the A2256 cluster relics reveals large-scale ordered magnetic fields which appear to trace the bright filaments in the relics. The polarization fraction across the relics ranges from 20% - 40% with the majority of the relics polarized above the 30% level. At the sensitivity of our current observations we place an upper limit of 20% on the polarization of the radio halo. Low frequency VLA observations of Abell 754 reveal extended, diffuse radio (halo) emission in the cluster core region as well as steep spectrum emission in the cluster periphery. The location, morphology, and spectral index of the peripheral emission are consistent with the properties of radio relics. The X-ray evidence of the ongoing mergers in both clusters, together with the polarization properties of A2256's radio relics supports recent suggestions of a merger-induced origin of the relic emission. Deciphering the complex radio properties of these clusters may thus provide the key to understanding the dynamical history of the systems.

In several merging clusters of galaxies so-called cluster radio relics have been observed. These are extended radio sources which do not seem to be associated with any radio galaxy. Two competing physical mechanisms to accelerate the radio emitting electrons have been proposed: (i) diffusive shock acceleration and (ii) adiabatic compression of fossil radio plasma by merger shock waves. Here the second scenario is investigated. We present detailed 3-dimensional magneto-hydrodynamical simulations of the passage of a radio plasma cocoon through a shock wave. Taking into account synchrotron, inverse Compton and adiabatic energy losses and gains we evolved the relativistic electron population to produce synthetic radio maps in Stokes I-, Q-, and U-polarisation. In the synthetic radio maps the electric polarisation vectors are mostly perpendicular to the filamentary radio structures.

Clusters of galaxies may contain cluster-wide, centrally located, diffuse radio sources, called halos. They have been found to show morphologies similar to those of the X-ray emission. To quantify this qualitative statement we performed a point-to-point comparison of the radio and the X-ray emission for four clusters of galaxies containing radio halos: Coma, Abell 2255, Abell 2319, Abell 2744. Our study leads to a linear relation between the radio and the X-ray surface brightness in two clusters, namely Abell 2255 and Abell 2744. In Coma and A2319 the radio and the X-ray brightnesses seem to be related with a sub-linear power law. Implications of these findings within simple radio halo formation models are briefly discussed.

Plasma containing relativistic particles appears in various forms in the inter galactic medium (IGM): As radio plasma released by active radio galaxies, as fossil radio plasma from former radio galaxies - so called `radio ghosts', as `cluster radio relics' in some clusters of galaxies, and as `cluster radio halos'. The impact of the different forms of radio plasma on the IGM and their use as diagnostic tools are briefly discussed.

Sensitive images of low-level, Mpc-sized radio cocoons offer new opportunities to probe large scale intergalactic gas flows outside clusters of galaxies. New radio images of high surface brightness sensitivity at strategically chosen wavelengths of the giant radio galaxy NGC 315 (Mack et al. 1997,1998) reveal significant asymmetries and particularities in the morphology, radio spectrum and polarization of the ejected radio plasma. We argue that the combination of these signatures provides a sensitive probe of an environmental shock wave. Analysis of optical redshifts in NGC 315 vicinity confirms its location to be near, or at a site of large-scale flow collisions in the 100 Mpc sized Pisces-Perseus Supercluster region. NGC 315 resides at the intersection of several galaxy filaments, and its radio plasma serves there as a `weather station' (Burns 1998) probing the flow of the elusive and previously invisible IGM gas. If our interpretation is correct, this is the first indication for a shock wave in flows caused by the cosmological large scale structure formation, which is located in a filament of galaxies. The possibility that the putative shock wave is a source of gamma-rays and ultra high energy cosmic rays is briefly discussed.

The isotropy in the ultra high energy cosmic ray (UHECR) flux observed by Yakutsk and AGASA experiments, is a very strong constraint to production and propagation models alike. Most of the scenarios proposed in the literature should produce a sizable anisotropy as either extragalactic luminous or dark matter is normally associated with the invoked particle sources. We explore the possibility that the magnetic fields in fossil cocoons of former radio galaxies -- so called {\it radio ghosts} -- are able to scatter UHECR in the intergalactic medium giving rise to the observed isotropy. We show, through numerical simulations, under which conditions this process can be operative and the magnitude of the effect. We further demonstrate, that if radio ghosts mix with the ambient medium, they might be able to produce the observed magnetic fields in clusters of galaxies. In the case of mixing, the UHECR isotropization would be even stronger than in our conservative estimates.

We give for a plasma with a history of several expansion and contraction phases an analytical model of the evolution of a contained relativistic electron population under synchrotron, inverse Compton and adiabatic energy losses or gains. This is applied to different scenarios for evolution of radio plasma inside the cocoons of radio galaxies, after the activity of the central engine has ceased. It is demonstrated that fossil radio plasma with an age of even up to 2 Gyr can be revived by compression in a shock wave of large-scale structure formation, caused during the merging events of galaxy clusters, or by the accretion onto galaxy clusters. We argue, that this is a highly plausible explanation for the observed cluster radio relics, which are the regions of diffuse radio emission found in clusters of galaxies, without any likely parent radio galaxy seen nearby. An implication of this model is the existence of a population of diffuse, ultra-steep spectrum, very low frequency radio sources located inside and possibly outside of clusters of galaxies, tracing the revival of aged fossil radio plasma by the shock waves associated with large-scale structure formation.

The strong activity of radio galaxies should have led to a nearly ubiquitous presence of fossil radio plasma in the denser regions of the inter-galactic medium as clusters, groups and filaments of galaxies. This fossil radio plasma can contain large quantities of relativistic particles (electrons and possibly protons) by magnetic confinement. These particles might be released and/or re-energized under environmental influences as turbulence and shock waves. Possible connections of such processes to the formation of the observed sources of diffuse radio emission in clusters of galaxies (the cluster radio halos and the cluster radio relics) are discussed.

A new model for the, so called, `cluster radio relics' is presented (Ensslin & Gopal-Krishna 2000). Fossil radio cocoons, resulting from the former activity of radio galaxies, should contain a low energy relativistic electron population and magnetic fields. Electrons with an age of even up to 2 Gyr can be re-accelerated adiabatically to radio emitting energies, if the fossil radio plasma gets compressed in an environmental shock wave. Such a wave can be caused by a merging event of galaxy clusters, or by the accretion onto galaxy clusters. An implication of this model is the existence of a population of diffuse, ultra-steep spectrum, very low-frequency radio sources located inside and possibly outside of clusters of galaxies, tracing the revival of aged fossil radio plasma by the shock waves associated with large-scale structure formation.

We investigate the possibility that radio halos of clusters of galaxies are caused by synchrotron emission of cosmic ray electrons (CRe), which were produced by cosmic ray protons (CRp) interacting hadronically with the intra-cluster medium (ICM) protons. We perform cosmological magneto-hydrodynamics (MHD) simulations to obtain a sample of ten magnetized galaxy clusters. They provide realistic models of the gas and magnetic field distribution, needed to predict the CRe production rates, their cooling, and their synchrotron emissivity. We assume a CRp population within the ICM with an energy density which has a constant ratio to thermal energy density. This ratio is adjusted in such a way that one of the simulated clusters reproduces the radio luminosity of the radio halo of the Coma cluster of galaxies. Our model exhibits the observed low degree of radio polarization and has a similar radial emission profile as the Coma cluster. We provide estimates for the expected gamma ray and neutrino flux. The necessary CRp/thermal energy ratio is 4 ... 14% ${(E_{p,min}/GeV)}^{-0.375}$ (for the range of magnetic field strengths suggested by Faraday measurements), where $E_{p,min}$ is the lower kinetic energy cutoff of the CRp with spectral index $\alpha_p \approx 2.375$. Assuming this ratio to be the same in the whole set of simulated clusters a $T_x-L_\nu$ relation is predicted which follows the observed relation well.

(Abridged) We investigate the spectral distortion of the cosmic microwave background (CMB) caused by relativistic plasma. Within the Thomson regime, an exact analytic expression for the photon scattering kernel of a momentum power-law electron distribution is given. The ultra-relativistic Sunyaev-Zeldovich (SZ) decrement measures the electron number and not the energy content. On the other hand, the relativistic SZ increment at higher frequencies depends strongly on the spectral shape of the electrons. We calculate the expected Comptonization due to the energy release of radio galaxies, which we estimate to be 3 10^66 erg Gpc^-3. We investigate Comptonization from (a) the part of the released energy which is thermalized and (b) the relativistic, remnant radio plasma, which may form a second, relativistic phase in the inter galactic medium, nearly unobservable for present day instruments (presence of so called 'radio ghosts'). We find a thermal Comptonization parameter due to (a) of y = 10^-6 and (b) an optical depth of relativistic electrons in old radio plasma of tau_rel < 10^-7. If a substantial fraction of the volume of clusters of galaxies is filled with such old radio plasma the SZ effect based determination of the Hubble constant is biased to lower values. Finally, it is shown that a supra-thermal population of electrons in the Coma cluster would produce a signature in the Wien-tail of the CMB, which is marginally detectable with a multifrequency measurement by the Planck satellite. Such a population is expected to exist, since its bremsstrahlung would explain Coma's recently reported high energy X-ray excess.

A short review on theoretical implications of non-thermal emission (radio, extreme ultraviolet, high energy X-ray) from the intra-cluster medium is given. The origin of cluster radio halos and cluster radio relics is discussed within the framework of a network of processes producing a non-thermal electron population. Emphasis is given to the role of old, remnant, presently invisible relativistic plasma released by former radio galaxies.

We report on an ongoing project to map a 15 degree x 15 degree region in Coma Berenices with the new 74 MHz receiver system of the VLA. The field is centered on the Coma cluster of galaxies. Coma Berenices has been observed in all VLA configurations, which allowed us to resolve several hundred point sources with 35 arcsecond resolution. We have successfully detected the diffuse emission from the Coma cluster radio halo, radio relic, and radio bridge. Difficulties and strategies of the observation and data reduction are discussed. Preliminary results are given.

We investigate the possibility that patches of old radio plasma (`radio ghosts') of former radio galaxies form a second distinct phase of the inter-galactic medium (IGM), not mixed with the thermal gas. The separation of this phase from the ambient gas and its resistance against eroding turbulent forces is given by magnetic fields, which are expected to be roughly in pressure equilibrium with the surrounding medium. Since patches of this plasma are largely invisible in the radio we use the term `radio ghost' to characterize their nature. Possibilities and difficulties of different detection strategies of ghosts are discussed. These involve radio emission, cosmic microwave background (CMB) and starlight Comptonization, and Faraday rotation. Re-activation of the electron population in shock waves of cosmological structure formation, which seems to lead to the cluster radio relic phenomena. We discuss the role radio ghosts can have: They are able to store relativistic particles for cosmological times, but are also able to release them under the influence of very strong turbulence. This might happen during a major merger event of clusters of galaxies. The released relativistic proton population could produce the observed radio halos of some cluster of galaxies via hadronic reactions with the background gas leading to the production of secondary electrons and positrons. Destroyed ghosts, mixed with the IGM can help to magnetize it. Finally, the strong field strength within ghosts should have a significant impact on the propagation of extragalactic high energy cosmic rays.

The excess in extreme-ultraviolet (EUV) radiation and high energy X-rays (HEX) of Coma may be modeled by power law spectra. Their origin could be inverse Compton (IC) emission by relativistic electrons. The scattered photon fields can either be the CMB or the starlight of the galaxies within the cluster. If the EUV excess is due to CMB-IC, a strongly inhomogeneous magnetized intra- cluster medium (ICM) is required, in which the density of the relativistic electrons is anticorrelated with the magnetic fields. This could be understood if the electrons were accelerated during a major merger event within the last 2 Gyr and cooled afterwards in the inhomogeneous fields. If the EUV excess is due to starlight-IC, a population of relativistic, very low energy electrons has to be present, with a high energy density. In order to survive Coulomb losses, these electrons have to be separated from the cluster gas by confining magnetic fields. This could be remnant radio plasma from the epoch of violent quasar activity. The observed narrow radial profile of the EUV excess emission is due to the narrow photon distribution for this model. Both models favor very inhomogeneous magnetic field and relativistic electron distributions. IC models for the HEX excess seem to be implausible: CMB-IC leads to a mainly unmagnetized ICM, in contradiction to Faraday rotation measurements. Starlight scattering electrons would overproduce EUV photons due to simultaneously CMB-IC. We propose that the observed HEX excess is due to bremsstrahlung of a small high energy power-law tail of the mainly thermal ICM electron distribution. Such a tail is expected due to in-situ acceleration by turbulence.

EGRET data on the Gamma ray emission from the inner Galaxy have shown a rather flat spectrum, extending to about 50 GeV. It is usually assumed that these gamma-rays arise from the interactions of cosmic ray nuclei with ambient matter. Cosmic Ray particles have been observed up to 3 10^20 eV, with many arguments suggesting, that up to about 3 10^18 eV they are of Galactic origin. Cosmic ray particles get injected by supernova explosions. Their injected spectrum is steepened by diffusive losses from the Galaxy to yield the observed spectrum. As cosmic ray particles roam around in the Galactic disk, they encounter molecular clouds and through p-p collisions produce gamma rays from pion decay. The flux and spectrum of these gamma rays is then a clear signature of cosmic rays throughout the Galaxy. Star formation activity peaks in the central region of the Galaxy, around the Galactic Center. Looking then at the gamma ray spectrum of the central region of our Galaxy yields clues as to where the cosmic ray particles interact, and with what spectrum. Using the FLUKA Monte-Carlo, we have modelled this spectrum, and find a best fit for a powerlaw spectral index of 2.34, rather close to the suggested injection spectrum for supernovae which explode into their own winds. This suggests that most cosmic ray interaction happens near the sources of injection; this is consistent with the spectrum of cosmic ray nuclei derived from spallation. One consequence is that cosmic ray heating and ionization should be strong in the Galactic Center region.

Simulations of structure formation in the Universe predict accretion shock waves at the boundaries of the large-scale structures as sheets, filaments, and clusters of galaxies. If magnetic fields are present at these shocks, particle acceleration should take place, and could contribute to the observed cosmic rays of high energies. When the radio plasma of an old invisible radio lobe is dragged into such a shock wave, the relativistic electron population will be reaccelerated and the plasma becomes radio-luminous again. Such tracers of the accretion shock waves are observed at the boundaries of some clusters of galaxies: the so-called cluster radio relics, which are large regions of diffuse radio emission, without any parent galaxy nearby. The observed properties of the cluster radio relics are naturally explained by accretion shock waves. Radio relics therefore give the first evidence for the existence of accretion shocks of the large-scale structure formation and they allow investigations of the shock properties.

We estimate the energy release of black hole formation to the intra-cluster medium of the Coma cluster and find that this is comparable to the present day energy content. Therefore an energetic and maybe hydrostatic influence is possible. Our calculations rely on the assumption of an universal black hole to galaxy mass ratio (more exactly: spheroidal mass component of the stellar population), for which there is growing evidence. On a cosmological scale, there is also an energy release of black hole formation comparable to what is expected to be present within the thermal gas, caused by the process of structure formation. This indicates an important dynamical influence, neglected by present day structure formation simulations. This estimate of cosmological black hole energy release is independent of the black hole to galaxy mass ratio, but consistent with its value.

We present evidence for the existence of shock waves caused by the formation of the large-scale structure. In some clusters of galaxies peripherally located sources of extended diffuse radio emission exist, the so-called cluster radio relics. They have steep radio spectra but no apparent cutoff, as old remnants of radio galaxies usually have. Therefore particle acceleration has to take place within them. We propose that shock structures of the cosmological large-scale matter flows are responsible for the acceleration of relativistic electrons: cluster accretion shocks and bow shocks of merger events. We develop a theory of radio plasma having traversed these shocks and compare it to observational data of nine radio relics (0038-096, 0917+75, 1140+203, 1253+275, 1712+64, 1706+78, 2006-56, 2010-57, 1401-33) and their host clusters (A85, A786, A1367, Coma, A2255, A2256, A3667, S753). The necessary accretion power, the spectral index of the radio spectrum, the acceleration efficiency of the shock, the diffusion coefficient in the post-shock region, and the predicted radio polarization in all of our examples fit into a coherent interpretation of the observational data. Since polarization measurements are available only for four sources, the predictions of our theory can be independently checked using other examples. The predicted values of the shock compression ratio, density and temperature of the infalling gas, magnetic field strength of the shocked and unshocked radio plasma are discussed within the frame of structure formation theory.

Relativistic electrons, seen in the large diffuse radio halo of the Coma cluster of galaxies, should scatter background photons to higher energies. We calculate the inverse Compton contributions from the microwave background, from the local radiation field of elliptical galaxies and from the thermal X-ray emission of the intra cluster medium to the different observed energy bands, and draw restrictions to the shape of the spectrum of the relativistic electron population. The expected electron spectra in different halo formation models are discussed, and signatures for a future distinction of these models are presented, tracing the injection processes and the influence of optical thickness of the radio halo for very low frequency emission due to synchrotron self absorption. Some of the upper limits on the electron spectrum translate into lower limits of the halo magnetic fields, since the observed synchrotron radiation is known. Despite uncertainties in the extrapolation of the radio spectrum to lower frequencies a lower limit to the central field strength of B_o > 0.3 muG seems to be sure. This is comparable to the value resulting from minimal energy arguments. If B_o < 1.2 muG, the recently reported extreme ultraviolet excess of Coma (Lieu et al. 1996) could result from inverse Compton scattered microwave photons. Formulae for transrelativistic Thomson scattering are given in the Appendix.

The masses of clusters of galaxies estimated by gravitational lensing exceed in many cases the mass estimates based on hydrostatic equilibrium. This may suggest the existence of nonthermal pressure. We ask if radio galaxies can heat and support the cluster gas with injected cosmic ray protons and magnetic field densities, which are permitted by Faraday rotation and gamma ray observations of clusters of galaxies. We conclude that they are powerful enough to do this within a cluster radius of roughly 1 Mpc. If present, nonthermal pressures could lead to a revised estimate of the ratio of baryonic mass to total mass, and the apparent baryonic overdensity in clusters would disappear. In consequence, $\Omega_{\rm cold}$, the clumping part of the cosmological density $\Omega_{o}$, would be larger than $0.4\,h_{50}^{-1/2}$.