Observations of Faraday rotation for extragalactic sources probe magnetic fields both inside and outside the Milky Way. Building on our earlier estimate of the Galactic contribution, we set out to ...estimate the extragalactic contributions. We discuss the problems involved; in particular, we point out that taking the difference between the observed values and the Galactic foreground reconstruction is not a good estimate for the extragalactic contributions. We point out a degeneracy between the contributions to the observed values due to extragalactic magnetic fields and observational noise and comment on the dangers of over-interpreting an estimate without taking into account its uncertainty information. To overcome these difficulties, we develop an extended reconstruction algorithm based on the assumption that the observational uncertainties are accurately described for a subset of the data, which can overcome the degeneracy with the extragalactic contributions. We present a probabilistic derivation of the algorithm and demonstrate its performance using a simulation, yielding a high quality reconstruction of the Galactic Faraday rotation foreground, a precise estimate of the typical extragalactic contribution, and a well-defined probabilistic description of the extragalactic contribution for each data point. We then apply this reconstruction technique to a catalog of Faraday rotation observations for extragalactic sources. The analysis is done for several different scenarios, for which we consider the error bars of different subsets of the data to accurately describe the observational uncertainties. By comparing the results, we argue that a split that singles out only data near the Galactic poles is the most robust approach. We find that the dispersion of extragalactic contributions to observed Faraday depths is most likely lower than 7 rad/m2, in agreement with earlier results, and that the extragalactic contribution to an individual data point is poorly constrained by the data in most cases.
We present resolve, a new algorithm for radio aperture synthesis imaging of extended and diffuse emission in total intensity. The algorithm is derived using Bayesian statistical inference techniques, ...estimating the surface brightness in the sky assuming a priori log-normal statistics. resolve estimates the measured sky brightness in total intensity, and the spatial correlation structure in the sky, which is used to guide the algorithm to an optimal reconstruction of extended and diffuse sources. During this process, the algorithm succeeds in deconvolving the effects of the radio interferometric point spread function. Additionally, resolve provides a map with an uncertainty estimate of the reconstructed surface brightness. Furthermore, with resolve we introduce a new, optimal visibility weighting scheme that can be viewed as an extension to robust weighting. In tests using simulated observations, the algorithm shows improved performance against two standard imaging approaches for extended sources, Multiscale-CLEAN and the Maximum Entropy Method.
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 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. Because cosmic magnetism is one of the key science projects of the new generation of radio interferometers, we have predicted the performance of our algorithm on mock data collected with these instruments. According to our tests, high-quality catalogs of a few thousands of sources should already enable us to investigate magnetic fields in the cosmic structure.
NIFTy (Numerical Information Field Theory) is a software package designed to enable the development of signal Inference algorithms that operate regardless of the underlying spatial grid and its ...resolution. Its object-oriented framework is written in Python, although It accesses libraries written in Cython, C++, and C for efficiency. NIFTY offers a toolkit that abstracts discretized representations of continuous spaces, fields in these spaces, and operators acting on fields Into classes. Thereby, the correct normalization of operations on fields is taken care of automatically without concerning the user. This allows for an abstract formulation and programming of inference algorithms, including those derived within Information field theory. Thus, NIFTY permits its user to rapidly prototype algorithms in 1D, and then apply the developed code in higher-dimensional settings of real world problems. The set of spaces on which NIFTy operates comprises point sets, n-dimensional regular grids, spherical spaces, their harmonic counterparts, and product spaces constructed as combinations of those. The functionality and diversity of the package is demonstrated by a Wiener filter code example that successfully runs without modification regardless of the space on which the inference problem is defined.
We aim to summarize the current state of knowledge regarding Galactic Faraday rotation in an all-sky map of the Galactic Faraday depth. For this we have assembled the most extensive catalog of ...Faraday rotation data of compact extragalactic polarized radio sources to date. In the map-making procedure we used a recently developed algorithm that reconstructs the map and the power spectrum of a statistically isotropic and homogeneous field while taking into account uncertainties in the noise statistics. This procedure is able to identify some rotation angles that are offset by an integer multiple of π. The resulting map can be seen as an improved version of earlier such maps and is made publicly available, along with a map of its uncertainty. For the angular power spectrum we find a power law behavior Cℓ ∝ ℓ-2.17 for a Faraday sky where an overall variance profile as a function of Galactic latitude has been removed, in agreement with earlier work. We show that this is in accordance with a 3D Fourier power spectrum P(k) ∝ k-2.17 of the underlying field neBr under simplifying geometrical and statistical assumptions.
We report the detection of diffuse radio emission which might be connected to a large-scale filament of the cosmic web covering a 8° × 8° area in the sky, likely associated with a z ≈ 0.1 overdensity ...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 intracluster 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 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.
Abstract
We observed the galaxy cluster CIZA J2242.8+5301 with the Sardinia Radio Telescope to provide new constraints on its spectral properties at high frequency. We conducted observations in three ...frequency bands centred at 1.4, 6.6 and 19 GHz, resulting in beam resolutions of 14, 2.9 and 1 arcmin, respectively. These single-dish data were also combined with archival interferometric observations at 1.4 and 1.7 GHz. From the combined images, we measured a flux density of S
1.4 GHz = (158.3 ± 9.6) mJy for the central radio halo and S
1.4 GHz = (126 ± 8) and (11.7 ± 0.7) mJy for the northern and the southern relics, respectively. After the spectral modelling of the discrete sources, we measured at 6.6 GHz S
6.6 GHz = (17.1 ± 1.2) and (0.6 ± 0.3) mJy for the northern and southern relics, respectively. Assuming simple diffusive shock acceleration, we interpret measurements of the northern relic with a continuous injection model represented by a broken power law. This yields an injection spectral index αinj = 0.7 ± 0.1 and a Mach number M = 3.3 ± 0.9, consistent with recent X-ray estimates. Unlike other studies of the same object, no significant steepening of the relic radio emission is seen in data up to 8.35 GHz. By fitting the southern relic spectrum with a simple power law (S
ν ∝ ν−α), we obtained a spectral index α ≈ 1.9 corresponding to a Mach number (M ≈ 1.8) in agreement with X-ray estimates. Finally, we evaluated the rotation measure of the northern relic at 6.6 GHz. These results provide new insights on the magnetic structure of the relic, but further observations are needed to clarify the nature of the observed Faraday rotation.
Abstract
The optimal source-finding strategy for linear polarization data is an unsolved problem, with many inhibitive factors imposed by the technically challenging nature of polarization ...observations. Such an algorithm is essential for Square Kilometre Array (SKA) pathfinder surveys, such as the Multifrequency Snapshot Sky Survey with the LOw Frequency ARray (LOFAR), as data volumes are significant enough to prohibit manual inspection. We present a new strategy of ‘Faraday Moments’ for source-finding in linear polarization with LOFAR, using the moments of the frequency-dependent full-Stokes data (i.e. the mean, standard deviation, skewness, and excess kurtosis). Through simulations of the sky, we find that moments can identify polarized sources with a high completeness: 98.5 per cent at a signal to noise of 5. While the method has low reliability, rotation measure (RM) synthesis can be applied per candidate source to filter out instrumental and spurious detections. This combined strategy will result in a complete and reliable catalogue of polarized sources that includes the full sensitivity of the observational bandwidth. We find that the technique can reduce the number of pixels on which RM Synthesis needs to be performed by a factor of ≈1 × 105 for source distributions anticipated with modern radio telescopes. Through tests on LOFAR data, we find that the technique works effectively in the presence of diffuse emission. Extensions of this method are directly applicable to other upcoming radio surveys such as the POlarization Sky Survey of the Universe's Magnetism with the Australia Square Kilometre Array Pathfinder, and the SKA itself.