Investigating cosmological weak lensing with the ray-bundle method

Using the ray-bundle method for calculating gravitational lens magnifications, we outline a method by which the magnification probability may be determined specifically in the weak lensing limit for cosmological models obtained from N -body simulations.
16 different models are investigated, which are variations on three broad classes of cold dark matter model: the standard model with  (Ω₀, λ ₀)=(1.0,0.0)  , the open model with  (Ω₀, λ ₀)=(0.3,0.0)  and the lambda model, which is a flat model with a cosmological constant  (Ω₀, λ ₀)=(0.3,0.7)  .
The effects of varying the Hubble parameter, H ₀, the power spectrum shape parameter, Γ, and the cluster mass normalization, σ ₈, are studied. It is shown that there is no signature of these parameters in the weak lensing magnification distributions. The magnification probability distributions are also shown to be independent of the numerical parameters such as the lens mass and simulation box size in the N -body simulations.     

Planetary microlensing at high magnification

Simulations of planetary microlensing at high magnification that were carried out on a cluster computer are presented. It was found that the perturbations owing to two-thirds of all planets occur in the time interval  −0.5 t _FWHM,0.5 t _FWHM  with respect to the peak of the microlensing light curve, where   t _FWHM  is typically ∼14 h. This implies that only this restricted portion of the light curve need be intensively monitored for planets – a very significant practical advantage. Nearly all planetary detections in high-magnification events will not involve caustic crossings. We discuss the issues involved in determining the planetary parameters in high magnification events. Earth-mass planets may be detected with 1-m class telescopes if their projected orbital radii lie within about 1.5–2.5 au. Giant planets are detectable over a much larger region. For multiplanet systems the perturbations caused by individual planets can be separated under certain conditions. The size of the source star needs to be determined independently, but the presence of spots on the source star is likely to be negligible, as is the effect of planetary motion during an event.     

A maximum-entropy method for reconstructing the projected mass distribution of gravitational lenses

The maximum-entropy method is applied to the problem of reconstructing the projected mass density of a galaxy cluster using its gravitational lensing effects on background galaxies. We demonstrate the method by reconstructing the mass distribution in a model cluster using simulated shear and magnification data to which Gaussian noise is added. The mass distribution is reconstructed directly and the inversion is regularized using the entropic prior for this positive additive distribution. For realistic noise levels, we find that the method faithfully reproduces the main features of the cluster mass distribution not only within the observed field but also slightly beyond it. We estimate the uncertainties in the reconstruction by calculating an analytic approximation to the covariance matrix of the reconstruction values of each pixel. This result is compared with error estimates derived from Monte Carlo simulations for different noise realizations and found to be in good agreement.     

IRAM observations of JVAS/CLASS gravitational lenses

We have searched for molecular absorption lines at millimetre wavelengths in 11 gravitational lens systems discovered in the JVAS/CLASS surveys of flat spectrum radio sources. Spectra of only one source 1030+074 were obtained in the 3-, 2- and 1.3-mm bands at the frequencies corresponding to common molecular transitions of CO and HCO⁺ as continuum emission was not found in any of the other sources. We calculated upper limits to the column density in molecular absorption for 1030+074, using an excitation temperature of 15 K, to be N _CO<6.3×10¹³ cm⁻² and N _HCO+<1.3×10¹¹ cm⁻² , equivalent to hydrogen column density of the order N _H<10¹⁸ cm⁻² , assuming standard molecular abundances. We also present the best upper limits of the continuum at the lower frequency for the other 10 gravitational lenses.     

Observations of time delays in gravitational lenses from intensity fluctuations: the coherence function

Measurements of the spectrum of the fluctuations of the output current of the quadratic detector of a telescope can be used to find unresolved astronomical gravitational lenses and determine time delays between their image components. These time delays can be used for astronomical studies. The spatial correlation coefficient of a source is an important parameter that quantifies the loss of contrast, caused by the extendedness of the source, in the spectral modulation of the intensity fluctuations. This work shows that the correlation coefficient must not be evaluated at the frequency of observation, but must instead be evaluated at the much lower beat frequencies of the spectrum of the fluctuations. This opens up a powerful, novel technique to find unresolved gravitational lenses and to study the lensing event and the source.     

