Detectability of extrasolar moons as gravitational microlenses
C. Liebig  J. Wambsganss
Astronomisches RechenInstitut, Zentrum
für Astronomie der Universität Heidelberg, Mönchhofstraße 1214, 69120
Heidelberg, Germany
Received 10 December 2009 / Accepted 25 April 2010
Abstract
We evaluate gravitational lensing as a technique for the detection of
extrasolar moons. Since 2004 gravitational microlensing has been
successfully applied as a detection method for extrasolar planets. In
principle, the method is sensitive to masses as low as an Earth mass or
even a fraction of it. Hence it seems natural to investigate the
microlensing effects of moons around extrasolar planets. We explore the
simplest conceivable triple lens system, containing one star, one
planet and one moon. From a microlensing point of view, this system can
be modelled as a particular triple with hierarchical mass ratios very
different from unity. Since the moon orbits the planet, the planetmoon
separation will be small compared to the distance between planet and
star. Such a configuration can lead to a complex interference of
caustics. We present detectability and detection limits by comparing
triplelens light curves to bestfit binary light curves as caused by a
doublelens system consisting of host star and planet  without moon.
We simulate magnification patterns covering a range of mass and
separation values using the inverse ray shooting technique. These
patterns are processed by analysing a large number of light curves and
fitting a binary case to each of them. A chisquared criterion is used
to quantify the detectability of the moon in a number of selected
triplelens scenarios. The results of our simulations indicate that it
is feasible to discover extrasolar moons via gravitational microlensing
through frequent and highly precise monitoring of anomalous Galactic
microlensing events with dwarf source stars.
Key words: gravitational lensing: micro  planets and satellites: general  methods: numerical  methods: statistical
1 Introduction
By now hundreds of extrasolar planets have been detected^{}. For all we know, none of the newly discovered extrasolar planets offers physical conditions permitting any form of life. But the search for planets potentially harbouring life and the search for indicators of habitability is ongoing. One of these indicators might be the presence of a large natural satellite  a moon  which stabilises the rotation axis of the planet and thereby the surface climate (Benn 2001). It has also been suggested that a large moon itself might be a good candidate for offering habitable conditions (Scharf 2006). In the solar system, most planets harbour moons. In fact, the moons in the solar system outnumber the planets by more than an order of magnitude. No moon has yet been detected around an extrasolar planet.
The majority of known exoplanets has been discovered through radial velocity measurements, with the first successful finding reported by Mayor & Queloz (1995). This method is not sensitive to satellites of those planets, because the stellar ``Doppler wobble'' is only affected by the orbital movement of the barycentre of a planet and its satellites, though higherorder effects could play a role eventually. Here we consider Galactic microlensing, which has led to the discovery of several relatively lowmass exoplanets since the first report of a successful detection by Bond et al. (2004), as a promising technique for the search for exomoons.
As early as 1999, it has been suggested that extrasolar moons might be detectable through transit observations (Sartoretti & Schneider 1999), either through direct observation of lunar occultation or through transit timing variations, as the moon and planet rotate around their common barycentre, causing time shifts of the transit ingress and egress (cf. also Holman & Murray 2005). In their simulations of spacebased gravitational microlensing Bennett & Rhie (2002) mention the possibility of discovering extrasolar moons similar to our own Moon. Later that year, Han & Han (2002) performed a detailed feasibility study whether microlensing offers the potential to discover an EarthMoon analogue, but concluded that finite source effects would probably be too severe to allow detections. Williams & Knacke (2004) published the quite original suggestion to look for spectral signatures of Earthsized moons in the absorption spectra of Jupitersized planets. Cabrera & Schneider (2007) proposed a sophisticated transit approach using ``mutual event phenomena'', i.e. photometric variation patterns due to different phases of occultation and light reflection of planet and satellite. Han (2008) undertook a new qualitative study of a number of triplelens microlensing constellations finding ``nonnegligible'' light curve signals to occur in the case of an Earthmass moon orbiting a 10 Earthmass planet, ``when the planetmoon separation is similar to or greater than the Einstein radius of the planet''. Lewis et al. (2008) analysed pulsar timeofarrival signals for lunar signatures. Kipping (2009a,b) refined and extended the transit timing models of exomoons to include transit duration variations, and Kipping et al. (2009) examined transit detectability of exomoons with Keplerclass photometry and concluded that in optimal cases moon detections down to should be possible.
We cover here several new aspects concerning the microlensing search for exomoons, extending the work of Han (2008). First, the detectability of lunar light curve perturbations is determined with a statistical significance test that does not need to rely on human judgement. Second, all parameters of the twodimensional threebody geometry, including the position angle of the moon with respect to the planetstar axis, are varied. Third, an unbiased extraction of light curves from the selected scenarios enables a tentative prediction of the occurrence rate of detectable lunar light curve signals. A more detailed account of this study is available as Liebig (2009).
The paper is structured as follows: in Sect. 2, we recall the fundamental equations of gravitational microlensing relevant to our work. Our method for quantifying the detection rates for extrasolar moons in selected lensing scenarios is presented in Sect. 3. In Sect. 4, we discuss the astrophysical implications of the input parameters of the simulations. Our results are presented in Sect. 5, together with a first interpretation and a discussion of potential problems of our method.
2 Basics of gravitational microlensing
The deflection of light by massive bodies is a consequence of the theory of general relativity (Einstein 1916) and has been experimentally verified since 1919 (Dyson et al. 1920), see Paczynski (1996) for an introduction to the field or Schneider et al. (2006) for a comprehensive review.
The typical scale of angular separations in gravitational
lensing is
the Einstein radius ,
the
angular radius of the ring of formally infinite image magnification
that appears when a source at a distance ,
a lens of mass Mat a distance ,
and the observer are perfectly aligned:
where G denotes the gravitational constant, c the speed of light. is the distance between lens plane and source plane; in the noncosmological distance scale of our galaxy holds true. With the Einstein radius also comes the characteristic time scale of transient gravitational lensing events, the Einstein time with the transverse velocity of the lens relative to the source.
Here we focus on the triple lens case with host star (S),
planet
(P) and moon (M), see
Fig. 1
for
illustration. The lens equation can be expressed using complex
coordinates, where
shall denote the angular source position and
the image
positions, cf. choice of notation in Witt
(1990)
and Gaudi et al. (1998).
stands for
the angular position of the
lensing body i. q_{ij}
is the mass ratio between lenses i and
j (
).
The lens equation gets the following
form, if the primary lens, the host star S, has
unit mass and is
placed in the origin of the lens plane,
This is a mapping from the lens plane to the source plane, which maps the images of a source star to its actual position in the source plane. As pointed out in Rhie (1997), and explicitly calculated in Rhie (2002), the triplelens equation is a tenthorder polynomial equation in .
In microlensing, the images cannot be resolved. Detectable is
only the
transient change in magnitude of the source star, when lens and source
star are in relative motion to each other. The total
magnification
of the base flux is obtained as the inverse of the determinant
of the Jacobian of the mapping Eq. (2),
Gravitational lensing changes the apparent solid angle of a source, not the surface brightness. The magnification is the ratio of the total solid angle of the images and the apparent solid angle of the unlensed source.
For a status report of the past, present and prospective future of planet searching via Galactic microlensing the recent white papers of 2008/2009 are a good source of reference (Gaudi et al. 2009; Beaulieu et al. 2008; Dominik et al. 2008; Bennett et al. 2009).
3 Method
The simplest gravitational lens system incorporating an extrasolar moon is a triplelens system consisting of the lensing star, a planet and a moon in orbit around that planet, as sketched in Fig. 1. Most likely, lunar effects will first show up as noticeable irregularities in light curves that have been initially observed and classified as light curves with planetary signatures. To measure the detectability of a given triplelens system among binary lenses, we have to determine whether the provided triplelens light curve differs statistically significantly from binarylens light curves.
Figure 1: Geometry of our triplelens scenario, not to scale. Five parameters have to be fixed: the mass ratios and , the angular separations in the lens plane and and the position angle of the moon . These parameters define uniquely the relative projected positions of the three bodies. 

