We show that the asymptotic scaling of the magnification volume cross section corresponding to an elliptic umbilic caustic surface is $\mu^{-2.5}$ in the two-image region and $\mu^{-2}$ in the four-image region, where $\mu$ is the total unsigned magnification.

In this second paper, we develop an analytical theory of quasi-equatorial lensing by Kerr black holes. In this setting we solve perturbatively our general lens equation with displacement given in Paper I, going beyond weak-deflection Kerr lensing to third order in our expansion parameter epsilon, which is the ratio of the angular gravitational radius to the angular Einstein radius. We obtain new formulas and results for the bending angle, image positions, image magnifications, total unsigned magnification, and centroid, all to third order in epsilon and including the displacement. New results on the time delay between images are also given to second order in epsilon, again including displacement. For all lensing observables we show that the displacement begins to appear only at second order in epsilon. When there is no spin, we obtain new results on the lensing observables for Schwarzschild lensing with displacement.

We develop a unified, analytic framework for gravitational lensing by Kerr black holes. In this first paper we present a new, general lens equation and magnification formula governing lensing by a compact object. Our lens equation assumes that the source and observer are in the asymptotically flat region and does not require a small angle approximation. Furthermore, it takes into account the displacement that occurs when the light ray's tangent lines at the source and observer do not meet on the lens plane. We then explore our lens equation in the case when the compact object is a Kerr black hole. Specifically, we give an explicit expression for the displacement when the observer is in the equatorial plane of the Kerr black hole as well as for the case of spherical symmetry.

We provide a geometric explanation for the existence of magnification relations for the A, D, E family of caustic singularities, which were established in recent work. In particular, it was shown that for families of general mappings between planes exhibiting any of these caustic singularities, and for any non-caustic target point, the total signed magnification of the corresponding pre-images vanishes. As an application to gravitational lensing, it was also shown that, independent of the choice of a lens model, the total signed magnification vanishes for a light source anywhere in the four-image region close to elliptic and hyperbolic umbilic caustics. This is a more global and higher-order analog of the well-known fold and cusp magnification relations. We now extend each of these mappings to weighted projective space, which is a compact orbifold, and show that magnification relations translate into a statement about the behavior of these extended mappings at infinity. This generalizes multi-dimensional residue techniques developed in previous work, and introduces weighted projective space as a new tool in the theory of caustic singularities and gravitational lensing.

In the final paper of this series, we extend our results on magnification invariants to the infinite family of A, D, E caustic singularities. We prove that for families of general mappings between planes exhibiting any caustic singularity of the A, D, E family, and for a point in the target space lying anywhere in the region giving rise to the maximum number of lensed images (real pre-images), the total signed magnification of the lensed images will always sum to zero. The proof is algebraic in nature and relies on the Euler trace formula.

The mathematical theory of gravitational lensing has revealed many generic and global properties. Beginning with multiple imaging, we review Morse-theoretic image counting formulas and lower bound results, and complex-algebraic upper bounds in the case of single and multiple lens planes. We discuss recent advances in the mathematics of stochastic lensing, discussing a general formula for the global expected number of minimum lensed images as well as asymptotic formulas for the probability densities of the microlensing random time delay functions, random lensing maps, and random shear, and an asymptotic expression for the global expected number of micro-minima. Multiple imaging in optical geometry and a spacetime setting are treated. We review global magnification relation results for model-dependent scenarios and cover recent developments on universal local magnification relations for higher order caustics.

Continuing our development of a mathematical theory of stochastic microlensing, we study the random shear and expected number of random lensed images of different types. In particular, we characterize the first three leading terms in the asymptotic expression of the joint probability density function (p.d.f.) of the random shear tensor at a general point in the lens plane due to point masses in the limit of an infinite number of stars. Up to this order, the p.d.f. depends on the magnitude of the shear tensor, the optical depth, and the mean number of stars through a combination of radial position and the stars' masses. As a consequence, the p.d.f.s of the shear components are seen to converge, in the limit of an infinite number of stars, to shifted Cauchy distributions, which shows that the shear components have heavy tails in that limit. The asymptotic p.d.f. of the shear magnitude in the limit of an infinite number of stars is also presented. Extending to general random distributions of the lenses, we employ the Kac-Rice formula and Morse theory to deduce general formulas for the expected total number of images and the expected number of saddle images. We further generalize these results by considering random sources defined on a countable compact covering of the light source plane. This is done to introduce the notion of {\it global} expected number of positive parity images due to a general lensing map. Applying the result to microlensing, we calculate the asymptotic global expected number of minimum images in the limit of an infinite number of stars, where the stars are uniformly distributed. This global expectation is bounded, while the global expected number of images and the global expected number of saddle images diverge as the order of the number of stars.

