D. H. F. M. Schnitzeler

A strong magnetic field around the supermassive black hole at the centre of the Galaxy

Earth’s nearest candidate supermassive black hole lies at the centre of the Milky Way. Its electromagnetic emission is thought to be powered by radiatively inefficient accretion of gas from its environment, which is a standard mode of energy supply for most galactic nuclei. X-ray measurements have already resolved a tenuous hot gas component from which the black hole can be fed. The magnetization of the gas, however, which is a crucial parameter determining the structure of the accretion flow, remains unknown. Strong magnetic fields can influence the dynamics of accretion, remove angular momentum from the infalling gas, expel matter through relativistic jets and lead to synchrotron emission such as that previously observed. Here we report multi-frequency radio measurements of a newly discovered pulsar close to the Galactic Centre and show that the pulsar’s unusually large Faraday rotation (the rotation of the plane of polarization of the emission in the presence of an external magnetic field) indicates that there is a dynamically important 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 radio to X-ray wavelengths—from the black hole.

An improved map of the galactic Faraday sky

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

Modelling the Galactic distribution of free electrons

An accurate picture of how free electrons are distributed throughout the Milky Way leads to more reliable distances for pulsars and more accurate maps of the magnetic field distribution in the Milky Way. In this paper we test eight models of the free electron distribution in the Milky Way that have been published previously, and we introduce four additional models that explore the parameter space of possible models further. These new models consist of a simple exponential thick-disc model, and updated versions of the models by Taylor & Cordes and Cordes & Lazio with more extended thick discs. The final model we introduce uses the observed Hα intensity as a proxy for the total electron column density, also known as the dispersion measure (DM). Since accurate maps of Hα intensity are now available, this final model can in theory outperform the other models. We use the latest available data sets of pulsars with accurate distances (through parallax measurements or association with globular clusters) to optimize the parameters in these models. In the process of fitting a new scale height for the thick disc in the model by Cordes & Lazio, we discuss why this thick disc cannot be replaced by the thick disc 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. We base our test on the ratios between the modelled and observed DMs, rather than on absolute deviations, and we identify systematic deviations between the modelled and observed DMs for the different models. For almost all models the ratio between the predicted and the observed DM cannot be described very well by a Gaussian distribution. We therefore calculate the deviations N between the modelled and observed DMs instead, and compare the cumulative distributions of N for the different models. Almost all models perform well, in that they predict DMs within a factor of 1.5–2 of the observed DMs for about 75 per cent 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 disc 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 disc 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 interstellar medium between us and the pulsar.

Estimating extragalactic Faraday rotation

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/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.

Pulsar spin-velocity alignment: Kinematic ages, birth periods and braking indices

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.

Infrared-faint radio sources: a new population of high-redshift radio galaxies

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.

Faraday tomography of the local interstellar medium with LOFAR: Galactic foregrounds towards IC 342

Magnetic fields pervade the interstellar medium (ISM), but are difficult to detect and characterize. 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 the associated 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 5° by 5° region centred on the nearby galaxy IC 342 (l = 138.2°,b = + 10.6°) using the LOFAR high-band antennae in the frequency range 115-178 MHz. We imaged this region at 4′.5x3′.84.5 × 3.8 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 a rich polarized structure, with up to 30 K of polarized emission at 150 MHz. We clearly detect two polarized features that extend over most of the field, but are clearly separated in Faraday depth. Simulations of the behaviour of the depolarization of Faraday-thick structures at such low frequencies show that such structures would be too strongly depolarized to explain the observations. These structures are therefore rejected as the source of the observed polarized features. Only Faraday thin structures will not be strongly depolarized at low frequencies; producing such structures requires localized variations in the ratio of synchrotron emissivity to Faraday depth per unit distance. Such variations can arise from several physical phenomena, such as a transition between regions of ionized and (mostly) neutral gas. We conclude that the observed polarized emission is Faraday thin, and propose that the emission originates from two mostly neutral clouds in the local ISM. Using maps of the local ISM to estimate distances to these clouds, we have modelled the Faraday rotation for this line of sight and estimated that the strength of the line of sight component of magnetic field of the local ISM for this direction varies between-0.86 and+0.12 μG (where positive is towards the Earth). We propose that this may be a useful method for mapping magnetic fields within the local ISM in all directions towards nearby neutral clouds.

Polarization signatures of unresolved radio sources

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.

Using SKA Rotation Measures to Reveal the Mysteries of the Magnetised Universe

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.

Rotation measure synthesis revisited

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.