The Chang–Refsdal lens revisited

This paper provides a complete theoretical treatment of the point-mass lens perturbed by constant external shear, often called the Chang–Refsdal lens. We show that simple invariants exist for the products of the (complex) positions of the four images, as well as moment sums of their signed magnifications. The image topographies and equations of the caustics and critical curves are also studied. We derive the fully analytic expressions for pre-caustics, which are the loci of non-critical points that map to the caustics under the lens mapping. They constitute boundaries of the region in the image domain that maps on to the interior of the caustics. The areas under the critical curves, caustics and pre-caustics are all evaluated, which enables us to calculate the mean magnification of the source within the caustics. Additionally, the exact analytic expression for the magnification distribution for the source in the triangular caustics is derived, as well as a useful approximate expression. Finally, we find that the Chang–Refsdal lens with additional convergence greater than unity (the 'overfocusing case') can exhibit third-order critical behaviour, if the 'reduced shear' is exactly equal to     , and that the number of images for N -point masses with non-zero constant shear cannot be greater than  5 N − 1  .     

Effects of microlensing on parameters of the images seen near the critical curves of gravitational lens galaxies

We investigate the effect of microlensing on parameters of the images of distant sources seen near the critical curves of complex gravitational lenses, which are represented as a sum of compact structures—microlenses (stars, star-like or planet-like bodies) and diffusely distributed matter (dust and gas clouds etc.). The observation of merging, cross-shaped, annular, or arc-shaped source images is an indication that the images are close to the critical curves of gravitational lenses. Our analysis and numerical solution have allowed us to determine the structures of the critical curves and caustics formed by macro-and microlenses, as well as to estimate the characteristic perturbations introduced by microlenses at their various positions relative to the critical curve of a regular gravitational lens. We show that, the closer are the microlenses to the critical curve, the larger is the discrepancy between our results and those obtained previously with standard (linearized) allowance for the effect of a regular gravitational lens.     

Adaptive contouring – an efficient way to calculate microlensing light curves of extended sources

The availability of a robust and efficient routine for calculating light curves of a finite source magnified due to bending of its light by the gravitational field of an intervening binary lens is essential for determining the characteristics of planets in such microlensing events, as well as for modelling stellar lens binaries and resolving the brightness profile of the source star. However, the presence of extended caustics, and the fact that the images of the source star cannot be determined analytically while their number depends on the source position (relative to the lens system), makes such a task difficult in general. Combining the advantages of several earlier approaches, an adaptive contouring algorithm is presented, which only relies on a small number of simple rules and operations on the adaptive search grid. By using the parametric representation of critical curves and caustics found by Erdl & Schneider, seed solutions to the adaptive grid are found, which ensures that no images or holes are missed.     

Search for exotic matter from gravitational microlensing observations of stars

We consider small-scale spheroidal clusters of weakly interacting massive particles in our Galaxy as non-compact gravitational microlenses and predict the appearance of caustics in the plane of a lensed source. The crossing of these caustics by a lensed star can produce a large variety of light curves, including some observed in actual microlensing events that have been interpreted as manifestations of binary gravitational lenses. We consider also observable effects during the gravitational microlensing of stars of non-zero angular size with a given brightness distribution across their disks by such an exotic objects as natural wormholes and objects whose space-time environment is described with the NUT metric. We demonstrate that, under certain conditions, the microlensing light curves, chromatic and polarizational effects due to the properties of the lens and the star disk brightness distributions can differ considerably from those observed for a Schwarzschild gravitational lens, so that their analysis can facilitate the identification of such objects.     

基于小波方法的平滑算法在引力透镜中的应用

基于离散小波变换(DWT)方法,提出了一种可用于计算三维数值模拟样本面密度的平滑算法.为检验方法的有效性,利用该算法研究了两组不同质量解析度的引力透镜数值模拟样本,样本采用了暗物质晕的等温椭球模型,使用蒙特卡罗方法生成.计算结果表明此算法能够在很高的精度上构建引力透镜模拟样本的面密度分布轮廓,由面密度计算出来的透镜的临界曲线和焦散曲线也能较好地和理论曲线吻合,结果是令人满意的.同时比较了三组不同的小波基的计算结果,包括Daub4,Daub6和B-spline 3th,给出了最优的选择.在不损失平滑效果的同时,此算法具有非常高的速度,非常适合于处理以后更高精度的N体数值模拟.     

