Search

The 21 cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation (EoR). This has led to the construction of low-frequency radio interferometric arrays, such as the Hydrogen Epoch of Reionization Array (HERA), aimed at systematically mapping this emission for the first time. Precision calibration, however, is a requirement in 21 cm radio observations. Due to the spatial compactness of HERA, the array is prone to the effects of mutual coupling, which inevitably lead to non-smooth calibration errors that contaminate the data. When unsmooth gains are used in calibration, intrinsically spectrally-smooth foreground emission begins to contaminate the data in a way that can prohibit a clean detection of the cosmological EoR signal. In this paper, we show that the effects of mutual coupling on calibration quality can be reduced by applying custom time-domain filters to the data prior to calibration. We find that more robust calibration solutions are derived when filtering in this way, which reduces the observed foreground power leakage. Specifically, we find a reduction of foreground power leakage by 2 orders of magnitude at k=0.5.

Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategies for mitigating mutual coupling. In this paper, we analyse 12 nights of data from the Hydrogen Epoch of Reionization Array and compare the data against simulations that include a computationally efficient and physically motivated semi-analytic treatment of mutual coupling. We find that simulated coupling features qualitatively agree with coupling features in the data; however, coupling features in the data are brighter than the simulated features, indicating the presence of additional coupling mechanisms not captured by our model. We explore the use of fringe-rate filters as mutual coupling mitigation tools and use our simulations to investigate the effects of mutual coupling on a simulated cosmological 21-cm power spectrum in a "worst case" scenario where the foregrounds are particularly bright. We find that mutual coupling contaminates a large portion of the "EoR Window", and the contamination is several orders-of-magnitude larger than our simulated cosmic signal across a wide range of cosmological Fourier modes. While our fiducial fringe-rate filtering strategy reduces mutual coupling by roughly a factor of 100 in power, a non-negligible amount of coupling cannot be excised with fringe-rate filters, so more sophisticated mitigation strategies are required., Comment: 19 pages, 12 figures, submitted to MNRAS

Cosmic Dawn, the onset of star formation in the early universe, can in principle be studied via the 21cm transition of neutral hydrogen, for which a sky-averaged absorption signal, redshifted to MHz frequencies, is predicted to be {\it O}(10-100)\,mK. Detection requires separation of the 21cm signal from bright chromatic foreground emission due to Galactic structure, and the characterisation of how it couples to instrumental response. In this work, we present characterisation of antenna gain patterns for the Large-aperture Experiment to detect the Dark Ages (LEDA) via simulations, assessing the effects of the antenna ground-plane geometries used, and measured soil properties. We then investigate the impact of beam pattern uncertainties on the reconstruction of a Gaussian absorption feature. Assuming the pattern is known and correcting for the chromaticity of the instrument, the foregrounds can be modelled with a log-polynomial, and the 21cm signal identified with high accuracy. However, uncertainties on the soil properties lead to \textperthousand\ changes in the chromaticity that can bias the signal recovery. The bias can be up to a factor of two in amplitude and up to few \% in the frequency location. These effects do not appear to be mitigated by larger ground planes, conversely gain patterns with larger ground planes exhibit more complex frequency structure, significantly compromising the parameter reconstruction. Our results, consistent with findings from other antenna design studies, emphasise the importance of chromatic response and suggest caution in assuming log-polynomial foreground models in global signal experiments., Comment: Accepted for publication in MNRAS

PySE is a Python software package for finding and measuring sources in radio telescope images. The software was designed to detect sources in the LOFAR telescope images, but can be used with images from other radio telescopes as well. We introduce the LOFAR Telescope, the context within which PySE was developed, the design of PySE, and describe how it is used. Detailed experiments on the validation and testing of PySE are then presented, along with results of performance testing. We discuss some of the current issues with the algorithms implemented in PySE and their inter- action with LOFAR images, concluding with the current status of PySE and its future development., Comment: 26 pages, 6 figures