We prove a theorem about magnification relations for all generic general caustic singularities up to codimension five: folds, cusps, swallowtail, elliptic umbilic, hyperbolic umbilic, butterfly, parabolic umbilic, wigwam, symbolic umbilic, 2nd elliptic umbilic, and 2nd hyperbolic umbilic. Specifically, we prove that for a generic family of general mappings between planes exhibiting any of these singularities, and for a point in the target lying anywhere in the region giving rise to the maximum number of real pre-images (lensed images), the total signed magnification of the pre-images will always sum to zero. The proof is algebraic in nature and makes repeated use of the Euler trace formula. We also prove a general algebraic result about polynomials, which we show yields an interesting corollary about Newton sums that in turn readily implies the Euler trace formula. The wide field imaging surveys slated to be conducted by the Large Synoptic Survey Telescope are expected to find observational evidence for many of these higher-order caustic singularities. Finally, since the results of the paper are for generic general mappings, not just generic lensing maps, the findings are expected to be applicable not only to gravitational lensing, but to any system in which these singularities appear.

Stochastic microlensing is a central tool in probing dark matter on galactic scales. From first principles, we initiate the development of a mathematical theory of stochastic microlensing. Beginning with the random time delay function and associated lensing map, we determine exact expressions for the mean and variance of these transformations. We characterize the exact p.d.f. of a normalized random time delay function at the origin, showing that it is a shifted gamma distribution, which also holds at leading order in the limit of a large number of point masses at a general point of the lens plane. For the large number of point masses limit, we also prove that the asymptotic p.d.f. of the random lensing map under a specified scaling converges to a bivariate normal distribution. We show analytically that the p.d.f. of the random scaled lensing map at leading order depends on the magnitude of the scaled bending angle due purely to point masses as well as demonstrate explicitly how this radial symmetry is broken at the next order. Interestingly, we found at leading order a formula linking the expectation and variance of the normalized random time delay function to the first Betti number of its domain. We also determine an asymptotic p.d.f. for the random bending angle vector and find an integral expression for the probability of a lens plane point being near a fixed point. Lastly, we show explicitly how the results are affected by location in the lens plane. The results of this paper are relevant to the theory of random fields and provide a platform for further generalizations as well as analytical limits for checking astrophysical studies of stochastic microlensing.

We prove that, independent of the choice of a lens model, the total signed magnification always sums to zero for a source anywhere in the four-image regions of swallowtail, elliptic umbilic, and hyperbolic umbilic caustics. This is a more global and higher-order analog of the well-known fold and cusp magnification relations, in which the total signed magnification in the two-image region of the fold, and the three-image region of the cusp, are both always zero. As an application, we construct a lensing observable for the hyperbolic umbilic magnification relation and compare it with the corresponding observables for the cusp and fold relations using a singular isothermal ellipsoidal lens. We demonstrate the greater generality of the hyperbolic umbilic magnification relation by showing how it applies to the fold image doublets and cusp image triplets, and extends to image configurations that are neither. We show that the results are applicable to the study of substructure on galactic scales using observed quadruple images of lensed quasars. The magnification relations are also proved for generic 1-parameter families of mappings between planes, extending their potential range of applicability beyond lensing.

A Kerr black hole with mass parameter m and angular momentum parameter a acting as a gravitational lens gives rise to two images in the weak field limit. We study the corresponding magnification relations, namely the signed and absolute magnification sums and the centroid up to post-Newtonian order. We show that there are post-Newtonian corrections to the total absolute magnification and centroid proportional to a/m, which is in contrast to the spherically symmetric case where such corrections vanish. Hence we also propose a new set of lensing observables for the two images involving these corrections, which should allow measuring a/m with gravitational lensing. In fact, the resolution capabilities needed to observe this for the Galactic black hole should in principle be accessible to current and near-future instrumentation. Since a/m >1 indicates a naked singularity, a most interesting application would be a test of the Cosmic Censorship conjecture. The technique used to derive the image properties is based on the degeneracy of the Kerr lens and a suitably displaced Schwarzschild lens at post-Newtonian order. A simple physical explanation for this degeneracy is also given.

