BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. Zaroubi

S. Zaroubi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.IMastro-ph.GAastro-phastro-ph.HEDatabasesphysics.ao-ph

Catalog footprint

What is connected

35works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2022arXiv

Statistical analysis of the causes of excess variance in the 21 cm signal power spectra obtained with the Low-Frequency Array

The detection of the 21 cm signal of neutral hydrogen from the Epoch of Reionization (EoR) is challenging due to bright foreground sources, radio frequency interference (RFI), the ionosphere, and instrumental effects. Even after correcting for these effects in the calibration step and applying foreground removal techniques, the remaining residuals in the observed 21 cm power spectra are still above the thermal noise, which is referred to as the "excess variance." We study potential causes of this excess variance based on 13 nights of data obtained with the Low-Frequency Array (LOFAR). We focused on the impact of gain errors, the sky model, and ionospheric effects on the excess variance by correlating the relevant parameters such as the gain variance over time or frequency, local sidereal time (LST), diffractive scale, and phase structure-function slope with the level of excess variance. Our analysis shows that excess variance has an LST dependence, which is related to the power from the sky. And the simulated Stokes I power spectra from bright sources and the excess variance show a similar progression over LST with the minimum power appearing at LST bin 6h to 9h. This LST dependence is also present in sky images of the residual Stokes I of the observations. In very-wide sky images, we demonstrate that the extra power comes exactly from the direction of bright and distant sources Cassiopeia A and Cygnus A with the array beam patterns. These results suggest that the level of excess variance in the 21 cm signal power spectra is related to sky effects and, hence, it depends on LST. In particular, very bright and distant sources such as Cassiopeia A and Cygnus A can dominate the effect. This is in line with earlier studies and offers a path forward toward a solution since the correlation between the sky-related effects and the excess variance is non-negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Probing Ionospheric Structures using the LOFAR radio telescope

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Galactic interstellar filaments as probed by LOFAR and Planck

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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$).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Initial deep LOFAR observations of Epoch of Reionization windows: I. The North Celestial Pole

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

LEDDB: LOFAR Epoch of Reionization Diagnostic Database

One of the key science projects of the Low-Frequency Array (LOFAR) is the detection of the cosmological signal coming from the Epoch of Reionization (EoR). Here we present the LOFAR EoR Diagnostic Database (LEDDB) that is used in the storage, management, processing and analysis of the LOFAR EoR observations. It stores referencing information of the observations and diagnostic parameters extracted from their calibration. This stored data is used to ease the pipeline processing, monitor the performance of the telescope and visualize the diagnostic parameters which facilitates the analysis of the several contamination effects on the signals. It is implemented with PostgreSQL and accessed through the psycopg2 python module. We have developed a very flexible query engine, which is used by a web user interface to access the database, and a very extensive set of tools for the visualization of the diagnostic parameters through all their multiple dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

LOFAR insights into the epoch of reionization from the cross power spectrum of 21cm emission and galaxies

Using a combination of N-body simulations, semi-analytic models and radiative transfer calculations, we have estimated the theoretical cross power spectrum between galaxies and the 21cm emission from neutral hydrogen during the epoch of reionization. In accordance with previous studies, we find that the 21cm emission is initially correlated with halos on large scales (> 30 Mpc), anti-correlated on intermediate (~ 5 Mpc), and uncorrelated on small (< 3 Mpc) scales. This picture quickly changes as reionization proceeds and the two fields become anti-correlated on large scales. The normalization of the cross power spectrum can be used to set constraints on the average neutral fraction in the intergalactic medium and its shape can be a tool to study the topology of reionization. When we apply a drop-out technique to select galaxies and add to the 21cm signal the noise expected from the LOFAR telescope, we find that while the normalization of the cross power spectrum remains a useful tool for probing reionization, its shape becomes too noisy to be informative. On the other hand, for a Lyalpha Emitter (LAE) survey both the normalization and the shape of the cross power spectrum are suitable probes of reionization. A closer look at a specific planned LAE observing program using Subaru Hyper-Suprime Cam reveals concerns about the strength of the 21cm signal at the planned redshifts. If the ionized fraction at z ~ 7 is lower that the one estimated here, then using the cross power spectrum may be a useful exercise given that at higher redshifts and neutral fractions it is able to distinguish between two toy models with different topologies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

LOFAR: The LOw-Frequency ARray

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Metals and ionizing photons from dwarf galaxies