The Large-Aperture Experiment to Detect the Dark Age (LEDA) was designed to detect the predicted O(100)mK sky-averaged absorption of the Cosmic Microwave Background by Hydrogen in the neutral pre- and intergalactic medium just after the cosmological Dark Age. The spectral signature would be associated with emergence of a diffuse Ly$\alpha$ background from starlight during 'Cosmic Dawn'. Recently, Bowman et al. (2018) have reported detection of this predicted absorption feature, with an unexpectedly large amplitude of 530 mK, centered at 78 MHz. Verification of this result by an independent experiment, such as LEDA, is pressing. In this paper, we detail design and characterization of the LEDA radiometer systems, and a first-generation pipeline that instantiates a signal path model. Sited at the Owens Valley Radio Observatory Long Wavelength Array, LEDA systems include the station correlator, five well-separated redundant dual polarization radiometers and backend electronics. The radiometers deliver a 30-85MHz band (16

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool., Comment: Accepted for publication in PASA. 27 pages, 10 figures, 5 tables

We present the highest-quality polarisation profiles to date of 16 non-recycled pulsars and four millisecond pulsars, observed below 200 MHz with the LOFAR high-band antennas. Based on the observed profiles, we perform an initial investigation of expected observational effects resulting from the propagation of polarised emission in the pulsar magnetosphere and the interstellar medium. The predictions of magnetospheric birefringence in pulsars have been tested using spectra of the pulse width and fractional polarisation from multifrequency data. The derived spectra offer only partial support for the expected effects of birefringence on the polarisation properties, with only about half of our sample being consistent with the model's predictions. It is noted that for some pulsars these measurements are contaminated by the effects of interstellar scattering. For a number of pulsars in our sample, we have observed significant variations in the amount of Faraday rotation as a function of pulse phase, which is possibly an artefact of scattering. These variations are typically two orders of magnitude smaller than that observed at 1400 MHz by Noutsos et al. (2009), for a different sample of southern pulsars. In this paper we present a possible explanation for the difference in magnitude of this effect between the two frequencies, based on scattering. Finally, we have estimated the magnetospheric emission heights of low-frequency radiation from four pulsars, based on the phase lags between the flux-density and the PA profiles, and the theoretical framework of Blaskiewicz, Cordes & Wasserman (1991). These estimates yielded heights of a few hundred km; at least for PSR B1133+16, this is consistent with emission heights derived based on radius-to-frequency mapping, but is up to a few times larger than the recent upper limit based on pulsar timing., Comment: 27 pages, 11 figures, 2 tables, A&A accepted

Context. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SKA, Comment: Published in A&A, 19 pages, 9 figures

Gravitational microlensing of planetary-mass objects (or "nanolensing", as it has been termed) can be used to probe the distribution of mass in a galaxy that is acting as a gravitational lens. Microlensing and nanolensing light curve fluctuations are indicative of the mass of the compact objects within the lens, but the size of the source is important, as large sources will smooth out a light curve. Numerical studies have been made in the past that investigate a range of sources sizes and masses in the lens. We extend that work in two ways - by generating high quality maps with over a billion small objects down to a mass of 2.5\times10-5M\odot, and by investigating the temporal properties and observability of the nanolensing events. The system studied is a mock quasar system similar to MG 0414+0534. We find that if variability of 0.1 mag in amplitude can be observed, a source size of ~ 0.1 Einstein Radius (ER) would be needed to see the effect of 2.5\times10-5M\odot masses, and larger, in the microlensing light curve. Our investigation into the temporal properties of nanolensing events finds that there are two scales of nanolensing that can be observed - one due to the crossing of nanolensing caustic bands, the other due to the crossing of nanolensing caustics themselves. The latter are very small, having crossing times of a few days, and requiring sources of size ~ 0.0001 ER to resolve. For sources of the size of an accretion disk, the nanolensing caustics are slightly smoothed-out, but can be observed on time scales of a few days. The crossing of caustic bands can be observed on times scales of about 3 months., Comment: 14 pages, 9 figures, 1 table, low-quality images

The kinematics and morphology of the broad emission line region (BELR) of quasars are the subject of significant debate. The two leading methods for constraining BELR properties are microlensing and reverberation mapping. Here we combine these two methods with a study of the microlensing behaviour of the BELR in Q2237+0305, as a change in continuum emission (a "flare") passes through it. Beginning with some generic models of the BELR - sphere, bicones, disk - we slice in velocity and time to produce brightness profiles of the BELR over the duration of the flare. These are numerically microlensed to determine whether microlensing of reverberation mapping provides new information about the properties of BELRs. We describe our method and show images of the models as they are flaring, and the unlensed and lensed spectra that are produced. Qualitative results and a discussion of the spectra are given in this paper, highlighting some effects that could be observed. Our conclusion is that the influence of microlensing, while not strong, can produce significant observable effects that will help in differentiating the properties of BELRs., Comment: 17 pages, 14 low resolution figures, 1 table, accepted for MNRAS. v2: Corrected velocities p16, 8 to 0.08, 9 to 0.09

Recent observations have provided strong evidence for the presence of an Electron Scattering Region (ESR) within the central regions of AGNs. This is responsible for reprocessing emission from the accretion disk into polarised radiation. The geometry of this scattering region is, however, poorly constrained. In this paper, we consider the influence of gravitational microlensing on polarised emission from the ESR in the quadruply imaged quasar, Q2237+0305, demonstrating how correlated features in the resultant light curve variations can determine both the size and orientation of the scattering region. This signal is due to differential magnification between perpendicularly polarised views of the ESR, and is clearest for a small ESR width and a large ESR radius. Cross- and auto-correlation measures appear to be independent of lens image shear and convergence parameters, making it ideal to investigate ESR features. As with many microlensing experiments, the time-scale for variability, being of order decades to centuries, is impractically long. However, with a polarization filter oriented appropriately with respect to the path that the quasar takes across the caustic structure, the ESR diameter and radius can be estimated from the auto- and cross-correlation of polarized light curves on much shorter time-scales., Comment: 11 pages, 12 figures, 1 table, accepted for MNRAS

To assess how future progress in gravitational microlensing computation at high optical depth will rely on both hardware and software solutions, we compare a direct inverse ray-shooting code implemented on a graphics processing unit (GPU) with both a widely-used hierarchical tree code on a single-core CPU, and a recent implementation of a parallel tree code suitable for a CPU-based cluster supercomputer. We examine the accuracy of the tree codes through comparison with a direct code over a much wider range of parameter space than has been feasible before. We demonstrate that all three codes present comparable accuracy, and choice of approach depends on considerations relating to the scale and nature of the microlensing problem under investigation. On current hardware, there is little difference in the processing speed of the single-core CPU tree code and the GPU direct code, however the recent plateau in single-core CPU speeds means the existing tree code is no longer able to take advantage of Moore's law-like increases in processing speed. Instead, we anticipate a rapid increase in GPU capabilities in the next few years, which is advantageous to the direct code. We suggest that progress in other areas of astrophysical computation may benefit from a transition to GPUs through the use of "brute force" algorithms, rather than attempting to port the current best solution directly to a GPU language -- for certain classes of problems, the simple implementation on GPUs may already be no worse than an optimised single-core CPU version., Comment: 11 pages, 4 figures, accepted for publication in New Astronomy