Magnification Cross Sections for the Elliptic Umbilic Caustic Surface
MMagnification Cross Sections for the Elliptic Umbilic Caustic Surface
Amir B. Aazami
Department of Mathematics & Computer Science,Clark University,Worcester, MA 01610; [email protected]
Charles R. Keeton
Department of Physics & Astronomy,Rutgers University, 136 Frelinghuysen Road,Piscataway, NJ 08854-8019; [email protected]
Arlie O. Petters
Departments of Mathematics and Physics, Duke University,Science Drive, Durham, NC 27708-0320; [email protected]
In gravitational lensing, magnification cross sections characterize the probability that a light sourcewill have magnification greater than some fixed value, which is useful in a variety of applications.The (area) cross section is known to scale as µ − for fold caustics and µ − . for cusp caustics. Weaim to extend these results to higher-order caustic singularities, focusing on the elliptic umbilic,which can be manifested in lensing systems with two or three galaxies. The elliptic umbilic has acaustic surface, and we show that the volume cross section scales as µ − . in the two-image regionand µ − in the four-image region, where µ is the total unsigned magnification. In both cases ourresults are supported both numerically and analytically. I. INTRODUCTION
Magnification cross sections are an important tool in gravitational lensing. Knowing the magnification cross sectionallows one to determine the probability that a light source will have magnification greater than some fixed value.This in turn gives information about the accuracy of cosmological models, since these predict different probabilitiesregarding source magnifications (see Schneider et al. 1992 [1], Kaiser 1992 [2], Bartelmann et al. 1998 [3], and Petterset al. 2001 [4, Chapter 13]). Knowing magnification cross sections is also important for observational programs thatuse lensing magnification to help detect extremely faint galaxies (see, e.g., Lotz et al. 2017 [5], and references therein).It is well known that for the so-called “fold” and “cusp” caustic singularities, the area cross sections scale asymp-totically as µ − and µ − . , respectively, where µ is the total unsigned magnification of a lensed source (for the cuspcaustic, it is assumed that the source lies in the one-image region locally). In the case of single-plane lensing, the foldscaling was determined by Blandford and Narayan 1986 [6], while the cusp scaling was determined by Mao 1992 [7]and Schneider and Weiss 1992 [8]. For multiple-plane lensing, the fold and cusp scalings were determined by Petterset al. [4, Chapter 13].In this paper we commence the study of magnification cross sections of “higher-order” caustic surfaces; in particular,we derive the (single-plane) asymptotic limit of the magnification cross section for the “elliptic umbilic” caustic surface,which is not a curve but rather a two-dimensional stable caustic surface in a three-dimensional parameter space (forprecise definitions, consult, e.g., Arnold 1973 [16], Callahan 1974 & 1977 [17, 18], Majthay 1985 [9], Arnold et al.1985 [10], Petters 1993 [19], [1, 4]); its magnification “cross section” is therefore a region with volume. As hypothesizedin Rusin et al. 2001 [11] and Blandford 2001 [12], the elliptic umbilic is likely to be manifested inside a triangle formedby three lensing galaxies, and to involve one positive-parity and three negative-parity images. As shown in Shin andEvans 2008 [13] and de Xivry and Marshall 2009 [14], elliptic umbilics can also appear in lensing by binary galaxies, ifthe binary separation is small enough. The asymptotic scaling of the elliptic umbilic volume cross section will dependon whether the source gives rise to two or four lensed images locally. We show that, in the two-image region, thisvolume cross section scales to leading order as µ − . , whereas in the four-image region, its leading order scales as µ − . In both cases, our results are supported numerically and analytically. In our derivation of µ − in the four-imageregion, we make use of a certain magnification relation that holds for higher-order caustic singularities, and the ellipticumbilic in particular, that was shown to hold in Aazami and Petters 2009 [15]. a r X i v : . [ a s t r o - ph . C O ] J u l FIG. 1: The elliptic umbilic caustic surface in three-dimensional parameter space { ( s , s , c ) } = R ; s , s are source planecoordinates and c ∈ R an additional parameter. When c = 0 the elliptic umbilic is the central point shown. A c -slice of thiscaustic is shown in Figure 2. I D ABC ← caustic curve FIG. 2: A generic c -slice of the elliptic umbilic caustic surface, which is the triangular-shaped solid curve. For level curves inthe two-image region, which is the (unbounded) region outside the closed caustic curve, the magnification volume cross section σ ( µ ) is, for µ large enough, the sum of the closed regions A, B, C . A source inside regions
A, B , or C has two lensed images andtotal unsigned magnification greater than µ . For level curves in the four-image region, which is inside the triangular-shapedsolid curve, σ ( µ ) is, for µ large enough, the region D inside the triangular-shaped caustic but outside the closed dashed curvewithin. A source in region D has four lensed images and total unsigned magnification greater than µ . Note that the actualcaustic and level “curve” are both surfaces in R ; what is shown here is a c -slice of this surface on the source plane. II. THE ELLIPTIC UMBILIC CAUSTIC SURFACE
Let ( x, y ) denote coordinates on the lens plane and ( s , s ) those on the source plane. Like all higher-order causticsurfaces, the elliptic umbilic caustic has parameters in addition to the two source plane coordinates. (Such higher-order parameters can, depending on the setting, be used to model the source redshift, radii of galaxies, ellipticities,distance along the line of sight, etc., of the lens system in question; see, e.g., [1, Chapter 8]). A gravitational lensingmap in the neighborhood of an elliptic umbilic critical point, as derived in [1, Chapter 5], takes the form η ( x, y ) = ( x − y , − xy + 4 cy ) = ( s , s ) , (1)where c ∈ R is a parameter in addition to the source plane coordinates s , s ; i.e., the elliptic umbilic caustic is asurface in the parameter space { ( s , s , c ) } = R , not a curve in the source plane S := { ( s , s ) } = R . Accordingly,the “magnification cross section” is a (three-dimensional) volume, and the asymptotic scaling is the leading order termin the limit as the magnification goes to infinity. Figure 1 shows the elliptic umbilic caustic surface, while Figure 2 − − − − − c = − . − − − − − c = − . − − − − − c = − . − − − − − c = − . − − − − − c = 0 . − − − − − c = 0 . − − − − − c = 0 . − − − − − c = 0 . − − − − − c = 1 .
00 0 . . . . . . . . . m a g n i fi c a t i o n FIG. 3: Level curves of the total unsigned magnification, for different values of c . shows a c -slice of it on the source plane S .If a source located at ( s , s ) on the source plane has a lensed image located at ( x, y ) on the lens plane, then themagnification ˜ µ of this lensed image is given by˜ µ ( x, y ) = 1det(Jac η )( x, y ) = 18 cx − x + y ) · (2)Now fix c ∈ R and µ >
0. Consider first the four-image region enclosed by the caustic curve in Figure 2. Let C µ denote the subset consisting of those source positions ( s , s ) in the four-image region with total unsigned magnificationequal to µ , where the unsigned magnification of each image ( x i , y i ) belonging to ( s , s ) is given by | ˜ µ ( x i , y i ) | in (2). C µ will consist of the closed dashed curve inside the caustic curve in Figure 2; any source inside the dashed region willhave total unsigned magnification less than µ , while a source in the region D will have total unsigned magnificationequal to µ . Likewise for the two-image region, which is the (unbounded) region outside the caustic curve: for µ largeenough, the enclosed regions A , B , and C comprise those sources whose two lensed images will have total unsignedmagnification greater than µ . (By symmetry, A , B , and C all have the same areas.)Whether in the two- or the four-image region, the asymptotic scaling there is determined as follows. The areaslabeled A , B , C , and D will in general be functions of µ and c ; denote any one of these, e.g., by A ( µ, c ). To obtain a volume section, denoted V ( µ ), we integrate V ( µ ) = (cid:90) c c A ( µ, c ) dc, (3)where c < c are arbitrary but of the same sign. Finally, we take the limit lim µ →∞ V ( µ ) and identify the leadingorder term; this is the asymptotic scaling of the elliptic umbilic magnification volume.We first examine the scalings numerically. For a given c -slice, we compute the total magnification on a grid inplane S , as shown in Figure 3. The grid is adaptive, meaning that more grid points are used in regions where the magnification, µ − − − − − − - i m a g e cr o sss e c t i o n , A ( µ , c ) c = 0 . c = 0 . c = 0 . magnification, µ − − − − - i m a g e cr o sss e c t i o n , V ( µ ) FIG. 4: (Left)
The colored points show the area cross section A ( µ, c ) in the two-image region, computed numerically fordifferent c -slices. The dotted lines show the scaling A ( µ, c ) ∝ µ − . . (Right) The points show the volume cross section V ( µ ) inthe two-image region, computed by integrating c over the range [0 . , . V ( µ ) ∝ µ − . . magnification, µ − − − - i m a g e cr o sss e c t i o n , A ( µ , c ) c = 0 . c = 0 . c = 0 . FIG. 5: The colored points show the area cross section A ( µ, c ) in the four-image region, computed numerically for different c -slices. The solid line shows A ( µ, c ) = π/ (2 µ ), in agreement with the result obtained analytically in (4) (with ˜ µ = µ/ magnification changes quickly and high resolution is needed to obtain accurate results. We use the grid to approximatethe area integral and compute A ( µ, c ). We then combine different c -slices to approximate the integral in (3) and obtain V ( µ ). Figure 4 shows examples of A ( µ, c ) and V ( µ ) for the two-image region, while Figure 5 shows A ( µ, c ) for thefour-image region.We support our numerical findings with the following analytical arguments. For the four-image region, there is infact a succinct argument that confirms the numerical result µ − , as follows. For any source in the four-image region(and assuming c (cid:54) = 0), the lensing map (1) has three lensed images with negative magnification and one lensed imagewith positive magnification, where the magnification of a lensed image is given by (2). Fix ˜ µ > S ˜ µ denotethe set of sources in the four-image region whose one positive-magnification image has magnification ˜ µ ; as usual, let µ denote the total unsigned magnification of this source. In fact S ˜ µ comprises the closed dashed curve inside the causticthat we saw in Figure 2. It is straightforward to compute that the area labeled D — that is, the area outside S ˜ µ butinside the caustic — is equal to 2 πc − (cid:16) πc − π µ (cid:17) = π µ · (4)Observe that this area is c -independent. This area is related to A ( µ, c ), the magnification (area) cross section, asfollows. Let (cid:126)s be a source in the four-image region whose total unsigned magnification is µ . Clearly (cid:126)s / ∈ S µ , since its FIG. 6: In each panel, the blue triangular-shaped curve is (a c -slice of) the elliptic umbilic caustic. In the leftmost figure, theblack curve comprises all source positions with at least one image with (unsigned) magnification ˜ µ , for some fixed value of ˜ µ .In the large magnification limit, the magnification of one of the images in the two-image region dominates that of the other,in which case the loops enclosing each cusp in the leftmost panel approach the regions A, B, C in Figure 2. The closed blackcurves in the middle and rightmost figures bound the areas
A, B , and C from above and below. It can be shown that thesetwo area bounds, which are functions of µ and c in general, both scale to leading order in µ as ∼ µ − . , where µ is the totalunsigned magnification of a source in the two-image region. one positive-magnification image must be less than µ . However, it is known that the four-image region of the ellipticumbilic, for any c (cid:54) = 0, satisfies the following magnification relation, (cid:88) i =1 ˜ µ i = 0 , with ˜ µ i is the signed magnification of lensed image i ; see [15] for a proof. (For clarity, we use the notation “˜ µ ”to denote the magnification belonging to an individual image, and reserve the notation “ µ ” to denote the totalunsigned magnification of a source.) It follows that (cid:126)s ∈ S µ/ , because the positive-magnification image must havemagnification µ/
2, since the other three magnifications are negative and must cancel it out. Thus the closed curve S µ/ is precisely the level curve of sources with total unsigned magnification µ . Observe that when ˜ µ = µ/ A ( µ, c ) = π/ (2 µ ), as in Figure 5. And thus the area cross section scales like (4). Figure 5 confirms this resultnumerically. Integrating (4) from c to c yields a volume cross section that clearly scales as V ( µ ) ∝ µ − .There remains, finally, the two-image region outside the caustic curve shown in Figure 2. Here, for µ > greater than µ comprise the regions labeled A , B , and C ,which three enclosed areas are equal by symmetry. In this case the area enclosed by them is not as easily derivableanalytically, due to the complicated nature of the intersection points of the dashed curves with the caustic. However,it is possible to bound the areas A , B , and C , from above and below, and to show that these lower and upper bounds,which are functions of µ and c , both scale to leading order in µ as ∼ µ − . , thereby supporting our numerical resultsin Figure 4. Figure 6 briefly describes the procedure, foregoing technical details. CONCLUSION
We have shown that the asymptotic scaling of the magnification volume cross section corresponding to an ellipticumbilic caustic surface is µ − . in the two-image region and µ − in the four-image region, where µ is the total unsignedmagnification. In both cases our results are supported both numerically and analytically. The goal was to extendto higher-order caustic singularities the well known (area) cross sections for fold ( µ − ) and cusp ( µ − . ) caustics, inparticular given that the elliptic umbilic caustic can be manifested in both binary and three-galaxy lensing systems. ACKNOWLEDGMENTS
Some of this work was done at the American Mathematical Society’s June, 2018, Mathematics Research Conferenceon the “Mathematics of Gravity and Light.” This material is based upon work supported in part by the NationalScience Foundation under Grant No. 1641020. [1] Schneider, P.; Ehlers, J.; Falco, E.
Gravitational Lenses ; Springer: Berlin, Germany, 1992.[2] Kaiser, N.
New Insights into the Universe ; Martinez, V.J., Portilla, M., Saez, S., Eds.; Springer: Berlin, Germany, 1992.[3] Bartelmann, M.; Huss, A.; Colberg, J.M.; Jenkins, A.; Pearce, F.R. Arc statistics with realistic cluster potentials IV.Clusters in different cosmologies.
Astron. Astrophys. , , 1–9.[4] Petters, A.O.; Levine, H.; Wambsganss, J. Singularity Theory and Gravitational Lensing ; Birkh¨auser: Basel, Switzerland,2001.[5] Lotz, J.M.; Koekemoer, A.; Coe, D.; Grogin, N.; Capak, P.; Mack, J.; Anderson, J.; Avila, R.; Barker, E.A.; Borncamp,D. The Frontier Fields: survey design and initial results.
Astrophys. J. , , 97.[6] Blandford, R.D.; Narayan, R. Fermat’s principle, caustics, and the classification of gravitational lens images. Astrophys.J. , , 568–582.[7] Mao, S. Gravitational microlensing by a single star plus external shear. Astrophys. J. , , 63–67.[8] Schneider, P.; Weiss, A. The gravitational lens equation near cusps. Astron. Astrophys. , , 1.[9] Majthay, A. Foundations of Catastrophe Theory ; Pitman: Wetherby, UK, 1985.[10] Arnold, V.I.; Gusein-Zade, S.M.; Varchenko, A.N.
Singularities of Differentiable Maps, Volume 1 ; Birkh¨auser: Basel,Switzerland, 1985.[11] Rusin, D.; Kochanek, C.S.; Norbury, M.; Falco, E.E.; Impey, C.D.; Lehr, J.; McLeod, B.A.; Rix, H.-W.; Keeton, C.R.;Muoz, J.A.; et al. B1359+154: A Six-Image Lens Produced by a z (cid:39) Astrophys. J. , , 594.[12] Blandford, R.D. The Future of Gravitational Optics. Publ. Astron. Soc. Pac. , , 789.[13] Shin, E.M.; Evans, N.W. Lensing by binary galaxies modelled as isothermal spheres. Mon. Not. R. Astron. Soc. , , 505–522.[14] Orban de Xivry, G.; Marshall, P. An atlas of predicted exotic gravitational lenses. Mon. Not. R. Astron. Soc. , ,2–20.[15] Aazami, A.B.; Petters, A.O. A universal magnification theorem for higher-order caustic singularities. J. Math. Phys. , , 032501.[16] Arnold, V.I. Normal forms of functions near degenerate critical points, the Weyl groups of A k , D k , E k and Lagrangiansingularities Func. Anal. Appl. , , 3–25.[17] Callahan, J. Singularities and plane maps. Am. Math. Mon. , , 211–240.[18] Callahan, J. Singularities and plane maps II: sketching catastrophes. Am. Math. Mon. , , 765–803.[19] Petters, A.O. Arnolds singularity theory and gravitational lensing. J. Math. Phys. ,33