Showing 1–50 of 78 results for author: Jelić, V

Faraday tomography of LoTSS-DR2 data: III. Revealing the Local Bubble and the complex of local interstellar clouds in the high-latitude inner Galaxy

Authors: Ana Erceg, Vibor Jelić, Marijke Haverkorn, Lovorka Gajović, Martin Hardcastle, Timothy W. Shimwell, Cyril Tasse

Abstract: In this work, we present a new mosaic created with the second release of LOFAR Two-Metre Sky Survey data (LoTSS-DR2), which probes polarised synchrotron emission in the high-latitude inner Galaxy. Our objective is to characterise the observed emission through multi-tracer analysis to better understand the volume and the structures that may be observed with LOFAR. Furthermore, we exploit Faraday de… ▽ More In this work, we present a new mosaic created with the second release of LOFAR Two-Metre Sky Survey data (LoTSS-DR2), which probes polarised synchrotron emission in the high-latitude inner Galaxy. Our objective is to characterise the observed emission through multi-tracer analysis to better understand the volume and the structures that may be observed with LOFAR. Furthermore, we exploit Faraday depth as a unique tool to probe the diffuse magnetised structure in the local ISM. We produced a mosaic Faraday cube of LoTSS-DR2 data by applying a rotation measure synthesis algorithm. From the cube, we constructed Faraday moment maps to characterise the nature of spectra. Additionally, we quantified the linear depolarisation canals using the Rolling Hough transform and used them to search for alignment with other data sets. Utilising LoTSS-DR2 observations alongside complementary data sets including Planck polarisation data, HI emission maps, and starlight polarisation measurements, we estimated the distance to the Faraday structures. The Faraday cube reveals an ordered structure across two-thirds of the observed area, whose orientation aligns well with that of both the HI filaments and the magnetic field. We estimate the minimum distance to the Faraday structures to be between 40 and 80 pc, which puts them in the vicinity of the Local Bubble wall. The emission is organised in a large gradient in Faraday depth whose origin we associate with the curved wall of the Local Bubble. Comparing our data with a model of the Local Bubble wall, we conclude that we might be probing a contribution of the medium inside the Local Bubble cavity as well, corresponding to the complex of local interstellar clouds. Moreover, we propose a toy model incorporating an ionised front of finite thickness into the Local Bubble wall, as a curved, cold neutral shell alone is insufficient to produce the observed gradient. △ Less

Submitted 9 July, 2024; v1 submitted 20 June, 2024; originally announced June 2024.

Comments: 16 pages, 26 figures, accepted for publication in A&A

LOFAR Deep Fields: Probing the sub-mJy regime of polarized extragalactic sources in ELAIS-N1. I. The catalog

Authors: S. Piras, C. Horellou, J. E. Conway, M. Thomasson, S. del Palacio, T. W. Shimwell, S. P. O'Sullivan, E. Carretti, I. Šnidaric, V. Jelic, B. Adebahr, A. Berger, P. N. Best, M. Brüggen, N. Herrera Ruiz, R. Paladino, I. Prandoni, J. Sabater, V. Vacca

Abstract: The aim of this study is to probe the sub-mJy polarized source population with LOFAR. We present the method used to stack LOFAR polarization datasets, the resulting catalog of polarized sources, and the derived polarized source counts. The ELAIS-N1 field was selected for a polarimetric study at 114.9-177.4 MHz. A total area of 25 deg2 was imaged at 6"- resolution in the Stokes Q and U parameters.… ▽ More The aim of this study is to probe the sub-mJy polarized source population with LOFAR. We present the method used to stack LOFAR polarization datasets, the resulting catalog of polarized sources, and the derived polarized source counts. The ELAIS-N1 field was selected for a polarimetric study at 114.9-177.4 MHz. A total area of 25 deg2 was imaged at 6"- resolution in the Stokes Q and U parameters. Alignment of polarization angles was done both in frequency and in Faraday space before stacking datasets from 19 eight-hour-long epochs. A search for polarized sources was carried out in the final, stacked dataset, and the properties of the detected sources were examined. The depolarization level of sources known to be polarized at 1.4 GHz was quantified. A one-sigma noise level of 19 μJy/beam was reached in the central part of the field after stacking. Twenty-five polarized sources were detected above 8σ, five of which had not been detected in polarization at any other radio frequencies before. Seven additional polarized components were found by lowering the threshold to 6σat positions corresponding to sources known to be polarized at 1.4 GHz. In two radio galaxies, polarization was detected from both radio lobes, so the final number of associated radio continuum sources is 31. The detected sources are weakly polarized, with a median degree of polarization of 1.75% for the sample of sources detected in polarized emission. The sources previously detected in polarization at 1.4 GHz are significantly depolarized at 150 MHz. The catalog is used to derive the polarized source counts at 150 MHz. This is the deepest and highest-resolution polarization study at 150 MHz to date. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 18 pages, 10 figures, accepted for publication in A&A

Faraday tomography of LoTSS-DR2 data: II. Multi-tracer analysis in the high-latitude outer Galaxy

Authors: Ana Erceg, Vibor Jelić, Marijke Haverkorn, Andrea Bracco, Lana Ceraj, Luka Turić, Juan D. Soler

Abstract: We conducted a follow-up study on the analysis of the LOw Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) mosaic in the high-latitude outer Galaxy presented in the first paper of this series. Here, we focus on the search for alignment between the magnetic field traced by dust, HI filaments, starlight optical linear polarisation, and linear depolarised structures (depolarisation canals) observ… ▽ More We conducted a follow-up study on the analysis of the LOw Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) mosaic in the high-latitude outer Galaxy presented in the first paper of this series. Here, we focus on the search for alignment between the magnetic field traced by dust, HI filaments, starlight optical linear polarisation, and linear depolarised structures (depolarisation canals) observed in low-frequency synchrotron polarisation. This alignment was previously found in several smaller fields observed with LOFAR, offering valuable insights into the nature of the interstellar medium and the 3D spatial distribution of the diffuse ionised medium. We aim to determine whether the alignment of the interstellar medium (ISM) phases observed through multiple tracers is a common occurrence or an exception. Additionally, in areas where depolarisation canals align with the magnetic field, we use starlight polarisation to constrain the distance to the structures associated with the observed canals. We employed the Rolling Hough Transform (RHT) and projected Rayleigh statistics (PRS) to identify and quantify the alignment between the different tracers. On the scale of the whole mosaic, we did not find any evidence of a universal alignment among the three tracers. However, in one particular area, the western region (Dec between $29^\circ$ and $70^\circ$ and RA between $\mathrm{7^h44^m}$ and $\mathrm{9^h20^m}$), we do find a significant alignment between the magnetic field, depolarisation canals, and HI filaments. Based on this alignment, we used the starlight polarisation of stars with known parallax distances to estimate that the minimum distance to the structures observed by LOFAR in this region lies within the range of 200 to 240 pc. We associate these structures with the edge of the Local Bubble. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 11 pages, 10 figures; accepted for publication in A&A; data products are available at https://data.fulir.irb.hr/islandora/object/irb:367

Journal ref: A&A, 687, A23 (2024)

Associating LOFAR Galactic Faraday structures with the warm neutral medium

Authors: F. Boulanger, C. Gry, E. B. Jenkins, A. Bracco, A. Erceg, V. Jelić, L. Turić