The Planck Surveyor mission and gravitational lenses

An extremely sensitive all-sky survey will be carried out in the millimetre/submillimetre waveband by the forthcoming ESA mission Planck Surveyor . The main scientific goal of the mission is to make very accurate measurements of the spatial power spectrum of primordial anisotropies in the cosmic microwave background radiation; however, hundreds of thousands of distant dusty galaxies and quasars will also be detected. These sources are much more likely to be gravitationally lensed by intervening galaxies compared with sources discovered in surveys in other wavebands. Here the number of lenses expected in the survey is estimated, and techniques for discriminating between lensed and unlensed sources are discussed. A practical strategy for this discrimination is presented, based on exploiting the remarkable sensitivity and resolving power of large ground-based millimetre/submillimetre-wave interferometer arrays. More than a thousand gravitational lenses could be detected: a sample that would be an extremely valuable resource in observational cosmology.     

Improving efficiency in radio surveys for gravitational lenses

Many lens surveys have hitherto used observations of large samples of background sources to select the small minority which are multiply imaged by lensing galaxies along the line of sight. Recently surveys such as SLACS and OLS have improved the efficiency of surveys by pre-selecting double-redshift systems from SDSS. We explore other ways to improve survey efficiency by optimum use of astrometric and morphological information in existing large-scale optical and radio surveys. The method exploits the small position differences between FIRST radio positions of lensed images and the SDSS lens galaxy positions, together with the marginal resolution of some larger gravitational lens systems by the FIRST beam. We present results of a small pilot study with the VLA and MERLIN, and discuss the desirable criteria for future surveys.     

Two new large-separation gravitational lenses from SDSS

We present discovery images, together with follow-up imaging and spectroscopy, of two large-separation gravitational lenses found by our survey for wide arcs [the CAmbridge Sloan Survey Of Wide ARcs in the skY (CASSOWARY)]. The survey exploits the multicolour photometry of the Sloan Digital Sky Survey to find multiple blue components around red galaxies. CASSOWARY 2 (or 'the Cheshire Cat') is composed of two massive early-type galaxies at   z = 0.426  and 0.432, respectively, lensing two background sources, the first a star-forming galaxy at   z = 0.97  and the second a high -redshift galaxy  ( z > 1.4)  . There are at least three images of the former source and probably four or more of the latter, arranged in two giant arcs. The mass enclosed within the larger arc of radius ∼11 arcsec is  ∼33 × 10¹² M_⊙  . CASSOWARY 3 comprises an arc of three bright images of a   z = 0.725  source, lensed by a foreground elliptical at   z = 0.274  . The radius of the arc is ∼4 arcsec and the enclosed mass is  ∼2.5 × 10¹² M_⊙  . Together with earlier discoveries like the Cosmic Horseshoe and the 8 o'clock Arc, these new systems, with separations intermediate between the arcsecond-separation lenses of typical strong galaxy lensing and arcminute-separation cluster lenses, probe the very high end of the galaxy mass function.     

Rotation in gravitational lenses

Gravitational lensing deflects light. A single lens deflector can only shear images, but cannot induce rotations. Multiple lens planes can induce rotations. Such rotations can be observed in quadruply imaged sources, and can be used to distinguish between two proposed solutions of the flux anomaly problem: substructures in lensing galaxies versus large-scale structure. We predict the expected amount of rotation due to large-scale structure in strong lensing systems, and show how this effect can be measured using ∼mas very long baseline interferometry astrometry of quadruple lenses with extended source structures. The magnitude of rotation is around 1°. The biggest theoretical uncertainty is the power spectrum of dark matter on very small scales. This procedure can potentially be turned around to measure the dark matter power spectrum on very small scales. We list the predicted rms rotation angles for several quadruple lenses with known lens and source redshifts.     

