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.     

Effects of a deformation of a star on the gravitational lensing

We study analytically a gravitational lens due to a deformed star, which is modelled by using a monopole and a quadrupole moment. Positions of the images are discussed for a source on the principal axis. We present explicit expressions for the lens equation for this gravitational lens as a single real 10th-order algebraic equation. Furthermore, we compute an expression for the caustics as a discriminant for the polynomial. Another simple parametric representation of the caustics is also presented in a more tractable form. A simple expression for the critical curves is obtained to clarify a topological feature of the critical curves; the curves are simply connected if and only if the distortion is sufficiently large.     

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.     

临界曲线附近的微引力透镜效应模拟

在光滑物质分布模型下,临界曲线是强引力透镜系统中像平面上一条放大率为无穷的线,而考虑少量离散质量的微透镜效应后,源平面上的放大率分布会出现复杂的结构,为暗物质成分的探测提供了一种有效途径.模拟临界曲线附近微透镜效应存在临界曲线上放大率无穷大和计算量巨大的困难.要达到所需的模拟精度,直接使用传统的光线追踪算法需要巨大的计算资源.为此发展了一个能实现海量计算的Graphics Processing Unit (GPU)并行方法来模拟临界曲线附近的微引力透镜效应.在型号为NVIDIA Tesla V100S PCIe 32 GB的GPU上,对于需要处理13000多个微透镜天体、发射1013量级光线的模拟,耗时在7000 s左右.在GPU并行的基础上,与直接的光线追踪算法相比,插值近似的引入使计算速度提升约两个数量级.利用该方法生成80个放大率分布图,并从中抽取800条光变曲线,进行了微焦散线数密度和峰值放大率的统计.     

Non-parametric inversion of gravitational lensing systems with few images using a multi-objective genetic algorithm

Galaxies acting as gravitational lenses are surrounded by, at most, a handful of images. This apparent paucity of information forces one to make the best possible use of what information is available to invert the lens system. In this paper, we explore the use of a genetic algorithm to invert in a non-parametric way strong lensing systems containing only a small number of images. Perhaps the most important conclusion of this paper is that it is possible to infer the mass distribution of such gravitational lens systems using a non-parametric technique. We show that including information about the null space (i.e. the region where no images are found) is prerequisite to avoid the prediction of a large number of spurious images, and to reliably reconstruct the lens mass density. While the total mass of the lens is usually constrained within a few per cent, the fidelity of the reconstruction of the lens mass distribution depends on the number and position of the images. The technique employed to include null space information can be extended in a straightforward way to add additional constraints, such as weak-lensing data or time-delay information.     

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.     

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 .     

Background-field effects in astrometric microlensing

We consider the motion of the image of a distant radiation source microlensed by one and many point masses in the presence of an external background field, which produces a shear in the lens equation. The trajectories of the image brightness centroid (BC) in a Chang-Refsdal lens are very sensitive to the magnitude of the shear and to the direction of motion of the microlenses, allowing these parameters to be unambiguously determined from the BC motion. Stochastic motions of the BC for an extended source are simulated in a problem with a large number of randomly arranged microlenses. We derive dependences of the root-mean-square variations in BC velocity on microlense density, on external shear, and on the continuous component of the mass density, which may be related to dark matter in extragalactic lenses.     

A new gravitational lens from the MUSCLES survey: ULAS J082016.1+081216

We present observations of a new double-image gravitational lens system, ULAS J082016.1+081216, of image separation 2.3 arcsec and high (∼6) flux ratio. The system is selected from the Sloan Digital Sky Survey (SDSS) spectroscopic quasar list using new high-quality images from the UKIRT (United Kingdom Infrared Telescope) Deep Sky Survey (UKIDSS). The lensed quasar has a source redshift of 2.024, and we identify the lens galaxy as a faint red object of redshift  0.803 ± 0.001  . Three other objects from the UKIDSS survey, selected in the same way, were found not to be lens systems. Together with the earlier lens found using this method, the SDSS–UKIDSS lenses have the potential to significantly increase the number of quasar lenses found in SDSS, to extend the survey to higher flux ratios and lower separations, and to give greater completeness which is important for statistical purposes.     

Colour change measurements of gravitational microlensing events by using the difference image analysis method

Detecting colour changes of a gravitational microlensing event induced by the limb-darkened extended source effect is important for obtaining useful information about both the lens and the source star. However, precise measurements of the colour changes are hampered by blending, which also causes colour changes of the event. In this paper, we show that although the colour change measured from the subtracted image by using the recently developed photometric method of the 'difference image analysis' (DIA) differs from the colour change measured by using the conventional method based on the extraction of the individual source stars' point spread functions, the curve of the colour changes (colour curve) constructed by using the DIA method enables one to obtain the same information about the lens and source star, but with significantly reduced uncertainties due to the absence of blending. We investigate the patterns of the DIA colour curves for both single lens and binary lens events by constructing colour change maps.     

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.     