We estimate the potential contribution of M < 10^9 Msun dwarf galaxies to the reionization and early metal-enrichment of the Milky Way environment, or circum-Galactic Medium. Our approach is to use the observed properties of ancient stars (> 12 Gyr old) measured in nearby dwarf galaxies to characterize the star-formation at high-z. We use a merger-tree model for the build-up of the Milky Way, which self-consistently accounts for feedback processes, and which is calibrated to match the present-day properties of the Galaxy and its dwarf satellites. We show that the high-z analogues of nearby dwarf galaxies can produce the bulk of ionizing radiation (>80%) required to reionize the Milky Way environment. Our fiducial model shows that the gaseous environment can be 50% reionized at z ~ 8 by galaxies with 10^7 Msun < M < 10^8 Msun. At later times, radiative feedback stops the star-formation in these small systems, and reionization is completed by more massive dwarf galaxies by z_rei = 6.4\pm 0.5. The metals ejected by supernova-driven outflows from M < 10^9 Msun dwarf galaxies almost uniformly fill the Milky Way environment by z ~ 5, enriching it to Z ~ 2 10^-2 Zsun. At z ~ 2 these early metals are still found to represent ~ 50% of the total mass of heavy elements in the circum-Galactic Medium.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radio Interferometric Calibration via Ordered-Subsets Algorithms: OS-LS and OS-SAGE calibrations

The main objective of this work is to accelerate the Maximum-Likelihood (ML) estimation procedure in radio interferometric calibration. We introduce the OS-LS and the OS-SAGE radio interferometric calibration methods, as a combination of the Ordered-Subsets (OS) method with the Least-Squares (LS) and Space Alternating Generalized Expectation maximization (SAGE) calibration techniques, respectively. The OS algorithm speeds up the ML estimation and achieves nearly the same level of accuracy of solutions as the one obtained by the non-OS methods. We apply the OS-LS and OS-SAGE calibration methods to simulated observations and show that these methods have a much higher convergence rate relative to the conventional LS and SAGE techniques. Moreover, the obtained results show that the OS-SAGE calibration technique has a superior performance compared to the OS-LS calibration method in the sense of achieving more accurate results while having significantly less computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The brightness and spatial distributions of terrestrial radio sources

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Imaging neutral hydrogen on large-scales during the Epoch of Reionization with LOFAR

The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (\sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it hard to disentangle the two on such large scales. The poor sensitivity on small scales on the one hand, and the foregrounds effect on large scales on the other hand, make direct imaging of the Epoch of Reionization of the Universe very difficult, and detection of the signal therefore is expected to be statistical.Despite these hurdles, in this paper we argue that for many reionization scenarios low resolution images could be obtained from the expected data. This is because at the later stages of the process one still finds very large pockets of neutral regions in the IGM, reflecting the clustering of the large-scale structure, which stays strong up to scales of \sim 120 comoving Mpc/h (\sim 1 degree). The coherence of the emission on those scales allows us to reach sufficient S/N (\sim 3) so as to obtain reionization 21 cm images. Such images will be extremely valuable for answering many cosmological questions but above all they will be a very powerful tool to test our control of the systematics in the data. The existence of this typical scale (\sim 120 comoving Mpc/h) also argues for designing future EoR experiments, e.g., with SKA, with a field of view of at least 4 degree.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planck intermediate results. VI: The dynamical structure of PLCKG214.6+37.0, a Planck discovered triple system of galaxy clusters

The survey of galaxy clusters performed by Planck through the Sunyaev-Zeldovich effect has already discovered many interesting objects, thanks to the whole coverage of the sky. One of the SZ candidates detected in the early months of the mission near to the signal to noise threshold, PLCKG214.6+37.0, was later revealed by XMM-Newton to be a triple system of galaxy clusters. We have further investigated this puzzling system with a multi-wavelength approach and we present here the results from a deep XMM-Newton re-observation. The characterisation of the physical properties of the three components has allowed us to build a template model to extract the total SZ signal of this system with Planck data. We partly reconciled the discrepancy between the expected SZ signal from X-rays and the observed one, which are now consistent at less than 1.2 sigma. We measured the redshift of the three components with the iron lines in the X-ray spectrum, and confirmed that the three clumps are likely part of the same supercluster structure. The analysis of the dynamical state of the three components, as well as the absence of detectable excess X-ray emission, suggest that we are witnessing the formation of a massive cluster at an early phase of interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planck intermediate results. X. Physics of the hot gas in the Coma cluster

We present an analysis of Planck satellite data on the Coma Cluster observed via the Sunyaev-Zeldovich effect. Planck is able, for the first time, to detect SZ emission up to r ~ 3 X R_500. We test previously proposed models for the pressure distribution in clusters against the azimuthally averaged data. We find that the Arnaud et al. universal pressure profile does not fit Coma, and that their pressure profile for merging systems provides a good fit of the data only at r<R_500: by r=2XR_500 it underestimates the observed y profile by a factor of ~2. This may indicate that at these larger radii either i) the cluster SZ emission is contaminated by unresolved SZ sources along the line of sight or ii) the pressure profile of Coma is higher at r>R_500 than the mean pressure profile predicted by the simulations. The Planck image shows significant local steepening of the y profile in two regions about half a degree to the west and to the south-east of the cluster centre. These features are consistent with the presence of shock fronts at these radii, and indeed the western feature was previously noticed in the ROSAT PSPC mosaic as well as in the radio. Using Planck y profiles extracted from corresponding sectors we find pressure jumps of 4.5+0.4-0.2 and 5.0+1.3-0.1 in the west and southeast, respectively. Assuming Rankine-Hugoniot pressure jump conditions, we deduce that the shock waves should propagate with Mach number M_w=2.03+0.09-0.04 and M_se=2.05+0.25-0.02 in the West and Southeast, respectively. Finally, we find that the y and radio-synchrotron signals are quasi-linearly correlated on Mpc scales with small intrinsic scatter. This implies either that the energy density of cosmic-ray electrons is relatively constant throughout the cluster, or that the magnetic fields fall off much more slowly with radius than previously thought.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Post-correlation filtering techniques for off-axis source and RFI removal

