Showing 1–25 of 25 results for author: Schnitzeler, D H F M

S-PASS/ATCA: a window on the magnetic universe in the southern hemisphere

Authors: D. H. F. M. Schnitzeler, E. Carretti, M. H. Wieringa, B. M. Gaensler, M. Haverkorn, S. Poppi

Abstract: 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 su… ▽ More 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. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 18 pages, 17 figures. The paper has been accepted for publication in MNRAS, and is in press. The catalogue can be downloaded by following this link: https://bit.ly/2E7jr3W

Diffuse polarized emission in the LOFAR Two-meter Sky Survey

Authors: C. L. Van Eck, M. Haverkorn, M. I. R. Alves, R. Beck, P. Best, E. Carretti, K. T. Chyży, T. Enßlin, J. S. Farnes, K. Ferrière, G. Heald, M. Iacobelli, V. Jelić, W. Reich, H. J. A. Röttgering, D. H. F. M. Schnitzeler

Abstract: 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… ▽ More 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 $μ$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. △ Less

Submitted 1 February, 2019; originally announced February 2019.

Comments: 23 pages, 16 figures, accepted by A&A. Figures have been degraded to comply with arXiv file size limits; full resolution figures will be available in the A&A version

Testing the accuracy of the ionospheric Faraday rotation corrections through LOFAR observations of bright northern pulsars

Authors: N. K. Porayko, A. Noutsos, C. Tiburzi, J. P. W. Verbiest, A. Horneffer, J. Künsemöller, S. Osłowski, M. Kramer, D. H. F. M. Schnitzeler, J. M. Anderson, M. Brüggen, J. -M. Grießmeier, M. Hoeft, D. J. Schwarz, M. Serylak, O. Wucknitz

Abstract: 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 a… ▽ More 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}$. △ Less

Submitted 4 December, 2018; originally announced December 2018.

Comments: 18 pages, 13 figures, 4 tables

Polarized point sources in the LOFAR Two-meter Sky Survey: A preliminary catalog

Authors: C. L. Van Eck, M. Haverkorn, M. I. R. Alves, R. Beck, P. Best, E. Carretti, K. T. Chyży, J. S. Farnes, K. Ferrière, M. J. Hardcastle, G. Heald, C. Horellou, M. Iacobelli, V. Jelić, D. D. Mulcahy, S. P. O'Sullivan, I. M. Polderman, W. Reich, C. J. Riseley, H. Röttgering, D. H. F. M. Schnitzeler, T. W. Shimwell, V. Vacca, J. Vink, G. J. White

Abstract: 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… ▽ More 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. △ Less

Submitted 13 January, 2018; originally announced January 2018.

Comments: 20 pages, 8 figures, 2 catalog tables (non-machine readable), accepted for publication in A&A

Journal ref: A&A 613, A58 (2018)

Finding a complex polarized signal in wide-band radio data

Authors: D. H. F. M. Schnitzeler

Abstract: 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 mo… ▽ More 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. △ Less

Submitted 27 October, 2017; originally announced October 2017.

Comments: 14 pages, 10 figures, 1 table. Accepted for publication in MNRAS

The Fan Region at 1.5 GHz. I: Polarized synchrotron emission extending beyond the Perseus Arm

Authors: A. S. Hill, T. L. Landecker, E. Carretti, K. Douglas, X. H. Sun, B. M. Gaensler, S. A. Mao, N. M. McClure-Griffiths, W. Reich, M. Wolleben, J. M. Dickey, A. D. Gray, M. Haverkorn, J. P. Leahy, D. H. F. M. Schnitzeler

Abstract: 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 l… ▽ More 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° < b < +10°; 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. △ Less

Submitted 7 February, 2017; originally announced February 2017.

Comments: Accepted to MNRAS

Journal ref: MNRAS 467 4631 (2017)

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

Authors: C. L. Van Eck, M. Haverkorn, M. I. R. Alves, R. Beck, A. G. de Bruyn, T. Enßlin, J. S. Farnes, K. Ferrière, G. Heald, C. Horellou, A. Horneffer, M. Iacobelli, V. Jelić, I. Martí-Vidal, D. D. Mulcahy, W. Reich, H. J. A. Röttgering, A. M. M Scaife, D. H. F. M. Schnitzeler, C. Sobey, S. S. Sridhar

Abstract: 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. Howev… ▽ More 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. △ Less

Submitted 2 December, 2016; originally announced December 2016.

Comments: 17 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 597, A98 (2017)

Finding a faint polarized signal in wide-band radio data

Authors: D. H. F. M. Schnitzeler, K. J. Lee