Open with DEXTER 
To clarify the terminology: We investigate light curves that show a deviation from the singlelens case due to the presence of the planetary caustic, i.e. which would be modelled as a starplusplanet system in a first approximation. We call lunar detectability the fraction of those light curves which display a significant deviation from a starplusplanet lens model due to the presence of the moon. We measure the difference between a given triplelens light curve (which is taken to be the ``true'' underlying light curve of the event) and its bestfit binarylens counterpart. The bestfit binarylens light curve is found by a leastsquare fit. We then employ statistics to see whether the triplelens light curve could be explained as a normal fluctuation within the error boundaries of the binarylens light curve. If this is not the case, the moon is considered detectable.
Detection and characterisation are two separate problems in the search for extrasolar planets or moons, though, and we do not make statements about the latter. When characterising an observed light curve with clear deviations from the binary model, it is still possible that ambiguous solutions  lunar and nonlunar  arise. This does not reflect on our results. We simply give the fraction of triplelens light curves significantly deviating from binarylens light curves, without exploring whether they can be uniquely characterised.
A qualitative impression of the lunar influence on the caustic structure can be gained from the magnification patterns in Fig. 2. Extracted example light curves are presented in Fig. 3. Going beyond the qualitative picture, in this paper we quantify the detectability of an extrasolar moon in gravitational microlensing light curves in selected scenarios.
Figure 2: Details of the analysed triplelens magnification maps for an example mass/separation scenario with , , and (the standard scenario as summarised in Table 1), showing the caustic configurations for twelve different lunar positions completing a full circular orbit in steps of 30 . The relative positions of star, planet and moon are sketched in the lower left of each panel (not to scale). The side lengths of the individual frames are . A darker shade of grey corresponds to a higher magnification. The straight black lines mark source trajectories, the resulting light curves are displayed in Fig. 3. 

Open with DEXTER 
Figure 3: Sample of triplelens light curves corresponding to the source trajectories depicted in Fig. 2. The solid light curves were extracted with an assumed solar source size, ; the thin, dashed lines show the same light curves with . The magnification scale is given in negative magnitude difference (mag), i.e. the unmagnified baseline flux of the source is 0. The time scale in units of the Einstein time t_{E}, assuming a uniform relative motion of source and lens. 

Open with DEXTER 
Figure 4: Magnification pattern displaying the planetary caustic of a triplelens system with mass ratios and , corresponding, e.g., to an Mdwarf, a Saturnmass planet and an Earthmass satellite of that planet. Darker shades correspond to areas of higher magnification. The third body shows its influence as the small perturbation on the bottom cusp of the planetary caustic. The position angle of the moon is , i.e. in the orientation of this map the moon is located above the planet. The complete set of parameters is given in Table 1. The side length of the pattern is in stellar Einstein radii. A source with radius (2 , 3 , 5 ) is indicated as a black disc (dashed circles) in the lower left. Taking the straight line that cuts the lunar perturbation of the caustic as the source trajectory of a solar sized source results in the light curve shown in Fig. 6 (solid line). Depending on the assumed source size, different light curves are obtained, see also Figs. 9a to c. 

Open with DEXTER 
3.1 Ray shooting
The triplelens Eq. (2) is analytically solvable, but Han & Han (2002) pointed out that numerical noise in the polynomial coefficients caused by limited computer precision was too high ( ) when solving the polynomial numerically for the very small mass ratios of moon and star ( ). To avoid this, we employ the inverse rayshooting technique, which has the further advantage of being able to account for finite source sizes and nonuniform source brightness profiles more easily. It also gives us the option of incorporating additional lenses (further planets or moons) without increasing the complexity of the calculations. This technique was developed by Schneider & Weiß (1986), Kayser et al. (1986) and Wambsganss (1990,1999) and already applied to planetary microlensing in Wambsganss (1997).
Inverse rayshooting means that light rays are traced from the observer back to the source plane. This is equivalent to tracing rays from the background source star to the observer plane. The influence of all masses in the lens plane on the light path is calculated. In the thin lens approximation, the deflection angle is just the sum of the deflection angles of every single lens. After deflection, all light rays are collected in pixel bins of the source plane. Thus a magnification pattern is produced (e.g. Fig. 4). Lengths in this map can be translated to angular separations or, assuming a constant relative velocity between source and lens, to time intervals. The number of collected light rays per pixel is proportional to the magnification of a background source with pixel size at the respective position in the source plane. The resolution in magnification of these numerically produced patterns is finite and depends on the total number of rays shot.
3.2 Light curve simulation
We obtain simulated microlensing light curves with potential traces of a moon by producing magnification maps (Fig. 4) of triplelens scenarios with mass ratios very different from unity. A light curve is then obtained as a onedimensional cut through the magnification pattern, convolved with the luminosity profile of a star. Only a few more assumptions are necessary to simulate realistic, in principle observable, light curves: angular source size, relative motion of lens and source, and lens mass. As a first approximation to the surface brightness profiles of stars we use a profile with radius and constant surface brightness (ignoring limbdarkening effects). The light curve is sampled at equidistant intervals. In Sect. 4.7 we discuss physical implications of the sampling frequency. In order to be able to make robust statistical statements about the detectability of a moon in a chosen setting, i.e. a certain magnification pattern, we analyse a grid of source trajectories that delivers about two hundred light curves, see also Fig. 5. To have an unbiased sample, the grid is chosen independently of lunar caustic features, but all trajectories are required to show pronounced deviations from the singlelens light curve due to the planet.
Figure 5: Illustration of how the magnification pattern of a given starplanetmoon configuration is evaluated. In order to have an unbiased statistical sample, the source trajectories are chosen independently of the lunar caustic features, though all light curves are required to pass through or close to the planetary caustic. This is realised through generating a grid of source trajectories that is only oriented at the planetary caustic. In the actual evaluation a 10 times denser grid than shown here was used. The magnification map has to be substantially larger than the caustic to ensure an adequate baseline for the light curve fitting. The side length of this pattern is , the position angle of the moon , all other parameters are chosen as listed in Table 1. Source sizes of 1 , 2 , 3 , and 5 are indicated in the lower left. 