Techniques to improve the data quality of interferometric radio observations are considered. Fundaments of fringe frequencies in the uv-plane are discussed and filters are used to attenuate radio-frequency interference (RFI) and off-axis sources. Several new applications of filters are introduced and tested. A low-pass filter in time and frequency direction on single baseline data is successfully used to lower the noise in the area of interest and to remove sidelobes coming from unmodelled off-axis sources and RFI. Related side effects of data integration, averaging and gridding are analysed, and shown to be able to cause ghosts and an increase in noise, especially when using long baselines or interferometric elements that have a large field of view. A novel projected fringe low-pass filter is shown to be potentially useful for first order source separation. Initial tests show that the filters can be several factors faster compared to common source separation techniques such as peeling and a variant of peeling that is currently being tested on LOFAR observations called "demixed peeling". Further testing is required to support the performance of the filters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Prospects for detecting the 21cm forest from the diffuse intergalactic medium with LOFAR

We discuss the feasibility of the detection of the 21cm forest in the diffuse IGM with the radio telescope LOFAR. The optical depth to the 21cm line has been derived using simulations of reionization which include detailed radiative transfer of ionizing photons. We find that the spectra from reionization models with similar total comoving hydrogen ionizing emissivity but different frequency distribution look remarkably similar. Thus, unless the reionization histories are very different from each other (e.g. a predominance of UV vs. x-ray heating) we do not expect to distinguish them by means of observations of the 21cm forest. Because the presence of a strong x-ray background would make the detection of 21cm line absorption impossible, the lack of absorption could be used as a probe of the presence/intensity of the x-ray background and the thermal history of the universe. Along a random line of sight LOFAR could detect a global suppression of the spectrum from z>12, when the IGM is still mostly neutral and cold, in contrast with the more well-defined, albeit broad, absorption features visible at lower redshift. Sharp, strong absorption features associated with rare, high density pockets of gas could be detected also at z~7 along preferential lines of sight.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The LOFAR radio environment

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lyman alpha emission from the first galaxies: Implications of UV backgrounds and the formation of molecules