Abstract: 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 de… ▽ More 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. △ Less

Submitted 23 October, 2017; v1 submitted 29 November, 2016; originally announced November 2016.

Comments: This is a resubmission of the original manuscript. The covariance matrix that we calculated in equation (17) was not correct, and several of the η in section 3.1 should be replaced with \hat η. We fixed these mistakes in an Erratum that has been accepted for publication in MNRAS, and in this new version of the manuscript that we uploaded to arXiv

Journal ref: Volume 466, 2017, Issue 1, p.378-391

Radio polarimetry of Galactic centre pulsars

Authors: D. H. F. M. Schnitzeler, R. P. Eatough, K. Ferrière, M. Kramer, K. J. Lee, A. Noutsos, R. M. Shannon

Abstract: 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… ▽ More 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. △ Less

Submitted 18 April, 2016; originally announced April 2016.

Comments: 7 pages, 3 figures. Accepted for publication in MNRAS

Tomography of the Galactic free electron density with the Square Kilometer Array

Authors: Maksim Greiner, Dominic H. F. M. Schnitzeler, Torsten A. Ensslin

Abstract: 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… ▽ More 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. △ Less

Submitted 10 December, 2015; originally announced December 2015.

Comments: 18 pages, 16 figures

Journal ref: A&A 590, A59 (2016)

Polarization signatures of unresolved radio sources

Authors: D. H. F. M. Schnitzeler, J. K. Banfield, K. J. Lee

Abstract: 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 Gauss… ▽ More 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. △ Less

Submitted 31 March, 2015; originally announced April 2015.

Comments: 20 pages, 15 figures. Accepted for publication in MNRAS

A radio-polarisation and rotation measure study of the Gum Nebula and its environment

Authors: C. R. Purcell, B. M. Gaensler, X. H. Sun, E. Carretti, G. Bernardi, M. Haverkorn, M. J. Kesteven, S. Poppi, D. H. F. M. Schnitzeler, L. Staveley-Smith

