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.     

Planetary microlensing signals from the orbital motion of the source star around the common barycentre

With several detections, the technique of gravitational microlensing has proven useful for studying planets that orbit stars at Galactic distances, and it can even be applied to detect planets in neighbouring galaxies. So far, planet detections by microlensing have been considered to result from a change in the bending of light and the resulting magnification caused by a planet around the foreground lens star. However, in complete analogy to the annual parallax effect caused by the revolution of the Earth around the Sun, the motion of the source star around the common barycentre with an orbiting planet can also lead to observable deviations in microlensing light curves that can provide evidence for the unseen companion. We discuss this effect in some detail and study the prospects of microlensing observations for revealing planets through this alternative detection channel. Given that small distances between lens and source star are favoured, and that the effect becomes nearly independent of the source distance, planets would remain detectable even if their host star is located outside the Milky Way with a sufficiently good photometry (exceeding present-day technology) being possible. From synthetic light curves arising from a Monte Carlo simulation, we find that the chances for such detections are not overwhelming and appear practically limited to the most massive planets (at least with current observational set-ups), but they are large enough for leaving the possibility that one or the other signal has already been observed. However, it may remain undetermined whether the planet actually orbits the source star or rather the lens star, which leaves us with an ambiguity not only with respect to its location, but also to its properties.     

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.     

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 .     

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.     

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.     

Variation of spot-induced anomalies in caustic-crossing binary microlensing event light curves

We investigate the pattern of anomalies in the light curves of caustic-crossing binary microlensing events induced by spot(s) on the lensed source star. To this end, we perform simulations of events with various models of spots. From these simulations we find that the spot-induced anomalies take various forms depending on the physical state of spots, which is characterized by the surface brightness contrast, the size, the number, the umbra/penumbra structure, the shape and the orientation with respect to the sweeping caustic. We also examine the feasibility of distinguishing the two possibly degenerate types of anomalies caused by a spot and a transiting planet and find that in many cases the degeneracy can be separated from the characteristic multiple deviation feature in the spot-induced anomaly pattern caused by the multiplicity of spots.     

Extended source effects in astrometric gravitational microlensing

Extended source size effects have been detected in photometric monitoring of gravitational microlensing events. We study similar effects in the centroid motion of an extended source lensed by a point mass. We show that the centroid motion of a source with uniform surface brightness can be obtained analytically. For a source with a circularly symmetric limb-darkening profile, the centroid motion can be expressed as a one-dimensional integral, which can be evaluated numerically. We find that when the impact parameter is comparable to the source radius, the centroid motion is significantly modified by the finite source size. In particular, when the impact parameter is smaller than the source radius, the trajectories become clover-leaf like. Such astrometric motions can be detected using space interferometers such as the Space Interferometry Mission . Such measurements offer exciting possibilities for determining stellar parameters, such as stellar radius, to excellent accuracy.     

Microlensing of an extended source by a power-law mass distribution

Microlensing promises to be a powerful tool for studying distant galaxies and quasars. As the data and models improve, there are systematic effects that need to be explored. Quasar continuum and broad-line regions may respond differently to microlensing due to their different sizes; to understand this effect, we study microlensing of finite sources by a mass function of stars. We find that microlensing is insensitive to the slope of the mass function but does depend on the mass range. For negative-parity images, diluting the stellar population with dark matter increases the magnification dispersion for small sources and decreases it for large sources. This implies that the quasar continuum and broad-line regions may experience very different microlensing in negative-parity lensed images. We confirm earlier conclusions that the surface brightness profile and geometry of the source have little effect on microlensing. Finally, we consider non-circular sources. We show that elliptical sources that are aligned with the direction of shear have larger magnification dispersions than sources with perpendicular alignment, an effect that becomes more prominent as the ellipticity increases. Elongated sources can lead to more rapid variability than circular sources, which raises the prospect of using microlensing to probe source shape.     

Probing subparsec structure in the Lyman α forest with gravitational microlensing

We present the results of microlens ray-tracing simulations showing the effect of absorbing material between a source quasar and a lensing galaxy in a gravitational lens system. We find that, in addition to brightness fluctuations due to microlensing, the strength of the absorption line relative to the continuum varies with time, with the properties of the variations depending on the structure of the absorbing material. We conclude that such variations will be measurable via ultraviolet spectroscopy of image A of the gravitationally lensed quasar Q2237+0305 if the Lyman α clouds between the quasar and the lensing galaxy possess structure on scales smaller than ∼0.1 pc. The time-scale for the variations is on the order of years to decades, although very short-term variability can occur. While the Lyman α lines may not be accessible at all wavelengths, this approach is applicable to any absorption system, including metal lines.     