Abstract: Faraday tomography observations with the Low Frequency Array (LOFAR) have unveiled a remarkable network of structures in polarized synchrotron emission at high Galactic latitudes. The observed correlation between LOFAR structures, dust polarization, and HI emission suggests a connection to the neutral interstellar medium (ISM). We investigated this relationship by estimating the rotation measure (… ▽ More Faraday tomography observations with the Low Frequency Array (LOFAR) have unveiled a remarkable network of structures in polarized synchrotron emission at high Galactic latitudes. The observed correlation between LOFAR structures, dust polarization, and HI emission suggests a connection to the neutral interstellar medium (ISM). We investigated this relationship by estimating the rotation measure (RM) of the warm neutral (partially ionized) medium (WNM) in the local ISM. Our work combines UV spectroscopy from FUSE and dust polarization observations from \Planck\ with LOFAR data. We derived electron column densities from UV absorption spectra toward nine background stars, within the field of published data from the LOFAR two-meter sky survey. The associated RMs were estimated using a local magnetic field model fitted to the dust polarization data of \Planck. A comparison with Faraday spectra at the position of the stars suggests that LOFAR structures delineate a slab of magnetized WNM and synchrotron emission, located ahead of the bulk of the warm ionized medium. This conclusion establishes an astrophysical framework for exploring the link between Faraday structures and the dynamics of the magnetized multiphase ISM. It will be possible to test it on a larger sample of stars when maps from the full northern sky survey of LOFAR become available. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A

Journal ref: A&A 687, A102 (2024)

Comparison of distance measurements to dust clouds using GRB X-ray halos and 3D dust extinction

Authors: Barbara Šiljeg, Željka Bošnjak, Vibor Jelić, Andrea Tiengo, Fabio Pintore, Andrea Bracco

Abstract: X-ray photons from energetic sources such as gamma-ray bursts (GRBs) can be scattered on dust clouds in the Milky Way, creating a time-evolving halo around the GRB position. X-ray observations of such halos allow the measurement of dust clouds distances in the Galaxy on which the scattering occurs. We present the first systematic comparison of the distances to scattering regions derived from GRB h… ▽ More X-ray photons from energetic sources such as gamma-ray bursts (GRBs) can be scattered on dust clouds in the Milky Way, creating a time-evolving halo around the GRB position. X-ray observations of such halos allow the measurement of dust clouds distances in the Galaxy on which the scattering occurs. We present the first systematic comparison of the distances to scattering regions derived from GRB halos with the 3D dust distribution derived from recently published optical-to-near infrared extinction maps. GRB halos were observed around 7 sources by the Swift XRT and the XMM-Newton EPIC instruments, namely GRB 031203, GRB 050713A, GRB 050724, GRB 061019, GRB 070129, GRB 160623A and GRB 221009A. We used four 3D extinction maps that exploit photometric data from different surveys and apply diverse algorithms for the 3D mapping of extinction, and compared the X-ray halo-derived distances with the local maxima in the 3D extinction density distribution. We found that in all GRBs we can find at least one local maximum in the 3D dust extinction map that is in agreement with the dust distance measured from X-ray rings. For GRBs with multiple X-ray rings, the dust distance measurements coincide with at least 3 maxima in the extinction map for GRB 160623A, and 5 maxima for GRB 221009A. The agreement of these independent distance measurements shows that the methods used to create dust extinction maps may potentially be optimized by the X-ray halo observations from GRBs. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 15 pages, 10 figures, accepted for publication in MNRAS

LOFAR Deep Fields: Probing faint Galactic polarised emission in ELAIS-N1

Authors: Iva Šnidarić, Vibor Jelić, Maaijke Mevius, Michiel Brentjens, Ana Erceg, Timothy W. Shimwell, Sara Piras, Cathy Horellou, Jose Sabater, Philip N. Best, Andrea Bracco, Lana Ceraj, Marijke Haverkorn, Shane P. O'Sullivan, Luka Turić, Valentina Vacca

Abstract: We present the first deep polarimetric study of Galactic synchrotron emission at low radio frequencies. Our study is based on 21 observations of the European Large Area Infrared Space Observatory Survey-North 1 (ELAIS-N1) field using the Low-Frequency Array (LOFAR) at frequencies from 114.9 to 177.4 MHz. These data are a part of the LOFAR Two-metre Sky Survey Deep Fields Data Release 1. We used ve… ▽ More We present the first deep polarimetric study of Galactic synchrotron emission at low radio frequencies. Our study is based on 21 observations of the European Large Area Infrared Space Observatory Survey-North 1 (ELAIS-N1) field using the Low-Frequency Array (LOFAR) at frequencies from 114.9 to 177.4 MHz. These data are a part of the LOFAR Two-metre Sky Survey Deep Fields Data Release 1. We used very low-resolution ($4.3'$) Stokes QU data cubes of this release. We applied rotation measure (RM) synthesis to decompose the distribution of polarised structures in Faraday depth, and cross-correlation RM synthesis to align different observations in Faraday depth. We stacked images of about 150 hours of the ELAIS-N1 observations to produce the deepest Faraday cube at low radio frequencies to date, tailored to studies of Galactic synchrotron emission and the intervening magneto-ionic interstellar medium. This Faraday cube covers $\sim36~{\rm deg^{2}}$ of the sky and has a noise of $27~{\rm μJy~PSF^{-1}~RMSF^{-1}}$ in polarised intensity. This is an improvement in noise by a factor of approximately the square root of the number of stacked data cubes ($\sim\sqrt{20}$), as expected, compared to the one in a single data cube based on five-to-eight-hour observations. We detect a faint component of diffuse polarised emission in the stacked cube, which was not detected previously. Additionally, we verify the reliability of the ionospheric Faraday rotation corrections estimated from the satellite-based total electron content measurements to be of $~\sim0.05~{\rm rad~m^{-2}}$. We also demonstrate that diffuse polarised emission itself can be used to account for the relative ionospheric Faraday rotation corrections with respect to a reference observation. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 15 pages, 15 figures, accepted for publication in A&A

Journal ref: A&A 674, A119 (2023)

The power of the rings: the GRB 221009A soft X-ray emission from its dust-scattering halo

Authors: Andrea Tiengo, Fabio Pintore, Beatrice Vaia, Simone Filippi, Andrea Sacchi, Paolo Esposito, Michela Rigoselli, Sandro Mereghetti, Ruben Salvaterra, Barbara Siljeg, Andrea Bracco, Zeljka Bosnjak, Vibor Jelic, Sergio Campana

Abstract: GRB 221009A is the brightest gamma-ray burst (GRB) ever detected and occurred at low Galactic latitude. Owing to this exceptional combination, its prompt X-ray emission could be detected for weeks in the form of expanding X-ray rings produced by scattering in Galactic dust clouds. We report on the analysis of 20 rings, generated by dust at distances ranging from 0.3 to 18.6 kpc, detected during tw… ▽ More GRB 221009A is the brightest gamma-ray burst (GRB) ever detected and occurred at low Galactic latitude. Owing to this exceptional combination, its prompt X-ray emission could be detected for weeks in the form of expanding X-ray rings produced by scattering in Galactic dust clouds. We report on the analysis of 20 rings, generated by dust at distances ranging from 0.3 to 18.6 kpc, detected during two XMM-Newton observations performed about 2 and 5 days after the GRB. By fitting the spectra of the rings with different models for the dust composition and grain size distribution, we reconstructed the spectrum of the GRB prompt emission in the 0.7-4 keV energy range as an absorbed power law with photon index 1-1.4 and absorption in the host galaxy nHz=(4.1-5.3)E21 cm-2. Taking into account the systematic uncertainties on the column density of dust contained in the clouds producing the rings, the 0.5-5 keV fluence of GRB 221009A can be constrained between 1E-3 and 7E-3 erg cm-2. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 14 pages, 7 figures, submitted to ApJL

Journal ref: The Astrophysical Journal Letters, 946, L30, 11 pp. (2023)

Structure in the Magnetic Field of the Milky Way Disk and Halo traced by Faraday Rotation

Authors: John M. Dickey, Jennifer West, Alec J. M. Thomson, T. L. Landecker, A. Bracco, E. Carretti, J. L. Han, A. S. Hill, Y. K. Ma, S. A. Mao, A. Ordog, Jo-Anne C. Brown, K. A. Douglas, A. Erceg, V. Jelic, R. Kothes, M. Wolleben

Abstract: Magnetic fields in the ionized medium of the disk and halo of the Milky Way impose Faraday rotation on linearly polarized radio emission. We compare two surveys mapping the Galactic Faraday rotation, one showing the rotation measures of extragalactic sources seen through the Galaxy (from Hutschenreuter et al 2022), and one showing the Faraday depth of the diffuse Galactic synchrotron emission from… ▽ More Magnetic fields in the ionized medium of the disk and halo of the Milky Way impose Faraday rotation on linearly polarized radio emission. We compare two surveys mapping the Galactic Faraday rotation, one showing the rotation measures of extragalactic sources seen through the Galaxy (from Hutschenreuter et al 2022), and one showing the Faraday depth of the diffuse Galactic synchrotron emission from the Global Magneto-Ionic Medium Survey. Comparing the two data sets in 5deg x 10deg bins shows good agreement at intermediate latitudes, 10 < |b| < 50 deg, and little correlation between them at lower and higher latitudes. Where they agree, both tracers show clear patterns as a function of Galactic longitude: in the Northern Hemisphere a strong sin(2 x longitude) pattern, and in the Southern hemisphere a sin(longitude + pi) pattern. Pulsars with height above or below the plane |z| > 300 pc show similar longitude dependence in their rotation measures. Nearby non-thermal structures show rotation measure shadows as does the Orion-Eridanus superbubble. We describe families of dynamo models that could explain the observed patterns in the two hemispheres. We suggest that a field reversal, known to cross the plane a few hundred pc inside the solar circle, could shift to positive z with increasing Galactic radius to explain the sin(2xlongitude) pattern in the Northern Hemisphere. Correlation shows that rotation measures from extragalactic sources are one to two times the corresponding rotation measure of the diffuse emission, implying Faraday complexity along some lines of sight, especially in the Southern hemisphere. △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: 37 pages, 26 figures, Ap. J. accepted

First look at the multiphase interstellar medium with synthetic observations of low-frequency Faraday tomography

Authors: A. Bracco, E. Ntormousi, V. Jelić, M. Padovani, B. Šiljeg, A. Erceg, L. Turić, L. Ceraj, I. Šnidarić

Abstract: Faraday tomography of radio polarimetric data below 200 MHz from LOFAR are providing us with a new perspective on the diffuse and magnetized interstellar medium (ISM). Of particular interest is the discovery of Faraday-rotated synchrotron polarization associated with neutral gas, as traced by atomic hydrogen (HI) and dust. Here we present the first in-depth numerical study of these LOFAR results.… ▽ More Faraday tomography of radio polarimetric data below 200 MHz from LOFAR are providing us with a new perspective on the diffuse and magnetized interstellar medium (ISM). Of particular interest is the discovery of Faraday-rotated synchrotron polarization associated with neutral gas, as traced by atomic hydrogen (HI) and dust. Here we present the first in-depth numerical study of these LOFAR results. We produce and analyze comprehensive synthetic observations of low-frequency synchrotron polarization from MHD simulations of colliding super shells in the multiphase ISM. We define five distinct gas phases over more than four orders of magnitude in gas temperature and density, ranging from hot, and warm, fully ionized gas to cold neutral medium. We focus on the contribution of each gas phase to synthetic observations of both rotation measure and synchrotron polarized intensity below 200 MHz. We also investigate the link between the latter and synthetic observations of optically thin HI gas. We find that, not only the fully ionized gas but also the warm partially ionized and neutral phases strongly contribute to the total rotation measure and polarized intensity. However, the contribution of each phase to the observables strongly depends on the choice of integration axis and the orientation of the mean magnetic field with respect to the shell-collision axis. Strong correlation between HI synthetic data and synchrotron polarized intensity, reminiscent of LOFAR results, is obtained with lines of sight perpendicular to the mean magnetic field direction. Our study suggests that multiphase modelling of MHD processes is needed in order to interpret observations of the radio sky at low frequency. This work is a first step toward understanding the complexity of low-frequency synchrotron emission that will be soon revolutionized by large-scale surveys with LOFAR and the SKA. △ Less

Submitted 15 April, 2022; v1 submitted 6 April, 2022; originally announced April 2022.

Comments: 17 pages, 20 figures, accepted by A&A, in press

Journal ref: A&A 663, A37 (2022)

Faraday tomography of LoTSS-DR2 data: I. Faraday moments in the high-latitude outer Galaxy and revealing Loop III in polarisation

Authors: Ana Erceg, Vibor Jelić, Marijke Haverkorn, Andrea Bracco, Timothy W. Shimwell, Cyril Tasse, John M. Dickey, Lana Ceraj, Alexander Drabent, Martin J. Hardcastle, Luka Turić

Abstract: Observations of synchrotron emission at low radio frequencies reveal a labyrinth of polarised Galactic structures. However, the explanation for the wealth of structures remains uncertain due to the complex interactions between the interstellar medium and the magnetic field. A multi-tracer approach to the analysis of large sky areas is needed. This paper aims to use polarimetric images from the LOF… ▽ More Observations of synchrotron emission at low radio frequencies reveal a labyrinth of polarised Galactic structures. However, the explanation for the wealth of structures remains uncertain due to the complex interactions between the interstellar medium and the magnetic field. A multi-tracer approach to the analysis of large sky areas is needed. This paper aims to use polarimetric images from the LOFAR Two metre Sky Survey (LoTSS) to produce the biggest mosaic of polarised emission in the northern sky at low radio frequencies (150 MHz) to date. The large area this mosaic covers allows for detailed morphological and statistical studies of polarised structures in the high-latitude outer Galaxy, including the well-known Loop III region. We produced a 3100 square degree Faraday tomographic cube using a rotation measure synthesis tool. We calculated the statistical moments of Faraday spectra and compared them with data sets at higher frequencies (1.4 GHz) and with a map of a rotation measure derived from extragalactic sources. The mosaic is dominated by polarised emission connected to Loop III. Additionally, the mosaic reveals an abundance of other morphological structures, mainly {narrow and extended} depolarisation canals, which are found to be ubiquitous. We find a correlation between the map of an extragalactic rotation measure and the LoTSS first Faraday moment image. The ratio of the two deviates from a simple model of a Burn slab (Burn 1966) along the line of sight, which highlights the high level of complexity in the magnetoionic medium that can be studied at these frequencies. △ Less

Submitted 4 July, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 20 pages, 25 figures, accepted for publication in A&A

Journal ref: A&A 663, A7 (2022)

The LOFAR Two-metre Sky Survey -- V. Second data release

Authors: T. W. Shimwell, M. J. Hardcastle, C. Tasse, P. N. Best, H. J. A. Röttgering, W. L. Williams, A. Botteon, A. Drabent, A. Mechev, A. Shulevski, R. J. van Weeren, L. Bester, M. Brüggen, G. Brunetti, J. R. Callingham, K. T. Chyży, J. E. Conway, T. J. Dijkema, K. Duncan, F. de Gasperin, C. L. Hale, M. Haverkorn, B. Hugo, N. Jackson, M. Mevius , et al. (81 additional authors not shown)

Abstract: In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for th… ▽ More In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

Multi-tracer analysis of straight depolarisation canals in the surroundings of the 3C 196 field

Authors: Luka Turić, Vibor Jelić, Rutger Jaspers, Marijke Haverkorn, Andrea Bracco, Ana Erceg, Lana Ceraj, Cameron van Eck, Saleem Zaroubi

Abstract: Faraday tomography of a field centred on the extragalactic point source 3C 196 with the LOw Frequency ARray (LOFAR) revealed an intertwined structure of diffuse polarised emission with straight depolarisation canals and tracers of the magnetized and multi-phase interstellar medium (ISM), such as dust and line emission from atomic hydrogen (HI). This study aims at extending the multi-tracer analysi… ▽ More Faraday tomography of a field centred on the extragalactic point source 3C 196 with the LOw Frequency ARray (LOFAR) revealed an intertwined structure of diffuse polarised emission with straight depolarisation canals and tracers of the magnetized and multi-phase interstellar medium (ISM), such as dust and line emission from atomic hydrogen (HI). This study aims at extending the multi-tracer analysis of LOFAR data to three additional fields in the surroundings of the 3C 196 field. For the first time, we study the three-dimensional structure of the LOFAR emission by determining the distance to the depolarisation canals. We use the Rolling Hough Transform to compare the orientation of the depolarisation canals with that of the filamentary structure seen in HI and, based on starlight and dust polarisation data, with that of the plane-of-the-sky magnetic field. Stellar parallaxes from $Gaia$ complement the starlight polarisation with the corresponding distances. Faraday tomography of the three fields shows a rich network of diffuse polarised emission at Faraday depths between $-10~{\rm rad~m^{-2}}$ and $+15~{\rm rad~m^{-2}}$. A complex system of straight depolarisation canals resembles that of the 3C 196 field. The depolarisation canals align both with the HI filaments and with the magnetic field probed by dust. The observed alignment suggests that an ordered magnetic field organises the multiphase ISM over a large area ($\sim$20$^{\circ}$). In one field, two groups of stars at distances below and above 200 pc, respectively, show distinct magnetic-field orientations. These are both comparable with the orientations of the depolarisation canals in the same field. We conclude that the depolarisation canals likely trace the same change of the magnetic field as probed by the stars, which corresponds to the edge of the Local Bubble. △ Less

Submitted 7 September, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: 14 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 654, A5 (2021)

Spectral index of synchrotron emission: insights from the diffuse and magnetised interstellar medium

Authors: Marco Padovani, Andrea Bracco, Vibor Jelić, Daniele Galli, Elena Bellomi

Abstract: The interpretation of Galactic synchrotron observations is complicated by the degeneracy between the strength of the magnetic field perpendicular to the line of sight (LOS), $B_\perp$, and the cosmic-ray electron (CRe) spectrum. Depending on the observing frequency, an energy-independent spectral energy slope $s$ for the CRe spectrum is usually assumed: $s=-2$ at frequencies below $\simeq$400 MHz… ▽ More The interpretation of Galactic synchrotron observations is complicated by the degeneracy between the strength of the magnetic field perpendicular to the line of sight (LOS), $B_\perp$, and the cosmic-ray electron (CRe) spectrum. Depending on the observing frequency, an energy-independent spectral energy slope $s$ for the CRe spectrum is usually assumed: $s=-2$ at frequencies below $\simeq$400 MHz and $s=-3$ at higher frequencies. Motivated by the high angular and spectral resolution of current facilities such as the LOw Frequency ARray (LOFAR) and future telescopes such as the Square Kilometre Array (SKA), we aim to understand the consequences of taking into account the energy-dependent CRe spectral energy slope on the analysis of the spatial variations of the brightness temperature spectral index, $β$, and on the estimate of the average value of $B_\perp$ along the LOS. We illustrate analytically and numerically the impact that different realisations of the CRe spectrum have on the interpretation of the spatial variation of $β$. We find that the common assumption of an energy-independent $s$ is valid only in special cases. We show that for typical magnetic field strengths of the diffuse interstellar medium ($\simeq$2$-$20 $μ$G), at frequencies of 0.1$-$10 GHz, the electrons that are mainly responsible for the synchrotron emission have energies in the range $\simeq$100 MeV$-$50 GeV. This is the energy range where the spectral slope, $s$, of CRe has its greatest variation. We also show that the polarisation fraction can be much smaller than the maximum value of $\simeq 70\%$ because the orientation of ${\bf B}_\perp$ varies across the telescope's beam and along the LOS. Finally, we present a look-up plot that can be used to estimate the average value of $B_\perp$ along the LOS from a set of values of $β$ measured at centimetre to metre wavelengths, for a given CRe spectrum. △ Less

Submitted 21 June, 2021; originally announced June 2021.

Comments: 12 pages, 11 figures, accepted by A&A

Journal ref: A&A 651, A116 (2021)

LOFAR Deep Fields: Probing a broader population of polarized radio galaxies in ELAIS-N1

Authors: N. Herrera Ruiz, S. P. O'Sullivan, V. Vacca, V. Jelić, B. Nikiel-Wroczyński, S. Bourke, J. Sabater, R. -J. Dettmar, G. Heald, C. Horellou, S. Piras, C. Sobey, T. W. Shimwell, C. Tasse, M. J. Hardcastle, R. Kondapally, K. T. Chyży, M. Iacobelli, P. N. Best, M. Brüggen, E. Carretti, I. Prandoni

Abstract: We present deep polarimetric observations of the European Large Area ISO Survey-North 1 (ELAIS-N1) field using the Low Frequency Array (LOFAR) at 114.9-177.4 MHz. The ELAIS-N1 field is part of the LOFAR Two-metre Sky Survey deep fields data release I. For six eight-hour observing epochs, we align the polarization angles and stack the 20"-resolution Stokes $Q$, $U$-parameter data cubes. This produc… ▽ More We present deep polarimetric observations of the European Large Area ISO Survey-North 1 (ELAIS-N1) field using the Low Frequency Array (LOFAR) at 114.9-177.4 MHz. The ELAIS-N1 field is part of the LOFAR Two-metre Sky Survey deep fields data release I. For six eight-hour observing epochs, we align the polarization angles and stack the 20"-resolution Stokes $Q$, $U$-parameter data cubes. This produces a 16 deg$^2$ image with 1$σ_{\rm QU}$ sensitivity of 26 $μ$Jy/beam in the central area. In this paper, we demonstrate the feasibility of the stacking technique, and we generate a catalog of polarized sources in ELAIS-N1 and their associated Faraday rotation measures (RMs). While in a single-epoch observation we detect three polarized sources, this number increases by a factor of about three when we consider the stacked data, with a total of ten sources. This yields a surface density of polarized sources of one per 1.6 deg$^2$. The Stokes $I$ images of three of the ten detected polarized sources have morphologies resembling those of FR I radio galaxies. This represents a greater fraction of this type of source than previously found, which suggests that more sensitive observations may help with their detection. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: This paper is part of the 1st data release of the LoTSS Deep Fields. Accepted for publication in A&A. 14 pages, 9 figures, 7 tables

Journal ref: A&A 648, A12 (2021)

The LOFAR Two Metre Sky Survey: Deep Fields. II. The ELAIS-N1 LOFAR deep field

Authors: J. Sabater, P. N. Best, C. Tasse, M. J. Hardcastle, T. W. Shimwell, D. Nisbet, V. Jelic, J. R. Callingham, H. J. A. Rottgering, M. Bonato, M. Bondi, B. Ciardi, R. K. Cochrane, M. J. Jarvis, R. Kondapally, L. V. E. Koopmans, S. P. O'Sullivan, I. Prandoni, D. J. Schwarz, D. J. B. Smith, L. Wang, W. L. Williams, S. Zaroubi

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) will cover the full northern sky and, additionally, aims to observe the LoTSS deep fields to a noise level of ~10 microJy/bm over several tens of square degrees in areas that have the most extensive ancillary data. This paper presents the ELAIS-N1 deep field, the deepest of the LoTSS deep fields to date. With an effective observing time of 163.7 hours, it rea… ▽ More The LOFAR Two-metre Sky Survey (LoTSS) will cover the full northern sky and, additionally, aims to observe the LoTSS deep fields to a noise level of ~10 microJy/bm over several tens of square degrees in areas that have the most extensive ancillary data. This paper presents the ELAIS-N1 deep field, the deepest of the LoTSS deep fields to date. With an effective observing time of 163.7 hours, it reaches a root mean square (RMS) noise level below 20 microJy/bm in the central region (and below 30 microJy/bm over 10 square degrees). The resolution is 6 arcsecs and 84862 radio sources were detected in the full area (68 sq. deg.) with 74127 sources in the highest quality area at less than 3 degrees from the pointing centre. The observation reaches a sky density of more than 5000 sources per sq. deg. in the central ~5 sq. deg. region. We present the calibration procedure, which addresses the special configuration of some observations and the extended bandwidth covered (115 to 177 MHz; central frequency 146.2 MHz) compared to standard LoTSS. We also describe the methods used to calibrate the flux density scale using cross-matching with sources detected by other radio surveys in the literature. We find the flux density uncertainty related to the flux density scale to be ~6.5%. By studying the variations of the flux density measurements between different epochs, we show that relative flux density calibration is reliable out to about a 3 degree radius, but that additional flux density uncertainty is present for all sources at about the 3 per cent level; this is likely to be associated with residual calibration errors, and is shown to be more significant in datasets with poorer ionosphere conditions. We also provide intra-band spectral indices, which can be useful to detect sources with unusual spectral properties. The final uncertainty in the flux densities is estimated to be ~10% for ELAIS-N1. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: Accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields. Electronic data catalogues will be made available on journal publication. 21 pages, 17 figures

Journal ref: A&A 648, A2 (2021)

The multiphase and magnetized neutral hydrogen seen by LOFAR

Authors: Andrea Bracco, Vibor Jelić, Antoine Marchal, Luka Turić, Ana Erceg, Marc-Antoine Miville-Deschênes, Elena Bellomi

Abstract: Faraday tomography of polarimetric observations at low frequency is a unique tool to study the structure of the magneto-ionic interstellar medium (ISM) based on Faraday depth. LOFAR data below 200 MHz revealed a plethora of features in polarization, whose origin remains unknown. Previous studies highlighted the remarkable association of such features to tracers of the magnetized-neutral ISM, as in… ▽ More Faraday tomography of polarimetric observations at low frequency is a unique tool to study the structure of the magneto-ionic interstellar medium (ISM) based on Faraday depth. LOFAR data below 200 MHz revealed a plethora of features in polarization, whose origin remains unknown. Previous studies highlighted the remarkable association of such features to tracers of the magnetized-neutral ISM, as interstellar dust and atomic hydrogen (HI). However, the physical conditions responsible for the correlation between magneto-ionic and neutral media are yet to be clarified. In this letter we investigate further the correlation between LOFAR data and HI observations at 21cm from the Effelsberg-Bonn HI Survey (EBHIS). We focus on the multiphase properties of the HI gas. We present the first statistical study on the morphological correlation between LOFAR tomographic data and the cold (CNM), luke-warm (LNM), and warm (WNM) HI phases, separately. We use the Regularized Optimization for Hyper-Spectral Analysis (ROHSA) approach to decompose the HI phases based on the Gaussian decomposition of the HI spectra. In at least two fields of view -- Fields 3C196 and A -- out of four -- Fields B and C -- we find a significant correlation between LOFAR and EBHIS data using the Histograms of Oriented Gradients (HOG). The absence of correlation in Fields B and C is caused by low signal-to-noise ratio in polarization. The observed HOG correlation in Fields 3C196 and A is associated with all HI phases and it is surprisingly dominant in the CNM and LNM phases. We discuss possible mechanisms that would explain the correlation between CNM, LNM, and WNM, with polarized emission at Faraday depths up to 10 rad m$^{-2}$. Our results show how the complex structure of the ionic medium seen by LOFAR data is tightly related to phase transition in the diffuse and magnetized neutral ISM traced by HI spectroscopic data. △ Less

Submitted 18 November, 2020; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 9 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 644, L3 (2020)

Comparing Foreground Removal Techniques for Recovery of the LOFAR-EoR 21cm Power Spectrum

Authors: Ian Hothi, Emma Chapman, Jonathan R. Pritchard, F. G. Mertens, L. V. E Koopmans, B. Ciardi, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, V. Jelić, S. Zaroubi

Abstract: We compare various foreground removal techniques that are being utilised to remove bright foregrounds in various experiments aiming to detect the redshifted 21cm signal of neutral hydrogen from the Epoch of Reionization. In this work, we test the performance of removal techniques (FastICA, GMCA, and GPR) on 10 nights of LOFAR data and investigate the possibility of recovering the latest upper limi… ▽ More We compare various foreground removal techniques that are being utilised to remove bright foregrounds in various experiments aiming to detect the redshifted 21cm signal of neutral hydrogen from the Epoch of Reionization. In this work, we test the performance of removal techniques (FastICA, GMCA, and GPR) on 10 nights of LOFAR data and investigate the possibility of recovering the latest upper limit on the 21cm signal. Interestingly, we find that GMCA and FastICA reproduce the most recent 2$σ$ upper limit of $Δ^2_{21} <$ (73)$^2$ mK$^2$ at $k=0.075~ h \mathrm{cMpc}^{-1}$, which resulted from the application of GPR. We also find that FastICA and GMCA begin to deviate from the noise-limit at \textit{k}-scales larger than $\sim 0.1 ~h \mathrm{cMpc}^{-1}$. We then replicate the data via simulations to see the source of FastICA and GMCA's limitations, by testing them against various instrumental effects. We find that no single instrumental effect, such as primary beam effects or mode-mixing, can explain the poorer recovery by FastICA and GMCA at larger \textit{k}-scales. We then test scale-independence of FastICA and GMCA, and find that lower \textit{k}-scales can be modelled by a smaller number of independent components. For larger scales ($k \gtrsim 0.1~h \mathrm{cMpc}^{-1}$), more independent components are needed to fit the foregrounds. We conclude that, the current usage of GPR by the LOFAR collaboration is the appropriate removal technique. It is both robust and less prone to overfitting, with future improvements to GPR's fitting optimisation to yield deeper limits. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: 15 pages, 11 figures, accepted by MNRAS on 02 November 2020

Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\boldsymbol{z \approx 9.1}$ from LOFAR

Authors: F. G. Mertens, M. Mevius, L. V. E Koopmans, A. R. Offringa, G. Mellema, S. Zaroubi, M. A. Brentjens, H. Gan, B. K. Gehlot, V. N. Pandey, A. M. Sardarabadi, H. K. Vedantham, S. Yatawatta, K. M. B. Asad, B. Ciardi, E. Chapman, S. Gazagnes, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, V. Jelić, R. Kooistra, R. Mondal, J. Schaye , et al. (1 additional authors not shown)

Abstract: A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due… ▽ More A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$σ$ upper limit of $Δ^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power. △ Less

Submitted 19 February, 2020; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: 27 pages, 14 figues, accepted in MNRAS (updated with reference to accompanying paper)

Constraining the intergalactic medium at $z\approx$ 9.1 using LOFAR Epoch of Reionization observations

Authors: R. Ghara, S. K. Giri, G. Mellema, B. Ciardi, S. Zaroubi, I. T. Iliev, L. V. E. Koopmans, E. Chapman, S. Gazagnes, B. K. Gehlot, A. Ghosh, V. Jelic, F. G. Mertens, R. Mondal, J. Schaye, M. B. Silva, K. M. B. Asad, R. Kooistra, M. Mevius, A. R. Offringa, V. N. Pandey, S. Yatawatta

Abstract: We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY an… ▽ More We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 22 pages, 15 Figures, 5 tables, Accepted for publication in MNRAS

Predictions for the 21cm-galaxy cross-power spectrum observable with SKA and future galaxy surveys

Authors: Dijana Vrbanec, Benedetta Ciardi, Vibor Jelic, Hannes Jensen, Ilian T. Iliev, Garrelt Mellema, Saleem Zaroubi

Abstract: In this paper we use radiative transfer + N-body simulations to explore the feasibility of measurements of cross-correlations between the 21cm field observed by the Square Kilometer Array (SKA) and high-z Lyman Alpha Emitters (LAEs) detected in galaxy surveys with the Subaru Hyper Supreme Cam (HSC), Subaru Prime Focus Spectrograph (PFS) and Wide Field Infrared Survey Telescope (WFIRST). 21cm-LAE c… ▽ More In this paper we use radiative transfer + N-body simulations to explore the feasibility of measurements of cross-correlations between the 21cm field observed by the Square Kilometer Array (SKA) and high-z Lyman Alpha Emitters (LAEs) detected in galaxy surveys with the Subaru Hyper Supreme Cam (HSC), Subaru Prime Focus Spectrograph (PFS) and Wide Field Infrared Survey Telescope (WFIRST). 21cm-LAE cross-correlations are in fact a powerful probe of the epoch of reionization as they are expected to provide precious information on the progress of reionization and the typical scale of ionized regions at different redshifts. The next generation observations with SKA will have a noise level much lower than those with its precursor radio facilities, introducing a significant improvement in the measurement of the cross-correlations. We find that an SKA-HSC/PFS observation will allow to investigate scales below ~10 Mpc/h and ~60 Mpc/h at z=7.3 and 6.6, respectively. WFIRST will allow to access also higher redshifts, as it is expected to observe spectroscopically ~900 LAEs per square degree and unit redshift in the range 7.5<z<8.5. Because of the reduction of the shot noise compared to HSC and PFS, observations with WFIRST will result in more precise cross-correlations and increased observable scales. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 8 pages, 8 figures, accepted for publication in MNRAS

Foregrounds and their mitigation

Authors: Emma Chapman, Vibor Jelić

Abstract: The low-frequency radio sky is dominated by the diffuse synchrotron emission of our Galaxy and extragalactic radio sources related to Active Galactic Nuclei and star-forming galaxies. This foreground emission is much brighter than the cosmological 21 cm emission from the Cosmic Dawn and Epoch of Reionization. Studying the physical properties of the foregrounds is therefore of fundamental importanc… ▽ More The low-frequency radio sky is dominated by the diffuse synchrotron emission of our Galaxy and extragalactic radio sources related to Active Galactic Nuclei and star-forming galaxies. This foreground emission is much brighter than the cosmological 21 cm emission from the Cosmic Dawn and Epoch of Reionization. Studying the physical properties of the foregrounds is therefore of fundamental importance for their mitigation in the cosmological 21 cm experiments. This chapter gives a comprehensive overview of the foregrounds and our current state-of-the-art knowledge about their mitigation. △ Less

Submitted 26 September, 2019; originally announced September 2019.

Comments: 55 pages, 15 figures, to appear as a book chapter in "The Cosmic 21-cm Revolution: Charting the first billion years of our Universe", Ed. Andrei Mesinger (Bristol: IOP Publishing Ltd.) AAS-IOP ebooks (http://www.iopscience.org/books/aas)

Peering into the Dark (Ages) with Low-Frequency Space Interferometers

Authors: Leon Koopmans, Rennan Barkana, Mark Bentum, Gianni Bernardi, Albert-Jan Boonstra, Judd Bowman, Jack Burns, Xuelei Chen, Abhirup Datta, Heino Falcke, Anastasia Fialkov, Bharat Gehlot, Leonid Gurvits, Vibor Jelić, Marc Klein-Wolt, Léon Koopmans, Joseph Lazio, Daan Meerburg, Garrelt Mellema, Florent Mertens, Andrei Mesinger, André Offringa, Jonathan Pritchard, Benoit Semelin, Ravi Subrahmanyan , et al. (6 additional authors not shown)

Abstract: Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at… ▽ More Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged] △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: 24 pages, 8 figures, 3 tables; White paper submitted to ESA Voyage 2050

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

The first power spectrum limit on the 21-cm signal of neutral hydrogen during the Cosmic Dawn at z=20-25 from LOFAR

Authors: B. K. Gehlot, F. G. Mertens, L. V. E. Koopmans, M. A. Brentjens, S. Zaroubi, B. Ciardi, A. Ghosh, M. Hatef, I. T. Iliev, V. Jelić, R. Kooistra, F. Krause, G. Mellema, M. Mevius, M. Mitra, A. R. Offringa, V. N. Pandey, A. M. Sardarabadi, J. Schaye, M. B. Silva, H. K. Vedantham, S. Yatawatta

Abstract: Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz fr… ▽ More Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz frequency range) using 14 hours of data obtained with the LOFAR-Low Band Antenna (LBA) array. We also demonstrate the application of a multiple pointing calibration technique to calibrate the LOFAR-LBA dual-pointing observations centred on the North Celestial Pole and the radio galaxy 3C220.3. We observe an unexplained excess of $\sim 30-50\%$ in Stokes $I$ noise compared to Stokes $V$ for the two observed fields, which decorrelates on $\gtrsim 12$ seconds and might have a physical origin. We show that enforcing smoothness of gain errors along frequency direction during calibration reduces the additional variance in Stokes $I$ compared Stokes $V$ introduced by the calibration on sub-band level. After subtraction of smooth foregrounds, we achieve a $2σ$ upper limit on the 21-cm power spectrum of $Δ_{21}^2 < (14561\,\text{mK})^2$ at $k\sim 0.038\,h\,\text{cMpc}^{-1}$ and $Δ_{21}^2 < (14886\,\text{mK})^2$ at $k\sim 0.038 \,h\,\text{cMpc}^{-1}$ for the 3C220 and NCP fields respectively and both upper limits are consistent with each other. The upper limits for the two fields are still dominated by systematics on most $k$ modes. △ Less