Open with DEXTER 
3.3 Fitting
For each triplelens scenario, we produce an additional magnification pattern of a corresponding binarylens system, where the mass of the moon is added to the planetary mass. As a first approximation, one can compare two light curves with identical source track parameters, see Fig. 6a. Numerically big differences can occur without a significant topological difference, because the additional third body not only introduces additional caustic features, it also slightly moves the location of the doublelens caustics.
Figure 6: Triplelens light curve (solid line) extracted from the magnification pattern in Fig. 4, compared with a light curve extracted from the corresponding binarylens magnification map (dashed line), with no third body and the planet mass increased by the previous moon mass. The difference of the two light curves is plotted as the dotted line. The central peak of the triplelens light curve is caused by lunar caustic perturbations, that cannot easily be reproduced with a binary lens. In a), identical source trajectory parameters are used to extract the binarylens light curve and large residuals remain that can be avoided with a different source trajectory. In b), the bestfit binarylens light curve is shown (dashed) in comparison with the triplelens light curve (solid). Here, the residuals can be attributed to the moon. If the deviation is significant, the moon is detectable in the triplelens light curve. A source with solar radius at a distance of 8 kpc is assumed. 

Open with DEXTER 
3.4 Light curve comparison
We want to use the properties of the distribution for a test of significance of deviation between the two simulated light curves of the triple lens and the corresponding binary lens.
Most commonly, the distribution is used to test the goodnessoffit of a model to experimental data. To compare two simulated curves instead, as in our case, one possible approach is to randomly generate artificial data around one of them and then consecutively fit the two theoretical curves to the artificial data with some free parameters. Two ``values'' Q^{2}_{1} and Q^{2}_{2}will result. Their difference, can then be calculated and a threshold value is chosen. The condition for reported detectability of the deviation is then . A detailed description of the algorithm and its application to microlensing was presented by Gaudi & Sackett (2000).
We decided to use a different approach here: In particular, we were led to consider other possible techniques by the wish to avoid ``feeding'' our knowledge of the data distribution to a random number generator in order to get an unbiased and random distributed sample of Q^{2}, when at the same time we have all the necessary information to calculate a much more representative value. We have developed a method to quantify significance of deviations between two theoretical functions that avoids the steps of data simulation and subsequent fitting.
Put simply, we calculate the mean of all ``values'' that would result, if data with a Gaussian scatter drawn from a triplelens light curve is fitted with a binary light curve. If the triplelens case and the bestfit binarylens case are indistinguishable, then Q^{2} will be distributed and will lie at or very close to , the mean of the probability density function. If the triplelens light curve and the binarylens light curve differ significantly, then Q^{2} is not distributed and in particular will lie outside the expected range for distributed random variables. In the latter case, the moon is detectable. This approach is presented in mathematical detail in Appendix A.
4 Choice of scenarios
This section presents the assumptions that we use for our simulations. We discuss the astrophysical parameter space that is available for simulations of a microlensing system consisting of star, planet and moon. By choosing the most probable or most pragmatic value for each of the parameters, we create a standard scenario, summarised in Table 1 and visualised in Fig. 10, that all other parameters are compared against during the analysis.
The observational search for microlensing events caused by extrasolar planets is carried out towards the Galactic bulge, currently limited to the fields of the widefield surveys OGLE (Udalski et al. 1992) and MOA (Muraki et al. 1999). This leads to some natural assumptions for the involved quantities. The source stars are typically close to the centre of the Galactic bulge at a distance of kpc, where the surface density of stars is very high. For the distance to the lens plane we adopt a value of kpc. We assume our primary lens mass to be an Mdwarf star with a mass of , because that is the most abundant type, cf. Fig. 5 of Dominik et al. (2008). Direct lens mass determination is only possible if additional observables, such as parallax, can be measured, cf. Gould (2009). We derive the corresponding Einstein radius of ms.
Five parameters describe our lens configuration, cf. Fig. 1: they are the mass ratios between planet and star and between moon and planet , the angular separations between planet and star and between moon and planet, and as the last parameter , the position angle of the moon with respect to the planetstar axis. These parameters are barely constrained by the physics of a threebody system, even if we do require mass ratios very different from unity and separations that allow for stable orbits.
The apparent source size and the sampling rate affect the shape of the simulated light curves. We have to define an observational uncertainty for the significance test in the final analysis of the light curves.
4.1 Mass ratio of planet and star q
We are interested in planets (and not binary stars). Accordingly, we want a small value for the mass ratio of planet and star . Our standard value will be 10^{3} which is the mass ratio of Jupiter and Sun, or a Saturnmass planet around an Mdwarf of . At a given projected separation between star and planet, this mass ratio determines the size of the planetary caustic, with a larger leading to a larger caustic. is varied to also examine scenarios with mass ratios corresponding to a Jovian mass around an Mdwarf and to a Saturn mass around the Sun. Summarised: for the mass ratio between planet and star we use the three values .
4.2 Mass ratio of moon and planet q
To be classified as a moon, the tertiary body must have a mass considerably smaller than the secondary. The standard case in our examination corresponds to the Moon/Earth mass ratio of . We are generous towards the higher mass end, and include a mass ratio of 10^{1} in our analysis, corresponding to the Charon/Pluto system. Both examples are singular in the solar system, but we argue that a more massive moon is more interesting, since it can effectively stabilise the obliquity of a planet, which is thought to be favouring the habitability of the planet (Benn 2001). At the low mass end of our analysis we examine . In Fig. 7, three different caustic interferences resulting from the three adopted mass ratios are shown. Summarised: For the mass ratio between moon and planet we use the three values .
4.3 Angular separation of planet and star
The physical separation of a planet and its host star, i.e. the semimajor axis , is not directly measurable in gravitational lensing. Only the angular separation can directly be inferred from an observed light curve, where is valid for a given physical separation at a lens distance . The angular separation of a binary will evoke a certain topology of caustics, illustrated e.g. in Fig. 1 of Cassan (2008). They gradually evolve from the close separation case with two small triangular caustics on the far side of the star and a small central caustic at the star position, to a large central caustic for the intermediate case, when the planet is situated near the stellar Einstein ring, . This is also called the ``resonant lensing'' case, as described in Wambsganss (1997). If the planet is moved further out, one obtains a small central caustic and a larger isolated, roughly kiteshaped planetary caustic that is elongated towards the primary mass in the beginning, but becomes more and more symmetric and diamondshaped if the planet is placed further outwards.
Figure 7: Caustic topology with different lunar mass ratios. From left to right, the lunar mass ratio increases from to 10^{1}. The side length of the displayed patterns is . All other parameters as given in Table 1, in particular and . Source sizes of to are indicated in the lower left. 

