In this paper we present a detailed analysis of the Faraday depth (FD) spectrum and its clean components obtained through the application of the commonly used technique of Faraday rotation measure synthesis to analyze spectro-polarimetric data. In order to directly compare the Faraday depth spectrum with physical properties of a magneto-ionic medium, we generated synthetic broad-bandwidth spectro-polarimetric observations from magnetohydrodynamic (MHD) simulations of a transonic, isothermal, compressible turbulent medium. We find that correlated magnetic field structures give rise to a combination of spiky, localized peaks at certain FD values, and broad structures in the FD spectrum. Although the majority of these spiky FD structures appear narrow, giving an impression of a Faraday thin medium, we show that they arise from strong synchrotron emissivity at that FD. Strong emissivity at a FD can arise because of both strong spatially-local polarized synchrotron emissivity at a FD or accumulation of weaker emissions along the distance through a medium that have Faraday depths within half the width of the rotation measure spread function. Such a complex Faraday depth spectrum is a natural consequence of MHD turbulence when the lines of sight pass through a few turbulent cells. This therefore complicates the convention of attributing narrow FD peaks to presence of a Faraday rotating medium along the line of sight. Our work shows that it is difficult to extract the FD along a line of sight from the Faraday depth spectrum using standard methods for a turbulent medium in which synchrotron emission and Faraday rotation occur simultaneously.

We present S-PASS/ATCA, the first wide-band radio polarimetry survey of compact sources in the southern sky. We describe how we selected targets for observations with the Australia Telescope Compact Array (ATCA) in the 16 cm band (1.3 - 3.1 GHz), our observing and calibration strategy, how we analysed the data, and how we tested the quality of the data. The data are made publicly available. The survey contains on average one source per five square degrees and has an angular resolution at 2.2 GHz of ~ 2'x1'. Sources with |RM|s > 150 rad m-2 are seen towards the Galactic plane and bright Hii regions, but are rare elsewhere on the sky. Sightlines that are separated by up to 3' show very similar RMs. Based on this observation, we argue that the Galactic foreground is the dominant contributor to RM, confirming previous results, and that the sources must have very simple distributions of Faraday-rotating and synchrotron-emitting media. Many sources that emit at a single RM have a spectral index in linear polarization that is (very) different from the spectral index in Stokes I. Analysing ratios of flux densities Q/I and U/I (to correct for spectral index effects) then leads to erroneous results. About 80 per cent of sightlines in our survey are dominated by emission at only one RM. Therefore, RMs that were determined previously from narrow-band observations at these frequencies are still safe to use.

Faraday tomography allows us to map diffuse polarized synchrotron emission from our Galaxy and use it to interpret the magnetic field in the interstellar medium (ISM). We have applied Faraday tomography to 60 observations from the LOFAR Two-meter Sky Survey (LOTSS) and produced a Faraday depth cube mosaic covering 568 square degrees at high Galactic latitudes, at 4.3' angular resolution and 1 rad m$^{-2}$ Faraday depth resolution, with a typical noise level of 50--100 $\mu$Jy per point spread function (PSF) per rotation measure spread function (RMSF) (40-80 mK RMSF$^{-1}$). While parts of the images are strongly affected by instrumental polarization, we observe diffuse polarized emission throughout most of the field, with typical brightness between 1 and 6 K RMSF$^{-1}$, and Faraday depths between $-7$ and +25 rad m$^{-2}$. We observed many new polarization features, some up to 15 degrees in length. These include two regions with very uniformly structured, linear gradients in the Faraday depth; we measured the steepness of these gradients as 2.6 and 13 rad m$^{-2}$ deg$^{-1}$. We also observed a relationship between one of the gradients and an HI filament in the local ISM. Other ISM tracers were also checked for correlations with our polarization data and none were found, but very little signal was seen in most tracers in this region. We conclude that the LOTSS data are very well suited for Faraday tomography, and that a full-scale survey with all the LOTSS data has the potential to reveal many new Galactic polarization features and map out diffuse Faraday depth structure across the entire northern hemisphere.

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

To constrain the large-scale magnetic field strengths in cosmologically distant galax- ies, we derive the probability distribution function of Faraday rotation measure (RM) when random lines of sight pass through a sample of disc galaxies, with axisymmetric large-scale magnetic fields. We find that the width of the RM distribution of the galaxy sample is directly related to the mean large-scale field strength of the galaxy population, provided the dispersion within the sample is lower than the mean value. In the absence of additional constraints on parameters describing the magneto-ionic medium of the intervening galaxies, and in the situation where RMs produced in the intervening galaxies have already been statistically isolated from other RM contributions along the lines of sight, our simple model of the magneto-ionic medium in disc galaxies suggests that the mean large-scale magnetic field of the population can be measured to within ~ 50% accuracy.