Submitted 20 July, 2019; v1 submitted 18 September, 2018; originally announced September 2018.

Comments: 19 pages, 16 figures, accepted for publication in MNRAS

Magnetically aligned straight depolarisation canals and the rolling Hough transform

Authors: Vibor Jelić, David Prelogović, Marijke Haverkorn, Jur Remeijn, Dora Klindžić

Abstract: Aims. We aim to characterize the properties of the straight depolarization canals detected in the Low Frequency Array (LOFAR) polarimetric observations of a field centered on the extragalactic source 3C 196. We also compare the canal orientations with magnet- ically aligned Hi filaments and the magnetic field probed by polarized dust emission. Methods. We used the rolling Hough transform (RHT) to… ▽ More Aims. We aim to characterize the properties of the straight depolarization canals detected in the Low Frequency Array (LOFAR) polarimetric observations of a field centered on the extragalactic source 3C 196. We also compare the canal orientations with magnet- ically aligned Hi filaments and the magnetic field probed by polarized dust emission. Methods. We used the rolling Hough transform (RHT) to identify and characterize the orientation of the straight depolarization canals in radio polarimetric data and the filaments in Hi data. Results. The majority of the straight depolarization canals and the Hi filaments are inclined by 10deg with respect to the Galactic plane and are aligned with the plane-of-sky magnetic field orientation probed by the Planck dust polarization data. The other distinct orientation, of 65deg with respect to the Galactic plane, is associated with the orientation of a bar-like structure observed in the 3C 196 field at 350 MHz. Conclusions. An alignment between three distinct tracers of the (local) interstellar medium (ISM) suggests that an ordered magnetic field plays a crucial role in confining different ISM phases. The majority of the straight depolarization canals are a result of a projection of the complicated 3D distribution of the ISM. The RHT analysis is a robust method for identifying and characterizing the straight depolarization canals observed in radio-polarimetric data. △ Less

