G. Bernardi

Foreground simulations for the LOFAR - Epoch of Reionization Experiment

Future high-redshift 21-cm experiments will suffer from a high degree of contamination, due both to astrophysical foregrounds and to non-astrophysical and instrumental effects. In order to reliably extract the cosmological signal from the observed data, it is essential to understand very well all data components and their influence on the extracted signal. Here we present simulated astrophysical foregrounds data cubes and discuss their possible statistical effects on the data. The foreground maps are produced assuming 5 degrees x 5 degrees windows that match those expected to be observed by the LOFAR epoch of reionization (EoR) key science project. We show that with the expected LOFAR-EoR sky and receiver noise levels, which amount to approximate to 52 mK at 150 MHz after 400 h of total observing time, a simple polynomial fit allows a statistical reconstruction of the signal. We also show that the polynomial fitting will work for maps with realistic yet idealized instrument response, i.e. a response that includes only a uniform uv coverage as a function of frequency and ignores many other uncertainties. Polarized Galactic synchrotron maps that include internal polarization and a number of Faraday screens along the line of sight are also simulated. The importance of these stems from the fact that the LOFAR instrument, in common with all current interferometric EoR experiments, has an instrumentally polarized response.

Post-correlation radio frequency interference classification methods

We describe and compare several post-correlation radio frequency interference (RFI) classification methods. As data sizes of observations grow with new and improved telescopes, the need for completely automated, robust methods for RFI mitigation is pressing. We investigated several classification methods and find that, for the data sets we used, the most accurate among them is the SumThreshold method. This is a new method formed from a combination of existing techniques, including a new way of thresholding. This iterative method estimates the astronomical signal by carrying out a surface fit in the time-frequency plane. With a theoretical accuracy of 95 per cent recognition and an approximately 0.1 per cent false probability rate in simple simulated cases, the method is in practice as good as the human eye in finding RFI. In addition, it is fast, robust, does not need a data model before it can be executed and works in almost all configurations with its default parameters. The method has been compared using simulated data with several other mitigation techniques, including one based upon the singular value decomposition of the time-frequency matrix, and has shown better results than the rest.

Fast large‐scale reionization simulations

We present an efficient method to generate large simulations of the epoch of reionization without the need for a full three-dimensional radiative transfer code. Large dark-matter-only simulations are post-processed to produce maps of the redshifted 21-cm emission from neutral hydrogen. Dark matter haloes are embedded with sources of radiation whose properties are either based on semi-analytical prescriptions or derived from hydrodynamical simulations. These sources could either be stars or power-law sources with varying spectral indices. Assuming spherical symmetry, ionized bubbles are created around these sources, whose radial ionized fraction and temperature profiles are derived from a catalogue of one-dimensional radiative transfer experiments. In case of overlap of these spheres, photons are conserved by redistributing them around the connected ionized regions corresponding to the spheres. The efficiency with which these maps are created allows us to span the large parameter space typically encountered in reionization simulations. We compare our results with other, more accurate, three-dimensional radiative transfer simulations and find excellent agreement for the redshifts and the spatial scales of interest to upcoming 21-cm experiments. We generate a contiguous observational cube spanning redshift 6 to 12 and use these simulations to study the differences in the reionization histories between stars and quasars. Finally, the signal is convolved with the Low Frequency Array (LOFAR) beam response and its effects are analysed and quantified. Statistics performed on this mock data set shed light on possible observational strategies for LOFAR.

Calibrating high-precision Faraday rotation measurements for LOFAR and the next generation of low-frequency radio telescopes

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

Detection and extraction of signals from the epoch of reionization using higher-order one-point statistics

Detecting redshifted 21-cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, and we suggest that using the correct statistics may tease out signatures of reionization. We generate mock data cubes simulating the output of the Low Frequency Array (LOFAR) EoR experiment. These cubes combine realistic models for Galactic and extragalactic foregrounds and the noise with three different simulations of the CS. We fit out the foregrounds, which are smooth in the frequency direction, to produce residual images in each frequency band. We denoise these images and study the skewness of the one-point distribution in the images as a function of frequency. We find that, under sufficiently optimistic assumptions, we can recover the main features of the redshift evolution of the skewness in the 21-cm signal. We argue that some of these features – such as a dip at the onset of reionization, followed by a rise towards its later stages – may be generic, and give us a promising route to a statistical detection of reionization.

High Galactic latitude polarized emission at 1.4 GHz and implications for cosmic microwave background observations