Deriving the Faraday rotation measure (RM) of quasar absorption line systems, which are tracers of high-redshift galaxies intervening background quasars, is a powerful tool for probing magnetic fields in distant galaxies. Statistically comparing the RM distributions of two quasar samples, with and without absorption line systems, allows one to infer magnetic field properties of the intervening galaxy population. Here, we have derived the analytical form of the probability distribution function (PDF) of RM produced by a single galaxy with an axisymmetric large-scale magnetic field. We then further determine the PDF of RM for one random sight line traversing each galaxy in a population with a large-scale magnetic field prescription. We find that the resulting PDF of RM is dominated by a Lorentzian with a width that is directly related to the mean axisymmetric large-scale field strength $\langle B_0 \rangle$ of the galaxy population if the dispersion of $B_0$ within the population is smaller than $\langle B_0 \rangle$. Provided that RMs produced by the intervening galaxies have been successfully isolated from other RM contributions along the line of sight, our simple model suggests that $\langle B_0 \rangle$ in galaxies probed by quasar absorption line systems can be measured within $\approx50$ per cent accuracy without additional constraints on the magneto-ionic medium properties of the galaxies. Finally, we discuss quasar sample selection criteria that are crucial to reliably interpret observations, and argue that within the limitations of the current database of absorption line systems, high-metallicity damped Lyman-$\alpha$ absorbers are best suited to study galactic dynamo action in distant disc galaxies.

The polarization properties of radio sources at very low frequencies (<200 MHz) have not been widely measured, but the new generation of low-frequency radio telescopes, including the Low Frequency Array (LOFAR: a Square Kilometre Array Low pathfinder), now gives us the opportunity to investigate these properties. In this paper, we report on the preliminary development of a data reduction pipeline to carry out polarization processing and Faraday tomography for data from the LOFAR Two-meter Sky Survey (LOTSS) and present the results of this pipeline from the LOTSS preliminary data release region (10h45m - 15h30m right ascension, 45 - 57 degrees declination, 570 square degrees). We have produced a catalog of 92 polarized radio sources at 150 MHz at 4.3 arcminute resolution and 1 mJy rms sensitivity, which is the largest catalog of polarized sources at such low frequencies. We estimate a lower limit to the polarized source surface density at 150 MHz, with our resolution and sensitivity, of 1 source per 6.2 square degrees. We find that our Faraday depth measurements are in agreement with previous measurements and have significantly smaller errors. Most of our sources show significant depolarization compared to 1.4 GHz, but there is a small population of sources with low depolarization indicating that their polarized emission is highly localized in Faraday depth. We predict that an extension of this work to the full LOTSS data would detect at least 3400 polarized sources using the same methods, and probably considerably more with improved data processing.

Polarised radio emission from PSR J1745-2900 has already been used to investigate the strength of the magnetic field in the Galactic Centre, close to Sagittarius A*. Here we report how persistent radio emission from this magnetar, for over four years since its discovery, has revealed large changes in the observed Faraday rotation measure, by up to 3500 rad m$^{-2}$ (a five per cent fractional change). From simultaneous analysis of the dispersion measure, we determine that these fluctuations are dominated by variations in the projected magnetic field, rather than the integrated free electron density, along the changing line of sight to the rapidly moving magnetar. From a structure function analysis of rotation measure variations, and a recent epoch of rapid change of rotation measure, we determine a minimum scale of magnetic fluctuations of size ~ 2 au at the Galactic Centre distance, inferring PSR J1745-2900 is just ~ 0.1 pc behind an additional scattering screen.

We present a new algorithm for fitting and classifying polarized radio sources, which is based on the QU fitting method introduced by O'Sullivan et al. and on our analysis of pulsars. Then we test this algorithm using Monte Carlo simulations of observations in the 16 cm band of the Australia Telescope Compact Array (1.3-3.1 GHz), to quantify how often the algorithm identifies the correct source model, how certain it is of this identification, and how the parameters of the injected and fitted models compare. In our analysis we consider the Akaike and Bayesian Information Criteria, and model averaging. For the observing setup we simulated, the Bayesian Information Criterion, without model averaging, is the best way for identifying the correct model and for estimating its parameters. Sources can only be identified correctly if their parameters lie inside a 'Goldilocks region': strong depolarization makes it impossible to detect sources that emit over a wide range in RM, whereas sources that emit over a narrow range in RM cannot be told apart from simpler sources or sources that emit at only one RM. We identify when emission at similar RMs is 'resolved', and quantify this in a way similar to the Rayleigh criterion in optics. Also, we identify pitfalls in RM synthesis that are avoided by QU fitting. Finally, we show how channel weights can be tweaked to produce apodized RM spectra, that observing time requirements in RM synthesis and QU fitting are the same, and we analyse when to stop RMClean.