The Lyman alpha line is a robust tracer of high redshift galaxies. We present estimates of Lyman alpha emission from a protogalactic halo illuminated by UV background radiation fields with various intensities. For this purpose, we performed cosmological hydrodynamics simulations with the adaptive mesh refinement code FLASH, including a detailed network for primordial chemistry,comprising the formation of primordial molecules, a multi-level model for the hydrogen atom as well as the photo-ionization and photo-dissociation processes in a UV background. We find that the presence of a background radiation field J_21 excites the emission of Lyman alpha photons, increasing the Lyman alpha luminosity up to two orders of magnitude. For a halo of \sim 10^10 M_sun, we find that a maximum flux of 5 \times 10^-15 erg cm^-2 s^-1 is obtained for J21 \times f_esc = 0.1, where f_esc is the escape fraction of the ionizing radiation. Depending on the environmental conditions, the flux may vary by three orders of magnitude. For J_21 \times f_esc > 0.1 the Lyman alpha luminosity decreases as the atomic hydrogen abundance becomes rather small. The fluxes derived here can be probed using Subaru and the upcoming James Webb Space Telescope. The emission of Lyman alpha photons is extended and comes from the envelope of the halo rather than its core. In the center of the halo, line trapping becomes effective above columns of 10^22 cm^-2 and suppresses the emission of Lyman alpha. In addition, cooling by primordial molecules may decrease the gas temperature in the central region, which further reduces Lyman alpha emission. In the central core, H_2 is photo-dissociated for a background flux of J_21 \geq 1000. For weaker radiation fields, i.e. J_21 < 0.1, H_2 and HD cooling are particularly strong in the center of the halo, leading to gas temperatures as low as \sim 100 K.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lyman alpha emission from the first galaxies: Signatures of accretion and infall in the presence of line trapping

The formation of the first galaxies is accompanied by large accretion flows and virialization shocks, during which the gas is shock-heated to temperatures of $\sim10^4$ K, leading to potentially strong fluxes in the Lyman alpha line. Indeed, a number of Lyman alpha blobs has been detected at high redshift. In this letter, we explore the origin of such Lyman alpha emission using cosmological hydrodynamical simulations that include a detailed model of atomic hydrogen as a multi-level atom and the effects of line trapping with the adaptive mesh refinement code FLASH. We see that baryons fall into the center of a halo through cold streams of gas, giving rise to a Lyman alpha luminosity of at least $\rm 10^{44} erg s^{-1}$ at $\rm z=4.7$, similar to observed Lyman alpha blobs. We find that a Lyman alpha flux of $\rm 5.0\times 10^{-17} erg cm^{-2} s^{-1}$ emerges from the envelope of the halo rather than its center, where the photons are efficiently trapped. Such emission can be probed in detail with the upcoming James Webb Space Telescope (JWST) and will constitute an important probe of gas infall and accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radio Interferometric Calibration Using The SAGE Algorithm

The aim of the new generation of radio synthesis arrays such as LOFAR and SKA is to achieve much higher sensitivity, resolution and frequency coverage than what is available now, especially at low frequencies. To accomplish this goal, the accuracy of the calibration techniques used is of considerable importance. Moreover, since these telescopes produce huge amounts of data, speed of convergence of calibration is a major bottleneck. The errors in calibration are due to system noise (sky and instrumental) as well as the estimation errors introduced by the calibration technique itself, which we call solver noise. We define solver noise as the distance between the optimal solution (the true value of the unknowns, uncorrupted by the system noise) and the solution obtained by calibration. We present the Space Alternating Generalized Expectation Maximization (SAGE) calibration technique, which is a modification of the Expectation Maximization algorithm, and compare its performance with the traditional Least Squares calibration based on the level of solver noise introduced by each technique. For this purpose, we develop statistical methods that use the calibrated solutions to estimate the level of solver noise. The SAGE calibration algorithm yields very promising results both in terms of accuracy and speed of convergence. The comparison approaches we adopt introduce a new framework for assessing the performance of different calibration schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A LOFAR RFI detection pipeline and its first results

Radio astronomy is entering a new era with new and future radio observatories such as the Low Frequency Array and the Square Kilometer Array. We describe in detail an automated flagging pipeline and evaluate its performance. With only a fraction of the computational cost of correlation and its use of the previously introduced SumThreshold method, it is found to be both fast and unrivalled in its high accuracy. The LOFAR radio environment is analysed with the help of this pipeline. The high time and spectral resolution of LOFAR have resulted in an observatory where only a few percent of the data is lost due to RFI.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Foregrounds for observations of the cosmological 21 cm line: II. Westerbork observations of the fields around 3C196 and the North Celestial Pole

