ISRN Astronomy and Astrophysics

Volume 2012, Article ID 906951, 21 pages

http://dx.doi.org/10.5402/2012/906951

## Analytical Methods in Gravitational Microlensing

Astronomical Observatory, Taras Shevchenko National University of Kyiv, Observatorna Street 3, Kiev 04053, Ukraine

Received 30 April 2012; Accepted 17 June 2012

Academic Editors: J. Gallimore and M. Sazhin

Copyright © 2012 V. I. Zhdanov et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We discuss analytical results dealing with photometric and astrometric gravitational microlensing. The first two sections concern approximation methods that allow us to get solutions of the general lens equation near fold caustics and cusp points up to any prescribed accuracy. Two methods of finding approximate solutions near the fold are worked out. The results are applied to derive new corrections to total amplifications of critical source images. Analytic expressions are obtained in case of the Gaussian, power-law, and limb-darkening extended source models; here we present the first nonzero corrections to the well-known linear caustic approximation. Possibilities to distinguish different source models in observations are discussed on the basis of statistical simulations of microlensed light curves. In the next section, we discuss astrometric microlensing effects in various cases of extended sources and extended lenses, including a simple model of weak statistical microlensing by extended dark matter clumps. Random walks of a distant source image microlensed by stochastic masses are estimated. We note that the bulk motion of foreground stars induces a small apparent rotation of the extragalactic reference frame. Compact analytical relations describing the statistics of such motions are presented.

#### 1. Introduction

Gravitational lensing theory contains a number of nice analytical results and beautiful theorems, which give a description of the observed effects on the level of quality and form a basis for quantitative calculations [1, 2]. This theory is based on the general relativity that describes the propagation of electromagnetic radiation in a curved space-time. Two effects are most important: the Einsteinian light bending in the gravitational field and the gravitational time delay. These effects are described by simple formulas that can be found in any textbook on relativistic gravity. However, their use in a particular gravitational lens system (GLS) can be rather complicated, and there is a considerable need for improvement of numerical methods and analytical studies. The present article deals mainly with some of the topics of the latter subject. We believe that there exists a considerable number of interesting problems concerning GLS that can be solved either without using extensive computations or where such computations can be reduced to a minimum.

One of the most spectacular effects in observations of GLS is a considerable brightness enhancement in one of the source images due to the microlensing [1–3]; this is often referred to as a high amplification event (HAE). In an extragalactic system, this is typically associated with the source crossing of the GLS caustic. It is well known that the only stable caustics of the two-dimensional mapping are folds and cusps [1, 3]. The crossing of the fold caustic is most probable (see, e.g., [4] in case of Q2237 + 0305), though the role of the cusps is also significant [1].

Below we discuss approximation methods that allow us to get solutions of the lens equation near caustics with any prescribed accuracy. We consider both the cases of a fold caustic and a cusp point using a general procedure of a power reduction. A special attention is paid to the case of the fold. Application of this procedure allows us to find an analytic form of the solution which is used further. In this view we discuss two methods of approximate solution of the lens equation [5, 6]. Both methods agree with each other in a common domain of validity; they are used to obtain first orders of expansion of the approximate solutions. The results are applied to derive the amplification of a point source near the fold and then to obtain the amplifications of extended sources (the Gaussian, power-law and limb-darkening source models). In order to obtain nonzero corrections to the total amplification of two critical images, it was necessary to take into consideration additional higher-order terms in the expansion of the lens mapping in comparison with earlier works (see, e.g., [7]). In Section 3.5, we discuss the possibility to distinguish different source models on the basis of computer simulations of microlensed light curves.

The following sections deal with astrometric microlensing. Here, we develop some of our earlier results [8–14]. We show that the trajectory of the image centroid depends strongly on the source size and on an external gravitational field (external shear); a microlensing by an extended lenses is also discussed including a simple “toy” model of microlensing by extended dark matter clumps. Then, we consider a statistics of weak astrometric microlensing of a distant source by a large number of foreground point masses in case of a small optical depth; this case corresponds to a weak microlensing by the Milky Way objects. Variable gravitational field of moving foreground masses induces small motions of remote sources’ images. We present compact analytical relations describing a statistics of these motions; it is different in case of a continuous (e.g., dark) matter and discrete masses (stars). This effect leads to a small apparent rotation of an extragalactic reference frame.

#### 2. Analytic Theory of Caustic Crossing Events

Investigation of the lens mapping near singular points [3] is an important part of the HAE studies. The lens equation near a fold may be expanded in powers of local coordinates; in the lowest order of this expansion the caustic is represented by a straight line; so this approximation is often referred to as “the linear caustic approximation.” In this approximation, the point source flux amplification depends on the distance to the caustic and contains two parameters [1]. In most cases, the linear caustic approximation is sufficient to treat the observed light curves over the range of HAEs at modern accuracy of the flux measurement. The need for a modification of this formula, for example, by taking into account the caustic curvature, is nevertheless being discussed for a long time [15–17]. One may hope that future improvement of the photometric accuracy will make it possible to obtain additional parameters of the lens mapping, which are connected to the mass distribution in a lensing galaxy. At the same time, consideration of the “postlinear” terms is sometimes sensible in order to explain even the presently available observational data [5, 6]. Note also that the corrections to the amplification on a macrolensing level were the subject of investigations in connection with the problem of “anomalous flux ratios” [7].

##### 2.1. Initial Equations and Notations

In this subsection, we recall some general notions of the gravitational lensing that may be found for example, in book [1]. We use the normalized lens equation in the form: is the lens mapping potential; this equation relates every point of the image plane to the point of the source plane. In the general case, there are several solutions of the lens equation (1) that represent images of one point source at ; we denote solution number with the index in parentheses.

In case of no continuous matter on the line of sight, the potential must be a harmonic function . Below, we will assume that this condition is fulfilled in the considered neighborhood of the critical point. We note, however, if we suppose that during HAE the continuous matter density is constant, the lens equation may be represented in the same form by a suitable renormalization of the variables.

The amplification of a separate image of a point source is , where is the Jacobian of the lens mapping. The critical curves of the lens mapping (1) are defined by equation ; they are mapped onto the caustics on the source plane. The stable critical points of two-dimensional mapping may be folds and cusps only, the folds being more probable in HAE. When a point source approaches the fold caustic from its convex side, two of its images approach the critical curve and their amplification tends to infinity; they disappear after the source crossing of the caustic. These two images are called critical. In case of the cusp point, there can be either three or one bright images near the critical curve depending on location of the source inside or outside of the caustic.

##### 2.2. Methods of Approximate Solution of the Lens Equation Near the Fold

Below, we outline general procedures that allow to find approximate solutions of the lens equation near folds and cusps, the main attention is being paid to the case of the folds. The standard consideration of the caustic crossing events deals with the Taylor expansion of the potential near some point of the critical curve in the image plane. Let this point be the coordinate origin and we suppose that (1) maps onto the coordinate origin on the source plane. We assume that the lens equation gives the coordinates of a point source as analytic functions of its image coordinates measured from the critical point . The problem is to construct the inverse dependence of the image position as functions of the source coordinates. Since the Jacobian is zero at the critical point, these functions cannot be analytic.

We consider two different methods. The first one deals with analytical expansions in powers of a small parameter, however, it results in nonanalytical functions of coordinates leading to nonintegrable terms in the amplification. The second method does not lead to such problems though it uses a somewhat more complicated representation of the solutions of the lens equation (containing roots of analytical functions). Both methods agree with each other in a common domain of validity; moreover, we use the second method to justify some expressions in the amplification formulas in terms of distributions to validate applications to the extended source models.

###### 2.2.1. Expansion of the Lens Mapping Near the Fold

We consider the folds in case of a harmonic potential, that is in case of no continuous matter on the line of sight. Nevertheless, we note that in a more general case we obtain analogous results [18].