Scalar field haloes as gravitational lenses

A non-topological soliton model with a repulsive scalar self-interaction of the Emden type provides a constant-density core, similarly as the empirical Burkert profile of dark matter (DM) haloes. As a further test, we derive the gravitational lens properties of our model, in particular, the demarcation curves between 'weak' and 'strong' lensing. Accordingly, strong lensing with typically three images is almost three times more probable for our solitonic model than for the Burkert fit. Moreover, some prospective consequences of a possible flattening of DM haloes are indicated.     

Study by MOA of extrasolar planets in gravitational microlensing events of high magnification

A search for extrasolar planets was carried out in three gravitational microlensing events of high magnification, MACHO  98–BLG–35  , MACHO  99–LMC–2  and OGLE  00–BUL–12  . Photometry was derived from observational images by the MOA and OGLE groups using an image subtraction technique. For MACHO  98–BLG–35  , additional photometry derived from the MPS and PLANET groups was included. Planetary modelling of the three events was carried out in a supercluster computing environment. The estimated probability for explaining the data on MACHO  98–BLG–35  without a planet is <1 per cent. The best planetary model has a planet of mass ∼(0.4–1.5)× M _Earth at a projected radius of either ∼1.5 or ∼2.3 au. We show how multiplanet models can be applied to the data. We calculate exclusion regions for the three events and find that Jupiter-mass planets can be excluded with projected radii from as wide as about 30 au to as close as around 0.5 au for MACHO  98–BLG–35  and OGLE  00–BUL–12  . For MACHO  99–LMC–2  , the exclusion region extends out to around 10 au and constitutes the first limit placed on a planetary companion to an extragalactic star. We derive a particularly high peak magnification of ∼160 for OGLE  00–BUL–12  . We discuss the detectability of planets with masses as low as Mercury in this and similar events.     

Strong gravitational lensing by spiral galaxies

We investigate gravitational lensing using a realistic model of disc galaxies. Most of the mass is contained in a large spherical isothermal dark matter halo, but the potential is modified significantly in the core by a gravitationally dominant exponential disc. The method used is adapted from a very general multilens ray-tracing technique developed by Mo¨ller. We investigate the effects of the disc-to-halo mass ratio, the disc scalelength, the disc inclination to the line of sight and the lens redshift on two strong-lensing cross-sections: the cross-section for multiple imaging and the cross-section for large magnifications, in excess of a factor of 10. We find that the multiple-imaging cross-section can be enhanced significantly by an almost edge-on Milky Way disc compared with a singular isothermal sphere (SIS) in individual cases; however, when averaged over all disc inclinations, the cross-section is only increased by about 50 per cent. These results are consistent with other recent work. The presence of a disc, however, increases the inclination-averaged high-magnification cross-section by an order of magnitude compared with a SIS. This result has important implications for magnification bias in future lens surveys, particularly those in the submillimetre waveband, where dust extinction in the lensing galaxy has no effect on the brightness of the images.     

Measuring the three-dimensional shear from simulation data, with applications to weak gravitational lensing

We have developed a new three-dimensional algorithm, based on the standard P³M method, for computing deflections resulting from weak gravitational lensing. We compare the results of this method with those of the two-dimensional planar approach, and rigorously outline the conditions under which the two approaches are equivalent. Our new algorithm uses a Fast Fourier Transform convolution method for speed, and has a variable softening feature to provide a realistic interpretation of the large-scale structure in a simulation. The output values of the code are compared with those from the Ewald summation method, which we describe and develop in detail. With an optimal choice of the high-frequency filtering in the Fourier convolution, the maximum errors, when using only a single particle, are about 7 per cent, with an rms error less than 2 per cent. For ensembles of particles, used in typical N -body simulations, the rms errors are typically 0.3 per cent. We describe how the output from the algorithm can be used to generate distributions of magnification, source ellipticity, shear and convergence for large-scale structure.