The bending angle of light is a central quantity in the theory of gravitational lensing. We develop an analytical perturbation framework for calculating the bending angle of light rays lensed by a Schwarzschild black hole. Using a perturbation parameter given in terms of the gravitational radius of the black hole and the light ray's impact parameter, we determine an invariant series for the strong-deflection bending angle that extends beyond the standard logarithmic deflection term used in the literature. In the process, we discovered an improvement to the standard logarithmic deflection term. Our perturbation framework is also used to derive as a consistency check, the recently found weak deflection bending angle series. We also reformulate the latter series in terms of a more natural invariant perturbation parameter, one that smoothly transitions between the weak and strong deflection series. We then compare our invariant strong deflection bending-angle series with the numerically integrated exact formal bending angle expression, and find less than 1% discrepancy for light rays as far out as twice the critical impact parameter. The paper concludes by showing that the strong and weak deflection bending angle series together provide an approximation that is within 1% of the exact bending angle value for light rays traversing anywhere between the photon sphere and infinity.

Braneworld gravity is a model that endows physical space with an extra dimension. In the type II Randall-Sundrum braneworld gravity model, the extra dimension modifies the spacetime geometry around black holes, and changes predictions for the formation and survival of primordial black holes. We develop a comprehensive analytical formalism for far-field black hole lensing in this model, using invariant quantities to compute all geometric optics lensing observables. We then make the first analysis of wave optics in braneworld lensing, working in the semi-classical limit. We show that wave optics offers the only realistic way to observe braneworld effects in black hole lensing. We point out that if primordial braneworld black holes exist, have mass M, and contribute a fraction f of the dark matter, then roughly 3e5 x f (M/1e-18 Msun)^(-1) of them lie within our Solar System. These objects, which we call "attolenses," would produce interference fringes in the energy spectra of gamma-ray bursts at energies ~100 (M/1e-18 Msun)^(-1) MeV (which will soon be accessible with the GLAST satellite). Primordial braneworld black holes spread throughout the universe could produce similar interference effects; the probability for "attolensing" may be non-negligible. If interference fringes were observed, the fringe spacing would yield a simple upper limit on M. Detection of a primordial black hole with M <~ 1e-19 Msun would challenge general relativity and favor the braneworld model. Further work on lensing tests of braneworld gravity must proceed into the physical optics regime, which awaits a description of the full spacetime geometry around braneworld black holes.

We study gravitational lensing by compact objects in gravity theories that can be written in a Post-Post-Newtonian (PPN) framework: i.e., the metric is static and spherically symmetric, and can be written as a Taylor series in m/r, where m is the gravitational radius of the compact object. Working invariantly, we compute corrections to standard weak-deflection lensing observables at first and second order in the ratio of the angular gravitational radius to the angular Einstein ring radius of the lens. We show that the first-order corrections to the total magnification and centroid position vanish universally for gravity theories that can be written in the PPN framework. This arises from some surprising, fundamental relations among the lensing observables in PPN gravity models. We derive these relations for the image positions, magnifications, and time delays. A deep consequence is that any violation of the universal relations would signal the need for a gravity model outside the PPN framework (provided that some basic assumptions hold). In practical terms, the relations will guide observational programs to test general relativity, modified gravity theories, and possibly the Cosmic Censorship conjecture. We use the new relations to identify lensing observables that are accessible to current or near-future technology, and to find combinations of observables that are most useful for probing the spacetime metric. We give explicit applications to the Galactic black hole, microlensing, and the binary pulsar J0737-3039.

We are developing a general, unified, and rigorous analytical framework for using gravitational lensing by compact objects to test different theories of gravity beyond the weak-deflection limit. In this paper we present the formalism for computing corrections to lensing observables for static, spherically symmetric gravity theories in which the corrections to the weak-deflection limit can be expanded as a Taylor series in one parameter, namely the gravitational radius of the lens object. We take care to derive coordinate-independent expressions and compute quantities that are directly observable. We compute first- and second-order corrections to the image positions, magnifications, and time delays. Interestingly, we find that the first-order corrections to the total magnification and centroid position vanish in all gravity theories that agree with general relativity in the weak-deflection limit, but they can remain nonzero in modified theories that disagree with GR in the weak-deflection limit. For the Reissner-Nordstrom metric and a related metric from heterotic string theory, our formalism reveals an intriguing connection between lensing observables and the condition for having a naked singularity, which could provide an observational method for testing the existence of such objects. We apply our formalism to the Galactic black hole and predict that the corrections to the image positions are at the level of 10 micro-arcseconds, while the correction to the time delay is a few hundredths of a second. These corrections would be measurable today if a pulsar were found to be lensed by the Galactic black hole; and they should be readily detectable with planned missions like MAXIM.

