High-resolution imaging of the cosmic mass distribution from gravitational lensing of pre-galactic H i

Low-frequency radio observations of neutral hydrogen during and before the epoch of cosmic re-ionization will provide ∼1000 quasi-independent source planes, each of precisely known redshift, if a resolution of ∼1 arcmin or better can be attained. These planes can be used to reconstruct the projected mass distribution of foreground material. Structure in these source planes is linear and Gaussian at high redshift  (30 < z < 300)  but is non-linear and non-Gaussian during re-ionization. At both epochs, significant power is expected down to subarcsecond scales. We demonstrate that this structure can, in principle, be used to make mass images with a formal signal-to-noise ratio (S/N) per pixel exceeding 10, even for pixels as small as an arcsecond. With an ideal telescope, both resolution and S/N can exceed those of even the most optimistic idealized mass maps from galaxy lensing by more than an order of magnitude. Individual dark haloes similar in mass to that of the Milky Way could be imaged with high S/N out to   z ∼ 10  . Even with a much less ambitious telescope, a wide-area survey of 21-cm lensing would provide very sensitive constraints on cosmological parameters, in particular on dark energy. These are up to 20 times tighter than the constraints obtainable from comparably sized, very deep surveys of galaxy lensing, although the best constraints come from combining data of the two types. Any radio telescope capable of mapping the 21-cm brightness temperature with good frequency resolution (∼0.05 MHz) over a band of width ≳10 MHz should be able to make mass maps of high quality. The planned Square Kilometre Array may be able to map the mass with moderate S/N down to arcminute scales, depending on the re-ionization history of the universe and the ability to subtract foreground sources.     

Precision of diffuse 21-cm lensing

We study the limits of accuracy for weak lensing maps of dark matter using diffuse 21-cm radiation from the pre-reionization epoch using simulations. We improve on previous 'optimal' quadratic lensing estimators by using shear and convergence instead of deflection angles. This is a generalization of the deflection estimator, and is more optimal for non-Gaussian sources. The cross-power spectrum of shear and convergence is an unbiased estimator of lensing power spectrum which does not require knowledge of the source four-point function. We find that non-Gaussianity provides a limit to the accuracy of weak lensing reconstruction, even if instrumental noise is reduced to zero. The best reconstruction result is equivalent to Gaussian sources with effective independent cell of side length  2.0  h ⁻¹ Mpc  . Using a source full map from z = 10 to 20, this limiting sensitivity allows mapping of dark matter at a signal-to-noise ratio greater than 1 out to l ≲ 6000, which is better than any other proposed technique for large-area weak lensing mapping.     

Noise properties of gravitational lens mass reconstruction

Gravitational lensing is potentially able to observe mass-selected haloes, and to measure the projected cluster mass function. An optimal mass selection requires a quantitative understanding of the noise behaviour in mass maps. This paper is an analysis of the noise properties in mass maps reconstructed from a maximum-likelihood method.
The first part of this work is the derivation of the noise power spectrum and the mass error bars as a straightforward extension of the Kaiser & Squires algorithm for the case of a correlated noise. Very good agreement is found between these calculations and the noise properties measured in the mass reconstructions limited to non-critical clusters of galaxies. It demonstrates that Kaiser & Squires and maximum-likelihood methods have similar noise properties and that the weak lensing approximation is valid for describing these properties .
In a second stage I show that the statistics of peaks in the noise follows accurately the peak statistics of a two-dimensional Gaussian random field (using the BBKS techniques) if the smoothing aperture contains enough galaxies. This analysis provides a full procedure for deriving the significance of any convergence peak as a function of its amplitude and profile.
I demonstrate that a detailed quantitative analysis of the structures in mass maps can be carried out, and that, to a very good approximation, a mass map is the sum of the lensing signal and known two-dimensional Gaussian random noise. A straightforward application is the measurement of the projected mass function in wide-field lensing surveys, down to small mass overdensities that are individually undetectable.     

Weak lensing from strong clustering