In the coming years a new insight into galaxy formation and the thermal history of the Universe is expected to come from the detection of the highly redshifted cosmological 21 cm line. The cosmological 21 cm line signal is buried under Galactic and extragalactic foregrounds which are likely to be a few orders of magnitude brighter. Strategies and techniques for effective subtraction of these foreground sources require a detailed knowledge of their structure in both intensity and polarization on the relevant angular scales of 1-30 arcmin. We present results from observations conducted with the Westerbork telescope in the 140-160 MHz range with 2 arcmin resolution in two fields located at intermediate Galactic latitude, centred around the bright quasar 3C196 and the North Celestial Pole. They were observed with the purpose of characterizing the foreground properties in sky areas where actual observations of the cosmological 21 cm line could be carried out. The polarization data were analysed through the rotation measure synthesis technique. We have computed total intensity and polarization angular power spectra. Total intensity maps were carefully calibrated, reaching a high dynamic range, 150000:1 in the case of the 3C196 field. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Post-correlation radio frequency interference classification methods

We describe and compare several post-correlation radio frequency interference classification methods. As data sizes of observations grow with new and improved telescopes, the need for completely automated, robust methods for radio frequency interference 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% recognition and an approximately 0.1% 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Impact of Lyman alpha Trapping on the Formation of Primordial Objects

Numerous cosmological simulations have been performed to study the formation of the first objects. We present the results of high resolution 3-D cosmological simulations of primordial objects formation using the adaptive mesh refinement code FLASH by including in an approximate manner the radiative transfer effects of Lyman alpha photons. We compare the results of a Lyman alpha trapping case inside gas clouds with atomic and molecular hydrogen cooling cases.The principal objective of this research is to follow the collapse of a zero metallicity halo with an effective equation of state (that accounts for the trapping) and to explore the fate of a halo in each of the three cases, specifically, the impact of thermodynamics on fragmentation of halos.Our results show that in the case of Lyman alpha trapping, fragmentation is halted and a massive object is formed at the center of a halo. The temperature of the gas remains well above $10^{4}$ K and the halo is not able to fragment to stellar masses. In the atomic cooling case, gas collapses into one or two massive clumps in contrast to the Lyman alpha trapping case. For the molecular hydrogen cooling case, gas cools efficiently and fragments.The formation of massive primordial objects is thus strongly dependent on the thermodynamics of the gas. A salient feature of our results is that for the formation of massive objects, e.g. intermediate mass black holes, feedback effects are not required to suppress $H_{2}$ cooling, as molecular hydrogen is collisionally dissociated at temperatures higher than $10^{4}$ K as a consequence of Lyman alpha trapping.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Foregrounds for observations of the cosmological 21 cm line: I. First Westerbork measurements of Galactic emission at 150 MHz in a low latitude field

We present the first results from a series of observations conducted with the Westerbork telescope in the 140--160 MHz range with a 2 arcmin resolution aimed at characterizing the properties of the foregrounds for epoch of reionization experiments. For the first time we have detected fluctuations in the Galactic diffuse emission on scales greater than 13 arcmin at 150 MHz, in the low Galactic latitude area known as Fan region. Those fluctuations have an $rms$ of 14 K. The total intensity power spectrum shows a power--law behaviour down to $\ell \sim 900$ with slope $β^I_\ell = -2.2 \pm 0.3$. The detection of diffuse emission at smaller angular scales is limited by residual point sources. We measured an $rms$ confusion noise of $\sim$3 mJy beam$^{-1}$. Diffuse polarized emission was also detected for the first time at this frequency. The polarized signal shows complex structure both spatially and along the line of sight. The polarization power spectrum shows a power--law behaviour down to $\ell \sim 2700$ with slope $β^P_\ell = -1.65 \pm 0.15$. The $rms$ of polarization fluctuations is 7.2 K on 4 arcmin scales. By extrapolating the measured spectrum of total intensity emission, we find a contamination on the cosmological signal of $δT= \sqrt{\ell (\ell+1) C^I_\ell / 2π} \sim 5.7$ K on 5 arcmin scales and a corresponding $rms$ value of $\sim$18.3 K at the same angular scale. The level of the polarization power spectrum is $δT \sim 3.3$ K on 5 arcmin scales. Given its exceptionally bright polarized signal, the Fan region is likely to represent an upper limit on the sky brightness at moderate and high Galactic latitude.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

On the Non-Gaussianity Observed in the COBE-DMR Sky Maps