A genetic algorithm for the non-parametric inversion of strong lensing systems

We present a non-parametric technique to infer the projected mass distribution of a gravitational lens system with multiple strong-lensed images. The technique involves a dynamic grid in the lens plane on which the mass distribution of the lens is approximated by a sum of basis functions, one per grid cell. We used the projected mass densities of Plummer spheres as basis functions. A genetic algorithm then determines the mass distribution of the lens by forcing images of a single source, projected back on to the source plane, to coincide as well as possible. Averaging several tens of solutions removes the random fluctuations that are introduced by the reproduction process of genomes in the genetic algorithm and highlights those features common to all solutions. Given the positions of the images and the redshifts of the sources and the lens, we show that the mass of a gravitational lens can be retrieved with an accuracy of a few percent and that, if the sources sufficiently cover the caustics, the mass distribution of the gravitational lens can also be reliably retrieved. A major advantage of the algorithm is that it makes full use of the information contained in the radial images, unlike methods that minimize the residuals of the lens equation, and is thus able to accurately reconstruct also the inner parts of the lens.     

The influence of central black holes on gravitational lenses

Recent observations indicate that many if not all galaxies host massive central black holes. In this paper we explore the influence of black holes on the lensing properties. We model the lens as an isothermal ellipsoid with a finite core radius plus a central black hole. We show that the presence of the black hole substantially changes the critical curves and caustics. If the black hole mass is above a critical value, then it will completely suppress the central images for all source positions. Realistic central black holes are likely to have masses below this critical value. Even in such subcritical cases, the black hole can suppress the central image when the source is inside a zone of influence, which depends on the core radius and black hole mass. In the subcritical cases, an additional image may be created by the black hole in some regions, which for some radio lenses may be detectable with high-resolution and large dynamic range VLBI maps. The presence of central black holes should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.     

A New Smoothing Algorithm and its Application in Gravitational Lensing_

Based on the DWT (discrete wavelet transform) method, we propose a new smoothing algorithm for computing surface densities from 3D numerical simulation samples. To check its effectiveness, we have applied this algorithm to two Monte-Carlo samples of gravitational lens simulation with different mass resolutions, generated from the isothermal ellipsoid model of dark matter halos. The calculated results indicate that this algorithm can reconstruct accurately the surface density distribution of the gravitational lens simulation sample, and that the lens caustics and critical curves derived from the surface densities agree well with the theoretical curves. We have compared the results calculated by using 3 different wavelet bases (Daub4, Daub6 and B-spline 3th), and identified the best one. Without sacrificing its smoothing capability, this algorithm has a very fast computing speed, suitable for later N-body numerical simulations, which require even higher resolutions.     

Gravitational lens under perturbations: symmetry of perturbing potentials with invariant caustics

When the gravitational lensing potential can be approximated by that of a circularly symmetric system affected by weak perturbations, it is found that the shape of the resulting (tangential) caustics is entirely specified by the local azimuthal behaviour of the affecting perturbations. This provides a common mathematical groundwork for understanding problems such as the close–wide  ( d ↔ d ⁻¹)  separation degeneracy of binary lens microlensing light curves and the shear–ellipticity degeneracy of quadruple image lens modelling.     

The anatomy of a quadruply imaged gravitational lens system

The key to using a strong gravitational lens system to measure the Hubble constant is to obtain an accurate model of the lens potential. In this paper, we investigate the properties of gravitational lens B1608+656, a quadruply imaged lens system with an extended source intensity distribution. Our analysis is valid for generic quadruply lensed systems. Limit curves and isophotal separatrices are defined for such systems, and we show that the isophotal separatrices must intersect at the critical curves and the satellite isophotes must be tangent to the limit curves. The most recent model of B1608+656 by Koopmans et al. satisfies these criteria for some, but not all, of the isophotal separatrices within the observational uncertainty. We study a non-parametric method of potential reconstruction proposed by Blandford, Surpi & Kundic and demonstrate that although the method works in principle and elucidates image formation, the initial potential only converges to the true model when it is within ∼1 per cent of the true model.     