For the harmonic potential near the fold, one can write (see, e.g., [1]) are the expansion coefficients. If the axis is directed toward the convexity of the caustic, then (at the fold points .

First, we will use a regular procedure [6] in order to construct solutions of (2) with a desired accuracy. This procedure is useful to study the light curve of the point source, which has a trajectory crossing the fold caustic under some nonzero angle. For , we substitute where and may be considered as a (small) parameter describing a vicinity to the caustic. In fact for small and for an appropriate choice of , we can work with . This is a formal substitution that makes easier operations with different orders of the expansion. On the other hand, for small , the value also may be considered as “small”; in this case, . Therefore, in the other treatment, after performing calculations, we can put and thus return to the initial variables.

Substitution of (3) into (2) yields where are analytical functions of all their variables. However, it should be emphasized that in fact we need finite orders of for some approximation order; therefore, we always deal with polynomials of a finite order in powers of .

###### 2.2.2. Direct Expansion of a Solution in Powers of

System (4) can be considered for fixed and variable , then the solution of (4) forms a curve which is useful to study a local behavior of critical image trajectories [19]. In this case, corresponds to a crossing of a caustic by a point source and may be considered as a time counted from the moment when two critical images appear. The results of [19] show that the solutions of (2) in terms of can be represented as the expansions in powers of . This can be seen directly for if we rewrite (4) in the form ready for iterations: , where either or is corresponding to different critical images.

Therefore, from the very beginning, we can look for the solutions in the form:

From zeroth iteration, we have The first iteration yields whence we obtain and so on.

Evidently, for every step of the iterative procedure, we obtain some polynomials of , but the dependence on is not analytical one. Moreover, the higher orders of approximation involve higher orders of leading to formally nonintegrable terms in the amplification factor (see below). This is a well-known phenomenon when we represent a nonanalytical expression (e.g., square root) as formal asymptotic expansion in powers of some parameter (e.g., expansion in powers of : ).

###### 2.2.3. Alternative Representation of Solutions

Now, we aim to obtain the other representation of the solution as an explicit algebraic function of analytical expansions in powers of coordinates of the source plane. For this purpose we use a power reduction which allows to obtain a simple approximate form of the lens equation on every step of the procedure.

The first equation of system (5) for sufficiently small can be solved iteratively with respect to yielding a solution in the form: In fact, we need finite orders of approximations, so only finite number of the iterations are involved. Then, after appropriate truncation, becomes a finite-order polynomial of all the arguments. Substitution of (9) into the second equation of system (4) leads to one equation with respect to one variable : with (again after some truncation) a polynomial of .

Now, we use the power reduction procedure so as to eliminate the terms containing with and to obtain a quadratic equation with respect to . We put and write the terms containing as Then, in the r.h.s. of (10), we can write ( !) where to be specific we have truncated the order of the expansion by some degree , coefficients are polynomial functions of . Thus, (10) can be written as . Note that the right-hand side is a linear function of .

Now, we solve (13) iteratively with respect to , and every step of the iteration procedure being combined with reduction (11) so as to exclude powers of larger than 1. In the order ~, we have In the next order, we must substitute (14) into the right-hand side of (13). This will yield the higher orders of ; these orders must be reduced making use of substitution (11) and the result of previous iterations. Thus, we again have an equation analogous to (14) with a linear function of in r.h.s. At the end, we obtain with analytical functions of their arguments.

At last, we return from to to obtain a quadratic equation: yielding a solution: The above approximation procedure can be fulfilled to yield approximate solutions with any degree of accuracy required. Note that in reality the calculations up to the order ~ involve a very limited number of coefficients after truncation of higher order terms in on every step of approximation. After derivation of , we find owing to (9).

From these considerations, it is easy to see that analytical structure of the solutions of the lens Equations (2) near the fold is where are analytical functions of . This allows us to seek for the solution in this form from the very beginning by substitution of (18) into the initial equations. This will be considered in more detail in Section 2.4.

##### 2.3. Approximation Procedure Near the Cusp

The coefficients of Taylor expansion in (2) depend on the caustic point near which the expansion is performed. If this point is a cusp, then and the above approximation methods of Sections 2.2.2 and 2.2.3 are no longer valid. In this case, the lens equation is

The lower order terms correspond to a familiar form of the lens equation near the cusp as described in [1, 21].

In order to perform a formal expansion, we write, instead of (3), Then the substitution into (19) yields, up to the first-order terms in ,

A detailed analysis of properties of this lens mapping in zeroth approximation (with respect to ) was carried out in [21, 22]. First-order corrections in a general case when the matter density in a vicinity of the line of sight is not zero were found in [23].

The solution can be sought by using the expansion in powers of . Here, we outline an alternative method using the power reduction like one described in Section 2.2.3, which allows determining an analytical structure of the solutions. First one must express by means of the other variables; this can be done iteratively using the first equation of system (21). Practically, we need a finite number of iterations. Then, we apply a power reduction to eliminate terms with so as to obtain a cubic equation for . We put On account of the second equation of system (21), after substitutions, we get an equation for of the form: Using an iterative procedure combined with (22), this enables us to obtain as follows: where are polynomial functions of and , the order of the polynomials depending upon the order of the approximation required. Thus, we come to the cubic equation with respect to : The roots can be obtained via the Cardano-Tartaglia formulas or the Viete trigonometric solution. Depending on sign of the discriminant of this equation, there are either three real roots (inside the cusp), or one real root (outside the cusp). The analytic properties of the solution are completely defined by the coefficients of (25), in particular, for , they are defined by and . It is important that for these coefficients are finite order polynomials of . This structure will be preserved in any order of approximation leading to (25).

##### 2.4. Amplification Near the Fold up to the Terms ~

###### 2.4.1. Solutions of the Lens Equation

Further, we deal with the results concerning the folds that follow from Sections 2.2.2 and 2.2.3.

The result of the expansions near the fold: as applied to the lens equation (2) in zeroth order is given by (7):

The next orders are as follows: where . The solutions up to this accuracy level have been obtained earlier [7, 19]. The contributions of this order are cancelled in calculations of the total amplification of two critical images. Therefore, to get a nontrivial correction to zero order amplification, we need the higher order approximations.

The second-order terms contain an expression, which is singular in :

Now, we proceed to second approach to construct approximate solutions of the lens equation in a vicinity of folds, which is free from such singularities. We look for two critical solutions (18) corresponding to different signs of . After substitution into the lens equation, we separate the terms containing integer and half-integer powers of , For example, where are analytical in . This yields a system of four independent equations for the variables ; this system can be reduced to a form convenient for iteration procedure [5]. It is important to note that at every iteration step we obtain an approximate solution in the form of finite order polynomials of , as it was stated at the end of Section 2.2.3. The solution of this system up to the terms is [5]

###### 2.4.2. Total Amplification of Critical Images of a Point Source Near the Fold

The solutions of the lens equation are then used to derive the Jacobians of the lens mapping (for both images near the critical curve). The value of yields the amplification of individual images. As we pointed out above, we need the total amplification of two critical images (the sum of two amplifications of the separate critical images). In the second-order approximation (using the expansion up to the terms ), this is where the constants , are expressed via the Taylor expansion coefficients from (1), and is the Heaviside step function. Note that is the caustic curvature at the origin which enters explicitly into the amplification formula. Parameters and are independent; explicit formulae for them may be found in [5, 6]. However, this is not needed when we use (32) for fitting the observational data, because these constants are whatever considered as free fitting parameters.

Formula (32) yields an effective approximation for the point source amplification near the coordinate origin provided that , and is not too small (see the term containing ). For a fixed source position, this can be satisfied always by an appropriate choice of the coordinate origin, so that the source will be situated almost on a normal to the tangent to the caustic.

If the source is on the caustic tangent or in the region between the caustic and the tangent, then formula (32) does not represent a good approximation to the point source amplification. Nevertheless, in case of an extended source, we will show that result (32) can be used to obtain approximations to the amplification of this source even as it intersects the caustic. However, to do this, we need to redefine correctly the convolution of (32) with a brightness distribution.

###### 2.4.3. Dark Matter on the Line of Sight

It is well known that the nonbaryonic dark matter (DM) dominates in galactic masses, though exact small-scale distribution of DM is a subject of studies. Simulations show that a complicated subhalo structure is possible [24]. In this connection, it is interesting to study how a *variable* DM density can affect the HAE characteristics.

Formulae (27), (32) have been obtained under an assumption that there is no continuous matter on the line of sight . It is easy to take into account an effect of a nonzero constant density by means of a rescaling of variables and change of coefficients in (2) [5].

For , the lens potential in (1) satisfies equation . The Taylor expansion analogous to (2) after substitution (3) takes on the form in the second approximation: Here, is the matter density at the origin, and the coefficients of expansion in (34) are expressed by means of derivatives of the lens potential . If is constant, then , , . Thus, system (34) contains four additional parameters as compared to the previously studied case . The solutions of system (34) have been studied in [18]; they are analogous to (27)–(29) and they have the same analytic structure. The formula for the amplification also preserves its form (32); the difference is only due to a change of explicit expressions for and in terms of new Taylor coefficients, which anyway cannot be determined from observations for realistic models of mass distributions. We, however, must make a reservation that here we suppose that the Taylor expansion in the lens equation is possible and it is effective (which presupposes that the continuous matter is sufficiently smooth); therefore, very inhomogeneous case of small objects with size is not involved in our consideration.

##### 2.5. Resume of Section 2

The main results of this section deal with solutions of the lens equation in the caustic region, namely, near the folds and cusps. We propose a general procedure of the power reduction that provides an approximate solution with a prescribed accuracy. Most attention is paid to the case of the fold caustic. We outlined two methods that enable us to obtain the critical solutions of the gravitational lens equation near a fold with any desired accuracy.

In order to obtain nontrivial corrections to near the fold obtained in the linear caustic approximation, the higher orders of the expansion of the lens equation must be taken into account as compared to works [7, 19]. The modified formula for contains 3 extra parameters in addition to those of the linear caustic approximation. We point out that any presence of a continuous (cf. dark) matter on the on the line of sight (with the same typical scales of Taylor expansion) does not change analytical structure and the number of fitting parameters in the formula for .

#### 3. Amplification of Extended Sources

##### 3.1. Preliminary Comments

In this section, we proceed to applications of the previous results to some extended source models typically used in fittings of the observed light curves. Interest to HAE in extragalactic GLSs is due to possibilities to study a brightness distributions in sources. This is especially interesting in connection with investigation of central regions of distant quasars. The idea, first proposed by Grieger et al. [25], uses an approximate formula of flux amplification during HAE, which contains a few fitting parameters. This makes possible some estimations of certain GLS characteristics, in particular, the source size [25]. For example, in case of the well-known GLS Q2237+0305 (Einstein Cross), several HAEs was observed [26–28] and the estimates of the source size have been obtained within different source models [16, 29–34].

A possibility do distinguish different source models is widely discussed elsewhere. Typical problems arising in determination of the source brightness distribution are as follows.

First, it is impossible to get complete information about the brightness distribution from the light curve observations. Information from separate HAE only, without making recourse to the whole light curve, and so forth, is still more limited. Observations provide only a one-dimensional luminosity profile of the source (integrated along the caustic); then, without using additional information, we cannot determine even the source size because we do not know the value and direction of the source velocity, ellipticity, and orientation with respect to the caustic.

Second, even for a circularly symmetric source and known normal velocity with respect to the fold caustic, determination of the luminosity profile is a kind of ill-posed mathematical problems: small variations of input data may lead to considerable changes of the solution. A standard way to mitigate this difficulty involves additional restrictions and/or using some simple “fiducial” models for brightness distribution. Some of these models are considered below in this paper. However, it should be clearly understood that the real picture of the central quasar region is more complicated that the simple brightness distributions of the following section; these models can be considered rather as reference ones. On the other hand, in view of the present-day accuracy of observations, it is difficult to distinguish even these simple source models on the basis of observational data. For example, the authors of [35] argue that the accretion disk can be modelled with any brightness profile (Gaussian, uniform, etc.), and this model will agree with the available data provided that an appropriate source size is chosen. On the other hand, a number of authors [4, 16, 33–39] discussed delicate questions concerning determination a fine quasar structure from HAE. For example, the authors of [36] wrote that the GLITP data [28] on Q2237+0305 admit only accretion disc models (see also [4, 39]). Obviously, the presence of an accretion disk in a central region of quasar is beyond any doubts, as well as the fact that the real appearance of the quasar core can be quite different from our simplified models. However, is it possible to prove the existence of the accretion disk in a concrete GLS *a posteriori*? This is an open question.

We note that since the work by Kochanek [37], followed by a number of authors [4, 35, 38–42], statistical methods dealing with the complete light curves of the GLS images have been developed. This approach is very attractive because it allows to take into account the whole aggregate of observational data on image brightness variations yielding estimates of the microlens masses and source model parameters. However, this treatment involves a large number of realizations of the microlensing field and requires a considerable computer time for such simulation. On the other hand, we must remember that (i) the source structure manifests itself only in the HAEs; far from the caustics the source looks like the point one and all the information about its structure is being lost; (ii) in reality, we have only one light curve, not a statistical ensemble, so probabilistic estimates of GLS parameters obtained in statistical simulations can differ from the real values in a concrete GLS. If we restrict ourselves to the HAE neighborhood, then part of the information is lost, but instead we use the most general model of the microlensing field described by a small set of coefficients in the lens mapping.

##### 3.2. From a Point Source to Extended One

After these reservations, we turn to the amplification of some extended source models.

Let be a surface brightness distribution of an extended source. If the source center is located at the point in the source plane, then the total microlensed flux from the source is where ; is the lens mapping. The result of using (35) obviously is equivalent to the result of the well-known ray-tracing method [1] (when the pixel sizes tend to zero).

An equivalent representation of this formula is where the point source amplification is the sum of amplifications of all the images.

Near a caustic, one can approximate , where is an amplification of all noncritical images that is supposed to be constant during HAE, and is the amplification of the critical images. Due to a relative motion of the lensing galaxy and the source (quasar), the flux is a function of time representing the lightcurve of some quasar image in GLS.

Formula (32) for contains the nonintegrable term~. Therefore, the question arises of how formula (32) can be used in situation when the extended source intersects a caustic and some part of the source is in the zone between the tangent and the caustic. In view of Sections 2.2.3 and 2.4, it is evident that the mentioned term is a result of the asymptotic expansion of the root in the approximate solution (18). Direct usage of the solution in the form (18) for calculation of the Jacobians of the lens mapping and then for the derivation of amplifications does not lead to any divergences and any nonintegrable terms in do not arise without using this expansion. Nevertheless, it is convenient to have a representation of in the form of expansion in powers of small parameter. Such an expansion can be fulfilled correctly after substitution of into integral (36). On this way, starting from form (18), it is easy to show that to define correctly, one must replace the term in (32) by the distribution (generalized function) [43]. We recall that the distribution of the variable is defined by the expression: for any test function .

After this redefinition, we have This formula can be used to correctly derive an approximate amplification of a sufficiently smooth extended source including the case where the source crosses the caustic.

##### 3.3. The Extended Source Models

###### 3.3.1. Gaussian and Power-Law Models

Below, we list most simple and commonly used brightness distributions of a source in GLS; without loss of generality, they are chosen to be normalized to 1: To compare different models of the brightness distribution, we have to use the same parameter that characterizes the size of an object. The r.m.s. size is often used: However, for slowly decreasing brightness profile (e.g., ), the r.m.s. size diverges. In case of the circularly symmetric sources the half-brightness radius is also widely used; it is defined by the relation:

In case of Gaussian source model, where stands for a size parameter; , .

Limb-darkening model (see, e.g., [44]): where and . Here, we assume . The half-brightness radius is . For fixed and , the brightness distribution (43) tends to the Gaussian one.

The models (43) and (42) describe a class of compact sources with fast brightness decrease. On the contrary, the power-law models [16, 33] describe a slow decrease at large : where is the power index, and is related to the r.m.s. radius as . The model (45) may be considered as an alternative to (43). The half-brightness radius of the source for this model is . Like (43), for fixed and , the brightness distribution (45) tends to the Gaussian one. For small , we have a “long-range” distribution; diverges for .

Linear combinations of different distributions (42), (43), (45) with different parameters yield rather a wide class of symmetric source models to fit any kind of data. On the other hand, (42) may be considered as a fiducial model to determine some parameters such as the source size and (43), (45) are useful in case when in addition we are interested in investigation of a brightness behavior at large . More physical models are considered in the next subsection.

###### 3.3.2. Accretion Disk Models

The accretion disk (AD) of Shakura-Sunyaev [45] has a more complicated profile. This model gives the energy density of the radiation from accretion disk around a nonrotating black hole as a function of radius whence the (normalized) brightness distribution is here being the radius of the inner edge of the accretion disk. For this AD model, the half-brightness radius is , . This formula describes a total brightness integrated over all radiation frequencies. The maximum brightness is at .

For a blackbody radiation, the temperature scales as , whence the specific intensity as a function of radius ~ for wavelength (AD1): here is a normalization factor. For the maximum disk temperature (at ), the peak value of intensity as a function of corresponds to and . Further, we adopt this value of for which , . Below in Section 3.5 we use these parameters with , so . Intensity (47) can be easily rewritten for any wavelength taking into account that .

Though intensities (46), (47) are quite different, their light curves in HAE can look very similar for an appropriate choice of parameters. The main feature that distinguishes the accretion disk models from the Gaussian one is the concavity of the light curve owing to the dark “hole” in the accretion disk center. We will return to this question in Section 3.5.

##### 3.4. Amplifications Including the Postlinear Corrections

###### 3.4.1. Amplification for the Gaussian Source

Formula (38) has been used [5] to derive the amplification of a Gaussian source with the brightness distribution (42), the limb-darkening source, and the power-law source (see the next subsections).

Further, we use the dimensionless coordinates of the source center and the functions: These functions can be expressed in terms of the confluent hypergeometric function or the parabolic cylinder function :

The substitution of (38) and (42) in (36) yields

Here, Note that the main term of (50) which corresponds to the linear caustic approximation was first obtained in the paper [20].

###### 3.4.2. Amplification for the Limb-Darkening Source

Analogous considerations allow us to obtain formulas for the amplification of extended sources for the limb-darkening (43) and power-law (45) brightness profiles; the results are represented analytically in terms of hypergeometric function [46].

Denote . We have for and for .

Also, for , we define

Then in case of the model with limb darkening (43), the critical images disappear when the source lies on the outer side of the caustic (i.e., for ). The amplification due to critical images takes on the form:

###### 3.4.3. Amplification for a Power-Law Source

The result for the amplification involves integrals: for , being the Beta-function.

We extend this to having in mind the definition of , so that Now, the amplification due to critical images takes on the form: The zeroth approximation to this formula has been derived in [16].

##### 3.5. Simulations of Microlensing Light Curves with Different Source Models

###### 3.5.1. Equations of Microlensing and Parameters of Simulations

Though we expect that a source structure reveals itself just during HAEs, it is instructive to compare whole light curves with different source models; this is fulfilled below. Nevertheless, we point out that below we mainly deal with an estimate for (see (62) below), which is just a characteristic of HAE. We assume the lens equation:
which describes microlensing by *N* extended masses (clumps) in absence of any external field; the surface mass density corresponding to (60) is
where the surface mass density of one clump is ; here stand for the size and mass of nth clump. In case of point masses ; however, in simulations of the point mass microlensing, we have chosen some small positive values of so as not to have problems with singularities. In Section 4 we also consider the case of a finite-size extended masses.

Here, we present results of straightforward calculations of a microlensed flux according to (35) to obtain the light curves for different realizations of the point microlenses positions. The parameters of numerical integration are presented in Table 1 along with the parameters of the microlensing field. The microlense positions were chosen in a random way with uniform distribution over the field. The length of trajectory has been taken sufficiently long so as to provide the caustic crossings. However, the trajectory in the simulations was situated far from the field boundaries, and the size of the field was chosen large enough to avoid boundary effects.

Here we compare the light curves for Gaussian, PL, LD, AD, and AD1 source models with brightness distributions from Section 3.3 The simulations were performed for a set of 100 random realizations of microlensing field with optical depth in order to have possibility to compare with some of the results of papers [47–50]. All source models had the same half-brightness radius, the light curves for all source models have been calculated for the same microlensing fields. All the microlenses were static, the total number of microlenses was 1470. The typical magnification pattern is shown on Figure 1; here the speed of source is so we can identify the source position as a function of time.

The simulations were carried out with half-brightness radius for all models; the power-law index was for the “long range” PL model; note that AD model also corresponds to this class of the power-law asymptotic dependence (with ). For LD model, we have chosen .

To compare the light curves for different source models and , we used the relative value: where is the amplification for th model along the trajectory of source linear motion; is an average over a number of realizations. Because we expect that maximum difference of brightness on lightcurves for different source models takes place during HAE, so typically is a characteristic of HAE.

From the “light curves” on Figure 2 we observe a significant difference between the “compact” (LD and Gaussian) and the “long-range” models. The long-range character of the latter reveals itself even on considerable distances from the caustics, where we expect that brightness of all the sources must have the same behavior as that of a point source. The differences between these two groups of models are essentially larger than the differences within each group (e.g., between Gaussian and LD models). This conclusion is confirmed by results of the statistical considerations over 100 realizations shown in Table 2 for half-brightness radius as the example.

We note that the results of comparison may depend on a complexity of the caustics involved into our consideration. One can expect, for example, that in case of complicated caustic crossing (such as crossing of the fold caustics in a vicinity of cusp, intersection of dense aggregations of caustics, among others) the light curve difference for different source models may be more significant. Our previous simulations involved *all* possible realizations of the microlensing field that may invoke the complicated caustic crossings. However, in case of a concrete GLS, we deal with a *single* light curve. One may ask: is it relevant to apply the results of statistical simulations to this single light curve, where we may have a kind of an observational selection. Therefore, it is necessary to check how this “complexity” affects the average value of . In this connection, we considered some modification of our statistical procedure: we have chosen (by eye) such realizations rather simple fold caustic crossings involved. However, the numerical results for simple caustic crossing events, such as shown as the example on Figure 3, appeared to be nearly the same as that of the Table 2. Though we must note that here the statistics have been worse because of a smaller number of the “simple” caustic realizations. So this conclusion requires further verification.

###### 3.5.2. Gaussian Fittings of Accretion Disk and Limb Darkening Models

Table 2 concerns a comparison of different models *with the same *. However, in reality, we do not know what source we must fit and one must check whether we can replace one model with a different one with some other source parameters to get better fitting.

We have fitted the limb-darkening and accretion disk model light curves with that of the Gaussian source of different radii. The half-brightness radius of LD and AD models was fixed: (this choice was determined by the AD1 model, see comments after formula (47)) and that of the Gaussian source was varied to get better fitting; here for PL model and for LD model. One example of the relative light curve differences according to (62) between LD, AD1 and the best fitted Gaussian models is shown in Figure 4. As we see from this figure, usage of an optimal size of the Gaussian model can diminish as compared to the average estimate from Table 2 dealing with the same for all the models.

To test this statement statistically, we carried out simulations with 20 magnification patterns having the “simple” caustic crossing. We found that fitting of AD light curves by the Gaussian model (with different ) yields ; this is roughly the same as in Table 2. Fittings of LD model light curves by the Gaussian model (with different ) yielding (less then half of Table 2); in this case, an improvement is noticeable. Therefore, the Gaussian source cannot reproduce all the models, though for some compact models the fitting results are rather good. Then, in the caustic crossing events different models may be distinguished on the accuracy level of modern photometric observations.

##### 3.6. Resume of Section 3

Here the formulas of Section 2 are applied to the Gaussian, power-law and limb darkening models of an extended source. The analytical formulas for these source models are obtained. We note that a fitting of the light curve of GLS Q2237+0305C [5, 6] showed that some of these higher-order corrections can be statistically significant even at the present accuracy level. This means that when we are looking for some fine effects in HAE due to the source size, a consistent treatment must also take into account possible effects of the above corrections.

In case of the accretion disk models, we have not succeeded to find compact “workable” expressions for , so this case has been treated numerically to compare possible observational differences of various source models.

Now we remind the main assumptions and the class of models used, which is necessary for understanding the results of the Section 3.5. The simulations use the GLS parameters similar to that of Q2237+030 [47, 51]. We consider the equal mass microlenses; we do not consider any mass distributions and/or populations of small (planetary) masses. Next, we consider most simple circular symmetric source models without effects of ellipticity, and so forth. The results of Section 3.5 can be summed up as follows. There are special classes of models: (i) accretion disk models having the dark region in the center and (ii) source models with a slow decrease of the (long-range) brightness distribution at large distances from the center. From our results, we see that (i) and (ii) show in HAEs a behavior which is different from that of the Gaussian and the limb-darkening models on a relative accuracy level of. There is then a possibility to distinguish such models on account of observations during HAEs. Outside HAE the surface brightness profile typically has little effect on microlensing (see [35, 52]).

We note that our conclusions involve concrete source models with a special choice of parameters (source sizes, microlensing optical depth, etc.) and a considerable work is to be fulfilled to obtain more general results. Also, our results are of a statistical nature, so our average estimate of may be different from analogous one in a separate light curve observation.

Note that adding a sufficient number of microlenses having a small (e.g., planetary) masses (cf. [51, 53–55]) must produce an additional caustics that in its turn can give rise to fine features of light curves during HAE. One may expect that this can mimic features due to brightness distribution over the source. This question deserves further investigation.

#### 4. Astrometric Gravitational Microlensing

##### 4.1. Preliminary Comments

Besides photometric observations of microlensed sources, additional information about the source and lens structure may come from astrometric effects. Gravitational fields of masses cause the bending of light rays; this effect has passed through numerous experimental tests within the solar system. This bending leads to brightness amplification of distant sources in GLSs. We say about the astrometric gravitational microlensing when we deal with shifts and motions of the distant source images due to the gravitational fields of foreground stellar-mass objects (stars, black holes, dark matter substructures, etc., either in the Milky Way or in the other galaxies).

Possibilities to detect an image motion of the remote source caused by the gravitational field of the Milky Way stars have been discussed for a long time [56–60]. The problem gained a new interest in connection with investigations of massive compact halo objects [61–63], see also [64–71] and references therein; for the extended source and extended lens effects see [11, 72–74]. The extragalactic topics are less studied [12, 71, 75, 76]. Nevertheless, the microarcsecond accuracy level appears to be sufficient to resolve the shape of a source image trajectory in extragalactic GLS in microlensing events. The corresponding image shifts can achieve the level ofas.

Astrometric measurements can provide a valuable information about mass distributions of microlenses, which can complement the photometric observations of GLS images [77, 78]. Observations of trajectories of the remote source images during the microlensing events, in addition to the photometric data, would give us a possibility to estimate the foreground field parameters and the optical depth of the continuous and discrete matter. This may be helpful in connection with the dark matter (DM) problem; in particular, this can be used to test the existence of the putative DM subhalos [79]. Though such positional effects are still beyond the accuracy of observations, the astrometric microlensing attracts a considerable attention during the last decade in connection with perspectives of spaceborn microarcsecond astrometry (the forthcoming Gaia mission, space-based VLBI missions, optical Space Interferometry Mission). It should be pointed out that, typically, the accuracy of astrometric positioning is essentially higher than the image resolution. The other positive feature of astrometric microlensing events is the slow dependence on the impact parameter and, hence, there is a larger probability of observing such events, once a required accuracy will be achieved.

There is a number of candidates for astrometric microlensing in our galaxy that can be detected in the near future [80–82]. On the other hand, Klioner [83] points out that unpredictable microlensing noise can spoil the determination of positions and proper motions of the objects resulting from future astrometric missions on a submicroarcsecond level. The gravitational image motion (GIM) due to microlensing effects may be comparable to proper motion of quasars [8, 9, 84–86].

##### 4.2. Extended Source and Extended Lens

Further, dealing with trajectories of the microlensed source images, we mean the trajectory of the source image centroid (SIC), that is, a weighted average over all elements of all the source images. The position of SIC can be defined as , is the lens mapping, is the center of the source in the source plane; is given by formula (35).

Having in mind that most of the astrometric microlensing effects are still outside observational possibilities, for astrometric purposes, it is quite reasonable to confine ourselves to the most simple source brightness profiles described in Section 3.3 We only note that for the microlensing of a star the limb darkening model is most appropriate, and structures of real extragalactic sources (e.g., quasars) are more complicated.

In this subsection, we consider microlensing by one mass in presence of an external shear. This can be the case of microlensing of a remote source by a small mass (planet) in presence of a background gravitational field of a star or microlensing by a star in the field of a putative DM clump.

First, we consider a weak microlensing system, that is, the linear distance from a microlensing mass to the line of sight is much larger than its Einstein radius. In this case, if a remote point source is microlensed by a foreground moving point mass, the trajectory of the source image is approximately the circle in the rest frame of the source [9, 66]. In case of strong lensing HAE (linear distance from the microlens to the line of sight is of the order of the Einstein radius), this circle will be distorted. In any case, the point source image is “repulsed” from the lensing point mass.

The trajectory qualitatively changes when the projection of the point mass onto the source plane crosses the extended source. In this case, the SIC motion trajectory typically has a three-leaf form [11, 72]. This has been first shown in case of a source with uniform surface brightness [72] and in case of the Gaussian source [11]. The reason for such behavior is that when the lens is projected onto the internal part of the source, the “repulsion” of SIC is changed by its “attraction” to the lens due to larger amplification near it. For sufficiently small values of impact parameter of the lens with respect to the extended source center, the SIC trajectory bends and passes through the source center forming a self-intersection.

The external field effects, which are essential in extragalactic GLS, distort considerably the trajectory of SIC [12, 13, 87]. Below, we present some simple examples of microlensing of an extended Gaussian source in presence of an external shear (Chang-Refsdal lens [88]). The configuration of the figures is as follows: the left hand panels show the SIC trajectories; the right-hand ones show corresponding dependencies of the amplification as functions of time. The SIC trajectories are shown in the rest frame of the source; in the right upper corner the source center trajectory is shown (schematically) with respect to the caustic. We show in figures the impact parameter () and the shear ().

Figure 5 shows the differences between the point-like source and the extended Gaussian one (42) in presence of an external shear. To compare astrometrical and photometrical effects for different source models, we represent in (Figures 6 and 7) the SIC trajectories and light curves of the Gaussian source and accretion disk (47) in band *V*, both with the same half-brightness radius of . The last one, Figure 8 shows the SIC trajectories in case of finite-size microlensing masses. We note that the SIC trajectory depends on the source trajectory and upon the source model more distinctly than the light curve.

The simplest situation is when the source track passes far from the caustic. In this case, the SIC trajectory is oval-like and light curve is continuous. However, when the source center track is close to the caustic, even if it does not cross it, the trajectory becomes much more complicated. In general, it can be self-intersected several times. When point source crosses the caustic, its SIC suffers a jump (Figure 5). This is explained by the emergence or by disappearance of pairs of images of infinite brightness, leading to sudden jumps of the weighted average of positions of the source elements. Analytical expressions for such jumps can be found in [13]. For an extended source, corresponding discontinuities of the trajectory are smoothed out but may be quite noticeable. When the small source passes near the cusp or crosses it, the total images brightness increases dramatically, but no jump of image motion occurs because of the existence of one more image of infinite (in case of a point source) brightness which does not disappear. This distinguishes the high amplification events corresponding to the fold and cusp caustic crossings.

The effect of high optical depth of continuous matter is demonstrated in Figure 7 to show possible observational signals of continuous DM near the line of the sight.

In papers [73, 89], astrometric and photometric effects of GLS with an extended microlens were considered. The microlensing effect caused by a compact body (e.g., a star) must be distinguished from the eclipse-like effect in case of a presence of quite dense but essentially continuous background (cf. “subhalos”) of DM. Such “subhalos” are supposed to be extended clumps of DM. The size of these clumps predicted by model is larger than typical Einstein radius (see, e.g., Diemand et al. [90]). Such an object can manifest itself in microlensing process acting together with another compact object. In this case, smooth-shaped or even quasielliptical caustics can exist instead of, or together with, cusped ones (Figure 7). When the source crosses the fold of the smooth caustic, two images of it disappear. As a result, after some peak of the image brightness, there may be swift diminution. Thus DM can cause the effect like “eclipse,” demonstrated on Figure 7. In an opposite case, when a subhalo have finite size smaller than the Einstein radius, finite-size lens models can be applied to describe the situation. As it was shown in [89], the finite-size opaque lens can cause “ears” on the trajectory when one of the images is obscured by the lens, and the convergence causes the trajectory disruption when the source crosses caustic. In the light curve, a finite lens size manifests itself in the magnification jumps when the image enters in the lens shadow. In case of an overcritical convergence, the magnification of the point source images becomes infinite when the source crosses the caustic, and the brightness decays swiftly when the source is inside the caustic.

We show here also some results obtained within the other “toy-example” for the microlensing by a subhalo of the size modelled by lens equation (60) with the only extended lens (). The results obtained within this model are shown on the Figure 9. As we can see here, the finite size of the subhalo causes an apparition of the additional quasielliptical caustic. Meanwhile, the SIC trajectories do not differ principally from the point-lens situation. However, here, we can see the dip on the lightcurve similar to that obtained in Zakharov and Sazhin lens model [89].

The effect of the finite-size lens manifest itself in smoother way if the source is extended (see Figure 8). In this case, it changes the SIC trajectory and the light curve has no jumps or gaps. The most prominent effect on the light curve is in the recession of the dip between two caustic crossings (i.e., when source is inside the caustic).

##### 4.3. Statistics of Weak Astrometric Microlensing

###### 4.3.1. Correlation Tensor of Image Shifts

In a general case of a complicated microlensing system, we need some statistical measure of image motion; possible characteristic is proposed below.

Let describes a deviation of the point source image from unperturbed position. In case of one microlens, this value typically varies as , where is a distance of the microlens from the line of sight. Because of this, we expect that in case of uniform distribution of many microlenses usual correlation function and diverge, and so we prefer to deal with a correlation tensor like that used in the turbulence theory [91]: This tensor is symmetric and it has a diagonal form in a system where the source moves along one of coordinate axes, for example, along -axis. So a determination of from observations may help to constrain a direction of the source motion.

We will illustrate calculation of in case of weak astrometric microlensing by *N* clumps having a uniform distribution over the lens plane. The system is described by the lens equation (60). We suppose that all the clumps are independent, and has a uniform distribution over the lens plane. To be precise, we assume that the spatial distribution function for every is
where is the area of domain . Further considerations deal with the limit of that will be blown up to cover the whole plane assuming
where we take into account that the surface per one clump may be different for clumps of different mass and size. Density goes to zero sufficiently fast as , so the clumps have a finite size and integral is convergent. To save the space, we do not consider a distribution of parameters assuming that these parameters are the same. This means that in the general case the results below must be averaged over possible .

According to (60), in case of weak microlensing the shift of the point source image is

For a straight line motion of a source, we put . Then, we have a diagonal form of with

Calculation of the integrals yields where

###### 4.3.2. Statistical Effect of Gravitational Dragging of a Distant Source Image by Point Masses

In this section, we consider statistical effects in the motion of the GLS images due to gravitational field of a large number of moving stars. Here, we consider a rarefied collection of lensing masses, when the typical distances between the source and the microlens are much larger than the corresponding Einstein radius. For brevity, we call this “weak” microlensing; this must be distinguished from the well-known extragalactic weak lensing, that is, image distortions caused by the cosmic shear. Unlike the previous sections, here we rule out strong microlensing events that are accompanied by a considerable brightness amplification.

In some situations it will be difficult to separate GIM from real proper motions, and this is important for accuracy of the fundamental reference frame based on extragalactic sources [83, 85, 86]. The gravitational field of point masses (which is essentially inhomogeneous) induces a stochastic motion of the source image in the reference frame of the source. On the other hand, smoothly distributed masses also induce a stochastic motion, but typical amplitude of the effect will be weaker. Indeed, we cannot take into account positions of all stars [83]; we only can work with their probability distribution leading to distribution of image motions of a source [8, 9]. The situation could be further complicated in presence of invisible DM clumps, if they exist on small scales.

In case of the point masses, it is possible to estimate the probability distribution of image motions (i.e., angular velocities) of the distant source image microlensed by collection of a randomly moving point masses [8]: where is a characteristic value depending upon the transverse velocity of point masses, their value, and their spatial density. It should be noted that the distribution (71) is obtained under the approximation of the weak microlensing [8], so that HAEs must be excluded.

The distribution (71) has a slowly decreasing tail yielding an infinite estimate of velocity squared . This tail is a formal consequence of pointness of the microlenses that makes possible very high GIMs. In fact this case is not covered by formula (71), which is valid for weak microlensing only. In case of smooth DM concentrations such velocities would be absent.

If a bulk motion of microlensing masses is present, the average requires a special consideration. It has been pointed out in [8, 9] that stochastic GIM of an extragalactic source induced by stellar motions is accompanied with a systematic component (not a statistical average) which depends upon bulk velocity of microlensing stars, including stars that are far away from the line of sight. The typical velocity of this component remains nonzero even for a homogeneous distribution of microlensing masses which move in certain direction (all masses have the same collective velocity component). This needs some explanation because the average position of the unmoving source image (averaged over a very long time) must be unchanged and, evidently, the statistical average of image velocity must be zero. However, when dealing with the extragalactic reference frame, in reality, we do not deal with a very large observational time and/or with a sufficiently large number of sources to obtain a statistical average. In this case, because the strong microlensing events yielding the fast image motions are very rare, they in practice are not observable unless the number of simultaneously observed sources is of the order of millions or more. This is well known from theory and observations of MACHOs. Then the statistical average over set of all events is not a good characteristic of GIM; the most probable value of GIM (which is different from zero) would be a better measure.

Otherwise, we can estimate probability for all necessary velocity intervals. Detailed approach on this way, however, involves consideration of strong lensing events with different extended source models. However, to have an order-of-magnitude estimate, we can avoid calculations of probability distributions, if we confine ourselves to the domain of weak microlensing events (“*W*”, for brevity) [9, 10]. Typical velocity which is a characteristic of such events is a statistical average over the domain *W*, with HAE being excluded. Note that these rare strong lensing events (*S*) essentially contribute to the usual statistical average velocity (i.e., average over all the space of events ). If we exclude *S*-events (this is easy to do in real observations because of strong brightness amplification of a source), then we deal with “truncated” average velocity over *W* events. This value is nonzero even in case of a homogeneous distribution of microlenses having a collective velocity. It is important to note that this value practically does not depend upon the exact definition of *S*-domain (i.e., upon a value of minimal impact parameter of microlensing masses with respect to the line of sight, which determines events “*S*”).

Note that the “truncated” average would better correspond to realistic observations when the typical velocity is estimated as an arithmetic mean over all reference sources used in extragalactic International Celestial Reference Frame (ICRF). In a near future, the number of these sources will not exceed several hundreds. Therefore, in fact, we deal with most probable dispositions of lensing masses which typically do not include the *S*-events: it is well known that the probability of the such events is of the order of per year in the Milky Way galaxy. Then, in ICRF considerations, it is reasonable to exclude such rare events that would hardly occur during the century.

In case of microlensing of distant sources by Milky Way stars, we may neglect cosmological curvature. Evidently, this does not contradict to assumption that the radiation source is at the infinity. Let unperturbed light ray move from the infinitely distant radiation source in negative direction of -axis of Cartesian coordinates {}, the observer being at the origin. Let position of microlensing point mass be , where is a two-dimensional vector in the transverse plane; that is, is the impact distance of the unperturbed ray with respect to the mass. In the case of weak microlensing, , . In this case, where is a two-dimensional vector describing the source image angular shift [9, 85, 92].

Further stands for velocity of the point microlens, is the velocity component parallel to the line of sight, and represents the transversal components. In virtue of (72), this leads to the source GIM that equals to (in radians per unit of time), where Here, is a function of the microlens position and its velocity .

Because of the smallness of the effect, the action of all Milky Way stars will be taken into account linearly. This will be performed by integrating (73). However, this requires justification. In fact, we are dealing with extended sources; this is not taken into account by (73). As we will see below, this may be important in our problem.

Thus, we return to (73), which describes contribution of a single star. We assume that a star at the point has velocity . Then, we consider sum over all stars in the Milky Way, which must further be averaged with the mass density . This enables us to pass to integration to yield an average GIM in domain : where we suppose that vanishes outside a bounded domain. Also, we suppose that is sufficiently smooth so as we can use approximation of weak microlensing inside the domain where .

Taking into account the explicit form (73), it is easy to see that in a general case considerable contribution may be due to stars at large impact distances from the line of sight. Also, one can show that the singularities in (74) for small are integrable and the integral is convergent. This allows us to avoid the question about exact value of the lower limit of in the definition of domain .

It is interesting to compare GIM (74) with an average over all events in case of a general continuous mass distribution . This distribution must satisfy the continuity equation:

In this case, (72) should be integrated with the density :

Note that this formula for can be considered either as the average over all events or simply as a shift in case of mass distribution .

The shift of remote source image in case of continuous mass distribution has been treated in [93]. These shifts appear to be of the order of radians, but this value cannot be observed from the solar system. The observable (at least in principle) value is a change of (76) with time: Formula (77) is, in fact, an average image velocity over all events: One might think that this must be the same as (74), but it is not true. The reason is that operations of averaging and differentiation here do not commute. To compare (74) with (77), one must use the continuity equation (75) and avoid singular points at in [10].

We have , where and are two-dimensional vectors, and the components of the first term (): are the same as the components of (74): . The components of the second term can be written as After some calculation on account of the Ostrogradsky-Gauss theorem, we have [10] where stands for the density of microlensing masses along the line of sight, is the transverse velocity of microlenses near the line of sight. Evidently, (81), differs from zero. For example, if an observer at rest is situated in the center of Milky Way (stationary mass distribution), we have (see (77)), but GIM does not equal to zero.

For an observer in the solar system this effect yields a nonzero GIM value leading to apparent rotation of extragalactic reference frame in the direction of the Milky Way rotation [10]. For sources near the Galactic plane, the effect amounts about arcseconds per year [10]. The effect will be even smaller for sources with higher Galactic latitudes. Therefore, in principle (if the corresponding accuracy be achieved!), the effect can be observed from the Solar system.

##### 4.4. Resume of Section 4

In this Section, we review main points of astrometric microlensing. We demonstrate effects due to a finite size of the source and extended lens as well; a special attention is paid to effects that may serve as signals of dark matter. Several examples of microlensing by one mass in presence of an external shear are given. This can be the case of microlensing of a remote source by a small mass (planet) in presence of a background gravitational field of a star or microlensing by a star in the field of a putative DM clump.

Then, we proceed to some statistical subjects of astrometric microlensing. We introduced the correlation tensor as a statistical characteristic of a stochastic image motion and applied this to a “toy” model of weak microlensing by dark matter clumps.

We have shown that gravitational field of foreground stars induces an image motion, which differs for a continuous and discrete mass density distribution. In the latter case, image performs random walks, as distinct from a regular motion in case of continuous matter, plus some additional motion due to the collective motion of the microlensing masses. This effect leads to apparent rotation of extragalactic reference frame in the direction of the Milky Way rotation. In principle, the effect can be observed, however, it is far beyond modern possibilities.

#### Acknowledgments

This work has been supported in part by “Cosmomicrophysics” Program of National Academy of Sciences of Ukraine. V. I. Zhdanov and V. M. Sliusar acknowledge the partial support of Swiss National Science Foundation within SCOPES Project IZ73Z0_128040.

#### References

- P. Schneider, J. Ehlers, and E. E. Falko,
*Gravitational Lenses*, Springer, New York, NY, USA, 1992. - A. F. Zakharov,
*Gravitational Lenses and Microlenses*, Yanus-K, Moscow, Russia, 1997. - A. O. Petters, H. Levine, and J. Wambsganss,
*Singularity Theory and Gravitational Lensing*, Birkhäuser, Boston, Mass, USA, 2001. - R. Gil-Merino, J. González-Cadelo, L. J. Goicoechea, V. N. Shalyapin, and G. F. Lewis, “Is there a caustic crossing in the lensed quasar Q2237$+$
0305 observational data record?”
*Monthly Notices of the Royal Astronomical Society*, vol. 371, no. 3, pp. 1478–1482, 2006. View at Publisher · View at Google Scholar · View at Scopus - A. N. Alexandrov and V. I. Zhdanov, “Asymptotic expansions and amplification of a gravitational lens near a fold caustic,”
*Monthly Notices of the Royal Astronomical Society*, vol. 417, pp. 541–554, 2011. View at Google Scholar - A. N. Alexandrov, V. I. Zhdanov, and E. V. Fedorova, “Asymptotic formulas for the magnification of a gravitational lens system near a fold caustic,”
*Astronomy Letters*, vol. 36, no. 5, pp. 329–337, 2010. View at Publisher · View at Google Scholar · View at Scopus - C. R. Keeton, B. S. Gaudi, and A. O. Petters, “Identifying lenses with small-scale structure. II. Fold lenses,”
*Astrophysical Journal*, vol. 635, no. 1, pp. 35–59, 2005. View at Publisher · View at Google Scholar · View at Scopus - V. I. Zhdanov and V. V. Zhdanova, “Analytical relations for time-dependent statistical microlensing,”
*Astronomy & Astrophysics*, vol. 299, pp. 321–325, 1995. View at Google Scholar - V. I. Zhdanov, “The general relativistic potential of astrometric studies at microarcsecond level,” in
*Astronomical and Astrophysical Objectives of Sub-Milliarcsecond Optical Astrometry*, E. Hog and P. K. Seidelmann, Eds., pp. 295–300, Kluwer, Dordrecht, The Netherlands, 1995. View at Google Scholar - V. I. Zhdanov, E. V. Fedorova, and A. N. Alexandrov, “Gravitational dragging of distant source images caused by Galaxy stars,”
*Kinematika i Fizika Nebesnykh Tel*, vol. 20, pp. 422–429, 2004 (Russian). View at Google Scholar - V. I. Zhdanov and S. A. Salata, “Motion of the image of a distant object microlensed by stars in a foreground galaxy,”
*Kinematics and Physics of Celestial Bodies*, vol. 14, pp. 156–161, 1998. View at Google Scholar - V. I. Zhdanov, S. A. Salata, and E. V. Fedorova, “Background-field effects in astrometric microlensing,”
*Astronomy Letters*, vol. 27, no. 9, pp. 562–567, 2001. View at Publisher · View at Google Scholar · View at Scopus - E. Fedorova, V. I. Zhdanov, and A. N. Alexandrov, “Motion of source image in Chang-Refsdal lens,”
*Journal of Physical Studies*, vol. 4, pp. 465–468, 2002. View at Google Scholar - A. N. Alexandrov, V. M. Sliusar, and V. I. Zhdanov, “Caustic crossing events and source models in gravitational lens systems,”
*Ukrainian Journal of Physics*, vol. 56, no. 4, pp. 389–400, 2011. View at Google Scholar · View at Scopus - C. J. Fluke and R. L. Webster, “Investigating the geometry of quasars with microlensing,”
*Monthly Notices of the Royal Astronomical Society*, vol. 302, no. 1, pp. 68–74, 1999. View at Google Scholar · View at Scopus - V. N. Shalyapin, “Caustic crossing in the gravitational lens Q2237$+$0305,”
*Astronomy Letters*, vol. 27, no. 3, pp. 150–155, 2001. View at Publisher · View at Google Scholar · View at Scopus - O. Pejcha and D. Heyrovský, “Extended-source effect and chromaticity in two-point-mass microlensing,”
*Astrophysical Journal Letters*, vol. 690, no. 2, pp. 1772–1796, 2009. View at Publisher · View at Google Scholar · View at Scopus - A. N. Alexandrov, S. M. Koval, and V. I. Zhdanov, “Asymptotic relations for high magnification events in presence of the dark matter,”
*Visnyk Kyivskogo Universytetu, Astronomiya*. In press. - A. N. Alexandrov, V. I. Zhdanov, and E. V. Fedorova, “Analytical relations for gravitational lens mapping in the vicinity of a critical curve,”
*Visnyk Kyivskogo Universytetu, Astronomiya*, vol. 39-40, pp. 52–59, 2003 (Ukrainian). View at Google Scholar - P. Schneider and A. Weiss, “A gravitational lens origin for AGN-variability? Consequences of micro-lensing,”
*Astronomy & Astrophysics*, vol. 171, pp. 49–65, 1987. View at Google Scholar - P. Schneider and A. Weiss, “The gravitational lens equation near cusps,”
*Astronomy & Astrophysics*, vol. 260, pp. 1–13, 1992. View at Google Scholar - B. S. Gaudi and A. O. Petters, “Gravitational microlensing near caustics. II. Cusps,”
*Astrophysical Journal Letters*, vol. 580, no. 1, pp. 468–489, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. B. Congdon, C. R. Keeton, and C. E. Nordgren, “Analytic relations for magnifications and time delays in gravitational lenses with fold and cusp configurations,”
*Monthly Notices of the Royal Astronomical Society*, vol. 389, no. 1, pp. 398–406, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. Diemand, M. Kuhlen, P. Madau et al., “Clumps and streams in the local dark matter distribution,”
*Nature*, vol. 454, no. 7205, pp. 735–738, 2008. View at Publisher · View at Google Scholar · View at Scopus - B. Grieger, R. Kayser, and S. Refsdal, “Gravitational micro-lensing as a clue to quasar structure,”
*Astronomy & Astrophysics*, vol. 194, pp. 54–64, 1988. View at Google Scholar - P. R. Woźniak, C. Alard, A. Udalski et al., “The optical gravitational lensing experiment monitoring of QSO 2237$+$0305,”
*Astrophysical Journal*, vol. 529, no. 1, pp. 88–92, 2000. View at Google Scholar · View at Scopus - A. Udalski, M. K. Szymański, M. Kubiak et al., “The optical gravitational lensing experiment. OGLE-III long term monitoring of the gravitational lens QSO 2237$+$0305,”
*Acta Astronomica*, vol. 56, no. 4, pp. 293–305, 2006. View at Google Scholar · View at Scopus - D. Alcalde, E. Mediavilla, O. Moreau et al., “QSO 2237$+$
0305 VR light curves from gravitational lenses international time project optical monitoring,”
*Astrophysical Journal Letters*, vol. 572, no. 2, pp. 729–734, 2002. View at Publisher · View at Google Scholar · View at Scopus - J. S. B. Wyithe, R. L. Webster, and E. L. Turner, “A measurement of the transverse velocity of Q2237$+$0305,”
*Monthly Notices of the Royal Astronomical Society*, vol. 309, no. 1, pp. 261–272, 1999. View at Google Scholar · View at Scopus - J. S. B. Wyithe, R. L. Webster, and E. L. Turner, “Interpretation of the OGLE Q2237$+$0305 microlensing light curve (1997–1999),”
*Monthly Notices of the Royal Astronomical Society*, vol. 318, no. 4, pp. 1120–1130, 2000. View at Google Scholar · View at Scopus - J. S. B. Wyithe, R. L. Webster, E. L. Turner, and D. J. Mortlock, “A gravitational microlensing determination of continuum source size in Q2237$+$0305,”
*Monthly Notices of the Royal Astronomical Society*, vol. 315, no. 1, pp. 62–68, 2000. View at Google Scholar · View at Scopus - A. Yonehara, “Evidence for a source size of less than 2000 AU in Quasar 2237$+$0305,”
*Astrophysical Journal Letters*, vol. 548, no. 2, pp. L127–L130, 2001. View at Publisher · View at Google Scholar · View at Scopus - V. N. Shalyapin, L. J. Goicoechea, D. Alcalde, E. Mediavilla, J. A. Muñoz, and R. Gil-Merino, “The nature and size of the optical continuum source in QSO 2237$+$0305,”
*Astrophysical Journal Letters*, vol. 579, no. 1, pp. 127–135, 2002. View at Publisher · View at Google Scholar · View at Scopus - M. B. Bogdanov and A. M. Cherepashchuk, “Reconstruction of the strip brightness distribution in a quasar accretion disk from gravitational microlensing data,”
*Astronomy Reports*, vol. 46, no. 8, pp. 626–633, 2002. View at Publisher · View at Google Scholar · View at Scopus - M. J. Mortonson, P. L. Schechter, and J. Wambsganss, “Size is everything: Universal features of quasar microlensing with extended sources,”
*Astrophysical Journal*, vol. 628, no. 2, pp. 594–603, 2005. View at Publisher · View at Google Scholar · View at Scopus - L. J. Goicoechea, D. Alcalde, E. Mediavilla, and J. A. Muñoz, “Determination of the properties of the central engine in microlensed QSOs,”
*Astronomy and Astrophysics*, vol. 397, no. 2, pp. 517–525, 2003. View at Publisher · View at Google Scholar · View at Scopus - C. S. Kochanek, “Quantitative interpretation of quasar microlensing light curves,”
*Astrophysical Journal*, vol. 605, no. 1, pp. 58–77, 2004. View at Publisher · View at Google Scholar · View at Scopus - V. G. Vakulik, R. E. Schild, G. V. Smirnov, V. N. Dudinov, and V. S. Tsvetkova, “Q2237$+$0305 source structure and dimensions from light-curve simulation,”
*Monthly Notices of the Royal Astronomical Society*, vol. 382, no. 2, pp. 819–825, 2007. View at Publisher · View at Google Scholar · View at Scopus - T. Anguita, R. W. Schmidt, E. L. Turner et al., “The multiple quasar Q2237$+$
0305 under a microlensing caustic,”
*Astronomy and Astrophysics*, vol. 480, no. 2, pp. 327–334, 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Poindexter, N. Morgan, and C. S. Kochanek, “The spatial structure of an accretion disk,”
*Astrophysical Journal Letters*, vol. 673, no. 1, pp. 34–38, 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Poindexter and C. S. Kochanek, “The transverse peculiar velocity of the Q2237$+$0305 lens galaxy and the mean mass of its stars,”
*Astrophysical Journal Letters*, vol. 712, no. 1, pp. 658–667, 2010. View at Publisher · View at Google Scholar · View at Scopus - S. Poindexter and C. S. Kochanek, “Microlensing evidence that a type 1 quasar is viewed face-on,”
*Astrophysical Journal Letters*, vol. 712, no. 1, pp. 668–673, 2010. View at Publisher · View at Google Scholar · View at Scopus - I. M. Gel'fand and G. E. Shilov,
*Generalized Functions*, vol. 1, Academic Press, New York, NY, USA, 1964. - M. Dominik, “Revealing stellar brightness profiles by means of microlensing fold caustics,”
*Monthly Notices of the Royal Astronomical Society*, vol. 353, no. 1, pp. 118–132, 2004. View at Publisher · View at Google Scholar · View at Scopus - N. I. Shakura and R. A. Sunyaev, “Black holes in binary systems. Observational appearance,”
*Astronomy & Astrophysics*, vol. 24, pp. 337–355, 1973. View at Google Scholar - H. Bateman and A. Erdélyi,
*Higher Transcendental Functions*, vol. 1, McGraw-Hill, New York, NY, USA, 1953. - R. Schmidt, R. L. Webster, and G. F. Lewis, “Weighing a galaxy bar in the lens Q2237 $+$0305,”
*Monthly Notices of the Royal Astronomical Society*, vol. 295, no. 2, pp. 488–496, 1998. View at Google Scholar · View at Scopus - X. Dai, G. Chartas, E. Agol, M. W. Bautz, and G. P. Garmire, “Chandra observations of QSO 2237$+$0305,”
*Astrophysical Journal Letters*, vol. 589, no. 1, pp. 100–110, 2003. View at Publisher · View at Google Scholar · View at Scopus - E. V. Fedorova, V. I. Zhdanov, C. Vignali, and G. G. C. Palumbo, “Q2237$+$
0305 in X-rays: spectra and variability with XMM-Newton,”
*Astronomy and Astrophysics*, vol. 490, no. 3, pp. 989–994, 2008. View at Publisher · View at Google Scholar · View at Scopus - V. Vakulik, R. Schild, V. Dudinov et al., “Observational determination of the time delays in gravitational lens system Q2237$+$0305,”
*Astronomy and Astrophysics*, vol. 447, no. 3, pp. 905–913, 2006. View at Publisher · View at Google Scholar · View at Scopus - J. S. B. Wyithe and E. L. Turner, “Determining the microlens mass function from quasar microlensing statistics,”
*Monthly Notices of the Royal Astronomical Society*, vol. 320, no. 1, pp. 21–30, 2001. View at Publisher · View at Google Scholar · View at Scopus - A. B. Congdon, C. R. Keeton, and S. J. Osmer, “Microlensing of an extended source by a power-law mass distribution,”
*Monthly Notices of the Royal Astronomical Society*, vol. 376, no. 1, pp. 263–272, 2007. View at Publisher · View at Google Scholar · View at Scopus - Y. Tsapras, K. Horne, S. Kane, and R. Carson, “Microlensing limits on numbers and orbits of extrasolar planets from the 1998-2000 OGLE events,”
*Monthly Notices of the Royal Astronomical Society*, vol. 343, no. 4, pp. 1131–1144, 2003. View at Publisher · View at Google Scholar · View at Scopus - S. Mao and B. Paczyński, “Gravitational microlensing by double stars and planetary systems,”
*Astrophysical Journal Letters*, vol. 374, no. 2, pp. L37–L40, 1991. View at Google Scholar · View at Scopus - M. Jaroszyński and B. Paczyński, “A possible planetary event OGLE-2002-BLG-055,”
*Acta Astronomica*, vol. 52, no. 4, pp. 361–367, 2002. View at Google Scholar · View at Scopus - F. Chollet, “A new method of measuring stellar masses,”
*Comptes Rendus de l'Académie des Sciences B*, vol. 288, pp. 163–165, 1979 (French). View at Google Scholar - M. Hosokawa, K. Ohnishi, and T. Fukushima, “Astrometric microlensing and degradation of reference frames,” in
*Proceedings of the 25th General Assembly on Highlights of Astronomy (IAU '03)*, O. Engvold, Ed., p. 602, Astronomical Society of the Pacific, San Francisco, Calif, USA, 2005. - M. Hosokawa, K. Ohnishi, T. Fukushima, and M. Takeuti, “Parallactic variation of gravitational lensing and measurement of stellar mass,”
*Astronomy & Astrophysics*, vol. 278, pp. L27–L30, 1993. View at Google Scholar - J. Kovalevsky, F. Mignard, and M. Froschle, “Space astrometry prospects and limitations,” in
*Proceedings of the International Astronomical Union (IAU '79)*, J. Kovalevsky and V. A. Brumberg, Eds., vol. 114, pp. 369–382, 1979. - K. Ohnishi, M. Hosokawa, and T. Fukushima, “Secular component of apparent proper motion of QSOs induced by gravitational lens of the galaxy,”
*ASP Conference Proceedings Series*, vol. 289, pp. 461–464, 2003. View at Google Scholar - E. Hog, I. D. Novikov, and A. G. Polnarev, “MACHO photometry and astrometry,”
*Astronomy & Astrophysics*, vol. 294, pp. 287–294, 1995. View at Google Scholar - M. Miyamoto and Y. Yoshii, “Astrometry for determining the MACHO mass and trajectory,”
*Astronomical Journal*, vol. 110, no. 3, pp. 1427–1432, 1995. View at Google Scholar · View at Scopus - M. A. Walker, “Microlensed image motions,”
*Astrophysical Journal Letters*, vol. 453, no. 1, pp. 37–39, 1995. View at Google Scholar · View at Scopus - V. Bozza, “Trajectories of the images in binary microlensing,”
*Astronomy and Astrophysics*, vol. 374, no. 1, pp. 13–27, 2001. View at Google Scholar · View at Scopus - D. P. Bennett, S. H. Rhie, A. C. Becker et al., “Discovery of a planet orbiting a binary star system from gravitational microlensing,”
*Nature*, vol. 402, no. 6757, pp. 57–59, 1999. View at Publisher · View at Google Scholar · View at Scopus - M. Dominik and K. C. Sahu, “Astrometric microlensing of stars,”
*Astrophysical Journal*, vol. 534, no. 1, pp. 213–226, 2000. View at Google Scholar · View at Scopus - A. Gould and C. Han, “Astrometric resolution of severely degenerate binary microlensing events,”
*Astrophysical Journal*, vol. 538, no. 2, pp. 653–656, 2000. View at Google Scholar · View at Scopus - C. Han, “On the astrometric behavior of binary microlensing events,”
*Monthly Notices of the Royal Astronomical Society*, vol. 325, pp. L1281–L1287, 2001. View at Google Scholar - C. Han, S. H. Park, and Y. S. Lee, “Distribution of caustic-crossing intervals for galactic binary-lens microlensing events,”
*Monthly Notices of the Royal Astronomical Society*, vol. 314, no. 1, pp. 59–64, 2000. View at Google Scholar · View at Scopus - C. Han, M. S. Chun, and K. Chang, “Astrometric properties of gravitational binary-microlens events and their applications,”
*Astrophysical Journal Letters*, vol. 526, no. 1, pp. 405–410, 1999. View at Google Scholar · View at Scopus - M. Treyer and J. Wambsganss, “Astrometric microlensing of quasars dependence on surface mass density and external shear,”
*Astronomy and Astrophysics*, vol. 416, no. 1, pp. 19–34, 2004. View at Publisher · View at Google Scholar · View at Scopus - S. Mao and H. J. Witt, “Extended source effects in astrometric gravitational microlensing,”
*Monthly Notices of the Royal Astronomical Society*, vol. 300, no. 4, pp. 1041–1046, 1998. View at Google Scholar · View at Scopus - C. H. Lee, S. Seitz, A. Riffeser, and R. Bender, “Finite-source and finite-lens effects in astrometric microlensing,”
*Monthly Notices of the Royal Astronomical Society*, vol. 407, no. 3, pp. 1597–1608, 2010. View at Publisher · View at Google Scholar · View at Scopus - V. I. Zhdanov, A. N. Alexandrov, and S. A. Salata, “Motion of images of microlensed extended sources: analytical relations and numerical estimates for moderate optical depths,”
*Kinematika i Fizika Nebesnykh Tel*, vol. 16, pp. 336–345, 2000 (Russian). View at Google Scholar - G. F. Lewis and R. A. Ibata, “Quasar image shifts resulting from gravitational microlensing,”
*Astrophysical Journal Letters*, vol. 501, no. 2, pp. 478–485, 1998. View at Publisher · View at Google Scholar · View at Scopus - S. A. Salata and V. I. Zhdanov, “Statistical astrometric microlensing of extended sources,”
*Astronomical Journal*, vol. 125, no. 3, pp. 1033–1037, 2003. View at Publisher · View at Google Scholar · View at Scopus - J. Wambsganss, “Gravitational microlensing,” in
*Gravitational Lensing: Strong, Weak, and Micro*, G. Meylan, P. North, and P. Jetzer, Eds., Saas-Fee Advanced Course 33, Part 4, pp. 453–540, Springer, Berlin, Germany, 2006. View at Google Scholar - J. Wambsganss, “Gravitational lensing in astronomy,”
*Living Reviews in Relativity*, vol. 1, article 12, 1998. View at Google Scholar - A. L. Erickcek and N. M. Law, “Astrometric microlensing by local dark matter subhalos,”
*Astrophysical Journal Letters*, vol. 729, no. 1, article 49, 2011. View at Publisher · View at Google Scholar · View at Scopus - S. Proft, M. Demleitner, and J. Wambsganss, “Prediction of astrometric microlensing events during the Gaia mission,”
*Astronomy & Astrophysics*, vol. 536, article A50, 11 pages, 2011. View at Google Scholar - M. V. Sazhin, O. S. Sazhina, and M. S. Pshirkov, “Apparent motions of quasars due to microlensing,”
*Astronomy Reports*, vol. 55, pp. 954–961, 2011. View at Google Scholar - V. A. Belokurov and N. W. Evans, “Astrometric microlensing with the GAIA satellite,”
*Monthly Notices of the Royal Astronomical Society*, vol. 331, no. 3, pp. 649–665, 2002. View at Publisher · View at Google Scholar · View at Scopus - S. A. Klioner, “A practical relativistic model for microarcsecond astrometry in space,”
*Astronomical Journal*, vol. 125, no. 3, pp. 1580–1597, 2003. View at Publisher · View at Google Scholar · View at Scopus - N. S. Kardashev, “Cosmological proper motion,”
*Astronomicheskii Zhurnal*, vol. 63, pp. 845–849, 1986. View at Google Scholar - M. V. Sazhin, “A fundamental limit to the accuracy of astrometric measurements,”
*Astronomy Letters*, vol. 22, no. 5, pp. 573–577, 1996. View at Google Scholar · View at Scopus - M. V. Sazhin, V. E. Zharov, A. V. Volynkin, and T. A. Kalinina, “Microarcsecond instability of the celestial reference frame,”
*Monthly Notices of the Royal Astronomical Society*, vol. 300, no. 1, pp. 287–291, 1998. View at Google Scholar · View at Scopus - E. Fedorova, “Binary gravitational microlensing of extragalactic sources,”
*Visnyk Kyivskogo Universytetu. Astronomiya*, vol. 45, pp. 33–39, 2009. View at Google Scholar - K. Chang and S. Refsdal, “Star disturbances in gravitational lens galaxies,”
*Astronomy & Astrophysics*, vol. 132, pp. 168–178, 1984. View at Google Scholar - A. F. Zakharov and M. V. Sazhin, “Non-compact astronomical objects as microlenses,”
*Astronomy & Astrophysics*, vol. 18, pp. 27–38, 1999. View at Google Scholar - J. Diemand, B. Moore, and J. Stadel, “Earth-mass dark-matter haloes as the first structures in the early Universe,”
*Nature*, vol. 433, no. 7024, pp. 389–391, 2005. View at Publisher · View at Google Scholar · View at Scopus - L. D. Landau and E. M. Lifshitz,
*Fluid Mechanics*, Pergamon, New York, NY, USA, 1959. - K. A. Piragas, V. I. Zhdanov, V. V. Zhdanova, and I. T. Zhuk, “Light propagation in a weak gravitational field of a stochastic system of pointlike sources,”
*Soviet Physics Journal*, vol. 29, no. 12, pp. 1019–1023, 1986. View at Publisher · View at Google Scholar · View at Scopus - A. A. Minakov and V. N. Shalyapin, “Effect of the gravitational field of the Galaxy on the apparent position, brightness and spatial density of remote radiation sources. I. A lens model of the Galaxy and deflection angles of rays,”
*Kinematics and Physics of Celestial Bodies*, vol. 6, pp. 49–59, 1990. View at Google Scholar