Open with DEXTER 
In analogy, a ``lunar resonant lensing zone'' can be defined, when measuring the angular separation of moon and planet in units of planetary Einstein radii . Viewing our solar system as a gravitational lens system at kpc and with kpc, Jupiter's Einstein radius in physical lengths would be 0.11 AU in the lens plane, that is almost 10 times the semimajor axis of Callisto. However, the smaller moons of Jupiter cover the region out to 2 Jovian Einstein radii. If there is a moon at an angular separation of the planet of , it will show its influence in any of the caustics topologies. But we focus our study on the wideseparation case for the following reasons: regarding the closeseparation triangular caustics, we argue that the probability to cross one of them is vanishingly small and, furthermore, the magnification substantially decreases as they move outwards from the star. The intermediate or resonant caustic is well ``visible'' because it is always located close to the peak of the singlelens curve. But all massive bodies of a given planetary system affect the central caustic with minor or major perturbations and deformations (Gaudi et al. 1998). There is no reason to expect that extrasolar systems are generally less densely populated by planets or moons than the solar system. Since we are looking for very small deviations, identifying the signature of the moon among multiple features caused by lowmass planets and possibly more moons would be increasingly complex. The wideseparation case is most favourable, because the planetary caustic is typically caused by a single planet and interactions due to the close presence of other planets are highly unlikely (Bozza 2000,1999). Any observed interferences can therefore be attributed to satellites of the planet. Planets with only one dominant moon are at least common in the solar system (Saturn & Titan, Neptune & Triton, Earth & Moon). Such a system would be the most straightforward to detect in real data.
This makes a strong case for concentrating on the wideseparation planetary caustic. As the standard we choose a separation of and also test a scenario with . For comparison: the maximum projected separation of Jupiter at the chosen distances and corresponds to 1.5 solar Einstein radii. Summarised: for the angular separation between planet and star we use the two values .
4.4 Angular separation of moon and planet
As the moon by definition orbits the planet, there is an upper
limit
to the semimajor axis of the moon .
It must not exceed the
distance between the planet and its inner Lagrange point. This
distance is called Hill radius and is calculated as
for circular orbits, with the semimajor axis of the planetary orbit. This can be translated to ``lensing parameters'', but the equation changes into the inequality
because for a physical distance , the apparent (i.e. projected) separation lies in the range depending on the inclination towards the line of sight. Thus, the Hill radius constraint cannot be regarded as a strict limit. Though, turning this argument around, a minimum region of secured stability exists for given mass ratios and angular planetstar separation . Bodies in prograde motion with an orbit below can have long term stability, for retrograde motion the limit is somewhat higher at (see Domingos et al. 2006, and references therein, particularly; Hunter 1967). The region of secured stability for our standard scenario has a radius of
which leads to AU in our setting. Our choice of is also motivated by a desire to have caustic interactions between lunar and planetary caustic, a resonance, because only under this condition the moon will give rise to light curve features that are distinctly different from those of a lowmass planet. If the moon is located at an angular separation from the planet of  depending on the position angle of the moon  roughly , it will influence the planetary caustic. For angular separations larger than , the lunar caustic barely influences the planetary caustic, and if it is detected, it will be detected as a distinct secondary caustic. Han (2008) already pointed out that the lunar signal does not vanish with an increased separation from the planet. But the lunar mass might be wrongly identified as a second, lowermass planet rather than a satellite of the first planet. For comparison, the major moons of the solar system are all at considerably smaller angular separations than 1.0 , if observed as a gravitational lens system from kpc with kpc, with our Moon being most favourable at 0.4 of Earth's Einstein radii at maximum projected separation. But as we have learnt through the discovery of hundreds of exoplanets, we do not have to expect solar system properties to be mirrored in exoplanetary systems.
We focus on configurations with strong interferences of planetary and lunar caustic, what one might call the ``resonant case of lunar lensing'', examples of which are displayed in Fig. 8. Therefore, our standard case will have an angular separation of moon and planet of one planetary Einstein radius. In total, we evaluate four settings: , in the chosen setting (Table 1), this corresponds to physical projected separations from 0.5 to 0.8 AU.
Figure 8: Magnification patterns with side lengths of showing the lunar perturbation of the planetary caustic for increasing angular separations of moon and planet of , while all other parameters are as in Table 1. The topology of interaction varies with the angular separation. Source sizes of to are indicated in the lower left. 

Open with DEXTER 
4.5 Position angle of moon
The fifth physical parameter for determining the magnification patterns is the position angle of the moon with respect to planetstar axis, as depicted in Fig. 1. It is the only parameter that is not fixed in our standard scenario. Instead, we evaluate this parameter in its entirety by varying it in steps of to complete a full circular orbit of the moon around the planet. By doing this, we are aiming at getting complete coverage of a selected mass/separation scenario. Figure 2 visualises how the position of the moon affects the planetary caustic. It is not to be expected that we will ever have an exactly frontal view of a perfectly circular orbit, but it serves well as a first approximation. Summarised: we choose 12 equally spaced values for the position angle of the moon relative to the planetstar axis, and the results for the various parameter sets are averaged over these 12 geometries.
4.6 Source size R
The source size strongly influences the amplitudes of the light curve fluctuations and the ``time resolution''. Larger sources blur out finer caustic structures, compare Fig. 9. For our source star assumptions, we have to take into account not only stellar abundances, but also the luminosity of a given stellar type. Main sequence dwarfs are abundant in the Galactic bulge, but they are faint, with apparent magnitudes of 15 mag to 25 mag at kpc. In the crowded microlensing fields, they are difficult to observe with satisfying photometry by groundbased telescopes, even if they are lensed and magnified. Giant stars, with apparent baseline magnitudes 13 mag to 17 mag, are more likely to allow precise photometric measurements from ground and, therefore, today's followup observations are biased towards these source stars.
The finite source size constitutes a serious limitation to the discovery of extrasolar moons. In fact, Han & Han (2002) stated that detecting satellite signals in the lensing light curves will be close to impossible, because the signals are smeared out by the severe finitesource effect. They tested various source sizes (and planetmoon separations) for an Earthmass planet with a Moonmass satellite. They find that even for a K0type source star ( ) any light curve modifications caused by the moon are washed out. Han (2008) increased lens masses and by assuming a solar source size , he finds that ``nonnegligible satellite signals occur'' in the light curves of planets of 10 to 300 Earth masses, ``when the planetmoon separation is similar to or greater than the Einstein radius of the planet'' and the moon has the mass of Earth. We use a solar sized source as our standard case, and present results for four more source radii: . In terms of stellar Einstein radii this corresponds to five different source size settings from to .
The ability to monitor smaller sources  dwarfs, lowmass main sequence stars  with good photometric accuracy is one of the advantages of spacebased observations. Hence a satellite mission is the ideal tool for this aspect of microlensing, routine detections of moons around planets can be expected with a spacebased monitoring program on a dedicated satellite.
Figure 9: Illustration of the increasing source radius, compare Fig. 6b: a) , b) , and c) , all other parameters as given in Table 1. The triplelens light curve (solid) is obtained from the magnification pattern in Fig. 4, the bestfit binarylens light curve (dashed) from the corresponding magnification pattern without a moon. The difference between the two curves is plotted as a dotted line. 