Submitted 18 June, 2018; originally announced June 2018.

Comments: 6 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 615, L3 (2018)

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)

Wide-field LOFAR-LBA power-spectra analyses: Impact of calibration, polarization leakage and ionosphere

Authors: B. K. Gehlot, L. V. E. Koopmans, A. G. de Bruyn, S. Zaroubi, M. A. Brentjens, K. M. B. Asad, M. Hatef, V. Jelic, M. Mevius, A. R. Offringa, V. N. Pandey, S. Yatawatta

Abstract: Contamination due to foregrounds (Galactic and Extra-galactic), calibration errors and ionospheric effects pose major challenges in detection of the cosmic 21 cm signal in various Epoch of Reionization (EoR) experiments. We present the results of a pilot study of a field centered on 3C196 using LOFAR Low Band (56-70 MHz) observations, where we quantify various wide field and calibration effects su… ▽ More Contamination due to foregrounds (Galactic and Extra-galactic), calibration errors and ionospheric effects pose major challenges in detection of the cosmic 21 cm signal in various Epoch of Reionization (EoR) experiments. We present the results of a pilot study of a field centered on 3C196 using LOFAR Low Band (56-70 MHz) observations, where we quantify various wide field and calibration effects such as gain errors, polarized foregrounds, and ionospheric effects. We observe a `pitchfork' structure in the 2D power spectrum of the polarized intensity in delay-baseline space, which leaks into the modes beyond the instrumental horizon (EoR/CD window). We show that this structure largely arises due to strong instrumental polarization leakage ($\sim30\%$) towards {Cas\,A} ($\sim21$ kJy at 81 MHz, brightest source in northern sky), which is far away from primary field of view. We measure an extremely small ionospheric diffractive scale ($r_{\text{diff}} \approx 430$ m at 60 MHz) towards {Cas\,A} resembling pure Kolmogorov turbulence compared to $r_{\text{diff}} \sim3 - 20$ km towards zenith at 150 MHz for typical ionospheric conditions. This is one of the smallest diffractive scales ever measured at these frequencies. Our work provides insights in understanding the nature of aforementioned effects and mitigating them in future Cosmic Dawn observations (e.g. with SKA-low and HERA) in the same frequency window. △ Less

Submitted 26 April, 2018; v1 submitted 22 September, 2017; originally announced September 2017.

Comments: 20 pages, 11 figures, accepted for publication in MNRAS

Polarization leakage in epoch of reionization windows: III. Wide-field effects of narrow-field arrays

Authors: K. M. B. Asad, L. V. E. Koopmans, V. Jelić, A. G. de Bruyn, V. N. Pandey, B. K. Gehlot

Abstract: Leakage of polarized Galactic diffuse emission into total intensity can potentially mimic the 21-cm signal coming from the epoch of reionization (EoR), as both of them might have fluctuating spectral structure. Although we are sensitive to the EoR signal only in small fields of view, chromatic sidelobes from further away can contaminate the inner region. Here, we explore the effects of leakage int… ▽ More Leakage of polarized Galactic diffuse emission into total intensity can potentially mimic the 21-cm signal coming from the epoch of reionization (EoR), as both of them might have fluctuating spectral structure. Although we are sensitive to the EoR signal only in small fields of view, chromatic sidelobes from further away can contaminate the inner region. Here, we explore the effects of leakage into the 'EoR window' of the cylindrically averaged power spectra (PS) within wide fields of view using both observation and simulation of the 3C196 and NCP fields, two observing fields of the LOFAR-EoR project. We present the polarization PS of two one-night observations of the two fields and find that the NCP field has higher fluctuations along frequency, and consequently exhibits more power at high-$k_\parallel$ that could potentially leak to Stokes $I$. Subsequently, we simulate LOFAR observations of Galactic diffuse polarized emission based on a model to assess what fraction of polarized power leaks into Stokes $I$ because of the primary beam. We find that the rms fractional leakage over the instrumental $k$-space is $0.35\%$ in the 3C196 field and $0.27\%$ in the NCP field, and it does not change significantly within the diameters of $15^\circ$, $9^\circ$ and $4^\circ$. Based on the observed PS and simulated fractional leakage, we show that a similar level of leakage into Stokes $I$ is expected in the 3C196 and NCP fields, and the leakage can be considered to be a bias in the PS. △ Less

Submitted 2 June, 2017; originally announced June 2017.

Comments: 13 pages, 8 figures, submitted to MNRAS

The VLA-COSMOS 3 GHz Large Project: Continuum data and source catalog release

Authors: V. Smolcic, M. Novak, M. Bondi, P. Ciliegi, K. P. Mooley, E. Schinnerer, G. Zamorani, F. Navarrete, S. Bourke, A. Karim, E. Vardoulaki, S. Leslie, J. Delhaize, C. L. Carilli, S. T. Myers, N. Baran, I. Delvecchio, O. Miettinen, J. Banfield, M. Balokovic, F. Bertoldi, P. Capak, D. A. Frail, G. Hallinan, H. Hao , et al. (15 additional authors not shown)

Abstract: We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations acr… ▽ More We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys. △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: 19 pages, 18 figures, 5 tables, accepted for publication in A&A

Upper limits on the 21-cm Epoch of Reionization power spectrum from one night with LOFAR

Authors: A. H. Patil, S. Yatawatta, L. V. E. Koopmans, A. G. de Bruyn, M. A. Brentjens, S. Zaroubi, K. M. B. Asad, M. Hatef, V. Jelic, M. Mevius, A. R. Offringa, V. N. Pandey, H. Vedantham, F. B. Abdalla, W. N. Brouw, E. Chapman, B. Ciardi, B. K. Gehlot, A. Ghosh, G. Harker, I. T. Iliev, K. Kakiichi, S. Majumdar, M. B. Silva, G. Mellema , et al. (3 additional authors not shown)

Abstract: We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $Δ^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-… ▽ More We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $Δ^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-$σ$) excess variance and a best 2-$σ$ upper limit of $Δ^2_{\rm 21} < (79.6 {\rm mK})^2$ at $k=0.053$$h$cMpc$^{-1}$ in the range $z=$9.6-10.6. The excess variance decreases when optimizing the smoothness of the direction- and frequency-dependent gain calibration, and with increasing the completeness of the sky model. It is likely caused by (i) residual side-lobe noise on calibration baselines, (ii) leverage due to non-linear effects, (iii) noise and ionosphere-induced gain errors, or a combination thereof. Further analyses of the excess variance will be discussed in forthcoming publications. △ Less

Submitted 28 February, 2017; originally announced February 2017.

Comments: 17 pages, 8 figures, accepted by ApJ

Recovering the HII region size statistics from 21-cm tomography

Authors: Koki Kakiichi, Suman Majumdar, Garrelt Mellema, Benedetta Ciardi, Keri L. Dixon, Ilian T. Iliev, Vibor Jelic, Leon V. E. Koopmans, Saleem Zaroubi, Philipp Busch

Abstract: We introduce a novel technique, called "granulometry", to characterize and recover the mean size and the size distribution of HII regions from 21-cm tomography. The technique is easy to implement, but places the previously not very well defined concept of morphology on a firm mathematical foundation. The size distribution of the cold spots in 21-cm tomography can be used as a direct tracer of the… ▽ More We introduce a novel technique, called "granulometry", to characterize and recover the mean size and the size distribution of HII regions from 21-cm tomography. The technique is easy to implement, but places the previously not very well defined concept of morphology on a firm mathematical foundation. The size distribution of the cold spots in 21-cm tomography can be used as a direct tracer of the underlying probability distribution of HII region sizes. We explore the capability of the method using large-scale reionization simulations and mock observational data cubes while considering capabilities of SKA1-low and a future extension to SKA2. We show that the technique allows the recovery of the HII region size distribution with a moderate signal-to-noise ratio from wide-field imaging ($\rm SNR\lesssim3$), for which the statistical uncertainty is sample variance dominated. We address the observational requirements on the angular resolution, the field-of-view, and the thermal noise limit for a successful measurement. To achieve a full scientific return from 21-cm tomography and to exploit a synergy with 21-cm power spectra, we suggest an observing strategy using wide-field imaging (several tens of square degrees) by an interferometric mosaicking/multi-beam observation with additional intermediate baselines (~2-4 km). △ Less

Submitted 20 October, 2017; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: 19 pages, 19 figures, the version accepted in MNRAS

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)

Probing Ionospheric Structures using the LOFAR radio telescope

Authors: M. Mevius, S. van der Tol, V. N. Pandey, H. K. Vedantham, M. A. Brentjens, A. G. de Bruyn, F. B. Abdalla, K. M. B. Asad, J. D. Bregman, W. N. Brouw, S. Bus, E. Chapman, B. Ciardi, E. R. Fernandez, A. Ghosh, G. Harker, I. T. Iliev, V. Jelić, S. Kazemi, L. V. E. Koopmans, J. E. Noordam, A. R. Offringa, A. H. Patil, R. J. van Weeren, S. Wijnholds , et al. (2 additional authors not shown)

Abstract: LOFAR is the LOw Frequency Radio interferometer ARray located at mid-latitude ($52^{\circ} 53'N$). Here, we present results on ionospheric structures derived from 29 LOFAR nighttime observations during the winters of 2012/2013 and 2013/2014. We show that LOFAR is able to determine differential ionospheric TEC values with an accuracy better than 1 mTECU over distances ranging between 1 and 100 km.… ▽ More LOFAR is the LOw Frequency Radio interferometer ARray located at mid-latitude ($52^{\circ} 53'N$). Here, we present results on ionospheric structures derived from 29 LOFAR nighttime observations during the winters of 2012/2013 and 2013/2014. We show that LOFAR is able to determine differential ionospheric TEC values with an accuracy better than 1 mTECU over distances ranging between 1 and 100 km. For all observations the power law behavior of the phase structure function is confirmed over a long range of baseline lengths, between $1$ and $80$ km, with a slope that is in general larger than the $5/3$ expected for pure Kolmogorov turbulence. The measured average slope is $1.89$ with a one standard deviation spread of $0.1$. The diffractive scale, i.e. the length scale where the phase variance is $1\, \mathrm{rad^2}$, is shown to be an easily obtained single number that represents the ionospheric quality of a radio interferometric observation. A small diffractive scale is equivalent to high phase variability over the field of view as well as a short time coherence of the signal, which limits calibration and imaging quality. For the studied observations the diffractive scales at $150$ MHz vary between $3.5$ and $30\,$ km. A diffractive scale above $5$ km, pertinent to about $90 \%$ of the observations, is considered sufficient for the high dynamic range imaging needed for the LOFAR Epoch of Reionization project. For most nights the ionospheric irregularities were anisotropic, with the structures being aligned with the Earth magnetic field in about $60\%$ of the observations. △ Less

Submitted 15 June, 2016; originally announced June 2016.

Comments: accepted for publication in Radio Science

Journal ref: Radio Sci. 51 (2016) 927-941

Systematic biases in low frequency radio interferometric data due to calibration: the LOFAR EoR case

Authors: Ajinkya H. Patil, Sarod Yatawatta, Saleem Zaroubi, Léon V. E. Koopmans, A. G. de Bruyn, Vibor Jelić, Benedetta Ciardi, Ilian T. Iliev, Maaijke Mevius, Vishambhar N. Pandey, Bharat K. Gehlot

Abstract: The redshifted 21 cm line of neutral hydrogen is a promising probe of the Epoch of Reionization (EoR). However, its detection requires a thorough understanding and control of the systematic errors. We study two systematic biases observed in the LOFAR EoR residual data after calibration and subtraction of bright discrete foreground sources. The first effect is a suppression in the diffuse foregroun… ▽ More The redshifted 21 cm line of neutral hydrogen is a promising probe of the Epoch of Reionization (EoR). However, its detection requires a thorough understanding and control of the systematic errors. We study two systematic biases observed in the LOFAR EoR residual data after calibration and subtraction of bright discrete foreground sources. The first effect is a suppression in the diffuse foregrounds, which could potentially mean a suppression of the 21 cm signal. The second effect is an excess of noise beyond the thermal noise. The excess noise shows fluctuations on small frequency scales, and hence it can not be easily removed by foreground removal or avoidance methods. Our analysis suggests that sidelobes of residual sources due to the chromatic point spread function and ionospheric scintillation can not be the dominant causes of the excess noise. Rather, both the suppression of diffuse foregrounds and the excess noise can occur due to calibration with an incomplete sky model containing predominantly bright discrete sources. We show that calibrating only on bright sources can cause suppression of other signals and introduce an excess noise in the data. The levels of the suppression and excess noise depend on the relative flux of sources which are not included in the model with respect to the flux of modeled sources. We discuss possible solutions such as using only long baselines to calibrate the interferometric gain solutions as well as simultaneous multi-frequency calibration along with their benefits and shortcomings. △ Less

Submitted 24 May, 2016; originally announced May 2016.

Comments: Submitted to MNRAS, 14 pages, 11 figures

Polarization leakage in epoch of reionization windows: II. Primary beam model and direction dependent calibration

Authors: K. M. B. Asad, L. V. E. Koopmans, V. Jelić, A. Ghosh, F. B. Abdalla, M. A. Brentjens, A. G. de Bruyn, B. Ciardi, B. K. Gehlot, I. T. Iliev, M. Mevius, V. N. Pandey, S. Yatawatta, S. Zaroubi

Abstract: Leakage of diffuse polarized emission into Stokes I caused by the polarized primary beam of the instrument might mimic the spectral structure of the 21-cm signal coming from the epoch of reionization (EoR) making their separation difficult. Therefore, understanding polarimetric performance of the antenna is crucial for a successful detection of the EoR signal. Here, we have calculated the accuracy… ▽ More Leakage of diffuse polarized emission into Stokes I caused by the polarized primary beam of the instrument might mimic the spectral structure of the 21-cm signal coming from the epoch of reionization (EoR) making their separation difficult. Therefore, understanding polarimetric performance of the antenna is crucial for a successful detection of the EoR signal. Here, we have calculated the accuracy of the nominal model beam of LOFAR in predicting the leakage from Stokes I to Q, U by comparing them with the corresponding leakage of compact sources actually observed in the 3C295 field. We have found that the model beam has errors of less than or equal to 10% on the predicted levels of leakage of ~1% within the field of view, i. e. if the leakage is taken out perfectly using this model the leakage will reduce to $10^{-3}$ of the Stokes I flux. If similar levels of accuracy can be obtained in removing leakage from Stokes Q, U to I, we can say, based on the results of our previous paper, that the removal of this leakage using this beam model would ensure that the leakage is well below the expected EoR signal in almost the whole instrumental k-space of the cylindrical power spectrum. We have also shown here that direction dependent calibration can remove instrumentally polarized compact sources, given an unpolarized sky model, very close to the local noise level. △ Less

Submitted 27 July, 2016; v1 submitted 15 April, 2016; originally announced April 2016.

Comments: 13 pages, 15 figures, accepted in MNRAS on 26 July 2016

Effects of the sources of reionization on 21-cm redshift-space distortions

Authors: Suman Majumdar, Hannes Jensen, Garrelt Mellema, Emma Chapman, Filipe B. Abdalla, Kai-Yan Lee, Ilian T. Iliev, Keri L. Dixon, Kanan K. Datta, Benedetta Ciardi, Elizabeth R. Fernandez, Vibor Jelić, Léon V. E. Koopmans, Saleem Zaroubi

Abstract: The observed 21-cm signal from the epoch of reionization will be distorted along the line-of-sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus p… ▽ More The observed 21-cm signal from the epoch of reionization will be distorted along the line-of-sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus potentially be used to extract information about the sources of reionization. In this paper, we study a collection of simulated reionization scenarios assuming different models for the sources of reionization. We show that the 21-cm anisotropy is best measured by the quadrupole moment of the power spectrum. We find that, unless the properties of the reionization sources are extreme in some way, the quadrupole moment evolves very predictably as a function of global neutral fraction. This predictability implies that redshift-space distortions are not a very sensitive tool for distinguishing between reionization sources. However, the quadrupole moment can be used as a model-independent probe for constraining the reionization history. We show that such measurements can be done to some extent by first-generation instruments such as LOFAR, while the SKA should be able to measure the reionization history using the quadrupole moment of the power spectrum to great accuracy. △ Less

Submitted 2 December, 2015; v1 submitted 24 September, 2015; originally announced September 2015.

Comments: 17 pages, 12 figures, 1 table. Accepted for publication in MNRAS. Replaced to match the accepted version

Journal ref: Monthly Notices of the Royal Astronomical Society 2016 456 (2): 2080-2094

Predictions for the 21cm-galaxy cross-power spectrum observable with LOFAR and Subaru

Authors: Dijana Vrbanec, Benedetta Ciardi, Vibor Jelić, Hannes Jensen, Saleem Zaroubi, Elizabeth R. Fernandez, Abhik Ghosh, Ilian T. Iliev, Koki Kakiichi, Léon V. E. Koopmans, Garrelt Mellema

Abstract: The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru… ▽ More The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru's Hyper Suprime Cam (HSC), we investigate the observability of such cross-power spectrum with these two instruments, which are both planning to observe the ELAIS-N1 field at z=6.6. In this paper we use N-body + radiative transfer (both for continuum and Lyalpha photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm-galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power spectrum expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum show clear anti-correlation on scales larger than ~ 60 h$^{-1}$ Mpc (corresponding to k ~ 0.1 h Mpc$^{-1}$), with levels of significance p=0.04 at z=6.6 and p=0.048 at z=7.3. On smaller scales, instead, the signal is completely contaminated. △ Less

Submitted 11 September, 2015; originally announced September 2015.

Journal ref: Monthly Notices of the Royal Astronomical Society 2016 457 (1): 666-675

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

Authors: G. H. Heald, R. F. Pizzo, E. Orrú, R. P. Breton, D. Carbone, C. Ferrari, M. J. Hardcastle, W. Jurusik, G. Macario, D. Mulcahy, D. Rafferty, A. Asgekar, M. Brentjens, R. A. Fallows, W. Frieswijk, M. C. Toribio, B. Adebahr, M. Arts, M. R. Bell, A. Bonafede, J. Bray, J. Broderick, T. Cantwell, P. Carroll, Y. Cendes , et al. (125 additional authors not shown)

Abstract: We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octave… ▽ More We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs. △ Less

Submitted 3 September, 2015; originally announced September 2015.

Comments: 23 pages, 19 figures. Accepted for publication in Astronomy & Astrophysics. MSSS Verification Field images and catalog data may be downloaded from http://vo.astron.nl

Galactic interstellar filaments as probed by LOFAR and Planck

Authors: S. Zaroubi, V. Jelić, A. G. de Bruyn, F. Boulanger, A. Bracco, R. Kooistra, M. I. R. Alves, M. A. Brentjens, K. Ferrière, T. Ghosh, L. V. E. Koopmans, F. Levrier, M. -A. Miville-Deschênes, L. Montier, V. N. Pandey, J. D. Soler

Abstract: Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structu… ▽ More Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structures along the line of sight, whereas the Planck dust polarization maps probe the orientation of the sky projected magnetic field component. Hence, no apparent correlation between the two is expected. Here we report a surprising, yet clear, correlation between the filamentary structures, detected with LOFAR, and the magnetic field orientation, probed by the Planck satellite. This finding points to a common, yet unclear, physical origin of the two measurements in this specific area in the sky. A number of follow-up multi- frequency studies are proposed to shed light on this unexpected finding. △ Less

Submitted 26 August, 2015; originally announced August 2015.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS Letters

Linear polarization structures in LOFAR observations of the interstellar medium in the 3C196 field

Authors: V. Jelić, A. G. de Bruyn, V. N. Pandey, M. Mevius, M. Haverkorn, M. A. Brentjens, L. V. E. Koopmans, S. Zaroubi, F. B. Abdalla, K. M. B. Asad, S. Bus, E. Chapman, B. Ciardi, E. R. Fernandez, A. Ghosh, G. Harker, I. T. Iliev, H. Jensen, S. Kazemi, G. Mellema, A. R. Offringa, A. H. Patil, H. K. Vedantham, S. Yatawatta

Abstract: This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. We have used the high band antennas (HBA) of LOFAR to image this region and Rotation Measure (RM) synthesis to unravel the distribution of polarized structures in Faraday depth. Th… ▽ More This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. We have used the high band antennas (HBA) of LOFAR to image this region and Rotation Measure (RM) synthesis to unravel the distribution of polarized structures in Faraday depth. The brightness temperature of the detected Galactic emission is $5-15~{\rm K}$ in polarized intensity and covers the range from -3 to +8 ${\rm rad~m^{-2}}$ in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of $+0.5~{\rm rad~m^{-2}}$ running from north to south, right through the centre of the field and parallel to the Galactic plane. There is also an interesting system of linear depolarization canals conspicuous in an image showing the peaks of Faraday spectra. We used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. For the first time, we see some common morphology in the RM cubes made at 150 and 350~{\rm MHz}. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies and previous LOFAR observations. Based on our results, we determined physical parameters of the ISM and proposed a simple model that may explain the observed distribution of the intervening magneto-ionic medium. The mean line-of-sight magnetic field component, $B_\parallel$, is determined to be $0.3\pm0.1~{\rm μG}$ and its spatial variation across the 3C196 field is $0.1~{\rm μG}$. The filamentary structure is probably an ionized filament in the ISM, located somewhere within the Local Bubble. This filamentary structure shows an excess in thermal electron density ($n_e B_\parallel>6.2~{\rm cm^{-3}μG}$) compared to its surroundings. △ Less

Submitted 21 September, 2015; v1 submitted 26 August, 2015; originally announced August 2015.

Comments: 16 pages, 10 figures, aceppted for publication in A&A

Journal ref: A&A 583, A137 (2015)

A Bayesian analysis of redshifted 21-cm HI signal and foregrounds: Simulations for LOFAR

Authors: Abhik Ghosh, L. V. E. Koopmans, Emma Chapman, Vibor Jelic

Abstract: Observations of the EoR with the 21-cm hyperfine emission of neutral hydrogen (HI) promise to open an entirely new window onto the formation of the first stars, galaxies and accreting black holes. In order to characterize the weak 21-cm signal, we need to develop imaging techniques which can reconstruct the extended emission very precisely. Here, we present an inversion technique for LOFAR baselin… ▽ More Observations of the EoR with the 21-cm hyperfine emission of neutral hydrogen (HI) promise to open an entirely new window onto the formation of the first stars, galaxies and accreting black holes. In order to characterize the weak 21-cm signal, we need to develop imaging techniques which can reconstruct the extended emission very precisely. Here, we present an inversion technique for LOFAR baselines at NCP, based on a Bayesian formalism with optimal spatial regularization, which is used to reconstruct the diffuse foreground map directly from the simulated visibility data. We notice the spatial regularization de-noises the images to a large extent, allowing one to recover the 21-cm power-spectrum over a considerable $k_{\perp}-k_{\para}$ space in the range of $0.03\,{\rm Mpc^{-1}}<k_{\perp}<0.19\,{\rm Mpc^{-1}}$ and $0.14\,{\rm Mpc^{-1}}<k_{\para}<0.35\,{\rm Mpc^{-1}}$ without subtracting the noise power-spectrum. We find that, in combination with using the GMCA, a non-parametric foreground removal technique, we can mostly recover the spherically average power-spectrum within $2σ$ statistical fluctuations for an input Gaussian random rms noise level of $60 \, {\rm mK}$ in the maps after 600 hrs of integration over a $10 \, {\rm MHz}$ bandwidth. △ Less

Submitted 16 June, 2015; originally announced June 2015.

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

The Cosmic Dawn and Epoch of Reionization with the Square Kilometre Array

Authors: L. V. E. Koopmans, J. Pritchard, G. Mellema, F. Abdalla, J. Aguirre, K. Ahn, R. Barkana, I. van Bemmel, G. Bernardi, A. Bonaldi, F. Briggs, A. G. de Bruyn, T. C. Chang, E. Chapman, X. Chen, B. Ciardi, K. K. Datta, P. Dayal, A. Ferrara, A. Fialkov, F. Fiore, K. Ichiki, I. T. Illiev, S. Inoue, V. Jelić , et al. (20 additional authors not shown)

Abstract: Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$α$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct… ▽ More Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$α$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged). △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015. PoS(AASKA14)001

Simulating the 21cm forest detectable with LOFAR and SKA in the spectra of high-z GRBs

Authors: B. Ciardi, S. Inoue, F. B. Abdalla, K. Asad, G. Bernardi, J. S. Bolton, M. Brentjens, A. G. de Bruyn, E. Chapman, S. Daiboo, E. R. Fernandez, A. Ghosh, L. Graziani, G. J. A. Harker, I. T. Iliev, V. Jelic, H. Jensen, S. Kazemi, L. V. E. Koopmans, O. Martinez, A. Maselli, G. Mellema, A. R. Offringa, V. N. Pandey, J. Schaye , et al. (4 additional authors not shown)

Abstract: We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the G… ▽ More We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9. △ Less

Submitted 20 July, 2015; v1 submitted 28 April, 2015; originally announced April 2015.

Comments: 5 pages, 4 figures, MNRAS, in press

Polarization leakage in Epoch of Reionization windows: I. LOFAR observations of the 3C196 field

Authors: K. M. B. Asad, L. V. E. Koopmans, V. Jelić, V. N. Pandey, A. Ghosh, F. B. Abdalla, G. Bernardi, M. A. Brentjens, A. G. de Bruyn, S. Bus, B. Ciardi, E. Chapman, S. Daiboo, E. R. Fernandez, G. Harker, I. T. Iliev, H. Jensen, O. Martinez-Rubi, G. Mellema, M. Mevius, A. R. Offringa, A. H. Patil, J. Schaye, R. M. Thomas, S. van der Tol , et al. (3 additional authors not shown)

Abstract: Detection of the 21-cm signal coming from the epoch of reionization (EoR) is challenging especially because, even after removing the foregrounds, the residual Stokes $I$ maps contain leakage from polarized emission that can mimic the signal. Here, we discuss the instrumental polarization of LOFAR and present realistic simulations of the leakages between Stokes parameters. From the LOFAR observatio… ▽ More Detection of the 21-cm signal coming from the epoch of reionization (EoR) is challenging especially because, even after removing the foregrounds, the residual Stokes $I$ maps contain leakage from polarized emission that can mimic the signal. Here, we discuss the instrumental polarization of LOFAR and present realistic simulations of the leakages between Stokes parameters. From the LOFAR observations of polarized emission in the 3C196 field, we have quantified the level of polarization leakage caused by the nominal model beam of LOFAR, and compared it with the EoR signal using power spectrum analysis. We found that at 134--166 MHz, within the central 4$^\circ$ of the field the $(Q,U)\rightarrow I$ leakage power is lower than the EoR signal at $k<0.3$ Mpc$^{-1}$. The leakage was found to be localized around a Faraday depth of 0, and the rms of the leakage as a fraction of the rms of the polarized emission was shown to vary between 0.2-0.3\%, both of which could be utilized in the removal of leakage. Moreover, we could define an `EoR window' in terms of the polarization leakage in the cylindrical power spectrum above the PSF-induced wedge and below $k_\parallel\sim 0.5$ Mpc$^{-1}$, and the window extended up to $k_\parallel\sim 1$ Mpc$^{-1}$ at all $k_\perp$ when 70\% of the leakage had been removed. These LOFAR results show that even a modest polarimetric calibration over a field of view of $\lesssim 4^\circ$ in the future arrays like SKA will ensure that the polarization leakage remains well below the expected EoR signal at the scales of 0.02-1 Mpc$^{-1}$. △ Less