We develop two algorithms, based on maximum likelihood (ML) inference, for estimating the parameters of polarized radio sources which emit at a single rotation measure (RM), e.g., pulsars. These algorithms incorporate the flux density spectrum of the source, either a power law or a scaled version of the Stokes I spectrum, and a variation in sensitivity across the observing band. We quantify the detection significance and measurement uncertainties in the fitted parameters, and we derive weighted versions of the RM synthesis algorithm which, under certain conditions, maximize the likelihood. We use Monte Carlo simulations to compare injected and recovered source parameters for a range of signal-to-noise ratios, investigate the quality of standard methods for estimating measurement uncertainties, and search for statistical biases. These simulations consider one frequency band each for the Australia Telescope Compact Array (ATCA), the Square Kilometre Array (SKA), and the Low Frequency Array (LOFAR). We find that results obtained for one frequency band cannot be easily generalized, and that methods which were developed in the past for correcting bias in individual frequency channels do not apply to wide-band data sets. The standard method for estimating the measurement uncertainty in RM is not accurate for sources with non-zero spectral indices. Furthermore, dividing Stokes Q and U by Stokes I to correct for spectral index effects, in combination with RM synthesis, does not maximize the likelihood.

The Fan Region is one of the dominant features in the polarized radio sky, long thought to be a local (distance < 500 pc) synchrotron feature. We present 1.3-1.8 GHz polarized radio continuum observations of the region from the Global Magneto-Ionic Medium Survey (GMIMS) and compare them to maps of Halpha and polarized radio continuum intensity from 0.408-353 GHz. The high-frequency (> 1 GHz) and low-frequency (< 600 MHz) emission have different morphologies, suggesting a different physical origin. Portions of the 1.5 GHz Fan Region emission are depolarized by about 30% by ionized gas structures in the Perseus Arm, indicating that this fraction of the emission originates >2 kpc away. We argue for the same conclusion based on the high polarization fraction at 1.5 GHz (about 40%). The Fan Region is offset with respect to the Galactic plane, covering -5{\deg} < b < +10{\deg}; we attribute this offset to the warp in the outer Galaxy. We discuss origins of the polarized emission, including the spiral Galactic magnetic field. This idea is a plausible contributing factor although no model to date readily reproduces all of the observations. We conclude that models of the Galactic magnetic field should account for the > 1 GHz emission from the Fan Region as a Galactic-scale, not purely local, feature.

The new generation of low-frequency radio telescopes, such as the Low Frequency Array (LOFAR: a Square Kilometre Array-low pathfinder), provides advancements in our capability of probing Galactic magnetism through low-frequency polarimetry. Maps of diffuse polarized radio emission and Faraday rotation can be used to infer properties of, and trace structure in, the magnetic fields in the ISM. However, to date very little of the sky has been probed at high angular and Faraday depth resolution. We observed a 5x5 degree region centred on the nearby galaxy IC342 using LOFAR in the frequency range 115-178 MHz at 4 arcmin resolution and performed Faraday tomography to detect foreground Galactic polarized synchrotron emission separated by Faraday depth (different amounts of Faraday rotation). Our Faraday depth cube shows rich polarized structure, with up to 30 K of polarized emission at 150 MHz. We detect two overlapping diffuse polarized features that are clearly separated in Faraday depth. Faraday-thick structures at such low frequencies would be too strongly depolarized to explain the observations and are therefore rejected. Only Faraday thin structures will not be strongly depolarized; producing such structures requires localized variations in the ratio of synchrotron emissivity to Faraday depth per unit distance, which can arise from several physical phenomena, such as a transition between regions of ionized and neutral gas. We conclude that the observed polarized emission is Faraday thin, and propose that the emission originates from two neutral clouds in the local ISM. We have modeled the Faraday rotation for this line of sight and estimated that the line of sight component of magnetic field of the local ISM for this direction varies between -0.86 and +0.12 uG. We propose that this may be a useful method for mapping magnetic fields within the local ISM.

The magneto-ionic structures of the interstellar medium of the Milky Way and the intergalactic medium are still poorly understood, especially at distances larger than a few kiloparsecs from the Sun. The three-dimensional (3D) structure of the Galactic magnetic field and electron density distribution may be probed through observations of radio pulsars, primarily owing to their compact nature, high velocities, and highly-polarized short-duration radio pulses. Phase 1 of the SKA, i.e. SKA1, will increase the known pulsar population by an order of magnitude, and the full SKA, i.e. SKA2, will discover pulsars in the most distant regions of our Galaxy. SKA1-VLBI will produce model-independent distances to a large number of pulsars, and wide-band polarization observations by SKA1-LOW and SKA1-MID will yield high precision dispersion measure, scattering measure, and rotation measure estimates along thousands of lines of sight. When combined, these observations will enable detailed tomography of the large-scale magneto-ionic structure of both the Galactic disk and the Galactic halo. Turbulence in the interstellar medium can be studied through the variations of these observables and the dynamic spectra of pulsar flux densities. SKA1-LOW and SKA1-MID will monitor interstellar weather and produce sensitive dynamic and secondary spectra of pulsar scintillation, which can be used to make speckle images of the ISM, study turbulence on scales between ~10^8 and ~10^13 m, and probe pulsar emission regions on scales down to $\sim$10 km. In addition, extragalactic pulsars or fast radio bursts to be discovered by SKA1 and SKA2 can be used to probe the electron density distribution and magnetic fields in the intergalactic medium beyond the Milky Way.