Open with DEXTER 
4.7 Sampling rate f
Typical exposure times for the microlensing light curves are 10 to 300 s at a 1.5 m telescope. The frequency of observations of an individual interesting event varies between once per 2 h to once every few minutes, for a single observing site only limited by exposure time and readout time of the instrument. One example of this is the peak coverage of MOA2007BLG400, see Dong et al. (2009), where frames were taken every six minutes. Higher observing frequency equals better coverage of the resulting light curve. A normal microlensing event is seen as a transient brightening that can last for a few days to weeks or even months. A planet will alter the light curve from fraction of a day to about a week; the effect of a lowmass planet or a moon will only last for a few hours or shorter, see e.g. Beaulieu et al. (2006). This duration is inversely proportional to the transverse velocity with which the lens moves relative to the line of sight or the relative proper motion between the background source and the lensing system. We fix this at a typical km s^{1} for kpc or .
For the simulations, a realistic, noncontinuous observing rate is mimicked by evaluating the simulated light curve at intervals that correspond to typical frequencies, but allowing for a constant sampling rate, which in practise is often prohibited by observing conditions. We sample the triplelens light curve in equally sized steps. With the assumed constant relative velocity, this translates to equal spacing in time. Aiming for a rate of , we choose a step size of . To see the effect of a decreased sampling rate, we have also examined the standard scenario (see Table 1) with sampling rates of factors of 1.5 and 2 longer. Summarised: we use three different sampling rates, min, 1/22 min, 1/30 min.
4.8 Assumed photometric uncertainty
The standard error of a data point in a microlensing light curve can vary significantly. Values between 5 and 100 mmag seem realistic, see, e.g., data plots in the recent event analyses of Gaudi et al. (2008) or Janczak et al. (2009). The photometric uncertainty directly enters the statistical evaluation of each light curve as described in Appendix A. We have drawn results for an (unrealistically) broad range from 0.5 mmag to 500 mmag. An error as small as mmag for a V=12.3star has recently been reached with highprecision photometry of an exoplanetary transit event at the Danish 1.54 m telescope at ESO La Silla (Southworth et al. 2009), which is one of the telescopes presently used for Galactic microlensing monitoring observations, but this low uncertainty cannot be transferred to standard microlensing observations, where crowded star fields make e.g. the use of defocussing techniques impossible. In our discussion, we regard only the more realistic range from to 100 mmag. We adopt an ideal scenario with a fixed value for all sampled points on the triplelens light curve and set mmag to be our standard for the comparison of different mass and separation scenarios. Summarised: we explore 5 values as photometric uncertainty: , 10, 20, 50, 100 mmag.
Table 1: Parameter values of our standard scenario.
Figure 10: Visualisation of the approximate lunar resonance zone (darkgrey ring), where planetary and lunar caustic overlap and interact. Also depicted is the planet lensing zone (light grey), which is the area of planet positions for which the planetary caustic(s) lie within the stellar Einstein radius. Distances and caustic inset are to scale. The physical scale depends not only on the masses of the lensing bodies, but also on the distances to the lens system and the background star. In our adopted standard scenario, cf. Table 1, the Einstein radius corresponds to 1.9 AU in the lens plane. This places the planet at a projected orbit radius of 2.5 AU, whereas the moon's orbit has a radius of 0.06 AU (corresponding to one planetary Einstein radius). The lunar resonance zone covers lunar orbit radii from 0.04 AU to 0.18 AU. 

Open with DEXTER 
5 Results: detectability of extrasolar moons in microlensing light curves
The numerical results of our work are presented and a first interpretation is made. We start by considering the result for a single magnification pattern in Sect. 5.1. We then present results that cover specific physical triplelens scenarios (Sect. 5.2). We also evaluate the parameter dependence of the detectability rate, by varying each parameter separately.
5.1 Evaluation of a single magnification pattern
As described in Sect. 3, all magnification maps are evaluated by taking a large sample of unbiased^{} light curves and then determining the fraction of light curves with significant lunar signals. The results for a single magnification map (here the pattern from Fig. 5), are displayed in Table 2.
Table 2: Detectability of the third mass in the magnification pattern of Fig. 5, for different assumed photometric uncertainties.
The percentages can be interpreted as the ``detection probability'' of the moon in a random observed light curve displaying a signature of the planetary caustic  ignoring the possible complication to fully characterise the signal. In this map, the detectability of the moon is about one third, provided we have an observational uncertainty of mmag. Obviously, with a larger photometric error, one is less sensitive to small deviations caused by the moon.5.2 Results for selected scenarios
In order to be able to make general statements about a given planetary/lunar system, we average the results for 12 magnification patterns representing a full lunar orbit with .
5.2.1 Variation of the photometric uncertainty
The scenario described in Table 1 is evaluated for different photometric uncertainties by analysing the 12 magnification maps shown in Fig. 2. The results are displayed in Table 3.
Table 3: Detectability of the moon in our standard scenario as a function of assumed photometric uncertainty.
For an photometric uncertainty of mmag, the result is that about 30% of the light curves show significant deviations from a bestfit binarylens light curve and display detectable signs of the moon. We recall again our requirements: the host planet of the satellite has been detected;
 a small source size, up to a few solar radii, is required;
 the moon is massive compared to satellites in the solar system;
 we have assumed light curves with a constant sampling rate of one frame every 15 min for about 50 to 70 h.