We investigate the effect of weak gravitational lensing in the limit of small angular scales where projected galaxy clustering is strongly non-linear. This is the regime likely to be probed by future weak lensing surveys. We use well-motivated hierarchical scaling arguments and the plane-parallel approximation to study multi-point statistical properties of the convergence field. These statistics can be used to compute the vertex amplitudes in tree models of hierarchical clustering; these can be compared with similar measurements from galaxy surveys, leading to a powerful probe of galaxy bias.     

Galaxy density profiles and shapes – II. Selection biases in strong lensing surveys

Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio and image separation distribution, with a flexible treatment of magnification bias to mimic different survey strategies. We present our results for two families of density profiles: cusped and deprojected Sérsic models. While we use fixed lens and source redshifts for most of the analysis, we show that the results are applicable to other redshift combinations, and we also explore the physics of how our results change for very different redshifts. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sérsic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.     

Weak gravitational flexion

Flexion is the significant third-order weak gravitational lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. Here we demonstrate how flexion measurements can be used to measure galaxy halo density profiles and large-scale structure on non-linear scales, via galaxy–galaxy lensing, dark matter mapping and cosmic flexion correlation functions. We describe the origin of gravitational flexion, and discuss its four components, two of which are first described here. We also introduce an efficient complex formalism for all orders of lensing distortion. We proceed to examine the flexion predictions for galaxy–galaxy lensing, examining isothermal sphere and Navarro–Frenk–White (NFW) profiles and both circularly symmetric and elliptical cases. We show that in combination with shear we can precisely measure galaxy masses and NFW halo concentrations. We also show how flexion measurements can be used to reconstruct mass maps in two-dimensional projection on the sky, and in three dimensions in combination with redshift data. Finally, we examine the predictions for cosmic flexion, including convergence–flexion cross-correlations, and we find that the signal is an effective probe of structure on non-linear scales.     

Halo masses for optically selected and for radio-loud AGN from clustering and galaxy–galaxy lensing

We compute two-point correlation functions and measure the shear signal due to galaxy–galaxy lensing for 80 000 optically identified and 5700 radio-loud active galactic nuclei (AGN) from Data Release 4 of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed   M _*  . Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive haloes than optical AGN: in our samples, their mean halo masses are  1.6 × 10¹³  and  8 × 10¹¹  h ⁻¹ M_⊙  , respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter haloes of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in haloes where gas cools quasi-statically.     

B0712+472: a new radio four-image gravitational lens

A new four-image gravitational lens system, B0712+472, has been discovered during the Cosmic Lens All-Sky Survey. This system consists of four flat-spectrum radio images that are also seen on a Hubble Space Telescope ( HST ) image, together with the lensing galaxy. We present MERLIN, VLA and VLBA maps and WHT spectra of the system as well as the HST images. The light distribution of the lensing galaxy is highly elongated and so too is the mass distribution deduced from modelling. We suggest a redshift of ∼1.33 for the lensed object; the lens redshift will require further investigation. The discovery of this new system further increases the ratio of four-image to two-image lens systems currently known, exacerbating problems of required ellipticity of matter distributions in lensing galaxies.     

Weak lensing of the CMB by large-scale structure

Several recent papers have studied lensing of the CMB by large-scale structures, which probes the projected matter distribution from z=10³ to z0. This interest is motivated in part by upcoming high resolution, high sensitivity CMB experiments, such as APEX/SZ, ACT, SPT or Planck, which should be sensitive to lensing. In this paper, we examine the reconstruction of the large-scale dark matter distribution from lensed CMB temperature anisotropies. We go beyond previous work in using numerical simulations to include higher order, non-Gaussian effects and find that the convergence and its power spectrum are biased, with the bias increasing with the angular resolution. We also study the contamination by the kinetic Sunyaev–Zel'dovich signal, which is spectrally indistinguishable from lensed CMB anisotropies, and find that it leads to an overestimate of the convergence. We finish by estimating the sensitivity of the previously cited experiments and find that all of them could detect the lensing effect, but would be biased at around the 10% level.     

S-Component Sources at 21 cm Wavelength in the Rising Phase of Cycle 23