Submitted 13 May, 2015; v1 submitted 5 March, 2015; originally announced March 2015.

Comments: 19 pages, 13 figures. Accepted for publication in MNRAS

Cosmic Dawn and Epoch of Reionization Foreground Removal with the SKA

Authors: Emma Chapman, Anna Bonaldi, Geraint Harker, Vibor Jelić, Filipe B. Abdalla, Gianni Bernardi, Jérôme Bobin, Fred Dulwich, Benjamin Mort, Mario Santos, Jean-Luc Starck

Abstract: The exceptional sensitivity of the SKA will allow observations of the Cosmic Dawn and Epoch of Reionization (CD/EoR) in unprecedented detail, both spectrally and spatially. This wealth of information is buried under Galactic and extragalactic foregrounds, which must be removed accurately and precisely in order to reveal the cosmological signal. This problem has been addressed already for the previ… ▽ More The exceptional sensitivity of the SKA will allow observations of the Cosmic Dawn and Epoch of Reionization (CD/EoR) in unprecedented detail, both spectrally and spatially. This wealth of information is buried under Galactic and extragalactic foregrounds, which must be removed accurately and precisely in order to reveal the cosmological signal. This problem has been addressed already for the previous generation of radio telescopes, but the application to SKA is different in many aspects. In this chapter we summarise the contributions to the field of foreground removal in the context of high redshift and high sensitivity 21-cm measurements. We use a state-of-the-art simulation of the SKA Phase 1 observations complete with cosmological signal, foregrounds and frequency-dependent instrumental effects to test both parametric and non-parametric foreground removal methods. We compare the recovered cosmological signal using several different statistics and explore one of the most exciting possibilities with the SKA --- imaging of the ionized bubbles. We find that with current methods it is possible to remove the foregrounds with great accuracy and to get impressive power spectra and images of the cosmological signal. The frequency-dependent PSF of the instrument complicates this recovery, so we resort to splitting the observation bandwidth into smaller segments, each of a common resolution. If the foregrounds are allowed a random variation from the smooth power law along the line of sight, methods exploiting the smoothness of foregrounds or a parametrization of their behaviour are challenged much more than non-parametric ones. However, we show that correction techniques can be implemented to restore the performances of parametric approaches, as long as the first-order approximation of a power law stands. △ Less