5.2.2 Changing the angular separation of planet and moon
Table 4 lines up the changing detectability for different values of the projected planetmoon separation , with all other parameters as in standard scenario (Table 1).
Table 4: Detectability of the moon as a function of the planetmoon separation for our standard scenario.
With the chosen range from to , we are within the regime of caustic interactions of planetary and lunar caustic. We find that the percentage of detectability is almost constant at a 30% level. Han (2008) defines regions of possible satellite detections as the region of angular moonplanet separations between the lower limit of one planetary Einstein radius and an upper limit at the projected Hillradius average. From our results, we conclude that detections with are well possible. We have not fully probed the lower limit region, but for the size of the lunar caustic perturbation decreases substantially.5.2.3 Variation of the moon mass
Table 5 shows the changing detectability for a varying mass ratio of lunar and planetary mass, with all other parameters as in our standard scenario (Table 1).
Table 5: Detectability of the moon depending on moonplanet mass ratio.
Over the range of the moonplanet mass ratio, , 10^{2}, and 10^{1} the detection rate increases from 2.4% to 85.6%, compare also Fig. 7. Not surprisingly: the more massive the moon, the easier it is to detect.5.2.4 Changing the planetary mass ratio
Table 6 lists the changing detectability for a varying mass ratio between planet and star .
Table 6: Dependence of lunar detectability as a function of the planetstar mass ratio, all other parameters as in our standard scenario.
We adjusted the parameters of the grid of source trajectories, so that the grid spacing scales with the total caustic size. For a decreasing mass of the planet, the apparent source size increases relative to the caustic size, so finer features will be blurred out more in the case of a smaller planet. Also the absolute change in magnification scales with the planetary mass. As expected, the detection efficiency is roughly proportional to the planetary caustic size.5.2.5 Different sampling rates
Table 7 presents the changing lunar detectability for different sampling rates.
Table 7: Detectability of the moon depending on the sampling rate of the simulated light curves.
A lower sampling frequency lowers the detection probability, but the effect is less pronounced than expected. This means it may be favourable to monitor a larger number of planetary microlensing events with high sampling frequency, rather than a very small number of them with ultrahigh sampling. Our assumed sampling rates can easily be met by followup observations, as they are presently performed for anomalous Galactic microlensing events.5.2.6 Source size variations
We analysed simulated light curves of the standard scenario with five different source star radii. They cover the range between a star the size of our Sun and a small giant star with at kpc.
Table 8: Detectability of the moon depending on the size of the source star, all other parameters as in our standard scenario.
A larger source blurs out all sharp causticcrossing features of a light curve, compare Fig. 9. We see from our results in Table 8 that the discovery of a moderately massive moon (see standard scenario assumptions in Table 1) is impossible for a source star size of or larger within our assumptions.
5.2.7 Increased separation between star and planet
For reasons given above in Sect. 4.3, we did not vary the separation between star and planet significantly, but stayed in the outer ranges of the planet lensing zone, corresponding to a welldetectable, wideseparation planetary caustic. For , we get a detectability of 33.3%, which is similar to our standard case ( ) with 30.6%.
6 Conclusion and outlook
In this work we provide probabilities for the detection of moons around extrasolar planets with gravitational lensing. We showed that massive extrasolar moons can principally be detected and identified via the technique of Galactic microlensing. From our results it can be concluded that the detection of an extrasolar moon  under favourable conditions  is within close reach of available observing technology. The unambiguous characterisation of observed lunar features however will be challenging. The examined lens scenarios model realistic triplelens configurations. An observing rate of about one frame per 15 min is desirable, which is high, but well within the range of what is now regularly performed in microlensing followup observations of anomalous events. Similarly, an observational uncertainty of about 20 mmag can be met with today's groundbased telescopes and photometric reduction techniques for sufficiently bright targets. Bright microlensing targets are mostly giant stars, which is an impediment to the detection of moons: Bulge giants with radii of order 10 or larger smooth out any lunar caustic feature. Therefore, in order to find moons, it is crucial to be able to perform precise photometry on small sources with angular sizes of the order of 10^{3} Einstein radii, corresponding to a few solar radii or smaller at a distance of 8 kpc. This means dwarf stars rather than giants need to be monitored in order to identify moons in the intervening planetary systems. Some improvement in resolution can be gained with the luckyimaging technique on mediumsized ground telescopes (Grundahl et al. 2009). Under very favourable circumstances, exomoons might already be detectable from ground. However, future spacebased telescopes, such as ESA's proposed Euclid^{} mission (cf. Réfrégier et al. 2010, Chapter 17) or the dedicated Microlensing Planet Finder (MPF) mission proposed to NASA (cf. Bennett et al. 2009), will surely have the potential to find extrasolar moons through gravitational lensing.
AcknowledgementsWe thank our referee, JeanPhilippe Beaulieu, for his constructive comments, which helped to improve the manuscript. C.L. would like to thank Sven Marnach for his helpful input on the significance test. This research has made use of NASA's Astrophysics Data System and the arXiv eprint service operated by Cornell University.
Appendix A Significance of deviation
This section explains the method that we use to determine the significance of deviation between the two simulated light curves of the triple lens and the bestfit binary lens model.
When we want to decide whether real observational data is better described with a triplelens or with a binarylens model, we can fit both models to the data and use the test to assess the significance of the deviation between the data and the models and pick the model with the better fit.
We do not have observational data, but numerically computed light curves. If we observed one of the simulated triplelens systems, we would expect the data points to have a Gaussian distribution around the theoretical values. We could produce artificial data by randomly scattering the simulated triplelens light curve data points around their theoretical values and then fit a binarylens model and the triplelens model to this artificial data set, as described in Sect. 3.4. This approach allows us to use the usual test, but is computationally intensive, because the models have to be fitted numerous times to get a statistically sound sample of values. In our approach, we directly calculate the value which is the expectation value of the above method.
We hypothesise, every simulated point on our triplelens light curve is in fact the mean of the distribution of a random variable X_{i} that is distributed according to the Gaussian probability density function f_{i}(x_{i}) with a standard deviation of . is the variance of the distribution^{}.
We now want to examine whether the X_{i}
could be described equally
well with a binarylens model .
We introduce the new
distributed^{} random variable
At this step, we could simulate data X_{i} in order to find a somewhat representative, randomly drawn value of Q^{2}, but instead we simply calculate what the mean value of all possible Q^{2} would be. We use the definitions, to find
Here, we keep in mind that, in general, . Using the parameters and of the distribution f_{i}(x_{i}), we reduce the equation by calculating
=  
=  
= 
We used the definitions of and and the property of the probability density function . In our case holds true without loss of accuracy, and indeed, we simplify for all i, that is we assume that the photometric uncertainty is the same for all data points, and argue that this does not pose a problem as long as is chosen to match the maximum photometric uncertainty. So we can further reduce to
This is the mean of all Q^{2} possibly resulting, when comparing the simulated binarylens light curve to randomly scattered triplelens light curve points. It shall be our measure of deviation. Fortunately, we can provide all remaining parameters:
 n
 is the number of degrees of freedom, equal to the number of compared data points.
 is the assumed standard deviation or photometric uncertainty of observations.
 equals the difference between the two compared light curves squared, which we already use for our least square fit.
Figure A.1: The probability density function plotted for illustration only. For increasing degrees of freedom the central limit theorem takes effect and the curve will very closely resemble a Gaussian distribution. If , i.e. if the probability P for a distributed random variable X to be larger or as large as the mean of Q^{2} is less than 0.01, we consider the deviation between the two compared curves to be significant. 