When the source in a 4-image gravitational lens lies close to a "fold" caustic, two of the lensed images lie close together. If the lens potential is smooth on the scale of the separation between the two close images, then the difference between their fluxes should approximately vanish, Rfold = (F1-F2)/(F1+F2) \approx 0. Violations of this "fold relation" in observed lenses are thought to indicate the presence of structure on scales smaller than the separation between the close images. We study the fold relation and find it to be more subtle and rich than was previously realized. The degree to which Rfold can differ from zero for realistic smooth lenses depends not only on the distance of the source from the caustic, but also on its location along the caustic, and on the angular structure of the lens potential. Thus, it is impossible to say from Rfold alone whether observed flux ratios are anomalous or not. Instead, we must consider the full distribution of Rfold values that can be obtained from smooth potentials that reproduce the separation between the two close images and the distance to the next nearest image. We analyze the generic and specific features of this distribution, and then show that 5 of the 12 known lenses with fold configurations have flux ratio anomalies: B0712+472, SDSS 0924+0219, PG 1115+080, B1555+375, and B1933+503. Combining this with our previous analysis revealing anomalies in 3 of the 4 known lenses with cusp configurations, we conclude that at least half (8/16) of all 4-image lenses that admit generic, local analyses exhibit flux ratio anomalies indicative of small-scale structure. The fold and cusp relations do not reveal the nature of the implied small-scale structure, but do provide the formal foundation for substructure studies. [Abridged]

The inability of standard models to explain the flux ratios in many 4-image gravitational lenses has been cited as evidence for significant small-scale structure in lens galaxies. That claim has generally relied on detailed lens modeling, so it is both model dependent and somewhat difficult to interpret. We present a more robust and generic method for identifying lenses with small-scale structure. For a close triplet of images associated with a source near a cusp caustic, the sum of the signed magnifications should approximately vanish. We derive realistic upper bounds on the sum, and argue that lenses with flux ratios that significiantly violate the bounds can be said to have structure in the lens potential on scales smaller than the image separation. Five observed lenses have such flux ratio ``anomalies'': B2045+265, 1RXS J1131-1231, and SDSS J0924+0219 have strong anomalies; B0712+472 has a strong anomaly at optical/near-IR wavelengths and a marginal anomaly at radio wavelengths; and RX J0911+0551 appears to have an anomaly, but this conclusion is subject to uncertainties about octopole modes in early-type galaxies. Analysis of the cusp relation does not identify the known anomaly in B1422+231, so methods that are more sophisticated (and less generic) than the cusp relation may be necessary to uncover flux ratio anomalies in some systems. Although these flux ratio anomalies might represent milli- or micro-lensing, we cannot identify the cause; we can only conclude that the lenses have significant structure in the potential on scales smaller than the separation between the images. Additional arguments must be invoked to specify the nature of this small-scale structure. [Abridged]

The standard treatment of gravitational lensing by a point mass lens M is based on a weak-field deflection angle a = 2/x, where x = (r c^2)/(2 G M) with r the distance of closest approach to the mass of a lensed light ray. It was shown that for a point mass lens, the total magnification and image centroid shift of a point source remain unchanged by relativistic corrections of second order in 1/x. This paper considers these issues analytically taking into account the relativistic images, under three standard lensing configuration assumptions, for a Schwarzschild black hole lens with background point and extended sources having arbitrary surface brightness profiles. We apply our results to the Galactic black hole for lensing scenarios where our assumptions hold. We show that a single factor characterizes the full relativistic correction to the weak-field image centroid and magnification. As the lens-source distance increases, the relativistic correction factor strictly decreases. In particular, we find that for point and extended sources about 10 pc behind the black hole (which is a distance significantly outside the tidal disruption radius of a sun-like source), the relativistic correction factor is miniscule, of order 10^{-14}. Therefore, for standard lensing configurations, any detectable relativistic corrections to microlensing by the Galactic black hole will most likely have to come from sources significantly closer to the black hole.