Effects of sherically-symmetric gravitational lenses produced by galaxies and clusters

Images are considered which are produced by a spherically symmetric gravitational lens from a source of small angular size. Simple analytic expressions are derived to describe the position, amplification, shape distortion and time delay for these images provided the mass distribution of the lens corresponds to the King model. Possible types of images are analysed (single and triple images, as well as degenerate double ones at caustics) which may be formed at different positions of the observer relative to the light source and lens. The focal length of galaxies of different types is evaluated on the basis of the obtained results and the conditions are determined at which these galaxies may considerably affect the images of distant sources.     

On the astrometric behaviour of binary microlensing events

Despite the suspected binarity for a significant fraction of Galactic lenses, the current photometric surveys detected binary microlensing events only for a small fraction of the total events. The detection efficiency is especially low for non-caustic crossing events, which comprise the majority of the binary lensing events, as a result of the absence of distinctive features in their light curves combined with small deviations from the standard light curve of a single point-mass event. In addition, even if they are detected, it will be difficult to determine the solution of the binary lens parameters owing to the severe degeneracy problem. In this paper, we investigate the properties of binary lensing event expected when they are astrometrically observed by using high-precision interferometers. For this, we construct vector field maps of excess centroid shifts, which represent the deviations of the binary lensing centroid shifts from those of a single lensing event as a function of source position. From the analysis of the maps, we find that the excess centroid shifts are substantial in a considerably large area around caustics. In addition, they have characteristic sizes and directions depending strongly on the source positions with respect to the caustics and the resulting trajectories of the light centroid (astrometric trajectories) have distinctive features, which can be distinguished from the deviations caused by other reasons. We classify the types of the deviations and investigate where they occur. Because of the strong dependence of the centroid shifts on the lens system geometry combined with the distinctive features in the observed astrometric trajectories, astrometric binary lensing observations will provide an important tool that can probe the properties of the Galactic binary lens population.     

The effect of a binary source companion on the astrometric microlensing behaviour

If gravitational microlensing occurs in a binary source system, both source components are magnified, and the resulting light curve deviates from the standard one of a single source event. However, in most cases only one source component is highly magnified and the other component (the companion) can be treated as a simple blending source: this is a blending approximation. In this paper we show that, unlike the light curves, the astrometric curves, representing the trajectories of the source image centroid, of an important fraction of binary source events will not be sufficiently well-modelled by the blending effect alone. This is because the centroid shift induced by the source companion endures to considerable distances from the lens. Therefore, in determining the lens parameters from astrometric curves to be measured by future high-precision astrometric instruments, it will be important to take the full effect of the source companion into consideration.     

Perturbation theory of N point-mass gravitational lens systems without symmetry: small mass-ratio approximation

This paper makes the first systematic attempt to determine using perturbation theory the positions of images by gravitational lensing due to arbitrary number of coplanar masses without any symmetry on a plane, as a function of lens and source parameters. We present a method of Taylor-series expansion to solve the lens equation under a small mass-ratio approximation. First, we investigate perturbative structures of a single-complex-variable polynomial, which has been commonly used. Perturbative roots are found. Some roots represent positions of lensed images, while the others are unphysical because they do not satisfy the lens equation. This is consistent with a fact that the degree of the polynomial, namely the number of zeros, exceeds the maximum number of lensed images if   N = 3  (or more). The theorem never tells which roots are physical (or unphysical). In this paper, unphysical ones are identified. Secondly, to avoid unphysical roots, we re-examine the lens equation. The advantage of our method is that it allows a systematic iterative analysis. We determine image positions for binary lens systems up to the third order in mass ratios and for arbitrary N point masses up to the second order. This clarifies the dependence on parameters. Thirdly, the number of the images that admit a small mass-ratio limit is less than the maximum number. It is suggested that positions of extra images could not be expressed as Maclaurin series in mass ratios. Magnifications are finally discussed.     