To study the strength and structure of the magnetic field in the Galactic centre (GC) we measured Faraday rotation of the radio emission of pulsars which are seen towards the GC. Three of these pulsars have the largest rotation measures (RMs) observed in any Galactic object with the exception of Sgr A*. Their large dispersion measures, RMs and the large RM variation between these pulsars and other known objects in the GC implies that the pulsars lie in the GC and are not merely seen in projection towards the GC. The large RMs of these pulsars indicate large line-of-sight magnetic field components between ~ 16-33 microgauss; combined with recent model predictions for the strength of the magnetic field in the GC this implies that the large-scale magnetic field has a very small inclination angle with respect to the plane of the sky (~ 12 degrees). Foreground objects like the Radio Arc or possibly an ablated, ionized halo around the molecular cloud G0.11-0.11 could contribute to the large RMs of two of the pulsars. If these pulsars lie behind the Radio Arc or G0.11-0.11 then this proves that low-scattering corridors with lengths >~ 100 pc must exist in the GC. This also suggests that future, sensitive observations will be able to detect additional pulsars in the GC. Finally, we show that the GC component in our most accurate electron density model oversimplifies structure in the GC.

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.

We know that magnetic fields are pervasive across all scales in the Universe and over all of cosmic time and yet our understanding of many of the properties of magnetic fields is still limited. We do not yet know when, where or how the first magnetic fields in the Universe were formed, nor do we fully understand their role in fundamental processes such as galaxy formation or cosmic ray acceleration or how they influence the evolution of astrophysical objects. The greatest challenge to addressing these issues has been a lack of deep, broad bandwidth polarimetric data over large areas of the sky. The Square Kilometre Array will radically improve this situation via an all-sky polarisation survey that delivers both high quality polarisation imaging in combination with observations of 7-14 million extragalactic rotation measures. Here we summarise how this survey will improve our understanding of a range of astrophysical phenomena on scales from individual Galactic objects to the cosmic web.

We investigate how the imprint of Faraday rotation on radio spectra can be used to determine the geometry of radio sources and the strength and structure of the surrounding magnetic fields. We model spectra of Stokes Q and U for frequencies between 200 MHz and 10 GHz for Faraday screens with large-scale or small-scale magnetic fields external to the source. These sources can be uniform or 2D Gaussians on the sky with transverse linear gradients in rotation measure (RM), or cylinders or spheroids with an azimuthal magnetic field. At high frequencies the spectra of all these models can be approximated by the spectrum of a Gaussian source; this is independent of whether the magnetic field is large-scale or small-scale. A sinc spectrum in polarized flux density is not a unique signature of a volume where synchrotron emission and Faraday rotation are mixed. A turbulent Faraday screen with a large field coherence length produces a spectrum which is similar to the spectrum of a partial coverage model. At low and intermediate frequencies, such a Faraday screen produces a significantly higher polarized signal than the depolarization model by Burn, as shown by a random walk model of the polarization vectors. We calculate RM spectra for four frequency windows. Sources are strongly depolarized at low frequencies, but RMs can be determined accurately if the sensitivity of the observations is sufficient. Finally, we show that RM spectra can be used to differentiate between turbulent foreground models and partial coverage models.

The Gum Nebula is 36 degree wide shell-like emission nebula at a distance of only 450 pc. It has been hypothesised to be an old supernova remnant, fossil HII region, wind-blown bubble, or combination of multiple objects. Here we investigate the magneto-ionic properties of the nebula using data from recent surveys: radio-continuum data from the NRAO VLA and S-band Parkes All Sky Surveys, and H-alpha data from the Southern H-Alpha Sky Survey Atlas. We model the upper part of the nebula as a spherical shell of ionised gas expanding into the ambient medium. We perform a maximum-likelihood Markov chain Monte-Carlo fit to the NVSS rotation measure data, using the H-halpha data to constrain average electron density in the shell $n_e$. Assuming a latitudinal background gradient in RM we find $n_e=1.3^{+0.4}_{-0.4} {\rm cm}^{-3}$, angular radius $\phi_{\rm outer}=22.7^{+0.1}_{-0.1} {\rm deg}$, shell thickness $dr=18.5^{+1.5}_{-1.4} {\rm pc}$, ambient magnetic field strength $B_0=3.9^{+4.9}_{-2.2} \mu{\rm G}$ and warm gas filling factor $f=0.3^{+0.3}_{-0.1}$. We constrain the local, small-scale (~260 pc) pitch-angle of the ordered Galactic magnetic field to $+7^{\circ}\lesssim\wp\lesssim+44^{\circ}$, which represents a significant deviation from the median field orientation on kiloparsec scales (~-7.2$^{\circ}$). The moderate compression factor $X=6.0\,^{+5.1}_{-2.5}$ at the edge of the H-alpha shell implies that the 'old supernova remnant' origin is unlikely. Our results support a model of the nebula as a HII region around a wind-blown bubble. Analysis of depolarisation in 2.3 GHz S-PASS data is consistent with this hypothesis and our best-fitting values agree well with previous studies of interstellar bubbles.

