ABSTRACT
Cosmic dawn, the onset of star formation in the early universe, can in principle be studied via the 21-cm transition of neutral hydrogen, for which a sky-averaged absorption signal, ...redshifted to MHz frequencies, is predicted to be O(10–100) mK. Detection requires separation of the 21-cm signal from bright chromatic foreground emission due to Galactic structure, and the characterization of how it couples to instrumental response. In this work, we present characterization 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 21-cm signal identified with high accuracy. However, uncertainties on the soil properties lead to percentage 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 per cent 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, emphasize the importance of chromatic response and suggest caution in assuming log-polynomial foreground models in global signal experiments.
ABSTRACT
Total-power radiometry with individual meter-wave antennas is a potentially effective way to study the Cosmic Dawn (z ∼ 20) through measurement of the sky brightness arising from the 21 cm ...transition of neutral hydrogen, provided this can be disentangled from much stronger Galactic and extra-galactic foregrounds. In the process, measured spectra of integrated sky brightness temperature can be used to quantify the foreground emission properties. In this work, we analyse a subset of data from the Large-aperture Experiment to Detect the Dark Age (LEDA) in the 50–87 MHz range and constrain the foreground spectral index β in the northern sky visible from mid-latitudes. We focus on two zenith-directed LEDA radiometers and study how estimates of β vary with local sidereal time (LST). We correct for the effect of gain pattern chromaticity and compare estimated absolute temperatures with simulations. We select a reference data set consisting of 14 d of observations in optimal conditions. Using this data set, we find, for one radiometer, that β varies from −2.55 at LST <6 h to a steeper −2.58 at LST ∼13 h, consistently with sky models and previous southern sky measurements. In the 13 − 24 h LST range, however, we find that β varies between −2.55 and −2.61 (data scatter ∼0.01). We observe a similar β versus LST trend for the second radiometer, although with slightly smaller |β| over the 24 h, in the −2.46 < β < −2.43 range (data scatter ∼ 0.02). Combining all data gathered during the extended campaign between mid-2018 and mid-2019, and focusing on the LST = 9−12.5 h range, we infer good instrument stability and find −2.56 < β < −2.50 with 0.09 < Δβ < 0.12.
LOFAR sparse image reconstruction Garsden, H.; Girard, J. N.; Starck, J. L. ...
Astronomy and astrophysics (Berlin),
03/2015, Letnik:
575, Številka:
A90
Journal Article
Recenzirano
Odprti dostop
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.
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, disc - 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.
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 source 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 × 10−5 M⊙, 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 a 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 × 10−5 M⊙ 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, and 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 disc, 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 time-scales of about three months.
Water masers have been observed in several high-redshift active galactic nuclei, including the gravitationally lensed quasar MG 0414+0534. This quasar is lensed into four images, and the water maser ...is detected in two of them. The broadening of the maser emission line and its velocity offset are consistent with a group of masers associated with a quasar jet. If the maser group is microlensed we can probe its structure and size by observing its microlensing behaviour over time. We present results of a high-resolution numerical analysis of microlensing of the maser in MG 0414+0534, using several physically motivated maser models covering a range of sizes and emission profiles. Time-varying spectra of the microlensed maser are generated, displayed and analysed, and the behaviour of the different models compared. The observed maser line in MG 0414+0534 is consistent with maser spots as in other quasar jets, provided substructure is de-magnified or currently lost in noise; otherwise smooth extended maser models are also candidates to generate the observed spectrum. Using measures of spectral variability we find that if the maser has small substructure of ∼0.002 pc then a variation of 0.12 mag in flux and 2.0 km s−1 in velocity centroid of the maser line could be observed within two decades. For the smallest maser model in this study a magnification of >35 is possible 22 per cent of the time, which is of significance in the search for other lensed masers.
Recent observations have provided strong evidence for the presence of an electron-scattering region (ESR) within the central regions of active galactic nuclei. This is responsible for reprocessing ...emission from the accretion disc into polarized radiation. The geometry of this scattering region is, however, poorly constrained. In this paper, we consider the influence of gravitational microlensing on polarized 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 polarized views of the ESR, and is clearest for a small ESR width and a large ESR radius. Cross-correlation and autocorrelation 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 autocorrelation and cross-correlation of polarized light curves on much shorter time-scales.
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.