In this paper we pursue the origin of the non-Gaussianity determined by a bispectrum analysis of the COBE-DMR 4-year sky maps. The robustness of the statistic is demonstrated by the rebinning of the data into 12 coordinate systems. By computing the bispectrum statistic as a function of various data partitions - by channel, frequency, and time interval, we show that the observed non-Gaussian signal is driven by the 53 GHz data. This frequency dependence strongly rejects the hypothesis that the signal is cosmological in origin. A jack-knife analysis of the coadded 53 and 90 GHz sky maps reveals those sky pixels to which the bispectrum statistic is particularly sensitive. We find that by removing data from the 53 GHz sky maps for periods of time during which a known systematic effect perturbs the 31 GHz channels, the amplitudes of the bispectrum coefficients become completely consistent with that expected for a Gaussian sky. We conclude that the non-Gaussian signal detected by the normalised bispectrum statistic in the publicly available DMR sky maps is due to a systematic artifact. The impact of removing the affected data on estimates of the normalisation of simple models of cosmological anisotropy is negligible.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

From Local Velocities to Microwave Background

The mass density field as extracted from peculiar velocities in our cosmological neighborhood is mapped back in time to the CMB in two ways. First, the density power spectrum ($P_k$) is translated into a temperature angular power spectrum of sub-degree resolution ($C_l$) and compared to observations. Second, the local density field is translated into a temperature map in a patch on the last-scattering surface of a distant observer. A likelihood analysis of the Mark III peculiar velocity data have constrained the range of parameters for $P_k$ within the family of COBE-normalized CDM models (Zaroubi et al 1996), favoring a slight tilt in the initial spectrum, $n<1$. The corresponding range of $C_l$'s is plotted against current observations, indicating that the CMB data can tighten the constraints further: only models with ``small'' tilt ($n\sim 0.9$) and ``high'' baryonic content ($Ω_b \sim 0.1$) could survive the two data sets simultaneously. The local mass density field that has been recovered from the velocities via a Wiener method is convolved with a Boltzmann calculation to recover $10'$ resolution temperature maps as viewed from different directions. The extent of the CMB patch and the amplitude of fluctuations depend on the choice of cosmological parameters, e.g., the local $100\hmpc$ sphere corresponds to $90'$ to $30'$ at the CMB for $Ω$ between 1 and 0 respectively. The phases of the temperature map are correlated with those of the density field, contrary to the contribution of the Sachs-Wolfe effect alone. This correlation suggests the possibility of an inverse reconstruction of the underlying density field from CMB data with interesting theoretical implications.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Large Scale Power Spectrum from Peculiar Velocities Via Likelihood Analysis

The power spectrum (PS) of {\it mass} density fluctuations, independent of ``biasing", is estimated from the Mark3 Catalog of Peculiar Velocities of galaxies using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the posterior probability of the model given the data. The method has been tested using detailed mock catalogs. We use generalized CDM models with and without COBE normalization. The robust result for all the models tested is $P(k) Ω^{1.2} = (4.1 \pm 0.7)x10^3 (Mpc h^{-1})^3$ at $k=0.1 h Mpc^{-1}$, with the peak constrained to the range $0.03 < k < 0.06 h Mpc^{-1}$. It is consistent with a direct computation of the PS (Kolatt \& Dekel 1996). When compared to galaxy-density surveys, the implied values for $β$ ($\equiv Ω^{0.6}/b$) are of order unity to within 25\%. A $Γ$-shape model, free of COBE normalization, is constrained by the velocity data to $Γ=0.5\pm 0.15$ and $σ_8 Ω^{0.6}=0.85\pm 0.1$. Within the family of COBE-normalized CDM models, the best tilted model ($Ω=1$, $n \approx 0.84 h_{50}^{-0.65}$) and the best open model ($n=1$, $Ω\approx 0.75 h_{50}^{-1.0}$) are more likely than the best $Λ$ model ($n=1$, $Λ=1-Ω$, $Ω\approx 0.70 h_{50}^{-1.2}$). The most likely CDM model with $Ω\leq 1$ is found to be of $Ω=1$, $h\approx 0.75$, and a {\it tilted} spectrum of $n= 0.8\pm 0.02$ with tensor fluctuations.The tightest constraint for the tilted-$Λ$ models is of the sort $Ωh_{50}^{1.2} n^ν= 0.7 \pm 0.08$, with $ν=3.8$ and $1.85$ with and without tensor fluctuations respectively.

0 citations0 reviewsNo code linkedOpen access