Open with DEXTER 
This probability is just the integral of the probability density function to the right of as illustrated in Fig. A.1. If this probability is very small, is outside the expected range for a distributed random variable Q^{2}. In that case, Q^{2} is obviously not distributed with n degrees of freedom, so we must conclude that there is a significant deviation between the triple lens model light curve and the binarylens model light curve.
We say we have a significant deviation between our two curves, if the mean value is so high that the probability Pfor any random variable X being larger or equally large is less than 1%. We interpret this to say, only if the probability for a given triple lens light curve with independent, normally distributed data points to be random fluctuation of the compared binarylens light curve is less than 1%, we consider it to be principally detectable.
There are two known sources for overestimating the detection rate coming with our method. First, systematic errors are not accounted for. A possible solution could be to add a further term to as in , where can be assumed to lie in the few percent region for real data. Secondly, we avoid refitting the binarylens model light curve and always compare to the one we got as bestfit to the exact (simulated) triplelens model light curve. This overestimation effect vanishes for and is negligible for a sufficiently large n. The high sampling frequency we use, ensures that the number of data points is always large enough (n > 250 in all cases).
While we cannot exactly quantify the uncertainty of the results, we know that they pose strict upper limits for the lunar detectability in the various scenarios. This knowledge would enable us to infer a tentative census of extrasolar moons, once the first successful microlensing detections have been achieved.
References
 Beaulieu, J.P., Albrow, M., Bennett, D. P., et al. 2006, Nature, 439, 437 [NASA ADS] [CrossRef] [PubMed] (In the text)
 Beaulieu, J.P., Kerins, E., Mao, S., et al. 2008, ESA white paper, submitted [arXiv:0808.0005v1] (In the text)
 Benn, C. R. 2001, Earth, Moon and Planets, 85, 61 (In the text)
 Bennett, D. P., & Rhie, S. H. 2002, ApJ, 574, 985 [NASA ADS] [CrossRef] (In the text)
 Bennett, D. P., Anderson, J., Beaulieu, J.P., et al. 2009, white paper, submitted [arXiv:0902.3000v1] (In the text)
 Bond, I. A., Udalski, A., Jaroszynski, M., et al. 2004, ApJ, 606, L155 [NASA ADS] [CrossRef] [MathSciNet] (In the text)
 Bozza, V. 1999, A&A, 348, 311 [NASA ADS] (In the text)
 Bozza, V. 2000, A&A, 355, 423 [NASA ADS] (In the text)
 Cabrera, J., & Schneider, J. 2007, A&A, 464, 1133 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Cassan, A. 2008, A&A, 491, 587 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Domingos, R. C., Winter, O. C., & Yokoyama, T. 2006, MNRAS, 373, 1227 [NASA ADS] [CrossRef] (In the text)
 Dominik, M., Jørgensen, U. G., Horne, K., et al. 2008, ESA white paper, submitted [arXiv:0808.0004v1] (In the text)
 Dong, S., Bond, I. A., Gould, A. P., et al. 2009, ApJ, 698, 1826 [NASA ADS] [CrossRef] (In the text)
 Dyson, F. W., Eddington, A. S., & Davidson, C. 1920, R. Soc. London Phil. Trans. Ser. A, 220, 291 [NASA ADS] [CrossRef] (In the text)
 Einstein, A. 1916, Ann. Phys., 354, 769 [NASA ADS] [CrossRef] (In the text)
 Gaudi, B. S., & Sackett, P. D. 2000, ApJ, 528, 56 [NASA ADS] [CrossRef] (In the text)
 Gaudi, B. S., Naber, R. M., & Sackett, P. D. 1998, ApJ, 502, L33 [NASA ADS] [CrossRef] (In the text)
 Gaudi, B. S., Bennett, D. P., Udalski, A., Gould, A., et al. 2008, Science, 319, 927 [NASA ADS] [CrossRef] [PubMed] (In the text)
 Gaudi, B. S., Beaulieu, J.P., Bennett, D. P., et al. 2009, white paper, [arXiv:0903.0880v1] (In the text)
 Gould, A. P. 2009, in ASP Conf. Ser., ed. K. Z. Stanek, 403, 86 (In the text)
 Grundahl, F., ChristensenDalsgaard, J., Kjeldsen, H., et al. 2009, The Stellar Observations Network Group  the Prototype [arXiv:0908.0436v1] (In the text)
 Han, C. 2008, ApJ, 684, 884 [NASA ADS] [CrossRef] (In the text)
 Han, C., & Han, W. 2002, ApJ, 580, 490 [NASA ADS] [CrossRef] (In the text)
 Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288 [NASA ADS] [CrossRef] [PubMed] (In the text)
 Hunter, R. B. 1967, MNRAS, 136, 245 [NASA ADS] (In the text)
 Janczak, J., Fukui, A., Dong, S., et al. 2010, ApJ, 711, 731 [NASA ADS] [CrossRef] (In the text)
 Kayser, R., Refsdal, S., & Stabell, R. 1986, A&A, 166, 36 [NASA ADS] (In the text)
 Kipping, D. M. 2009a, MNRAS, 392, 181 [NASA ADS] [CrossRef] (In the text)
 Kipping, D. M. 2009b, MNRAS, 396, 1797 [NASA ADS] [CrossRef] (In the text)
 Kipping, D. M., Fossey, S. J., & Campanella, G. 2009, MNRAS, 400, 398 [NASA ADS] [CrossRef] (In the text)
 Lewis, K. M., Sackett, P. D., & Mardling, R. A. 2008, ApJ, 685, L153 [NASA ADS] [CrossRef] (In the text)
 Liebig, C. 2009, Diploma thesis, Universität Heidelberg (In the text)
 Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [NASA ADS] [CrossRef] (In the text)
 Muraki, Y., Sumi, T., Abe, F., et al. 1999, Progress of Theoretical Physics Supplement, 133, 233 [NASA ADS] [CrossRef] (In the text)
 Paczynski, B. 1996, Ann. Rev. A&A, 34, 419 [NASA ADS] [CrossRef] (In the text)
 Réfrégier, A., Amara, A., Kitching, T. D., et al. 2010, Euclid Imaging Consortium Science Book [arXiv:1001.0061v1] (In the text)
 Rhie, S. H. 1997, ApJ, 484, 63 [NASA ADS] [CrossRef] (In the text)
 Rhie, S. H. 2002, How Cumbersome is a Tenth Order Polynomial? The Case of Gravitational Triple Lens Equation [arXiv:astroph/0202294v1] (In the text)
 Sartoretti, P., & Schneider, J. 1999, A&AS, 134, 553 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] (In the text)
 Scharf, C. A. 2006, ApJ, 648, 1196 [NASA ADS] [CrossRef] (In the text)
 Schneider, P., & Weiß, A. 1986, A&A, 164, 237 [NASA ADS] (In the text)
 Schneider, P., Kochanek, C. S., & Wambsganss, J. 2006, Gravitational Lensing: Strong, Weak and Micro, SaasFee Advanced Courses 33 (Berlin, Heidelberg: SpringerVerlag) (In the text)
 Southworth, J., Hinse, T. C., Jørgensen, U. G., et al. 2009, MNRAS, 396, 1023 [NASA ADS] [CrossRef] (In the text)
 Udalski, A., Szymanski, M., Ka▯uzny, J., Kubiak, M., & Mateo, M. 1992, Acta Astron., 42, 253 [NASA ADS] (In the text)
 Wambsganss, J. 1990, Ph.D. Thesis, LudwigMaximiliansUniversität München (In the text)
 Wambsganss, J. 1997, MNRAS, 284, 172 [NASA ADS] (In the text)
 Wambsganss, J. 1999, J. Comput. Appl. Math., 109, 353 [NASA ADS] [CrossRef] (In the text)
 Williams, D. M., & Knacke, R. F. 2004, Astrobiology, 4, 400 [NASA ADS] [CrossRef] [PubMed] (In the text)
 Witt, H. J. 1990, A&A, 236, 311 [NASA ADS] (In the text)