The Synthesis Radio Telescope at the Dominion Radio Astrophysical Observatory is being used to make a long-term study of the relationship between sources of the slowly-varying component at 21-cm wavelength and the active regions with which they are associated. Although the arc-minute angular resolution of the radio telescope is insufficient to map individual active regions, the unique capability of this radio telescope to map the whole solar disc in a single operation makes it useful for identifying and measuring sources of the slowly-varying component, which at this wavelength are largely due to free–free thermal emission from plasma concentrations trapped in active region magnetic fields. The brightness temperatures of the emission are higher in large, complex active regions, suggesting the plasma slab producing the emission is thicker and/or denser in such regions. We find also that the relationship between the brightness temperature and the size and complexity of the host region is a function of the phase of the 10–13-year solar activity cycle.     

Neutral hydrogen gas distribution and kinematics in NGC 7479

We have investigated the effects of a bar and an asymmetric spiral structure on the neutral hydrogen distribution and kinematics in the strongly barred spiral galaxy NGC 7479. The strongest 21-cm line emission at 1-kpc resolution comes from the western spiral arm which appears to be slightly inclined to the plane of the main disc. In contrast, the area within the radius of the bar is devoid of 21-cm line emission. The radio continuum emission at 21 cm follows the bar dust lanes, but beyond 3 kpc from the nucleus the radio continuum emission has a peculiar morphology, unlike that of optical and near-infrared images. We did not detect any low surface brightness gas-rich companions in the near neighbourhood of NGC 7479. This leads us to propose that the strong western spiral arm was created in a recent minor merger.     

QSO lensing magnification associated with galaxy groups

We simulated both the matter and light (galaxy) distributions in a wedge of the Universe and calculated the gravitational lensing magnification caused by the mass along the line-of-sight of galaxies and galaxy groups identified in sky surveys. A large volume redshift cone containing cold dark matter particles mimics the expected cosmological matter distribution in a flat universe with low matter density and a cosmological constant. We generate a mock galaxy catalogue from the matter distribution and identify thousands of galaxy groups in the luminous sky projection. We calculate the expected magnification around galaxies and galaxy groups and then the induced quasi-stellar object (QSO)–lens angular correlation due to magnification bias. This correlation is observable and can be used both to estimate the average mass of the lens population and to make cosmological inferences. We also use analytical calculations and various analyses to compare the observational results with theoretical expectations for the cross-correlation between faint QSOs from the 2dF Survey and nearby galaxies and groups from the Automated Plate Measurement and Sloan Digital Sky Survey Early Data Release. The observed QSO–lens anticorrelations are stronger than the predictions for the cosmological model used. This suggests that there could be unknown systematic errors in the observations and data reduction, or that the model used is not adequate. If the observed signal is assumed to be solely due to gravitational lensing, then the lensing is stronger than expected, due to more massive galactic structures or more efficient lensing than simulated.     

The dynamical state of RX J1347.5−1145 from a combined strong lensing and X-ray analysis

We perform a combined X-ray and strong lensing analysis of RX J1347.5−1145, one of the most luminous galaxy clusters at X-ray wavelengths. We show that evidence from strong lensing alone, based on published Very Large Telescope (VLT) and new Hubble Space Telescope ( HST ) data, strongly argues in favour of a complex structure. The analysis takes into account arc positions, shapes and orientations, and is done thoroughly in the image plane. The cluster inner regions are well fitted by a bimodal mass distribution, with a total projected mass of   M _tot= (9.9 ± 0.3) × 10¹⁴ M_⊙  h ⁻¹  within a radius of 360 kpc  h ⁻¹ (1.5 arcmin). Such a complex structure could be a signature of a recent major merger as further supported by X-ray data. A temperature map of the cluster, based on deep Chandra observations, reveals a hot front located between the first main component and an X-ray emitting south-eastern subclump. The map also unveils a filament of cold gas in the innermost regions of the cluster, most probably a cooling wake caused by the motion of the cD inside the cool core region. A merger scenario in the plane of the sky between two dark matter subclumps is consistent with both our lensing and X-ray analyses, and can explain previous discrepancies with mass estimates based on the virial theorem.     

The existence and detection of optically dark galaxies by 21-cm surveys