Submitted 19 January, 2015; originally announced January 2015.

Comments: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015

SKA - EoR correlations and cross-correlations: kSZ, radio galaxies, and NIR background

Authors: Vibor Jelic, Benedetta Ciardi, Elizabeth Fernandez, Hiroyuki Tashiro, Dijana Vrbanec

Abstract: The Universe's Cosmic Dawn (CD) and Epoch of Reionization (EoR) can be studied using a number of observational probes that provide complementary or corroborating information. Each of these probes suffers from its own systematic and statistical uncertainties. It is therefore useful to consider the mutual information that these data sets contain. In this paper, we discuss a potential of cross-correl… ▽ More The Universe's Cosmic Dawn (CD) and Epoch of Reionization (EoR) can be studied using a number of observational probes that provide complementary or corroborating information. Each of these probes suffers from its own systematic and statistical uncertainties. It is therefore useful to consider the mutual information that these data sets contain. In this paper, we discuss a potential of cross-correlations between the SKA cosmological 21 cm data with: (i) the kinetic Sunyaev- Zel'dovich (kSZ) effect in the CMB data; (ii) the galaxy surveys; and (iii) near infrared (NIR) backgrounds. △ Less

Submitted 17 January, 2015; originally announced January 2015.

Comments: Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015

The Effect of Foreground Mitigation Strategy on EoR Window Recovery