Footnotes
 ... detected^{}
 exoplanet.eu
 ... unbiased^{}
 I.e. independent of lunar features, but required to pass close to or through the planetary caustic.
 ... Euclid^{}
 http://sci.esa.int/euclid
 ... distribution^{}
 We recall the definitions of mean and
variance. For the random variable X with probability density
function f(x), they are in general
= =
where we also introduce the notation for the expectation value of X.  ...distributed^{}
 The distribution can be applied
to the sum of n independent, normally distributed random
variables Z_{i} with the mean
and the
standard deviation
for all i,
See any standard book with an introduction to calculus of probability.  ... Library^{}
 http://www.gnu.org/software/gsl/
All Tables
Table 1: Parameter values of our standard scenario.
Table 2: Detectability of the third mass in the magnification pattern of Fig. 5, for different assumed photometric uncertainties.
Table 3: Detectability of the moon in our standard scenario as a function of assumed photometric uncertainty.
Table 4: Detectability of the moon as a function of the planetmoon separation for our standard scenario.
Table 5: Detectability of the moon depending on moonplanet mass ratio.
Table 6: Dependence of lunar detectability as a function of the planetstar mass ratio, all other parameters as in our standard scenario.
Table 7: Detectability of the moon depending on the sampling rate of the simulated light curves.
Table 8: Detectability of the moon depending on the size of the source star, all other parameters as in our standard scenario.
All Figures
Figure 1: Geometry of our triplelens scenario, not to scale. Five parameters have to be fixed: the mass ratios and , the angular separations in the lens plane and and the position angle of the moon . These parameters define uniquely the relative projected positions of the three bodies. 

Open with DEXTER  
In the text 
Figure 2: Details of the analysed triplelens magnification maps for an example mass/separation scenario with , , and (the standard scenario as summarised in Table 1), showing the caustic configurations for twelve different lunar positions completing a full circular orbit in steps of 30 . The relative positions of star, planet and moon are sketched in the lower left of each panel (not to scale). The side lengths of the individual frames are . A darker shade of grey corresponds to a higher magnification. The straight black lines mark source trajectories, the resulting light curves are displayed in Fig. 3. 

Open with DEXTER  
In the text 
Figure 3: Sample of triplelens light curves corresponding to the source trajectories depicted in Fig. 2. The solid light curves were extracted with an assumed solar source size, ; the thin, dashed lines show the same light curves with . The magnification scale is given in negative magnitude difference (mag), i.e. the unmagnified baseline flux of the source is 0. The time scale in units of the Einstein time t_{E}, assuming a uniform relative motion of source and lens. 

Open with DEXTER  
In the text 
Figure 4: Magnification pattern displaying the planetary caustic of a triplelens system with mass ratios and , corresponding, e.g., to an Mdwarf, a Saturnmass planet and an Earthmass satellite of that planet. Darker shades correspond to areas of higher magnification. The third body shows its influence as the small perturbation on the bottom cusp of the planetary caustic. The position angle of the moon is , i.e. in the orientation of this map the moon is located above the planet. The complete set of parameters is given in Table 1. The side length of the pattern is in stellar Einstein radii. A source with radius (2 , 3 , 5 ) is indicated as a black disc (dashed circles) in the lower left. Taking the straight line that cuts the lunar perturbation of the caustic as the source trajectory of a solar sized source results in the light curve shown in Fig. 6 (solid line). Depending on the assumed source size, different light curves are obtained, see also Figs. 9a to c. 

Open with DEXTER  
In the text 
Figure 5: Illustration of how the magnification pattern of a given starplanetmoon configuration is evaluated. In order to have an unbiased statistical sample, the source trajectories are chosen independently of the lunar caustic features, though all light curves are required to pass through or close to the planetary caustic. This is realised through generating a grid of source trajectories that is only oriented at the planetary caustic. In the actual evaluation a 10 times denser grid than shown here was used. The magnification map has to be substantially larger than the caustic to ensure an adequate baseline for the light curve fitting. The side length of this pattern is , the position angle of the moon , all other parameters are chosen as listed in Table 1. Source sizes of 1 , 2 , 3 , and 5 are indicated in the lower left. 

Open with DEXTER  
In the text 
Figure 6: Triplelens light curve (solid line) extracted from the magnification pattern in Fig. 4, compared with a light curve extracted from the corresponding binarylens magnification map (dashed line), with no third body and the planet mass increased by the previous moon mass. The difference of the two light curves is plotted as the dotted line. The central peak of the triplelens light curve is caused by lunar caustic perturbations, that cannot easily be reproduced with a binary lens. In a), identical source trajectory parameters are used to extract the binarylens light curve and large residuals remain that can be avoided with a different source trajectory. In b), the bestfit binarylens light curve is shown (dashed) in comparison with the triplelens light curve (solid). Here, the residuals can be attributed to the moon. If the deviation is significant, the moon is detectable in the triplelens light curve. A source with solar radius at a distance of 8 kpc is assumed. 

Open with DEXTER  
In the text 
Figure 7: Caustic topology with different lunar mass ratios. From left to right, the lunar mass ratio increases from to 10^{1}. The side length of the displayed patterns is . All other parameters as given in Table 1, in particular and . Source sizes of to are indicated in the lower left. 

Open with DEXTER  
In the text 
Figure 8: Magnification patterns with side lengths of showing the lunar perturbation of the planetary caustic for increasing angular separations of moon and planet of , while all other parameters are as in Table 1. The topology of interaction varies with the angular separation. Source sizes of to are indicated in the lower left. 

Open with DEXTER  
In the text 
Figure 9: Illustration of the increasing source radius, compare Fig. 6b: a) , b) , and c) , all other parameters as given in Table 1. The triplelens light curve (solid) is obtained from the magnification pattern in Fig. 4, the bestfit binarylens light curve (dashed) from the corresponding magnification pattern without a moon. The difference between the two curves is plotted as a dotted line. 

Open with DEXTER  
In the text 
Figure 10: Visualisation of the approximate lunar resonance zone (darkgrey ring), where planetary and lunar caustic overlap and interact. Also depicted is the planet lensing zone (light grey), which is the area of planet positions for which the planetary caustic(s) lie within the stellar Einstein radius. Distances and caustic inset are to scale. The physical scale depends not only on the masses of the lensing bodies, but also on the distances to the lens system and the background star. In our adopted standard scenario, cf. Table 1, the Einstein radius corresponds to 1.9 AU in the lens plane. This places the planet at a projected orbit radius of 2.5 AU, whereas the moon's orbit has a radius of 0.06 AU (corresponding to one planetary Einstein radius). The lunar resonance zone covers lunar orbit radii from 0.04 AU to 0.18 AU. 

Open with DEXTER  
In the text 
Figure A.1: The probability density function plotted for illustration only. For increasing degrees of freedom the central limit theorem takes effect and the curve will very closely resemble a Gaussian distribution. If , i.e. if the probability P for a distributed random variable X to be larger or as large as the mean of Q^{2} is less than 0.01, we consider the deviation between the two compared curves to be significant. 

Open with DEXTER  
In the text 
Copyright ESO 2010