Astronomical wide-field imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new ...wide-field interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependences of casa's w-projection and our new imager are analysed and analytical functions are derived that describe the required computing cost for both imagers. On data from the Murchison Widefield Array, we find our new method to be an order of magnitude faster than w-projection, as well as being capable of full-sky imaging at full resolution and with correct polarization correction. We predict the computing costs for several other arrays and estimate that our imager is a factor of 2–12 faster, depending on the array configuration. We estimate the computing cost for imaging the low-frequency Square Kilometre Array observations to be 60 PetaFLOPS with current techniques. We find that combining w-stacking with the w-snapshot algorithm does not significantly improve computing requirements over pure w-stacking. The source code of our new imager is publicly released.
Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield ...Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of Δ2 ≤ 3.9 × 103 mK2 at k = 0.20 h Mpc−1 and z = 7 using 21 hr of data, improving previous MWA limits by almost an order of magnitude.
Abstract
We compute the spherically averaged power spectrum from four seasons of data obtained for the Epoch of Reionization (EoR) project observed with the Murchison Widefield Array (MWA). We ...measure the EoR power spectrum over k = 0.07–3.0 h Mpc−1 at redshifts $z$ = 6.5–8.7. The largest aggregation of 110 h on EoR0 high band (3340 observations), yields a lowest measurement of (43 mK)2 = 1.8 × 103 mK2 at k = 0.14 h Mpc−1 and $z$ = 6.5 (2σ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013–2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13 591 2-min snapshots (453 h), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12 569 observations (419 h). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8935 observations (298 h of data) to form integrated power spectra over the different fields, pointings, and redshift ranges.
We present a survey for transient and variable sources, on time-scales from 28 s to ∼1 yr, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no ...transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of <4.1 × 10−7 deg−2 for the shortest time-scale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest time-scale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of <82 FRBs per sky per day for dispersion measures <700 pc cm−3. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities ≳0.5 Jy and identify no sources with significant variability in their light curves. However, we note that substantial further work is required to fully characterize both the short-term and low-level variability within this field.
21 cm epoch of reionization (EoR) observations promise to transform our understanding of galaxy formation, but these observations are impossible without unprecedented levels of instrument ...calibration. We present end-to-end simulations of a full EoR power spectrum (PS) analysis including all of the major components of a real data processing pipeline: models of astrophysical foregrounds and EoR signal, frequency-dependent instrument effects, sky-based antenna calibration, and the full PS analysis. This study reveals that traditional sky-based per-frequency antenna calibration can only be implemented in EoR measurement analyses if the calibration model is unrealistically accurate. For reasonable levels of catalogue completeness, the calibration introduces contamination in otherwise foreground-free PS modes, precluding a PS measurement. We explore the origin of this contamination and potential mitigation techniques. We show that there is a strong joint constraint on the precision of the calibration catalogue and the inherent spectral smoothness of antennas, and that this has significant implications for the instrumental design of the SKA (Square Kilometre Array) and other future EoR observatories.
ABSTRACT
Observations in the lowest Murchison Widefield Array (MWA) band between 75 and 100 MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at ...redshift ∼13–17. Using 15 h of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction-independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction-dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 h of the cleanest data. Using these data, we obtain 2σ upper limits on the 21 cm power spectrum in the range of $0.1~ h~{\mathrm{ Mpc}}^{-1}\lessapprox k \lessapprox 1 ~ ~h~{\mathrm{ Mpc}}^{-1}$ and at z = 14.2, 15.2, and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $k=0.14 ~h~{\mathrm{ Mpc}}^{-1}$ and at z = 15.2 with a possibility of a few per cent of signal loss due to direction-independent calibration.
The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and ...redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hr of data targeting one of the MWA's epoch of reionization (EoR) fields observed in 2016. We use both improved analysis techniques presented in Barry et al. and several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8, and 7.1 based on our deep analysis on this 40 hr data set. These limits span a range in k-space of 0.18 h Mpc−1 < k < 1.6 h Mpc−1, with a lowest measurement of Δ2 ≤ 2.39 × 103 mK2 at k = 0.59 h Mpc−1 and z = 6.5.
GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results ...from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3 mm and 2 mm) and the 40m Yebes telescope (1.3 cm and 7 mm) to determine the fractional abundances of CO, HCO
, HCN, CS, SO, HCS
, and N
H
in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from A
~ 3 to ~20 mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of 1-5×10
cm
; and ii) the dense phase, located at A
> 10 mag, with molecular hydrogen densities >10
cm
. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C
/C/CO transition zone (A
~ 3 mag) where C/H ~ 8×10
and C/O~0.8-1, until the beginning of the dense phase at A
~ 10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H ~ (0.4 - 2.2) ×10
, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H ~8×10
). Based on our chemical modeling, we constrain the value of
to ~ (0.5 - 1.8) ×10
s
in the translucent cloud.