Abstract: 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-alph… ▽ More 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 $φ_{\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} μ{\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. △ Less

Submitted 22 February, 2015; originally announced February 2015.

Comments: 33 pages, 16 figures. Accepted by The Astrophysical Journal

Rotation measure synthesis revisited

Authors: D. H. F. M. Schnitzeler, K. J. Lee

Abstract: 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 g… ▽ More 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. △ Less

Submitted 5 November, 2014; originally announced November 2014.

Comments: Accepted for publication in MNRAS. 6 pages, 4 figures

Estimating extragalactic Faraday rotation

Authors: Niels Oppermann, Henrik Junklewitz, Maksim Greiner, Torsten A. Enßlin, Takuya Akahori, Ettore Carretti, Bryan M. Gaensler, Ariel Goobar, Lisa Harvey-Smith, Melanie Johnston-Hollitt, Luke Pratley, Dominic H. F. M. Schnitzeler, Jeroen M. Stil, Valentina Vacca

Abstract: (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 recon… ▽ More (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. △ Less

Submitted 16 February, 2015; v1 submitted 14 April, 2014; originally announced April 2014.

Comments: 20 + 6 pages, 19 figures; minor changes after bug-fix; version accepted for publication by A&A; results are available at http://www.mpa-garching.mpg.de/ift/faraday/

Journal ref: A&A 575, A118 (2015)

Radio galaxies and their magnetic fields out to z <= 3

Authors: J. K. Banfield, D. H. F. M. Schnitzeler, S. J. George, R. P. Norris, T. H. Jarrett, A. R. Taylor, J. M. Stil

Abstract: 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… ▽ More 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. △ Less

Submitted 16 July, 2014; v1 submitted 6 April, 2014; originally announced April 2014.

Comments: 12 pages, 8 figures, accepted for publication into MNRAS

Infrared-Faint Radio Sources: A New Population of High-redshift Radio Galaxies

Authors: Jordan D. Collier, Julie K. Banfield, Ray P. Norris, Dominic H. F. M. Schnitzeler, Amy E. Kimball, Miroslav D. Filipović, Tom H. Jarrett, Carol J. Lonsdale, Nick F. H. Tothill

Abstract: 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.… ▽ More 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. △ Less

Submitted 24 February, 2014; v1 submitted 3 December, 2013; originally announced December 2013.

Comments: 23 pages, 17 figures. Submitted to MNRAS

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

Authors: R. P. Eatough, H. Falcke, R. Karuppusamy, K. J. Lee, D. J. Champion, E. F. Keane, G. Desvignes, D. H. F. M. Schnitzeler, L. G. Spitler, M. Kramer, B. Klein, C. Bassa, G. C. Bower, A. Brunthaler, I. Cognard, A. T. Deller, P. B. Demorest, P. C. C. Freire, A. Kraus, A. G. Lyne, A. Noutsos, B. Stappers, N. Wex

Abstract: 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 o… ▽ More 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. △ Less

Submitted 14 August, 2013; originally announced August 2013.

Comments: Advance online publication in Nature, 14/08/13

Modelling the Galactic distribution of free electrons

Authors: D. H. F. M. Schnitzeler

Abstract: 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 m… ▽ More 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) △ Less

Submitted 15 August, 2012; originally announced August 2012.

Comments: 16 pages, 10 figures, 5 tables. Accepted for publication in the Monthly Notices of the RAS by the Royal Astronomical Society and Blackwell Publishing

Complex Faraday depth structure of Active Galactic Nuclei as revealed by broadband radio polarimetry

Authors: S. P. O'Sullivan, S. Brown, T. Robishaw, D. H. F. M. Schnitzeler, N. M. McClure-Griffiths, I. J. Feain, A. R. Taylor, B. M. Gaensler, T. L. Landecker, L. Harvey-Smith, E. Carretti

Abstract: 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… ▽ More 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. △ Less

Submitted 10 February, 2012; v1 submitted 16 January, 2012; originally announced January 2012.

Comments: Accepted for publication in MNRAS. 19 pages, 15 figures, 4 tables

An improved map of the Galactic Faraday sky

Authors: N. Oppermann, H. Junklewitz, G. Robbers, M. R. Bell, T. A. Enßlin, A. Bonafede, R. Braun, J. C. Brown, T. E. Clarke, I. J. Feain, B. M. Gaensler, A. Hammond, L. Harvey-Smith, G. Heald, M. Johnston-Hollitt, U. Klein, P. P. Kronberg, S. A. Mao, N. M. McClure-Griffiths, S. P. O'Sullivan, L. Pratley, T. Robishaw, S. Roy, D. H. F. M. Schnitzeler, C. Sotomayor-Beltran , et al. (6 additional authors not shown)

Abstract: 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 statist… ▽ More 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. △ Less

Submitted 5 March, 2012; v1 submitted 26 November, 2011; originally announced November 2011.

Comments: 16 pages, 11 figures. Update in one data catalog. All results are available at http://www.mpa-garching.mpg.de/ift/faraday/

The latitude dependence of the rotation measures of NVSS sources

Authors: D. H. F. M. Schnitzeler

Abstract: 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… ▽ More 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. △ Less

Submitted 2 November, 2010; originally announced November 2010.

Comments: 6 pages, 3 figures, 1 table. Published in MNRAS Letters; the definitive version is available at wileyonlinelibrary.com, http://onlinelibrary.wiley.com/doi/10.1111/j.1745-3933.2010.00957.x/pdf

Journal ref: Monthly Notices Letters - 2010 volume: 409 Issue: 1 Page: L97

Measurement of a Magnetic Field in a Leading Arm High Velocity Cloud

Authors: N. M. McClure-Griffiths, G. J. Madsen, B. M. Gaensler, D. McConnell, D. H. F. M. Schnitzeler

Abstract: 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 ha… ▽ More 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-α 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 μ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. △ Less

Submitted 29 September, 2010; originally announced September 2010.

Comments: 17 pages, 3 figures

WSRT Faraday tomography of the Galactic ISM at λ~ 0.86 m

Authors: D. H. F. M. Schnitzeler, P. Katgert, A. G. de Bruyn

Abstract: 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… ▽ More 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. △ Less

Submitted 23 October, 2008; originally announced October 2008.

Comments: 12 pages, 11 figures, 2 tables. Accepted for publication in A&A

WSRT Faraday tomography of the Galactic ISM at λ\sim 0.86 m

Authors: D. H. F. M. Schnitzeler, P. Katgert, A. G. de Bruyn

Abstract: 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 angl… ▽ More 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. △ Less

Submitted 18 June, 2007; originally announced June 2007.

Comments: 5 pages, 5 figures. Accepted for publication as a Letter to the Editor in A&A

The WENSS and Dwingeloo surveys and the Galactic magnetic field

Authors: D. H. F. M. Schnitzeler, P. Katgert, M. Haverkorn, A. G. de Bruyn

Abstract: 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… ▽ More 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. △ Less

Submitted 19 October, 2006; originally announced October 2006.

Comments: accepted for publication in A&A; 16 pages, 9 figures, 2 tables, 1 online figure