Authors: Emma Chapman, Saleem Zaroubi, Filipe Abdalla, Fred Dulwich, Vibor Jelić, Benjamin Mort

Abstract: The removal of the Galactic and extragalactic foregrounds remains a major challenge for those wishing to make a detection of the Epoch of Reionization 21-cm signal. Multiple methods of modelling these foregrounds with varying levels of assumption have been trialled and shown promising recoveries on simulated data. Recently however there has been increased discussion of using the expected shape of… ▽ More The removal of the Galactic and extragalactic foregrounds remains a major challenge for those wishing to make a detection of the Epoch of Reionization 21-cm signal. Multiple methods of modelling these foregrounds with varying levels of assumption have been trialled and shown promising recoveries on simulated data. Recently however there has been increased discussion of using the expected shape of the foregrounds in Fourier space to define an EoR window free of foreground contamination. By carrying out analysis within this window only, one can avoid the foregrounds and any statistical bias they might introduce by instead removing these foregrounds. In this paper we discuss the advantages and disadvantages of both foreground removal and foreground avoidance. We create a series of simulations with noise levels in line with both current and future experiments and compare the recovered statistical cosmological signal from foreground avoidance and a simplified, frequency independent foreground removal model. We find that while, for current generation experiments, foreground avoidance enables a better recovery at $k_{perp} > 0.6 \mathrm{Mpc}^{-1}$, foreground removal is able to recover significantly more signal at small $k_{los}$ for both current and future experiments. We also relax the assumption that the foregrounds are smooth by introducing a Gaussian random factor along the line-of-sight and then also spatially. We find that both methods perform well for foreground models with line-of-sight and spatial variations around $0.1\%$ however at levels larger than this foregrounds removal shows a greater signal recovery. △ Less

