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.

Observations of several gravitationally microlensed quasars in X-rays revealed variations in the profile of the iron K$\alpha$ line in the course of microlensing events. We explore the effect by simulating a microlensing caustic crossing a spatially resolved model of emission from a thin accretion disk around a Kerr black hole. We demonstrate the sequence of spectral changes during the event, in particular the appearance of additional peaks and edges in the line profile due to microlensing. We trace the origin of these features to points on the disk, at which the total energy shift ($g$-factor) contours are tangent to the caustic. Contours tangent from the inner side of the caustic generate peaks, while those tangent from its outer side generate edges. We derive analytical shapes of the generated features and map the peak strength as a function of position of the tangent point on the disk. Since the features are determined by the positional geometry of the caustic relative to the $g$-factor contours, the same type of behavior can be expected in a much broader range of emission models. The sequence of line profile changes thus serves as a sensitive probe of the geometry and physics of the innermost region of the quasar accretion disk.

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.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Recently Sumi et al. (2011) reported evidence for a large population of planetary-mass objects (PMOs) that are either unbound or orbit host stars in orbits > 10 AU. Their result was deduced from the statistical distribution of durations of gravitational microlensing events observed by the MOA collaboration during 2006 and 2007. Here we study the feasibility of measuring the mass of an individual PMO through microlensing by examining a particular event, MOA-2011-BLG-274. This event was unusual as the duration was short, the magnification high, the source-size effect large and the angular Einstein radius small. Also, it was intensively monitored from widely separated locations under clear skies at low air masses. Choi et al. (2012) concluded that the lens of the event may have been a PMO but they did not attempt a measurement of its mass. We report here a re-analysis of the event using re-reduced data. We confirm the results of Choi et al. and attempt a measurement of the mass and distance of the lens using the terrestrial parallax effect. Evidence for terrestrial parallax is found at a 3 sigma level of confidence. The best fit to the data yields the mass and distance of the lens as 0.80 +/- 0.30 M_J and 0.80 +/- 0.25 kpc respectively. We exclude a host star to the lens out to a separation ~ 40 AU. Drawing on our analysis of MOA-2011-BLG-274 we propose observational strategies for future microlensing surveys to yield sharper results on PMOs including those down to super-Earth mass.

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

We present the discovery and mass measurement of the cold, low-mass planet MOA-2009-BLG-266Lb, made with the gravitational microlensing method. This planet has a mass of m_p = 10.4 +- 1.7 Earth masses and orbits a star of mass M_* = 0.56 +- 0.09 Solar masses at a semi-major axis of a = 3.2 (+1.9 -0.5) AU and an orbital period of P = 7.6 (+7.7 -1.5} yrs. The planet and host star mass measurements are enabled by the measurement of the microlensing parallax effect, which is seen primarily in the light curve distortion due to the orbital motion of the Earth. But, the analysis also demonstrates the capability to measure microlensing parallax with the Deep Impact (or EPOXI) spacecraft in a Heliocentric orbit. The planet mass and orbital distance are similar to predictions for the critical core mass needed to accrete a substantial gaseous envelope, and thus may indicate that this planet is a "failed" gas giant. This and future microlensing detections will test planet formation theory predictions regarding the prevalence and masses of such planets.

Aims: We present a detailed analysis of OGLE 2004-BLG-482, a relatively high-magnification single-lens microlensing event which exhibits clear extended-source effects. These events are relatively rare, but they potentially contain unique information on the stellar atmosphere properties of their source star, as shown in this study. Methods: Our dense photometric coverage of the overall light curve and a proper microlensing modelling allow us to derive measurements of the OGLE 2004-BLG-482 source star's linear limb-darkening coefficients in three bands, including standard Johnson-Cousins I and R, as well as in a broad clear filter. In particular, we discuss in detail the problems of multi-band and multi-site modelling on the expected precision of our results. We also obtained high-resolution UVES spectra as part of a ToO programme at ESO VLT from which we derive the source star's precise fundamental parameters. Results: From the high-resolution UVES spectra, we find that OGLE 2004-BLG-482's source star is a red giant of MK type a bit later than M3, with Teff = 3667 +/- 150 K, log g = 2.1 +/- 1.0 and an assumed solar metallicity. This is confirmed by an OGLE calibrated colour-magnitude diagram. We then obtain from a detailed microlensing modelling of the light curve linear limb-darkening coefficients that we compare to model-atmosphere predictions available in the literature, and find a very good agreement for the I and R bands. In addition, we perform a similar analysis using an alternative description of limb darkening based on a principal component analysis of ATLAS limb-darkening profiles, and also find a very good agreement between measurements and model predictions.

