Since the first observation of triple-lens gravitational microlensing in 2006, analyses of six more events have been published by the end of 2018. In three events the lens was a star with two planets; four involved a binary star with a planet. Other possible triple lenses, such as triple stars or stars with a planet with a moon, are yet to be detected. The analysis of triple-lens events is hindered by the lack of understanding of the diversity of their caustics and critical curves. We present a method for identifying the full range of critical curves for a triple lens with a given combination of masses in an arbitrary spatial configuration. We compute their boundaries in parameter space, identify the critical-curve topologies in the partitioned regions, and evaluate their probabilities of occurrence. We demonstrate the analysis on three triple-lens models. For three equal masses the computed boundaries divide the parameter space into 39 regions yielding nine different critical-curve topologies. The other models include a binary star with a planet, and a hierarchical star--planet--moon combination of masses. Both have the same set of 11 topologies, including new ones with doubly nested critical-curve loops. The number of lensing regimes thus depends on the combination of masses -- unlike in the double lens, which has the same three regimes for any mass ratio. The presented approach is suitable for further investigations, such as studies of the changes occurring in nonstatic lens configurations due to orbital motion of the components or other parallax-type effects.

Among the 25 planetary systems detected up to now by gravitational microlensing, there are two cases of a star with two planets, and two cases of a binary star with a planet. Other, yet undetected types of triple lenses include triple stars or stars with a planet with a moon. The analysis and interpretation of such events is hindered by the lack of understanding of essential characteristics of triple lenses, such as their critical curves and caustics. We present here analytical and numerical methods for mapping the critical-curve topology and caustic cusp number in the parameter space of $n$-point-mass lenses. We apply the methods to the analysis of four symmetric triple-lens models, and obtain altogether 9 different critical-curve topologies and 32 caustic structures. While these results include various generic types, they represent just a subset of all possible triple-lens critical curves and caustics. Using the analyzed models, we demonstrate interesting features of triple lenses that do not occur in two-point-mass lenses. We show an example of a lens that cannot be described by the Chang-Refsdal model in the wide limit. In the close limit we demonstrate unusual structures of primary and secondary caustic loops, and explain the conditions for their occurrence. In the planetary limit we find that the presence of a planet may lead to a whole sequence of additional caustic metamorphoses. We show that a pair of planets may change the structure of the primary caustic even when placed far from their resonant position at the Einstein radius.

The interpretation of gravitational microlensing events caused by planetary systems or multiple stars is based on the n-point-mass lens model. The first planets detected by microlensing were well described by the two-point-mass model of a star with one planet. By the end of 2014, four events involving three-point-mass lenses had been announced. Two of the lenses were stars with two planetary companions each; two were binary stars with a planet orbiting one component. While the two-point-mass model is well understood, the same cannot be said for lenses with three or more components. Even the range of possible critical-curve topologies and caustic geometries of the three-point-mass lens remains unknown. In this paper we provide new tools for mapping the critical-curve topology and caustic cusp number in the parameter space of n-point-mass lenses. We perform our analysis in the image plane of the lens. We show that all contours of the Jacobian are critical curves of re-scaled versions of the lens configuration. Utilizing this property further, we introduce the cusp curve to identify cusp-image positions on all contours simultaneously. In order to track cusp-number changes in caustic metamorphoses, we define the morph curve, which pinpoints the positions of metamorphosis-point images along the cusp curve. We demonstrate the usage of both curves on simple two- and three-point-mass lens examples. For the three simplest caustic metamorphoses we illustrate the local structure of the image and source planes.