A small mass,large image separation gravitational lens model

It is shown that gravitational lensing could produce images of comparable brightness which are separated by an angular distance larger than the angular size of the Einstein ring of the lens distribution. Hence, specific lens configurations allow given image separations to be derived with significantly less mass than a standard single lens model. A specific example is supplied, a simple case of two properly positioned lenses acting instead of one. Observational consequences are discussed that would result if such a model was used to explain the candidate gravitational lens systems Hazard 1146+111 B, C and PKS 1145–071. Such a lens distribution might, in many cases, be verifiable with VLBI techniques.     

Asymptotic formulas for the magnification of a gravitational lens system near a fold caustic

An approximate formula for the magnification of a point source near a fold caustic obtained in the first linear caustic approximation is widely used in the theory of gravitational lens systems. Here, this formula is refined to include the post-linear terms that have been found both for a point source and for an extended Gaussian source in the absence of continuous matter on the line of sight. The formulas are reduced to a form containing three additional parameters; the derivation of nontrivial corrections requires including the expansion terms in the lens equation up to the fourth order. The modified formula for an extended source is used to analyze strong microlensing events in the gravitational lens system Q2237+0305 (the Einstein Cross). For such an event on the light curve of image C (1999, OGLE data), the corrections found are statistically significant.     

An investigation of gravitational lens determinations of H0 in quintessence cosmologies

There is growing evidence that the majority of the energy density of the Universe is not baryonic or dark matter, but rather it resides in an exotic component with negative pressure. The nature of this 'quintessence' influences our view of the Universe, modifying angular diameter and luminosity distances. Here, we examine the influence of a quintessence component upon gravitational lens time-delays. As well as a static quintessence component, an evolving equation of state is also considered. It is found that the equation of state of the quintessence component and its evolution influence the value of the Hubble constant derived from gravitational lenses. However, the differences between evolving and non-evolving cosmologies are relatively small. We undertake a suite of Monte Carlo simulations to examine the potential constraints that can be placed on the universal equation of state from the monitoring of gravitational lens systems, and demonstrate that at least an order of magnitude more lenses than currently known will have to be discovered and analysed to accurately probe any quintessence component.     

A simple model for lensing by black holes in galactic nuclei

The lensing properties of the Plummer model with a central point mass and external shear are derived, including the image multiplicities, critical curves and caustics. This provides a simple model for a flattened galaxy with a central supermassive black hole. For the Plummer model with black hole, the maximum number of images is four, provided the black hole mass is less than an upper bound which is calculated analytically. This introduces a method to constrain black hole masses by counting images, thus applicable at cosmological distance. With shear, the maximum number of images is six and we illustrate the occurrence of an astroid caustic and two metamorphoses.     

Gravitational lens effect of degenerate neutrino celestial bodies

We took the degenerate neutrino celestial body (NCO) to be transparent. We solved the geodesic equation for the light rays emitted by sources inside and outside the body. We then derived the expression for the amplication factor K by a transparent gravitational lens, and evaluated K and the redshift under the metric of the NCO. Our calculations show that, for light travelling in certain directions from sources in certain positions, the brightness may be amplified several hundred times. We point out that some quasars may be galactic nuclei inside NCOs.     

Chromaticity of gravitational microlensing events

In this paper, we investigate the colour changes of gravitational microlensing events caused by the two different mechanisms of differential amplification for a limb-darkened extended source and blending. From this investigation, we find that the colour changes of limb-darkened extended source events (colour curves) have dramatically different characteristics depending on whether the lens transits the source star or not. We show that for a source transit event, the lens proper motion can be determined by simply measuring the turning time of the colour curve instead of fitting the overall colour or light curves. We also find that even for a very small fraction of blended light, the colour changes induced by blending are equivalent to those induced by limb darkening, causing serious distortion in the observed colour curve. Therefore, to obtain useful information about the lens and source star from the colour curve of an event, it will be essential to correct for blending. We discuss various methods of blending correction .