We analyse the polarized emission at 1.4 GHz in a 3° x 3° area at high Galactic latitude (b ∼ -40°). The region, centred in (a = 5 h , δ = -49°), was observed with the Australia Telescope Compact Array (ATCA) radio-interferometer, whose 3-30 arcmin angular sensitivity range allows the study of scales appropriate for cosmic microwave background polarization (CMBP) investigations. The angular behaviour of the diffuse emission is analysed through the E- and B-mode angular power spectra. These follow a power law C X l lβ X with slopes β E = -1.97 ± 0.08 and β B = -1.98 ± 0.07. The emission is found to be approximately a factor 25 fainter than in Galactic plane regions. The comparison of the power spectra with other surveys indicates that this area is intermediate between strong and negligible Faraday rotation effects. A similar conclusion can be reached by analysing both the frequency and Galactic latitude behaviours of the diffuse Galactic emission of the 408-1411 MHz Leiden survey data. We present an analysis of the Faraday rotation effects on the polarized power spectra and find that the observed power spectra can be enhanced by a transfer of power from large to small angular scales. The extrapolation of the spectra to 32 and 90 GHz of the cosmic microwave background (CMB) window suggests that Galactic synchrotron emission leaves the CMBP E-mode uncontaminated at 32 GHz. The level of the contamination at 90 GHz is expected to be more than 4 orders of magnitude below the CMBP spectrum. Extrapolating to the relevant angular scales, this region also appears adequate for investigation of the CMBP B-modes for models with tensor-to-scalar fluctuation power ratio T/S ≥ 0.01. We also identify polarized point sources in the field, providing a nine object list, which is complete down to the polarized flux limit of S p lim = 2 mJy.

The Sky Polarization Observatory

Abstract The Sky Polarization Observatory (SPOrt) is an ASI-funded experiment specifically designed to measure the sky polarization at 22, 32 and 90 GHz, which was selected in 1997 by ESA to be flown on the International Space Station. Starting in 2006 and for at least 18 months, it will be taking direct and simultaneous measurements of the Stokes parameters Q and U at 660 sky pixels, with FWHM=7°. Due to development efforts over the past few years, the design specifications have been significantly improved with respect to the first proposal. Here we present an up-to-date description of the instrument, which now warrants a pixel sensitivity of 1.7 μK for the polarization of the cosmic background radiation, assuming two years of observations. We discuss SPOrt scientific goals in the light of WMAP results, in particular in connection with the emerging double-reionization cosmological scenario.

Probing atmospheric electric fields in thunderstorms through radio emission from cosmic-ray induced air showers

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

THE SPECTRAL VARIABILITY OF THE GHZ-PEAKED SPECTRUM RADIO SOURCE PKS 1718-649 AND A COMPARISON OF ABSORPTION MODELS

Using the new wideband capabilities of the Australia Telescope Compact Array (ATCA), we obtain spectra for PKS 1718-649, a well-known gigahertz-peaked spectrum radio source. The observations, between approximately 1 and 10 GHz over three epochs spanning approximately 21 months, reveal variability both above the spectral peak at ~3 GHz and below the peak. The combination of the low and high frequency variability cannot be easily explained using a single absorption mechanism, such as free-free absorption or synchrotron self-absorption. We find that the PKS 1718-649 spectrum and its variability are best explained by variations in the free-free optical depth on our line-of-sight to the radio source at low frequencies (below the spectral peak) and the adiabatic expansion of the radio source itself at high frequencies (above the spectral peak). The optical depth variations are found to be plausible when X-ray continuum absorption variability seen in samples of Active Galactic Nuclei is considered. We find that the cause of the peaked spectrum in PKS 1718-649 is most likely due to free-free absorption. In agreement with previous studies, we find that the spectrum at each epoch of observation is best fit by a free-free absorption model characterised by a power-law distribution of free-free absorbing clouds. This agreement is extended to frequencies below the 1 GHz lower limit of the ATCA by considering new observations with Parkes at 725 MHz and 199 MHz observations with the newly operational Murchison Widefield Array. These lower frequency observations argue against families of absorption models (both free-free and synchrotron self-absorption) that are based on simple homogenous structures.

The Parkes Galactic Meridian Survey: observations and CMB polarization foreground analysis

We present observations and cosmic microwave background (CMB) foreground analysis of the Parkes Galactic Meridian Survey, an investigation of the Galactic latitude behaviour of the polarized synchrotron emission at 2.3 GHz with the Parkes Radio Telescope. The survey consists of a 5 degrees wide strip along the Galactic meridian l = 254 degrees extending from the Galactic plane to the South Galactic pole. We identify three zones distinguished by polarized emission properties: the disc, the halo and a transition region connecting them. The halo section lies at latitudes vertical bar b vertical bar > 40 degrees and has weak and smooth polarized emission mostly at large scale with steep angular power spectra of median slope beta(med) similar to -2.6. The disc region covers the latitudes vertical bar b vertical bar <20 degrees and has a brighter, more complex emission dominated by the small scales with flatter spectra of median slope beta(med) = -1.8. The transition region has steep spectra as in the halo, but the emission increases towards the Galactic plane from halo to disc levels. The change of slope and emission structure at b similar to -20 degrees is sudden, indicating a sharp disc-halo transition. The whole halo section is just one environment extended over 50 degrees with very low emission which, once scaled to 70 GHz, is equivalent to the CMB B-mode emission for a tensor-toscalar perturbation power ratio r(halo) = (3.3 +/- 0.4) x 10(-3). Applying a conservative cleaning procedure, we estimate an r detection limit of delta r similar to 2 x 10(-3) at 70 GHz (3 sigma confidence limit) and, assuming a dust polarization fraction of <12 per cent, dr similar to 1 x 10(-2) at 150 GHz. The 150-GHz limit matches the goals of planned sub-orbital experiments, which can therefore be conducted at this high frequency. The 70-GHz limit is close to the goal of proposed next-generation space missions, which thus might not strictly require space-based platforms.