Gravitational microlensing is not only a successful tool for discovering distant exoplanets, but it also enables characterization of the lens and source stars involved in the lensing event. In high magnification events, the lens caustic may cross over the source disk, which allows a determination of the angular size of the source and additionally a measurement of its limb darkening. When such extended-source effects appear close to maximum magnification, the resulting light curve differs from the characteristic Paczynski point-source curve. The exact shape of the light curve close to the peak depends on the limb darkening of the source. Dense photometric coverage permits measurement of the respective limb-darkening coefficients. In the case of microlensing event OGLE 2008-BLG-290, the K giant source star reached a peak magnification of about 100. Thirteen different telescopes have covered this event in eight different photometric bands. Subsequent light-curve analysis yielded measurements of linear limb-darkening coefficients of the source in six photometric bands. The best-measured coefficients lead to an estimate of the source effective temperature of about 4700 +100-200 K. However, the photometric estimate from colour-magnitude diagrams favours a cooler temperature of 4200 +-100 K. As the limb-darkening measurements, at least in the CTIO/SMARTS2 V and I bands, are among the most accurate obtained, the above disagreement needs to be understood. A solution is proposed, which may apply to previous events where such a discrepancy also appeared.

We explore the sensitivity of two-point-mass gravitational microlensing to the extended nature of the source star, as well as the related sensitivity to its limb darkening. We demonstrate that the sensitive region, usually considered to be limited to a source-diameter-wide band along the caustic, is strongly expanded near cusps, most prominently along their outer axis. In the case of multi-component caustics, facing cusps may form a region with a non-negligible extended-source effect spanning the gap between them. We demonstrate that for smaller sources the size of the sensitive region extending from a cusp measured in units of source radii increases, scaling as the inverse cube root of the radius. We study the extent of different sensitivity contours and show that for a microlensed Galactic bulge giant the probability of encountering at least a 1% extended-source effect is higher than the probability of caustic crossing by 40-60% when averaged over a typical range of lens-component separations, with the actual value depending on the mass ratio of the components. We derive analytical expressions for the extended-source effect and chromaticity for a source positioned off the caustic. These formulae are more generally applicable to any gravitational lens with a sufficiently small source. Using exactly computed amplifications we test the often used linear-fold caustic approximation and show that it may lead to errors on the level of a few percent even in near-ideal caustic-crossing events. Finally, we discuss several interesting cases of observed binary and planetary microlensing events and point out the importance of our results for the measurement of stellar limb darkening from microlensing light curves.

Thirteen exo-planets have been discovered using the gravitational microlensing technique (out of which 7 have been published). These planets already demonstrate that super-Earths (with mass up to ~10 Earth masses) beyond the snow line are common and multiple planet systems are not rare. In this White Paper we introduce the basic concepts of the gravitational microlensing technique, summarise the current mode of discovery and outline future steps towards a complete census of planets including Earth-mass planets. In the near-term (over the next 5 years) we advocate a strategy of automated follow-up with existing and upgraded telescopes which will significantly increase the current planet detection efficiency. In the medium 5-10 year term, we envision an international network of wide-field 2m class telescopes to discover Earth-mass and free-floating exo-planets. In the long (10-15 year) term, we strongly advocate a space microlensing telescope which, when combined with Kepler, will provide a complete census of planets down to Earth mass at almost all separations. Such a survey could be undertaken as a science programme on Euclid, a dark energy probe with a wide-field imager which has been proposed to ESA's Cosmic Vision Programme.

We explore the sensitivity of limb darkening coefficients computed from stellar atmosphere models to different least-squares fitting methods. We demonstrate that conventional methods are strongly biased to fitting the stellar limb. Our suggested method of fitting by minimizing the radially integrated squared residual yields improved fits with better flux conservation. The differences of the obtained coefficients from commonly used values are observationally significant. We show that the new values are in better agreement with solar limb darkening measurements as well as with coefficients reported from analyses of eclipsing binary light curves.