(abridged) 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. We vary our assumptions about the data, showing that the dispersion of extragalactic contributions to observed Faraday depths is most likely lower than 7 rad/m^2, in agreement with earlier results, and that the extragalactic contribution to an individual data point is poorly constrained by the data in most cases.

Faraday rotation of polarised background sources is a unique probe of astrophysical magnetic fields in a diverse range of foreground objects. However, to understand the properties of the polarised sources themselves and of depolarising phenomena along the line of sight, we need to complement Faraday rotation data with polarisation observations over very broad bandwidths. Just as it is impossible to properly image a complex source with limited u-v coverage, we can only meaningfully understand the magneto-ionic properties of polarised sources if we have excellent coverage in $\lambda^2$-space. We here propose a set of broadband polarisation surveys with the Square Kilometre Array, which will provide a singular set of scientific insights on the ways in which galaxies and their environments have evolved over cosmic time.

Magnetic fields in the Milky Way are present on a wide variety of sizes and strengths, influencing many processes in the Galactic ecosystem such as star formation, gas dynamics, jets, and evolution of supernova remnants or pulsar wind nebulae. Observation methods are complex and indirect; the most used of these are a grid of rotation measures of unresolved polarized extragalactic sources, and broadband polarimetry of diffuse emission. Current studies of magnetic fields in the Milky Way reveal a global spiral magnetic field with a significant turbulent component; the limited sample of magnetic field measurements in discrete objects such as supernova remnants and HII regions shows a wide variety in field configurations; a few detections of magnetic fields in Young Stellar Object jets have been published; and the magnetic field structure in the Galactic Center is still under debate. The SKA will unravel the 3D structure and configurations of magnetic fields in the Milky Way on sub-parsec to galaxy scales, including field structure in the Galactic Center. The global configuration of the Milky Way disk magnetic field, probed through pulsar RMs, will resolve controversy about reversals in the Galactic plane. Characteristics of interstellar turbulence can be determined from the grid of background RMs. We expect to learn to understand magnetic field structures in protostellar jets, supernova remnants, and other discrete sources, due to the vast increase in sample sizes possible with the SKA. This knowledge of magnetic fields in the Milky Way will not only be crucial in understanding of the evolution and interaction of Galactic structures, but will also help to define and remove Galactic foregrounds for a multitude of extragalactic and cosmological studies.

We re-formulate rotation measure (RM) synthesis for data sets with discrete frequency channels and an arbitrary channel response function. The most commonly used version of the formalism by Brentjens & De Bruyn assumes a top-hat response function in wavelength squared, while real data sets can often be approximated better with a top-hat in frequency. We simulate mock data sets for various source geometries, using a top-hat response function in frequency, and we compare the quality of the RM spectra that are found with both formalisms. We include the response function of the simulated data to calculate exact RM spectra using our formalism. We show that the formalism by Brentjens & De Bruyn produces accurate results even if depolarization at the lowest frequencies in the observing band is severe. If RMs are large, our formalism reconstructs the emitted signal more accurately, with a higher amplitude and (in most cases) a narrower RM spread function. Our formalism can also detect sources with larger (absolute) RMs for a given sensitivity level of the observations.