Submitted 22 January, 2016; v1 submitted 20 August, 2014; originally announced August 2014.

Comments: 14 pages, 10 figures, accepted by MNRAS

Journal ref: MNRAS, Volume 458, Issue 3, p.2928-2939, May 2016

The VLA-COSMOS Survey: V. 324 MHz continuum observations

Authors: Vernesa Smolcic, Paolo Ciliegi, Vibor Jelic, Marco Bondi, Eva Schinnerer, Chris L. Carilli, Dominik A. Riechers, Mara Salvato, Alen Brkovic, Peter Capak, Olivier Ilbert, Alexander Karim, Henry McCracken, Nick Z. Scoville

Abstract: We present 90 cm VLA imaging of the COSMOS field, comprising a circular area of 3.14 square degrees at 8.0"x6.0" angular resolution with an average rms of 0.5 mJy/beam. The extracted catalog contains 182 sources (down to 5.5sigma), 30 of which are multi-component sources. Using Monte Carlo artificial source simulations we derive the completeness of the catalog, and we show that our 90 cm source co… ▽ More We present 90 cm VLA imaging of the COSMOS field, comprising a circular area of 3.14 square degrees at 8.0"x6.0" angular resolution with an average rms of 0.5 mJy/beam. The extracted catalog contains 182 sources (down to 5.5sigma), 30 of which are multi-component sources. Using Monte Carlo artificial source simulations we derive the completeness of the catalog, and we show that our 90 cm source counts agree very well with those from previous studies. Using X-ray, NUV-NIR and radio COSMOS data to investigate the population mix of our 90 cm radio sample, we find that our sample is dominated by active galactic nuclei (AGN). The average 90-20 cm spectral index (S_nu~nu**alpha, where S_nu is the flux density at frequency nu, and alpha the spectral index) of our 90 cm selected sources is -0.70, with an interquartile range of -0.90 to -0.53. Only a few ultra-steep-spectrum sources are present in our sample, consistent with results in the literature for similar fields. Our data do not show clear steepening of the spectral index with redshift. Nevertheless, our sample suggests that sources with spectral indices steeper than -1 all lie at z>1, in agreement with the idea that ultra-steep-spectrum radio sources may trace intermediate-redshift galaxies (z>1). △ Less

Submitted 17 July, 2014; originally announced July 2014.

Comments: 10 pages, 12 figures, accepted for publication in MNRAS

Lunar occultation of the diffuse radio sky: LOFAR measurements between 35 and 80 MHz

Authors: H. K. Vedantham, L. V. E. Koopmans, A. G. de Bruyn, S. J. Wijnholds, M. Brentjens, F. B. Abdalla, K. M. B. Asad, G. Bernardi, S. Bus, E. Chapman, B. Ciardi, S. Daiboo, E. R. Fernandez, A. Ghosh, G. Harker, V. Jelic, H. Jensen, S. Kazemi, P. Lambropoulos, O. Martinez-Rubi, G. Mellema, M. Mevius, A. R. Offringa, V. N. Pandey, A. H. Patil , et al. (69 additional authors not shown)

Abstract: We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<ν<80$ MHz since it is `colder' than the diffuse Galactic background i… ▽ More We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<ν<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$). △ Less

Submitted 16 July, 2014; originally announced July 2014.

Comments: 15 pages, 11 figures, 1 table

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

Authors: V. Jelic, A. G. de Bruyn, M. Mevius, F. B. Abdalla, K. M. B. Asad, G. Bernardi, M. A. Brentjens, S. Bus, E. Chapman, B. Ciardi, S. Daiboo, E. R. Fernandez, A. Ghosh, G. Harker, H. Jensen, S. Kazemi, L. V. E. Koopmans, P. Labropoulos, O. Martinez-Rubi, G. Mellema, A. R. Offringa, V. N. Pandey, A. H. Patil, R. M. Thomas, H. K. Vedantham , et al. (84 additional authors not shown)

Abstract: This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes.… ▽ More This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%. △ Less

Submitted 8 July, 2014; originally announced July 2014.

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

Journal ref: A&A 568, A101 (2014)