Very little attention has been paid to the properties of optical wavefronts and caustic surfaces due to gravitational lensing. Yet the wavefront-based point of view is natural and provides insights into the nature of the caustic surfaces on a gravitationally lensed lightcone. We derive analytically the basic equations governing the wavefronts, lightcones, caustics on wavefronts, and caustic surfaces on lightcones in the context of weak-field, thin-screen gravitational lensing. These equations are all related to the potential of the lens. In the process, we also show that the standard single-plane gravitational lensing map extends to a new mapping, which we call a wavefront lensing map. Unlike the standard lensing map, the Jacobian matrix of a wavefront lensing map is not symmetric. Our formulas are then applied to caustic ``surfing.'' By surfing a caustic surface, a space-borne telescope can be fixed on a gravitationally lensed source to obtain an observation of the source at very high magnification over an extended time period, revealing structure about the source that could not otherwise be resolved. Using our analytical expressions for caustic sheets, we present a scheme for surfing a caustic sheet of a lensed source in rectilinear motion. Detailed illustrations are also presented of the possible types of wavefronts and caustic sheets due to nonsingular and singular elliptical potentials, and singular isothermal spheres, including an example of caustic surfing for a singular elliptical potential lens.

We present a rigorous, detailed study of the generic, quantitative properties of gravitational microlensing near cusp catastrophes. We derive explicit formulas for the total magnification and centroid of the images created for sources outside, on, and inside the cusped caustic. We obtain results on how the image magnifications scale with respect to separation from the cusped caustic for arbitrary source positions. Along the axis of symmetry of the cusp, the total magnification $\mu$ scales as $\mu \propto u^{-1}$, where $u$ is the distance of the source from the cusp, whereas perpendicular to this axis, $\mu \propto u^{-2/3}$. When the source passes through a point $\bu_0$ on a fold arc abutting the cusp, the image centroid has a jump discontinuity; we present a formula for the size of the jump in terms of the local derivatives of the lens potential. The total magnifications for a small extended source located both on, and perpendicular to, the axis of symmetry are also derived, for both uniform and limb darkened surface brightness profiles. We find that the difference in magnification between a finite and point source is $< 5%$ for separations of $> 2.5$ source radii from the cusp point, while the effect of limb-darkening is $< 1%$ in the same range. Our predictions for the astrometric and photometric behavior of both pointlike and finite sources passing near a cusp are illustrated and verified using numerical simulations of the cusp-crossing Galactic binary-lens event MACHO-1997-BUL-28. We discuss several possible applications of our results to Galactic binary lenses and quasar microlensing.

We study the local behavior of gravitational lensing near fold catastrophes. Using a generic form for the lensing map near a fold, we determine the observable properties of the lensed images, focusing on the case when the individual images are unresolved, i.e., microlensing. Allowing for images not associated with the fold, we derive analytic expressions for the photometric and astrometric behavior near a generic fold caustic. We show how this form reduces to the more familiar linear caustic, which lenses a nearby source into two images which have equal magnification, opposite parity, and are equidistant from the critical curve. In this case, the simplicity and high degree of symmetry allows for the derivation of semi-analytic expressions for the photometric and astrometric deviations in the presence of finite sources with arbitrary surface brightness profiles. We use our results to derive some basic properties of astrometric microlensing near folds, in particular we predict for finite sources with uniform and limb darkening profiles, the detailed shape of the astrometric curve as the source crosses a fold. We find that the astrometric effects of limb darkening will be difficult to detect with the currently planned accuracy of the Space Interferometry Mission. We verify our results by numerically calculating the expected astrometric shift for the photometrically well-covered Galactic binary lensing event OGLE-1999-BUL-23, finding excellent agreement with our analytic expressions. Our results can be applied to any lensing system with fold caustics, including Galactic binary lenses and quasar microlensing.

The generic, qualitative, local behavior of center-of-light curves near folds and cusps are studied. The results apply to any finite number of lens planes.

We study the critical curves, caustics, and multiple imaging due to one of the simplest many-body gravitational lens configurations: equal-mass point masses on the vertices of a regular polygon. Some examples of the critical curves and caustics are also displayed. We pose the problem of determining the maximum number of lensed images due to regular-polygon and general point mass configurations. Our numerical simulations suggest a maximum that is linear, rather than quadratic, in the number of point masses.