We present polarisation properties at $1.4\,$GHz of two separate extragalactic source populations: passive quiescent galaxies and luminous quasar-like galaxies. We use data from the {\it Wide-Field Infrared Survey Explorer} data to determine the host galaxy population of the polarised extragalactic radio sources. The quiescent galaxies have higher percentage polarisation, smaller radio linear size, and $1.4\,$GHz luminosity of $6\times10^{21}<L_{\rm 1.4}<7\times10^{25}\,$W Hz$^{-1}$, while the quasar-like galaxies have smaller percentage polarisation, larger radio linear size at radio wavelengths, and a $1.4\,$GHz luminosity of $9\times10^{23}<L_{\rm 1.4}<7\times10^{28}\,$W Hz$^{-1}$, suggesting that the environment of the quasar-like galaxies is responsible for the lower percentage polarisation. Our results confirm previous studies that found an inverse correlation between percentage polarisation and total flux density at $1.4\,$GHz. We suggest that the population change between the polarised extragalactic radio sources is the origin of this inverse correlation and suggest a cosmic evolution of the space density of quiescent galaxies. Finally, we find that the extragalactic contributions to the rotation measures (RMs) of the nearby passive galaxies and the distant quasar-like galaxies are different. After accounting for the RM contributions by cosmological large-scale structure and intervening Mg\,{II} absorbers we show that the distribution of intrinsic RMs of the distant quasar-like sources is at most four times as wide as the RM distribution of the nearby quiescent galaxies, if the distribution of intrinsic RMs of the WISE-Star sources itself is at least several rad m$^{-2}$ wide.

We present a sample of 1317 Infrared-Faint Radio Sources (IFRSs) that, for the first time, are reliably detected in the infrared, generated by cross-correlating the Wide-Field Infrared Survey Explorer (WISE) all-sky survey with major radio surveys. Our IFRSs are brighter in both radio and infrared than the first generation IFRSs that were undetected in the infrared by the Spitzer Space Telescope. We present the first spectroscopic redshifts of IFRSs, and find that all but one of the IFRSs with spectroscopy has z > 2. We also report the first X-ray counterparts of IFRSs, and present an analysis of radio spectra and polarization, and show that they include Gigahertz-Peaked Spectrum, Compact Steep Spectrum, and Ultra-Steep Spectrum sources. These results, together with their WISE infrared colours and radio morphologies, imply that our sample of IFRSs represents a population of radio-loud Active Galactic Nuclei at z > 2. We conclude that our sample consists of lower-redshift counterparts of the extreme first generation IFRSs, suggesting that the fainter IFRSs are at even higher redshift.

A VLA Sky Survey of the extragalactic sky at S band (2-4 GHz) with polarization information can uniquely probe the magneto-ionic medium in a wide range of astrophysical environments over cosmic time. For a shallow all-sky survey, we expect to detect over 4 million sources in total intensity $>$ 0.45 mJy beam$^{-1}$ and over 2.2$\times$10$^5$ sources in polarized intensity. With these new observations, we expect to discover new classes of polarized radio sources in very turbulent astrophysical environments and those with extreme values of Faraday depth. Moreover, by determining reliable Faraday depths and by modeling depolarization effects, we can derive properties of the magneto-ionic medium associated with AGNs, absorption line systems and galaxies, addressing the following unresolved questions: (1) What is the covering fraction, the degree of turbulence and the origin of absorption line systems? (2) What is the thermal content in AGNs and radio galaxies? (3) How do AGNs and galaxies evolve over cosmic time? (4) What causes the increase in percentage polarization with decreasing flux densities at the low flux density end of the polarized source count? (5) What is the growth rate of large-scale magnetic fields in galaxies?

The centre of our Milky Way harbours the closest candidate for a supermassive black hole. The source is thought to be powered by radiatively inefficient accretion of gas from its environment. This form of accretion is a standard mode of energy supply for most galactic nuclei. X-ray measurements have already resolved a tenuous hot gas component from which it can be fed. However, the magnetization of the gas, a crucial parameter determining the structure of the accretion flow, remains unknown. Strong magnetic fields can influence the dynamics of the accretion, remove angular momentum from the infalling gas, expel matter through relativistic jets and lead to the observed synchrotron emission. Here we report multi-frequency measurements with several radio telescopes of a newly discovered pulsar close to the Galactic Centre and show that its unusually large Faraday rotation indicates a dynamically relevant magnetic field near the black hole. If this field is accreted down to the event horizon it provides enough magnetic flux to explain the observed emission from the black hole, from radio to X-rays.

This paper presents a detailed investigation of the dependence of pulsar spin-velocity alignment, which has been observed for a sample of 58 pulsars, on pulsar age. At first, our study considers only pulsar characteristic ages, resulting in no change in the degree of correlation as a function of age, up to at least 100 Myr. Subsequently, we consider a more reliable estimate of pulsar age, the kinematic age, assuming that pulsars are born near the Galactic plane. We derive kinematic ages for 52 pulsars, based on the measured pulsar proper motions and positions, by modelling the trajectory of the pulsars in a Galactic potential. The sample of 52 pulsar kinematic ages constitutes the largest number of independently estimated pulsar ages to date. Using only the 33 most reliable kinematic ages from our simulations, we revisit the evolution of spin- velocity alignment, this time as a function of kinematic age. We find that the strong correlation seen in young pulsars is completely smeared out for pulsars with kinematic ages above 10 Myr, a length of time beyond which we expect the gravitational pull of the Galaxy to have a significant effect on the directions of pulsar velocities. In the discussion, we investigate the impact of large distance uncertainties on the reliability of the calculated kinematic ages. Furthermore, we present a detailed investigation of the implications of our revised pulsar ages for the braking-index and birth-period distributions. Finally, we discuss the predictions of various SN-kick mechanisms and their compatibility with our results.