Analysis of the Q2237+0305 light-curve variability with regularization technique

The microlensing high-amplification events in the light curves of the gravitationally lensed quasar Q2237+0305 observed by the OGLE group and GLITP collaboration are analysed. The significant brightness amplification in the A and C components in 1999 observational season are considered under the assumption of the fold caustic crossing. Under this assumption we applied the model-independent method based on regularization technique for one-dimensional profile restoration of the quasar accretion disc brightness distribution. The recovered brightness distribution of the source seems to obey the standard model for the accretion disc. The estimated effective radius of the quasar emitting region is in agreement with the previous studies.     

Another channel to detect close-in binary companions via gravitational microlensing

Gaudi & Gould showed that close companions of remote binary systems can be efficiently detected by using gravitational microlensing via the deviations in the lensing light curves induced by the existence of the lens companions. In this paper, we introduce another channel to detect faint close-in binary companions by using microlensing. This method utilizes a caustic-crossing binary lens event with a source also composed of binary stars, where the companion is a faint star. Detection of the companion is possible because the flux of the companion can be highly amplified when it crosses the lens caustic. The detection is facilitated since the companion is more amplified than the primary because it, in general, has a smaller size than the primary, and thus experiences less finite source effect. The method is an extension of the previous one suggested to detect close-in giant planets by Graff & Gaudi and Lewis & Ibata and further developed by Ashton & Lewis. From the simulations of realistic Galactic bulge events, we find that companions of K-type main-sequence or brighter stars can be efficiently detected from the current type of microlensing follow-up observations by using the proposed method. We also find that compared with the method of detecting lens companions for which the efficiency drops significantly for binaries with separations ≲0.2 of the angular Einstein ring radius, θ _E, the proposed method has an important advantage of being able to detect companions with substantially smaller separations down to ∼     .     

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.     

Gravitational microlensing of gamma-ray bursts at medium optical depth

Gravitational lensing of a gamma-ray burst (GRB) by a single point mass will produce a second, delayed signal. Several authors have discussed using microlensed GRBs to probe a possible cosmological population of compact objects. We analyse a closely related phenomenon: the effect of microlensing by low to medium optical depth in compact objects on the averaged observed light curve of a sample of GRBs. We discuss the cumulative measured flux as a function of time resulting from delays caused by microlensing by cosmological compact objects. The time-scale and curvature of this function describe unique values for the compact object mass and optical depth. For GRBs with durations larger than the detector resolution, limits could be placed on the mass and optical depth of cosmological compact objects. The method does not rely on the separation of lensed bursts from those that are spatially coincident.     

The distribution of microlensed light-curve derivatives: the relationship between stellar proper motions and transverse velocity

We present a method for computing the probability distribution of microlensed light-curve derivatives both in the case of a static lens with a transverse velocity, and in the case of microlensing that is produced through stellar proper motions. The distributions are closely related in form, and can be considered equivalent after appropriate scaling of the input transverse velocity. The comparison of the distributions in this manner provides a consistent way to consider the relative contribution to microlensing (both large and small fluctuations) of the two classes of motion, a problem that is otherwise an extremely expensive computational exercise. We find that the relative contribution of stellar proper motions to the microlensing rate is independent of the mass function assumed for the microlenses, but is a function of optical depth and shear. We find that stellar proper motions produce a higher overall microlensing rate than a transverse velocity of the same magnitude. This effect becomes more pronounced at higher optical depth. With the introduction of shear, the relative rates of microlensing become dependent on the direction of the transverse velocity. This may have important consequences in the case of quadruply lensed quasars such as Q2237+0305, where the alignment of the shear vector with the source trajectory varies between images.     

Identification of bright lenses from the astrometric observations of gravitational microlensing events

