ABSTRACT
A new upper limit on the 21 cm signal power spectrum at a redshift of z ≈ 9.1 is presented, based on 141 h 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 $\Delta ^2_{21} \lt (73)^2\, \mathrm{mK^2}$ at $k = 0.075\, \mathrm{h\, cMpc^{-1}}$ is found, an improvement by a factor ≈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 radiofrequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power.
ABSTRACT
The ARCADE2 and LWA1 experiments have claimed an excess over the cosmic microwave background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio ...background is between 0.1 per cent and 22 per cent of the CMB at 1.42 GHz, it could explain the tentative EDGES low-band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen. We use the upper limit on the 21-cm signal from the Epoch of Reionization (z = 9.1) based on 141 h of observations with LOFAR to evaluate the contribution of the high-redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming haloes, we find (at 95 per cent CL) that the cosmological radio background can be at most 9.6 per cent of the CMB at 1.42 GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of 0.1–9.6 per cent over the CMB (at 1.42 GHz) is still allowed and could explain the EDGES low-band detection. We also constrain the thermal and ionization state of the gas at z = 9.1, and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES high-band data, LOFAR data constrain scenarios with inefficient X-ray sources, and cases where the Universe was ionized by stars in massive haloes only.
ABSTRACT
We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift ≈ 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 ≳0.13 and a distribution of the ionized regions with a characteristic size ≳ 8 h−1 comoving megaparsec (Mpc) and a full width at half-maximum (FWHM) ≳16 h−1 Mpc is ruled out. For an IGM with a uniform spin temperature TS ≳ 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 ≲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 ≲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.
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 frequency range) using 14 h 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{{\ \rm per\ cent}}$ in Stokes / noise compared to Stokes V for the two observed fields, which decorrelates on ≳12 s 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 $\Delta _{21}^2 \lt (14561\, \text{mK})^2$ at $k\sim 0.038\, h\, \text{cMpc}^{-1}$ and $\Delta _{21}^2 \lt (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.
We present the first limits on the Epoch of Reionization 21 cm H i power spectra, in the redshift range z = 7.9-10.6, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total, 13.0 hr ...of data were used from observations centered on the North Celestial Pole. After subtraction of the sky model and the noise bias, we detect a non-zero (1- ) excess variance and a best 2- upper limit of 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 nonlinear effects, (iii) noise and ionosphere-induced gain errors, or a combination thereof. Further analyses of the excess variance will be discussed in forthcoming publications.
The distances to individual wind-driven bubbles such as Planetary Nebulae (PNe) can be determined using expansion parallaxes: the angular expansion velocity in the sky is compared to the radial ...velocity of gas measured spectroscopically. Since the one is a pattern velocity, and the other a matter velocity, these are not necessarily the same. Using the jump conditions for both shocks and ionization fronts, I show that for typical PNe the pattern velocity is 20 to 30% larger than the material velocity, and the derived distances are therefore typically 20 to 30% too low. I present some corrected distances and suggest approaches to be used when deriving distances using expansion parallaxes.
Aims. This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C 196 field, one of the primary fields of the LOFAR-Epoch of ...Reionization key science project. Methods. 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. Results. The brightness temperature of the detected Galactic emission is 5−15 K in polarized intensity and covers the range from –3 to +8 rad m-2 in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of +0.5 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 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. Conclusions. The mean line-of-sight magnetic field component, B∥, is determined to be 0.3 ± 0.1 μG and its spatial variation across the 3C 196 field is 0.1 μ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 (neB∥> 6.2 cm-3μG) compared to its surroundings.
ABSTRACT
We investigate systematic effects in direction-dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key ...Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of 50–250 λ. We show that suppression of faint signals can effectively be avoided by calibrating these short baselines using only the longer baselines. However, this approach causes an excess variance on the short baselines due to small gain errors induced by overfitting during calibration. We apply a regularized expectation–maximization algorithm with consensus optimization (sagecal-co) to real data with simulated signals to show that overfitting can be largely mitigated by penalising spectrally non-smooth gain solutions during calibration. This reduces the excess power with about a factor of 4 in the simulations. Our results agree with earlier theoretical analysis of this bias-variance trade off and support the gain-calibration approach to the LOFAR 21-cm signal data.
We present the results of radiation-magnetohydrodynamic simulations of the formation and expansion of H ii regions and their surrounding photodissociation regions (PDRs) in turbulent, magnetized, ...molecular clouds on scales of up to 4 pc. We include the effects of ionizing and non-ionizing ultraviolet radiation and X-rays from population synthesis models of young star clusters. For all our simulations we find that the H ii region expansion reduces the disordered component of the magnetic field, imposing a large-scale order on the field around its border, with the field in the neutral gas tending to lie along the ionization front, while the field in the ionized gas tends to be perpendicular to the front. The highest pressure-compressed neutral and molecular gas is driven towards approximate equipartition between thermal, magnetic and turbulent energy densities, whereas lower pressure neutral/molecular gas bifurcates into, on the one hand, quiescent, magnetically dominated regions and, on the other hand, turbulent, demagnetized regions. The ionized gas shows approximate equipartition between thermal and turbulent energy densities, but with magnetic energy densities that are 1-3 orders of magnitude lower. A high velocity dispersion (∼8 km s−1) is maintained in the ionized gas throughout our simulations, despite the mean expansion velocity being significantly lower. The magnetic field does not significantly brake the large-scale H ii region expansion on the length and time-scales accessible to our simulations, but it does tend to suppress the smallest scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the H ii region. However, the relative luminosity of ionizing and non-ionizing radiation has a much larger influence than the presence or absence of the magnetic field. When the star cluster radiation field is relatively soft (as in the case of a lower mass cluster, containing an earliest spectral type of B0.5), then fragmentation is less vigorous and a thick, relatively smooth PDR forms.
We present the first large-scale radiative transfer simulations of cosmic reionization, in a simulation volume of (100 h
−1 Mpc)3. This is more than a two orders of magnitude improvement over ...previous simulations. We achieve this by combining the results from extremely large, cosmological, N-body simulations with a new, fast and efficient code for 3D radiative transfer, c
2-ray, which we have recently developed. These simulations allow us to do the first numerical studies of the large-scale structure of reionization which at the same time, and crucially, properly take account of the dwarf galaxy ionizing sources which are primarily responsible for reionization. In our realization, reionization starts around z∼ 21, and final overlap occurs by z∼ 11. The resulting electron-scattering optical depth is in good agreement with the first-year Wilkinson Microwave Anisotropy Probe (WMAP) polarization data. We show that reionization clearly proceeded in an inside-out fashion, with the high-density regions being ionized earlier, on average, than the voids. Ionization histories of smaller-size (5-10 comoving Mpc) subregions exabit a large scatter about the mean and do not describe the global reionization history well. This is true even when these subregions are at the mean density of the universe, which shows that small-box simulations of reionization have little predictive power for the evolution of the mean ionized fraction. The minimum reliable volume size for such predictions is ∼30 Mpc. We derive the power spectra of the neutral, ionized and total gas density fields and show that there is a significant boost of the density fluctuations in both the neutral and the ionized components relative to the total at arcmin and larger scales. We find two populations of H ii regions according to their size, numerous, mid-sized (∼10-Mpc) regions and a few, rare, very large regions tens of Mpc in size. Thus, local overlap on fairly large scales of tens of Mpc is reached by z∼ 13, when our volume is only about 50 per cent ionized, and well before the global overlap. We derive the statistical distributions of the ionized fraction and ionized gas density at various scales and for the first time show that both distributions are clearly non-Gaussian. All these quantities are critical for predicting and interpreting the observational signals from reionization from a variety of observations like 21-cm emission, Lyα emitter statistics, Gunn-Peterson optical depth and small-scale cosmic microwave background secondary anisotropies due to patchy reionization.