In this paper we test 8 models of the free electron distribution in the Milky Way that have been published previously, and we introduce 4 additional models that explore the parameter space of possible models further. These new models consist of a simple exponential thick disk model, and updated versions of the models by Taylor & Cordes and Cordes & Lazio with more extended thick disks. The final model we introduce uses the observed H-alpha intensity as a proxy for the total electron column density, also known as the dispersion measure (DM). We use the latest available data sets of pulsars with accurate distances (through parallax measurements or association with globular clusters) to optimise the parameters in these models. In the process of fitting a new scale height for the thick disk in the model by Cordes & Lazio we discuss why this thick disk cannot be replaced by the thick disk that Gaensler et al. advocated in a recent paper. In the second part of our paper we test how well the different models can predict the DMs of these pulsars at known distances. Almost all models perform well, in that they predict DMs within a factor of 1.5-2 of the observed DMs for about 75% of the lines of sight. This is somewhat surprising since the models we tested range from very simple models that only contain a single exponential thick disk to very complex models like the model by Cordes & Lazio. We show that the model by Taylor & Cordes that we updated with a more extended thick disk consistently performs better than the other models we tested. Finally, we analyse which sightlines have DMs that prove difficult to predict by most models, which indicates the presence of local features in the ISM between us and the pulsar. (abridged)

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 use 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 pi. 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 with a power law index of -2.14 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) proportional to k^-2.14 of the underlying field n_e times B_r under simplifying geometrical and statistical assumptions.

We present a detailed study of the Faraday depth structure of four bright (> 1 Jy), strongly polarized, unresolved, radio-loud quasars. The Australia Telescope Compact Array (ATCA) was used to observe these sources with 2 GHz of instantaneous bandwidth from 1.1 to 3.1 GHz. This allowed us to spectrally resolve the polarization structure of spatially unresolved radio sources, and by fitting various Faraday rotation models to the data, we conclusively demonstrate that two of the sources cannot be described by a simple rotation measure (RM) component modified by depolarization from a foreground Faraday screen. Our results have important implications for using background extragalactic radio sources as probes of the Galactic and intergalactic magneto-ionic media as we show how RM estimations from narrow-bandwidth observations can give erroneous results in the presence of multiple interfering Faraday components. We postulate that the additional RM components arise from polarized structure in the compact inner regions of the radio source itself and not from polarized emission from Galactic or intergalactic foreground regions. We further suggest that this may contribute significantly to any RM time-variability seen in RM studies on these angular scales. Follow-up, high-sensitivity VLBI observations of these sources will directly test our predictions.

Here we describe the Compact Array Broadband Backend (CABB) and present first results obtained with the upgraded Australia Telescope Compact Array (ATCA). The 16-fold increase in observing bandwidth, from 2 x 128 MHz to 2 x 2048 MHz, high bit sampling, and addition of 16 zoom windows (each divided into a further 2048 channels) provide major improvements for all ATCA observations. The benefits of the new system are: (1) hugely increased radio continuum and polarization sensitivity as well as image fidelity, (2) substantially improved capability to search for and map emission and absorption lines over large velocity ranges, (3) simultaneous multi-line and continuum observations, (4) increased sensitivity, survey speed and dynamic range due to high-bit sampling, and (5) high velocity resolution, while maintaining full polarization output. The new CABB system encourages all observers to make use of both spectral line and continuum data to achieve their full potential. Given the dramatic increase of the ATCA capabilities in all bands (ranging from 1.1 to 105 GHz) CABB enables scientific projects that were not feasible before the upgrade, such as simultaneous observations of multiple spectral lines, on-the-fly mapping, fast follow-up of radio transients (e.g., the radio afterglow of new supernovae) and maser observations at high velocity resolution and full polarization. The first science results presented here include wide-band spectra, high dynamic-range images, and polarization measurements, highlighting the increased capability and discovery potential of the ATCA.