When a source star is gravitationally microlensed by a dark lens, the centroid of the source star image is displaced relative to the position of the unlensed source star, with an elliptical trajectory. Recently, routine astrometric follow-up measurements of these source star image centroid shifts by using high-precision interferometers have been proposed to measure the lens proper motion, which can resolve the lens parameter degeneracy in the photometrically determined Einstein time-scale. When an event is caused by a bright lens, on the other hand, the astrometric shift is affected by the light from the lens, but one cannot identify the existence of the bright lens from the observed trajectory because the resulting trajectory of the bright lens event is also an ellipse. As results, lensing parameters determined from the trajectory differ from those of a dark lens event, causing an incorrect identification of the lens population. In this paper, we show that although the shape and size of the astrometric centroid shift trajectory are changed because of the bright lens, the angular speed of centroid shifts around the apparent position of the unlensed source star is not affected by the lens brightness. Therefore, one can identify the existence of a bright lens and determine its brightness by comparing the lens parameters determined from the 'angular speed curve' with those determined from the trajectory of observed centroid shifts. Once the lens brightness is determined, one can correct for the lens proper motion. As the proposed method provides information about both the lens brightness (dark or bright) and the corrected values of the physical parameters of the lens, one can constrain the nature of massive compact halo objects (MACHOs) significantly better.     

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.     

Pixel lensing as a way to detect extrasolar planets in M31

We study the possibility to detect extrasolar planets in M31 through pixel-lensing observations. Using a Monte Carlo approach, we select the physical parameters of the binary lens system, a star hosting a planet, and we calculate the pixel-lensing light curve taking into account the finite source effects. Indeed, their inclusion is crucial since the sources in M31 microlensing events are mainly giant stars. Light curves with detectable planetary features are selected by looking for significant deviations from the corresponding Paczyński shapes. We find that the time-scale of planetary deviations in light curves increase (up to 3–4 d) as the source size increases. This means that only few exposures per day, depending also on the required accuracy, may be sufficient to reveal in the light curve a planetary companion. Although the mean planet mass for the selected events is about     , even small mass planets  ( M _P < 20 M_⊕)  can cause significant deviations, at least in the observations with large telescopes. However, even in the former case, the probability to find detectable planetary features in pixel-lensing light curves is at most a few per cent of the detectable events, and therefore many events have to be collected in order to detect an extrasolar planet in M31. Our analysis also supports the claim that the anomaly found in the candidate event PA-99-N2 towards M31 can be explained by a companion object orbiting the lens star.     

Analytic relations between the observed gravitational microlensing parameters with and without the effect of blending

When a microlensing light curve is contaminated by blended light from unresolved stars near the line of sight to the lensed star, the light curve shape and corresponding parametrization for the event will differ from the values expected when the event is not affected by blending. As a result, blending makes it difficult to identify the major lens population and to estimate the amount of lensing matter. In order to estimate the effect of blending on the result of lensing experiments, it is, therefore, essential to know how the observed lensing parameters change depending on the fraction of blended light. Previously, the changed lensing parameters were obtained with a statistical method that not only required a large amount of computation time but also was prone to uncertainty. In this paper, we derive analytic relations between the lensing parameters with and without the effect of blending. By using these relations, we investigate the dependence of the observed lensing parameters on the amount of blended light, the impact parameter and the threshold amplification for event detection.     

Parameter degeneracies and (un)predictability of gravitational microlensing events

Some of the difficulties in determining the underlying physical properties that are relevant for observed anomalies in microlensing light curves, such as the mass and separation of extrasolar planets orbiting the lens star, or the relative source–lens parallax, are already anchored in factors that limit the amount of information available from ordinary microlensing events and in the way these are being parametrized. Moreover, a real-time detection of deviations from an ordinary light curve while these are still in progress can only be done against a known model of the latter, and such is also required for properly prioritizing ongoing events for monitoring in order to maximize scientific returns. Despite the fact that ordinary microlensing light curves are described by an analytic function that only involves a handful of parameters, modelling these is far less trivial than one might be tempted to think. A well-known degeneracy for small impacts, and another one for the initial rise of an event, makes an interprediction of different phases impossible, while in order to determine a complete set of model parameters, the fundamental characteristics of all these phases need to be properly assessed. While it is found that the wing of the light curve provides valuable information about the time-scale that absorbs the physical properties, the peak flux of the event can be meaningfully predicted only after about a third of the total magnification has been reached. Parametrizations based on observable features not only ease modelling by bringing the covariance matrix close to diagonal form, but also allow good predictions of the measured flux without the need to determine all parameters accurately. Campaigns intending to infer planet populations from observed microlensing events need to invest some fraction of the available time into acquiring data that allow to properly determine the magnification function.