One explanation for the disparity between cold dark matter (CDM) predictions of galaxy numbers and observations could be that there are numerous dark galaxies in the Universe. These galaxies may still contain baryons, but no stars, and may be detectable in the 21-cm line of atomic hydrogen. The results of surveys for such objects, and simulations that do/do not predict their existence, are controversial. In this paper, we use an analytical model of galaxy formation, consistent with CDM, to first show that dark galaxies are certainly a prediction of the model. Secondly, we show that objects like VIRGOHI21, a dark galaxy candidate recently discovered by us, while rare are predicted by the model. Thirdly, we show that previous 'blind' H  i surveys have placed few constraints on the existence of dark galaxies. This is because they have either lacked the sensitivity and/or velocity resolution or have not had the required detailed optical follow up. We look forward to new 21-cm blind surveys [Arecibo Legacy Fast ALFA (ALFALFA) survey and Arecibo Galactic Environments Survey (AGES)] using the Arecibo multibeam instrument which should find large numbers of dark galaxies if they exist.     

Galaxy density profiles and shapes – I. Simulation pipeline for lensing by realistic galaxy models

Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, which we aim to investigate in a consistent way. In this first paper in a series, we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii.     

High resolution maps of two southern radio galaxies 1331-09 and 1417-20

Two southern radio galaxies 1331-09 and 1417-20 have been observed with the Fleurs synthesis radio telescope at 1415 MHz. The deduced geometrical and physical parameters show that 1331-09 is a giant radio galaxy having two outer components extending 1.9 Mpc and a central component coinciding with its optical counterpart. Many similarities of its features to CygA suggest that hot spots structure at the outer edge of its components should be expected in a higher resolution map. The source 1417-20 is a rather typical radio galaxy with halo-nucleus structure.     

Tracking the shadows through GMRT

The structures of faint high redshift galaxies cannot be observed directly. But if a luminous quasar is located farther along their line of sight, high resolution absorption lines offer a valuable and reliable probe to their structure. GMRT is suited to monitor the absorption spectra, if the redshifted neutral hydrogen or OH doublet fall in one of the windows of the telescope. We present the OH doublet absorption spectra for the system B0218+357, taken at GMRT this year at resolution of approx. 9.5 km/sec with an rms noise of the order of 1 mJy. Based on our study of the OH doublet and 21cm neutral hydrogen line we infer that, in the lensing spiral galaxy of B0218 + 357, neutral hydrogen and OH coexist in tenous clouds and there is possibly a hole in the central part of the galaxy. In contrast, the gas is seen in high density clouds in the lens in an otherwise similar system PKS1830-211.     

Galaxy–dark matter correlations applied to galaxy–galaxy lensing: predictions from the semi-analytic galaxy formation models

We use semi-analytic models of galaxy formation combined with high-resolution N -body simulations to make predictions for galaxy–dark matter correlations and apply them to galaxy–galaxy lensing. We analyse cross-power spectra between the dark matter and different galaxy samples selected by luminosity, colour or star formation rate. We compare the predictions with the recent detection by the Sloan Digital Sky Survey (SDSS). We show that the correlation amplitude and the mean tangential shear depend strongly on the luminosity of the sample on scales below 1  h ⁻¹ Mpc, reflecting the correlation between the galaxy luminosity and the halo mass. The cross-correlation cannot, however, be used to infer the halo profile directly because different halo masses dominate on different scales and because not all galaxies are at the centres of the corresponding haloes. We compute the redshift evolution of the cross-correlation amplitude and compare it with those of galaxies and dark matter. We also compute the galaxy–dark matter correlation coefficient and show that it is close to unity on scales above 1  h ⁻¹ Mpc for all considered galaxy types. This would allow one to extract the bias and the dark matter power spectrum on large scales from the galaxy and galaxy–dark matter correlations.     

A complete infrared Einstein ring in the gravitational lens system B1938 + 666

We report the discovery, using NICMOS on the Hubble Space Telescope , of an arcsec-diameter Einstein ring in the gravitational lens system B1938+666. The lensing galaxy is also detected, and is most likely an early-type galaxy. Modelling of the ring is presented and compared with the radio structure from MERLIN maps. We show that the Einstein ring is consistent with the gravitational lensing of an extended infrared component, centred between the two radio components.