In this Letter I use the variation of the spread in rotation measure (RM) with Galactic latitude to separate the Galactic from the extragalactic contributions to RM. This is possible since the latter does not depend on Galactic latitude. As input data I use RMs from the catalogue by Taylor, Stil, and Sunstrum, supplemented with published values for the spread in RM (`sigmaRM') in specific regions on the sky. I test 4 models of the free electron column density (which I will abbreviate to `DMinf') of the Milky Way, and the best model builds up DMinf on a characteristic scale of a few kpc from the Sun. sigmaRM correlates well with DMinf. The measured sigmaRM can be modelled as a Galactic contribution, consisting of a term sigmaRM,MW that is amplified at smaller Galactic latitudes as 1/sin|b|, in a similar way to DMinf, and an extragalactic contribution, sigmaRM,EG, that is independent of latitude. This model is sensitive to the relative magnitudes of sigmaRM,MW and sigmaRM,EG, and the best fit is produced by sigmaRM,MW approx. 8 rad/m^2 and sigmaRM,EG approx. 6 rad/m^2. The 4 published values for sigmaRM as a function of latitude suggest an even larger sigmaRM,MW contribution and a smaller sigmaRM,EG. This result from the NVSS RMs and published sigmaRM shows that the Galactic contribution dominates structure in RM on scales between about 1degr -- 10degr on the sky. I work out which factors contribute to the variation of sigmaRM with Galactic latitude, and show that the sigmaRM,EG I derived is an upper limit. Furthermore, to explain the modelled sigmaRM,MW requires that structure in <B||> has a 1-sigma spread <~ 0.4 microG.

Using a recent catalogue of extragalactic Faraday rotation derived from the NRAO VLA Sky Survey we have found an agreement between Faraday rotation structure and the HI emission structure of a High Velocity Cloud (HVC) associated with the Leading Arm of the Magellanic System. We suggest that this morphological agreement is indicative of Faraday rotation through the HVC. Under this assumption we have used 48 rotation measures through the HVC, together with estimates of the electron column density from H-\alpha\ measurements and QSO absorption lines to estimate a strength for the line-of-sight component of the coherent magnetic field in the HVC of <B_{||}> > 6 {\rm \mu G}$. A coherent magnetic field of this strength is more than sufficient to dynamically stabilize the cloud against ram pressure stripping by the Milky Way halo and may also provide thermal insulation for the cold cloud. We estimate an upper limit to the ratio of random to coherent magnetic field of $B_{r}/B_{||} < 0.8$, which suggests that the random field does not dominate over the coherent field as it does in the Magellanic Clouds from which this HVC likely originates.

We investigate the properties of the Galactic ISM by applying Faraday tomography to a radio polarization data set in the direction of the Galactic anti-centre. We address the problem of missing large-scale structure in our data, and show that this does not play an important role for the results we present. The main peak of the Faraday depth spectra in our data set is not measurably resolved for about 8% of the lines of sight. An unresolved peak indicates a separation between the regions with Faraday rotation and synchrotron emission. However, cosmic rays pervade the ISM, and synchrotron emission would therefore also be produced where there is Faraday rotation. We suggest that the orientation of the magnetic field can separate the two effects. By modelling the thermal electron contribution to the Faraday depth, we map the strength of the magnetic field component along the line of sight. Polarized point sources in our data set have rotation measures that are comparable to the Faraday depths of the diffuse emission in our data. Our Faraday depth maps show narrow canals of low polarized intensity. We conclude that depolarization over the telescope beam produces at least some of these canals. Finally, we investigate the properties of one conspicuous region in this data set and argue that it is created by a decrease in line-of-sight depolarization compared to its surroundings.

We investigate the distribution and properties of Faraday rotating and synchrotron emitting regions in the Galactic ISM in the direction of the Galactic anti-centre. We apply Faraday tomography to a radio polarization dataset that we obtained with the WSRT. We developed a new method to calculate a linear fit to periodic data, which we use to determine rotation measures from our polarization angle data. From simulations of a Faraday screen + noise we could determine how compatible the data are with Faraday screens. An unexpectedly large fraction of 14% of the lines-of-sight in our dataset show an unresolved main component in the Faraday depth spectrum. For lines-of-sight with a single unresolved component we demonstrate that a Faraday screen in front of a synchrotron emitting region that contains a turbulent magnetic field component can explain the data.

In this paper we use the diffuse Galactic synchrotron emission to study the properties of the Galactic magnetic field. We combined data from the low-frequency WENSS survey with single-dish observations carried out by Brouw and Spoelstra to (partially) fill in the central gap in the (u,v) plane that is missing from the interferometer data. The small bandwidth of the WENSS data meant that we could not determine rotation measures (RM) directly. Instead we interpreted the polarization angle gradients that we derived as gradients in RM, and to do this we found a new way to efficiently and reliably fit linear gradients to periodic data. RM are available for the single-dish observations, and we found that we have to rescale these single-dish RM to match the RM gradients that we derive from the combined WENSS/single-dish dataset with the single-dish only RM gradients. We tentatively show that the difference in beamsize between the datasets can be responsible for this. We interpret the scaled-up RM we find in terms of a simple toy model, and by combining our results with those by Haverkorn et al. (2004) we reconstruct the full 3D magnetic field vector along a number of lines-of-sight. For these lines-of-sight we derive the properties of the magnetic field component perpendicular to the Galactic plane.