ULAS J234311.93-005034.0: a gravitational lens system selected from UKIDSS and SDSS

We report the discovery of a new gravitational lens system. This object, ULAS J234311.93-005034.0, is the first to be selected by using the new UKIRT Infrared Deep Sky Survey (UKIDSS), together with the Sloan Digital Sky Survey (SDSS). The ULAS J234311.93-005034.0 system contains a quasar at redshift 0.788 which is doubly imaged, with separation 1.4 arcsec. The two quasar images have the same redshift and similar, though not identical, spectra. The lensing galaxy is detected by subtracting point spread functions from R -band images taken with the Keck telescope. The lensing galaxy can also be detected by subtracting the spectra of the A and B images, since more of the galaxy light is likely to be present in the latter. No redshift is determined from the galaxy, although the shape of its spectrum suggests a redshift of about 0.3. The object's lens status is secure, due to the identification of two objects with the same redshift together with a lensing galaxy. Our imaging suggests that the lens is found in a cluster environment, in which candidate arc-like structures, that require confirmation, are visible in the vicinity. Further discoveries of lenses from the UKIDSS survey are likely as part of this programme, due to the depth of UKIDSS and its generally good seeing conditions.     

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.     

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

基于离散小波变换(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.     

General polytropic spheres as gravitational lenses

We propose hydrostatic polytropic spheres governed by the Lane-Emden equation (LEE) of index $n$ as a novel set of physical models for axially averaged gravitational lenses anywhere in the Universe, alternative to the familiar singular isothermal sphere (SIS) and the Navarro–Frenk–White (NFW) profile, as such general polytropic spheres are conceptually simple, versatile in representing a series of equations of state, and able to address both the inner core and cusp features. As LEE is nonlinear, there exist several distinct classes of LEE solutions to serve as physical lens models. With a few scaling parameters, the complete problem can be readily reconstructed with full physical dimensions. A given mass density profile satisfying LEE produces lensing effects that are solely determined by a dimensionless parameter $q$ which contains geometric and kinematic information about the source-lens-observer system. The lens mapping and tangential shear or distortion profile are derived, first analytically for special cases and then asymptotically at the outskirts or near the edge of the lens. Numerical procedures for calculating full lensing profiles of a general lens are developed. Our results include the analytical “singular polytropic sphere” (SPS) profile which generalizes the SIS model and may outperform the latter in modeling dark matter halos among others. We further point out that dynamic models of general polytropic spheres in self-similar evolution can serve as several broad classes of gravitational lenses and produce time-dependent lensing effects slow or fast depending on the pertinent time scales. Astrophysical sources that can be lensed include electromagnetic wave sources in the entire frequency band, gravitational wave sources in the entire frequency band, gravitons even possibly with finite masses, neutrino sources of three different types, neutron sources, and ultra high energy cosmic rays (UHECRs) of electrically charged particles which can also interact with magnetic fields. We discuss and elabrate applications to dark matter halos, hypermassive black holes and supermassive black holes in the entire Universe including the early Universe, magnetized supermassive stars, static and dynamically evolving spherical and cylindrical lenses in contexts of astrophysics and cosmology.     

Astrometric microlensing: a channel to detect multiple lens systems

If a source star is gravitationally microlensed by a multiple lens system, the resulting light curve can have significant deviations from the standard form of a single lens event. The chance of producing significant deviations becomes important when the separations between the component lenses are equivalent to the combined angular Einstein ring radius of the system. For multiple lens systems composed of more than two lenses, however, this condition is difficult to meet because the orbits of such systems are unstable. Even if events are caused by a multiple lens system with stable orbits where a pair of lenses are closely located and the other component (a third body) has a wide separation from the pair, identifying the lens multiplicity photometrically will be difficult because the event will be identified by either a binary lens event caused by the close pair of lenses or a single lens event caused by the third body. In this paper, we show that if a seemingly binary lens event is followed up astrometrically using future high-precision interferometers, the existence of an additional third body can be identified via a repeating event. We show that the signatures of third bodies can be unambiguously identified from the characteristic distortions they make in the centroid shift trajectories. We also show that owing to the long-range astrometric effect of third bodies, the detection efficiency will be considerable even for third bodies with large separations from their close lens pairs.     

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.     

Submillimetre observations of luminous z>4 radio-quiet quasars and the contribution of AGN to the submm source population

We present polarization observations of the gravitational lens system B1422+231 made at 8.4 GHz using the VLBA and the 100-m telescope at Effelsberg. All four images of the quasar show structure on the milliarcsec scale. The three bright images show tangential stretching as expected from lens models. Some basic properties of gravitational lensing are exhibited by this system. The surface brightness of images A and B are the same and the parity reversal expected in image B is revealed, for the first time, by polarization observations. There is a large differential Faraday rotation between images A and B.     

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.