We present an analysis of OGLE 2004-BLG-254, a high-magnification and relatively short duration microlensing event in which the source star, a Bulge K3-giant, has been spatially resolved by a point-like lens. We have obtained dense photometric coverage of the event light curve with OGLE and PLANET telescopes, as well as a high signal-to-noise ratio spectrum taken while the source was still magnified by 20, using the UVES/VLT spectrograph. Our dense coverage of this event allows us to measure limb darkening of the source star in the I and R bands. We also compare previous measurements of linear limb-darkening coefficients involving GK-giant stars with predictions from ATLAS atmosphere models. We discuss the case of K-giants and find a disagreement between limb-darkening measurements and model predictions, which may be caused by the inadequacy of the linear limb-darkening law.

We study the influence of general lens and source velocities on the gravitational deflection of light by single and two-point-mass microlenses with general axis orientation. We demonstrate that in all cases the lens equation preserves its form exactly. However, its parameters -- the Einstein radius and the binary-lens separation -- are influenced by the lens velocity. In Galactic microlensing settings the velocity mainly affects the inferred separation for wide binary-star or star+planet microlenses oriented close to the line of sight. We briefly discuss the case of lenses moving with highly relativistic velocities.

Observations of microlensing transit events can be used to measure the limb darkening of the lensed star. We discuss the advantages and drawbacks of several microlensing light curve inversion methods. The method of choice in this work is inversion by means of decomposition of the stellar surface brightness profile. We construct an ideal basis by principal component analysis of brightness profiles obtained from model atmosphere calculations. Limb darkening approximations using such a basis are superior to those using standard power-law limb darkening laws. We perform a full analysis of simulated single-lens microlensing transit events including a detailed error analysis of the method. In realistic events with a low impact parameter the brightness profile of the source can be recovered with a relative accuracy of 2% from the center of the source disk to 0.9 of the disk radius. We show that in the particular case of the observed MACHO Alert 95-30 event the intrinsic complex variability of the lensed red giant hinders efforts to recover its surface features.

During microlensing events with a small impact parameter, the amplification of the source flux is sensitive to the surface brightness distribution of the source star. Such events provide a means for studying the surface structure of target stars in the ongoing microlensing surveys, most efficiently for giants in the Galactic bulge. In this work we demonstrate the sensitivity of point-mass microlensing to small spots with radii $r_s\lesssim0.2$ source radii. We compute the amplification deviation from the light curve of a spotless source and explore its dependence on lensing and spot parameters. During source-transit events spots can cause deviations larger than 2%, and thus be in principle detectable. Maximum relative deviation usually occurs when the lens directly crosses the spot. Its numerical value for a dark spot with sufficient contrast is found to be roughly equal to the fractional radius of the spot, i.e., up to 20% in this study. Spots can also be efficiently detected by the changes in sensitive spectral lines during the event. Notably, the presence of a spot can mimic the effect of a low-mass companion of the lens in some events.

Gravitational microlensing provides a new technique for studying the surfaces of distant stars. Microlensing events are detected in real time and can be followed up with precision photometry and spectroscopy. This method is particularly adequate for studying red giants in the Galactic bulge. Recently we developed an efficient method capable of computing the lensing effect for thousands of frequencies in a high-resolution stellar spectrum. Here we demonstrate the effects of microlensing on synthesized optical spectra of red giant model atmospheres. We show that different properties of the stellar surface can be recovered from time-dependent photometry and spectroscopy of a point-mass microlensing event with a small impact parameter. In this study we concentrate on center-to-limb variation of spectral features. Measuring such variations can reveal the depth structure of the atmosphere of the source star.

We present an efficient method for computing lightcurves of an elliptical source which is microlensed by a point mass. The amplification of an extended source involves a two-dimensional integral over its surface brightness distribution. We show that for a general surface brightness profile with an elliptical symmetry, this integral can be reduced to one dimension. We derive analytical results for the entire lightcurve in the limit of low (e.g. planetary) mass lenses, and for the wings of all microlensing lightcurves in general. In both cases, the lightcurve carries information about deviations of the source from elliptical symmetry, e.g. due to spots. The method is used to find the amplification of a circular red giant photosphere and an inclined accretion disk. We demonstrate that microlensing of an emission line from a disk can be used to infer the disk velocity structure and surface brightness profile.