Issue 
A&A
Volume 543, July 2012



Article Number  A63  
Number of page(s)  12  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201118140  
Published online  27 June 2012 
Minimizing followup for spacebased transit surveys using full lightcurve analysis^{⋆}
Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
email: nefs@strw.leidenuniv.nl
Received: 23 September 2011
Accepted: 8 May 2012
Context. One of the biggest challenges facing large transit surveys is the elimination of falsepositives from the vast number of transit candidates. A large amount of expensive followup time is spent on verifying the nature of these systems.
Aims. We investigate to what extent information from the lightcurves can identify blend scenarios and eliminate them as planet candidates, to significantly decrease the amount of followup observing time required to identify the true exoplanet systems.
Methods. If a lightcurve has a sufficiently high signaltonoise ratio, a distinction can be made between the lightcurve of a stellar binary blended with a third star and the lightcurve of a transiting exoplanet system. We first simulate lightcurves of stellar blends and transiting planet systems to determine what signaltonoise level is required to make the distinction between blended and nonblended systems as function of transit depth and impact parameter. Subsequently we test our method on real data from the first IRa01 field observed by the CoRoT satellite, concentrating on the 51 candidates already identified by the CoRoT team.
Results. Our simulations show that blend scenarios can be constrained for transiting systems at low impact parameters. At high impact parameter, blended and nonblended systems are indistinguishable from each other because they both produce Vshaped transits. About 70% of the planet candidates in the CoRoT IRa01 field are best fit with an impact parameter of b > 0.85, while less than 15% are expected in this range considering random orbital inclinations. By applying a cut at b < 0.85, meaning that ~15% of the potential planet population would be missed, the candidate sample decreases from 41 to 11. The lightcurves of 6 of those are best fit with such low host star densities that the planettostar size ratii imply unrealistic planet radii of R > 2 R_{Jup}. Two of the five remaining systems, CoRoT1b and CoRoT4b, have been identified as planets by the CoRoT team, for which the lightcurves alone rule out blended light at 14% (2σ) and 31% (2σ). One system possesses a Mdwarf secondary, one a candidate Neptune.
Conclusions. We show that in the first CoRoT field, IRa01, 85% of the planet candidates can be rejected from the lightcurves alone, if a cut in impact parameter of b < 0.85 is applied, at the cost of a <15% loss in planet yield. We propose to use this method on the Kepler database to study the fraction of real planets and to potentially increase the efficiency of followup.
Key words: planets and satellites: detection / planetary systems / planets and satellites: fundamental parameters
Appendix A is available in electronic form at http://www.aanda.org
© ESO, 2012
1. Introduction
With the CoRoT and Kepler space observatories in full swing (Baglin et al. 2006; Borucki et al. 2003), which both deliver thousands of lightcurves with unprecedented photometric precision and cadence, we have moved into an exciting new era of exoplanet research. Now, the characterisation of small, possibly rocky planets has finally become a realistic prospective (e.g. Corot7b, Leger et al. 2009; Kepler10b, Batalha et al. 2011). One of the biggest challenges is to seperate real planets from the significant fraction of (astrophysical) falsepositives that can mimic a genuine transit signal (e.g. Batalha et al. 2010). Groundbased transit surveys have revealed that stellar eclipsing binaries (EBs) blended with light from a third star are the main source of contamination (e.g. Udalski et al. 2002). Also, for superEarth planet candidates blends with a background transiting Jupitersized planet system can be important. In these systems the eclipse depth, shape and ellipsoidal light variations of an EB are diluted by the effects of chance alignment of a foreground or background star or associated companion inside a photometric aperture set by either the pixel scale or the point spread function. In addition, light from a third star in the photometric aperture can bias the fitted parameters of a planet transit system. High resolution, high signaltonoise spectra are normally required to exclude binary scenarios by excluding their large radial velocity or bisector variations, a process that can be very timeconsuming.
Stellar blends are common in spacebased transit surveys as apertures are relatively large (e.g. 19″ × 21″ for CoRoT), and target fields are crowded since the number of target stars is maximized in this way. To weed out falsepositives, the CoRoT team relies on an extensive groundbased followup campaign for onoff photometry to identify the transited star in the CoRoT aperture (Deeg et al. 2009) and high resolution imaging observations to identify possible stars that dilute the lightcurve of a planet candidate. Even so, many candidates remain unresolved and defy easy characterisation after such a campaign. Kepler uses its unique astrometric precision to minimise the number of blends, which can be identified by a position shift of the flux centroid during transit, but will still require enormous groundbased efforts on the remaining ~1200 candidates (e.g. Borucki et al. 2011). Together with the new influx of planet candidates from current surveys, possible future missions (such as PLATO; e.g. Catala et al. 2011) and groundbased efforts to hunt for planets around lowmass stars, the telescope demand for full followup may grow enormously. Therefore, any new technique or strategy that can eliminate even a moderate fraction of all candidates from the discovery lightcurves, prior to followup, is extremely valuable.
In this paper we investigate to what extent information from the lightcurves themselves can identify blend scenarios and eliminate them as planet candidates and on the other hand rule out blend scenarios in the case of true planet systems. Our key motivation is that the lightcurves of blended systems can not be perfectly fit by pure transit models and neither can genuine transits be fit by blended light models. In Sect. 2 we introduce our lightcurve fitting procedure and in Sect. 3 we apply it to simulated data of a transiting hot Jupiter and superEarth. While such a procedure provides a natural tool to distinguish blends from genuine planetary systems by lightcurve fitting, it breaks down for transits with high impact parameters. We therefore only consider transiting systems with impact parameter b < 0.85, loosing potentially ~15% of the planet catch, but significantly decreasing (by an order of magnitude) the required amount of followup observations. In Sect. 4 we apply our method to the candidates of the CoRoT IRa01 field, whose candidates are almost completely characterised through an extensive followup campaign, and discuss the results in Sect. 5.
2. Method
2.1. Transit fitting
Several methods have been presented in the literature to identify blended systems and to select the best planet candidates. Seager & MallenOrnellas (2003) proposed a diagnostic that involves fitting a trapezoid to the transit lightcurve to obtain estimates for the transit parameters and subsequently identify the best candidates. In this paper we use a method very similar to that used by Snellen et al. (2009) to reject blend scenarios for the transiting hot Jupiter OGLE2TRL9. It involves leastsquare fitting of a lightcurve using the standard transit parameters (see below) plus an additional parameter representing the extra light from a third light source. If the fit is significantly better with extra light, the lightcurve is from a blended system. If this is not the case, an upper limit to the third light fraction can be set to a degree depending on the signaltonoise of the data. This procedure is in essence similar to Blender, which is used by the Kepler team (e.g. Torres et al. 2011). However, Blender simulates physical systems involving so many parameters that it is impractical to run on a large number of candidates. Here we are not interested in the true nature of the second object (whether it is a background, foreground or physically related star), just in its possible influence on the transit lightcurve.
We assume at this point that lightcurves with obvious signs of the presence of a stellar binary, such as ellipsoidal light variations and/or secondary eclipses, have been excluded from the candidate list. Note that a useful upper limit to the amount of ellipsoidal light variation, and the likelihood of a genuine planetary secondary, can be obtained by taking a Fourier transform of the data with the transit signal removed. We therefore do not require EBOP (Popper & Etzel 1981) to model the complex binary effects in the lightcurve, but rather utilize an IDL routine that incorporates the analytical transit model of Mandel & Agol (2002; M&A). Our system simply consists of a secondary transiting a primary with possible additional light from a tertiary.
2.2. Transit parameters
We treat the transit midtime T_{0} and the orbital period P as fixed parameters, resulting from the candidate selection process. For extra simplicity we keep the limb darkening parameters fixed at the tabulated solar values for CoRoT white light, assuming quadratic parameters (a,b) = (0.44, 0.23) from Sing et al. (2010). Although this gives a small bias (<0.06 in impact parameter) for primary stars of different stellar type, the method is not meant for precise planet characterization and does not influence the characterization of potential blended and nonblended systems. Our transit model has three free parameters; the ratio of secondary over primary radii (R_{2}/R_{1}), the impact parameter of the transit b, which is the smallest projected distance of the centre of the secondary to that of the primary in units of R_{1}, and the density of the primary star ρ_{1}. This density can be converted to the scaled orbital radius (a/R_{1}), assuming that M_{1} ≫ M_{2}, through (1)The relative projected distances z between secondary and primary are computed from the input orbital phases φ, (2)
together with (R_{2}/R_{1}), these are used as input to a custommade IDL program, incorporating the routine from M&A, that computes the theoretical models. We introduce light to this transit system by adding the blended light fraction k, (3)
where F_{eclipse} is the original transit lightcurve. We then devise the following chisquare statistic to compare the lightcurve to the data F_{obs,i} with uncertainty σ_{obs,i},
Note that we assume circular orbits. This has no influence on the characterization of blended and nonblended systems, but it does affect the derived host star density, and is therefore important for the estimate of the radius of the secondary object. This is further discussed in Sect. 5.
2.3. MCMC
To obtain the bestmatching system parameters, we use a Monte Carlo Markov chain (MCMC) χ^{2} optimisation technique (e.g. Tegmark et al. 1998) to map out the probability distribution for each lightcurve parameter. The MCMC is found to be a more robust technique to obtain a global parameter solution in multiparameter space than (downhill) gridbased methods, due to the resolution inefficiency of the latter (e.g. Serra et al. 2011). In the MCMC algorithm, the parameters p_{i} are perturbed by an amount drawn from a normal distribution according to: , where f is the jump function and σ_{p} the standard deviation of the sampling distribution for each p. Subsequently χ^{2} is recalculated for these perturbed parameters and a Gaussian likelyhood ℒ ∝ exp( − χ^{2}/2) is determined. These random jumps in parameter space are accepted or rejected according to the MetropolisHastings rule (Metropolis et al. 1953; Hastings 1970). If the perturbed parameter set has a higher likelyhood ℒ′ than its progenitor, it will be accepted as a new chain point, otherwise it will be accepted with a probability of ℒ′/ℒ. We run the algorithm many times to build up a ’chain’ of parameter values and tweak σ_{p} and f such that ~ 40% of the jumps are accepted. After creating multiple chains from different starting conditions, we check proper model convergence and mixing of the individual chains using the Gelman and Rubin R statistic (Gelman & Rubin 1992). To save time, first k is set to zero at the minimum χ^{2} determined with MCMC analysis. Subsequently k is increased in small steps (but always kept fixed during the MCMC) with the previously found parameters as starting values. In this way the parameter values (adopting the median of the distribution) and the uncertainties in the parameters are determined as function of k in an efficient way.
Fig. 1
Simulated lightcurve for a transiting exoplanet system consisting of a hot Jupiter in a 2.5 day orbit around a solar type star with impact parameter b = 0.2 (black dots). The solid curves show diluted binary models with bestfit parameters determined by MCMC, for blended light fraction k = [0.2, 0.5, 0.8, 0.95]. 
3. Tests on synthetic lightcurves
In this section, we test our method on synthetic lightcurves to determine the required precision to detect or exclude third light in a particular transit system. We perform these simulations for two candidate systems: (i) a hot Jupiter orbiting a solar type star and (ii) a superEarth around a similar host.
3.1. Transiting hot Jupiter
We simulated a set of transit lightcurves for a hot Jupiter with R_{2} = 1 R_{Jup} and P = 2.5 days, orbiting a star with a solar density, for a range of impact parameters. The lightcurve for an impact parameter of b = 0.2 is shown in Fig. 1. As explained in the previous section, our method finds the best fit for a range in blended light fraction k. Of course, in this simulation a perfect fit is obtained for k = 0. As can be seen in Fig. 1, an increasingly worse fit is obtained for increasing k, most obviously seen by comparing the k = 0.95 model to the synthetic data. This latter model fit assumes that 95% of the light is from a third object, meaning that the unblended transit is actually a factor 20 deeper, hence 20% deep instead of 1%. It implies that R_{2}/R_{1} ~ 0.45, resulting in a much longer duration transit unless it is grazing. This results in the bestfitting k = 0.95 model being much more Vshaped than the synthetic lightcurve of the planet. We can now convert the differences between the synthetic lightcurves and model fits to χ^{2} values for each combination of b and k by assigning uncertainties to the synthetic data. In this way we can determine what photometric precision is required to exclude a certain blended light fraction in the lightcurves as a function of b. Figure 2 shows the precision per 5 min of data required to exclude a blended light fraction k at a 3σ level in a system with an impact parameter of b = 0.2, 0.5, 0.8, and 0.95. The required precision becomes more stringent for lower values of k and higher values of b. For b = 0.2, 80% blended light (k = 0.8) can be excluded in a lightcurve with a precision of only 2 × 10^{3} per 5 min, while for b = 0.8, 20% of blended light can only be rejected if the lightcurve has a precision of 4 × 10^{5} per 5 min.
Fig. 2
The photometric precision per 5 min required to exclude a blended light fraction k at 3σ for a hot Jupiter around a solar type host star (R_{2}/R_{1} = 0.1 and ρ_{ ∗ } = ρ_{⊙}), as a function of the system parameters b and k. The four solid curves are for impact parameters b = 0.2, 0.5, 0.8, and 0.95. The upper and lower horizontal dotted lines indicate the range in precision for objects in the IRa01 field, determined by Aigrain et al. (2009). 
3.2. Transiting superEarth
We performed also tests on a superEarth with R_{2} = 2.5 R_{ ⊕ } orbiting a sunlike star, following the same procedure as described above. Since the transit itself is a factor ~ 16 more shallow than for a Jupitersize planet, the level of precision required to reject blend scenarios is also significantly higher, as can be seen in Fig. 3. Note however that even for a blended light fraction of k = 0.95, the radius of transiting object R_{2} is still in the Jupitersize regime. Hence only if the blended light fraction is very high, k > 0.95, can an eclipsing binary mimic a superEarth transit.
Fig. 3
As for Fig. 2, but for a 2.5 R_{Earth} SuperEarth planet around a solar type host in a 2.5 day orbital. We can exclude 80% blended light at the 3σ level at a moderate impact parameter of b = 0.5. The horizontal dashed line refers to the precision reached in the discovery lightcurve of CoRoT 7b, the first rocky SuperEarth planet (Leger et al. 2009). 
4. Tests on candidates in the CoRoT IRa01 field
4.1. The data set
In this section we test our method on real data, using the lightcurves of the candidates selected by the CoRoT team from CoRoT field IRa01 (Carpano et al. 2009). In this first field targeted by CoRoT, 3898 bright stars were observed in chromatic mode (with a blue, green and red channel) and another 5974 in a single monochromatic “white” band in a 66 day staring run towards the Galactic anticenter. From the 50 initial candidates, a subsample of 29 promising targets received extensive followup as discussed in Moutou et al. (2009). Two of these have so far been identified as genuine planets: CoRoT1b, a low density R_{p} = 1.49 R_{Jup} transiting hot Jupiter around a G0V host (Barge et al. 2008) and CoRoT4b, a R_{p} = 1.19 R_{Jup} hot Jupiter around a F8V host (Aigrain et al. 2008). Seventeen additional systems were solved using the photometric and spectroscopic followup observations (Moutou et al. 2009). We choose to test our method on the 45 bright candidates with more than one transit observed, using the publicly available N2level data.
4.2. Precleaning of the lightcurves
We first combine the multicolor lightcurves into one single “white lightcurve” for each candidate under the assumption that the CoRoT analysis teams did not detect any significant variation of eclipse depth with wavelength, which would already have been a clear sign of blending effects. We first clip each lightcurve by removing outliers at the 5σ level. These outliers are mostly associated with the epochs at which the satellite passes the South Atlantic Anommaly (SAA) or moves in/out of the Earth’s shadow. We then iteratively refine the midtimes T_{0} and the orbital period P using the Kweevan Woerden method (Kwee & van Woerden 1956) and crosscorrelation with a theoretical transit model (e.g. Rauer et al. 2009). Individual transit events that show temporary jumps in flux, caused by the impact of energetic particles (mainly protons) onto the CCD (“hot pixels”), are excluded from our analysis. For 16 out of the initial 50 CoRoT IRa01 candidates (32%) we had to remove one or more transits from the lightcurve that were affected by such particle hits. Each individual lightcurve was then phasefolded around every transit. To normalise the data, we fit either a first order polynomial in a small range in phase ( ± 0.1 from midtransit) around each transit or a higher order polynomial (order n = 13) in a larger phase range (typically ± 0.4 in phase), depending on which approach delivers the lowest rms in and out of eclipse and the least red noise (Pont et al. 2006).
Fig. 4
Fourier diagram of an example noise spectrum prior to lightcurve cleaning. Amplitude of the bestfitting sine curve on the vertical axis is plotted against frequency. Peaks around frequencies of 1.0 and ~14 are due to remaining systematics related to the satellite orbit and Earth’s rotational period. 
Figure 4 shows a typical example of the dominant frequencies still remaining after the polynomial fit. For most objects, distinct peaks exist around periods of 103 min and at 24 h. We identify these peaks with remaining systematics, related to the satellite’s orbit and Earth’s rotational period, caused by ingress and egress of the spacecraft from Earth’s shadow, variations in gravity and magnetic field and changes in the levels of thermal and reflected light from the Earth (e.g. Aigrain et al. 2009). By folding the outofeclipse data onto the dominant frequencies of the Fourier diagram, we then fit a sinusoidal function to the remaining systematics, followed by median averaging over all transits. We subsequently binned the lightcurves and assign errors, according to the standard deviation divided by the square root of the number of points in each bin.
4.3. Fitting the lightcurves
Each lightcurve is first fitted with the method explained in Sect. 2, assuming k = 0, yielding the starting parameter sets (R_{2}/R_{1},b,ρ_{ ∗ }) for our blend analysis. In Fig. 5 we show the resulting MCMC distribution of impact parameter b versus R_{2}/R_{1} for all the 45 candidates. CoRoT WinIDs (a shortcut of the CoRoT run identification number, e.g. IRa01E12046) for each candidate are indicated, with yellow for the two confirmed planets CoRoT1b and CoRoT4b, in blue those candidates that have been confirmed to be blended systems by Moutou et al. (2009), and in red unsettled cases (either suspected early type stars with only few or very broad spectral lines for further radial velocity followup observations with HARPS or confirmed genuine EBs with nonplanetary secondary masses). As can be seen, a large fraction of the candidates are, assuming no blended light, best fitted with a very high (often larger than unity) impact parameter. This is even more clear in the distribution of fitted impact parameters as shown in Fig. 6. For 32 out of the 45 (~70%) candidates b > 0.85, while from geometric arguments it is expected that ~15% of planets would be found at such a high impact parameter. Assuming that all eleven candidates at b < 0.85 are nonblended systems only ~1.6 objects are expected at b > 0.85. Since our tests in Sect. 3 have shown that it is very difficult to distinguish blends from nonblended systems at high impact parameters due to their Vshaped lightcurves, we apply a cut in the candidate list at b < 0.85, knowing that we will potentially remove only a small fraction of the planet yield, in the case of the CoRoTa01 field <0.3 planets. From this it can be seen that it is highly likely that all candidates with b > 0.85 are blended and/or grazing eclipsing binaries. For the eleven remaining candidates we used the transit parameters from the k = 0 model to refit the lightcurve with an increasing value of k, as outlined in Sect. 3. In this way we redetermine the best fit solution and χ^{2} as a function of k. As an example we show the best fit transit models for a range of k and the χ^{2} as function of k for candidate E14617 in Fig. 7. As can be seen, the lightcurve can only be well fitted by models with a low k. E.g. the χ^{2} of the best fitting k = 0.5 model is ~40% higher than that for k = 0. The 2 sigma upper limit for the fraction of blended light (Δχ^{2}) is k = 0.20. We performed this same analysis for all eleven remaining candidates for which the χ^{2} versus k plots are shown in the Appendix, together with their best fit lightcurves. None of these candidates are better fitted by a high k model than a low k model, indicating that all blended systems have moved out of the remaining sample since they are all fitted with a high impact parameter. For six objects a significant fraction of blended light can be excluded from the lightcurve alone, including CoRoT1b and CoRoT4b. It would therefore not have been necessary to check whether the variability in these candidates came from the target star or not and the followup could have immediately concentrated on radial velocity measurements.
Fig. 5
The MCMC solution for all IRa01 candidates in (b,(R_{p}/R_{ ∗ })) space. Note the strong parameter degeneracy at high b. Yellow CoRoT WinIDs are the two confirmed planets CoRoT1b and CoRoT4b, blue objects are confirmed blends from the followup work presented in Moutou et al. (2009), and red sources are either unsettled cases or confirmed genuine binaries with nonplanetary secondary masses from the radial velocity variations. 
Fig. 6
Distribution of fitted impact parameter of the CoRoT IRa01 candidates. The distribution is strongly peaked around b = 1.0, indicating a significant population of (blended) EB contaminants. For a genuine planet distribution we would expect a flat histogram that falls off at high impact parameter. The dotted vertical line indicates the b = 0.85 cutoff we have proposed in this paper. 
Fig. 7
Panel a): the reduced χ^{2} as function of blended light fraction k. The horizontal dashed line shows the 2,3 and 5σ rejection criteria. Our lightcurve model directly indicates an early type main sequence stellar host, with a 2σ upper limit for blended light of k ~ 20%. The low stellar density implies a large secondary radius, rejecting the planet hypothesis. Panel b): best fitting EB models with blended light fraction k = [0.2,0.5,0.90,0.95], clearly showing that solutions with low k are favoured. Note that an orbit with an eccentricity of e = 0.5, orientated in the right way, could increase the estimated stellar density to that of the Sun, and decrease R_{2} to 2 R_{Jup}. This ambiguity can be easily removed by taking a single spectrum of the star, resolving its spectral type. 
All parameters of the remaining candidates are shown in Table 1. An additional cut in the candidate list is made using a combination of the best fit mean stellar densities ρ_{ ∗ } and R_{2}/R_{1}, as shown in Fig. 8. Six of the candidates have host stars with densities corresponding to Astars, resulting in unrealisticly large secondary radii of > 2.0 R_{Jup}. Note that there is currently no consensus on the upper limit of planet size, meaning that by setting a hard limit on planet radius we may exclude very large or bloated (hot) Jupiters. However, there are currently only 4 out of 219 transiting exoplanets reported with radii larger than 1.8 R_{Jup} (http://www.exoplanet.eu). Also, the probability that the secondary is a midtype Mdwarf rather than a genuine planet increases when considering larger radii. This results in a remaining planet candidate sample of 5 objects instead of the original 45 using arguments based on the lightcurve alone. These five objects have been marked with filled symbols in Fig. 8. Details on each system are discussed in Appendix A.
5. Discussion
In this paper we investigated to what extent we can use the high signaltonoise lightcurves of spacebased exoplanet transit surveys to identify blended light scenarios, and eliminate them as planet candidates. We concentrated on the 51 exoplanet candidates from the first CoRoT IRa01 field (Carpano et al. 2009). About 70% of the 51 planet candidates in the CoRoT IRa01 field are best fit with an impact parameter of b > 0.85, which at face value already indicates that the candidate distribution is strongly contaminated by blended and/or grazing systems. We find that by cutting a candidate sample such that those objects with high impact parameter are removed, at the cost of losing a small fraction of potential planets, a significant reduction in required followup observations can be achieved. Of all candidates, only 5 remain in the final sample of which two are genuine planet systems, one is a low mass transiting M dwarf and one is a candidate Neptune.
The Vshaped lightcurves of neargrazing planet systems are strongly degenerate with blended eclipsing binary systems and can therefore not be distinguished from each other. How many planets are potentially missed by invoking the cut in impact parameter? Of the known transiting exoplanets, ~ 6% has an impact parameter larger than 0.85 and ~16% an impact parameter of more than 0.75^{1}. The cumulative probability of a particular transit at a given impact parameter greater or equal to a cutoff value b_{X} and transit depth ΔF is given by: (5)Note that this expression is different from the equation presented in Seager & MallenOrnellas (2003), because the maximum impact parameter in their formula is determined by the grazing condition b_{max} = 1 − R_{2}/R_{1}, yielding a minus sign in Eq. (5). For a 1 R_{Jup} planet around a solar type star ~ 22% would potentially be missed by setting the cut in impact parameter (~6% according to Seager and MallenOrnellas). However, extremely grazing systems will be very shallow and of short duration and will therefore provide very limited physical information. For example, a grazing 1 R_{Jup} with impact parameter b = 1.05, will show a transit with a duration of 30% and only 20% of the depth of a transit with b = 0. Therefore, the actual planet loss fraction will be closer to the predictions of Seager & MallenOrnellas (2003), i.e. <15%.
The candidate sample that survives the impact parameter cut.
Fig. 8
The R_{2}/R_{1} size ratio versus the log of the stellar density for the CoRoT candidates in the IRa01 field, assuming k = 0. The dotted lines mark the densities of A to M type mainsequence hosts. The five filled dots are the candidates that survive both our cuts in impact parameter and secondary size. The two confirmed transiting hot Jupiters CoRoT1b and CoRoT4b are shown as yellow filled dots. Open circles are the candidates we have excluded using our cuts. Blue circles indicate sources which have been identified as blended EBs by the CoRoT team followup, and red circles are either unsettled cases in the CoRoT followup or systems identified as genuine EBs through their radial velocities. The four solid curves indicate R_{2} = [0.10,0.50,1.0,2.0] R_{Jup}, assuming the main sequence massradius relation of Cox (2000) for the primary. 
In this paper we have made the assumption of circular orbits, but radial velocity surveys teach us that such an assumption is not valid for longer periods (e.g. Butler et al. 2006). In addition, Barnes (2007) shows that a planet with an eccentric orbit is more likely to transit by a factor of (1 − e^{2})^{1} than a planet in a circular orbit with the same semimajor axis. A significant population of transiting exoplanets with an eccentric orbit is therefore expected for long duration spacebased surveys. Because the planet orbital velocity varies from to between periastron and apastron in an eccentric orbit, transit duration can vary as function of e and ω (the angle of pericenter). This leads to a wrong fit of the host star density (e.g. Kipping 2010a; Tingley et al. 2011), therefore directly affecting our estimate of the secondary radius R_{2}. We therefore can not reliable make the planettostar ratio versus host star density cut in the eccentric orbit case for longer period planets (P > 3.0 days). Fortunately, the fitted impact parameter, R_{2}/R_{1} and blended light fraction k are not affected by an eccentric orbit. This means that we can still first apply a cut in impact parameter b < 0.85 and remove likely blends. To subsequently determine the real host star density it is sufficient to take a single highresolution spectrum to determine ρ_{1} and estimate R_{2}. Using this spectroscopically determined density an upper limit to ecos(ω) can be set. One particular case in our sample is CoRoT4617 with an orbital period of P = 19.76 days. Assuming a circular orbit, the host star is estimated to have a density only ~ 4% of that of the Sun, in accordance with an early Bstar. This would imply that the radius of the secondary object has R_{2} ~ 8 R_{Jup}. However, an orbit with an eccentricity of e = 0.5, orientated in the right way, could increase the estimated stellar density to that of the Sun, and decrease R_{2} to 2 R_{Jup}. This ambiguity can be easily removed by taking a single spectrum of the star, resolving its spectral type.
The method presented here is designed to remove falsepositives in candidate lists through the identification of blended light. We do not assign a likelihood of planetary nature to the remaining candidates, meaning that we do not assess whether these are genuine planet systems, we just removed those systems which are not (except for a small fraction of collateral damage). However, it is anyway interesting to link blended light fractions to the population of random background eclipsing binaries. Assuming that 1:300 of field stars are eclipsing binaries (Devor et al. 2008), and 1:1000 stars have a transiting hot Jupiter, we require an average of 0.3 background stars within the PSF, and within the magnitude range set by the limit of blended light, to have an equal likelihood for the two scenarios, and to end up with half of the remaining objects as eclipsing binaries. For a typical magnitude (V = 14.0) for the candidate star, taking into account the size of the CoRoT PSF, this is reached at a Δmag = ~1.5, corresponding to k = 0.8. For 8 out of the 10 remaining targets this level of blended light is excluded at the > 3σ level. For a typical limit of k < 0.6, the background eclipsing binary can at most be 0.5 magnitudes fainter than the target star, making this scenario a factor ~5 less likely. Do note however that this does not take into account physical triple systems, for which radial velocity followup is required to exclude them. Recently, the Kepler team have announced the discovery of ~1200 planet candidates (Burucki et al. 2011). We propose to use the method presented here on this candidate list, to identify clear blend systems to obtain a better estimate on the fraction of planet systems in this sample.
6. Conclusions
In this paper we have investigated to what extent information from lightcurves of a spacebased exoplanet transit survey can identify blended light scenarios and eliminate them as planet candidates, to significantly decrease the required amount of followup time. If a lightcurve has sufficiently high signaltonoise, a distinction can be made between a blended eclipsing binary and a transiting exoplanet. We first have simulated lightcurves of stellar blends and transiting planet systems to determine the required signaltonoise as a function of impact parameter and transit depth. Our simulations show that blend scenarios can be distinguished from transiting systems at low impact parameter. At high impact parameter, blended and nonblended systems both produce Vshaped transits and are indistinguishable from each other. We have subsequently tested our method on real data from the first IRa01 field of CoRoT, concentrating on the 51 candidates already identified by the CoRoT
team (Carpano et al. 2009). We show that 70% of the planet candidates in the CoRoT IRa01 field are best fit with an impact parameter of b > 0.85, whereas ~15% are expected assuming random orbital orientations. By applying a cut at b < 0.85, meaning that ~15% of the potential planet population would be missed, the candidate sample decreases from 41 to 11. The lightcurves of 6 of those are best fit with such a low host star density that the planettostar size ratio implies an unrealistic planet radius of R_{2} > 2 R_{Jup}. From the remaining five, two systems, CoRoT1b and CoRoT4b, have been identified by the CoRoT team as planets, for which the lightcurves alone rule out blended light at a 25%(2σ) and 40%(2σ). One other candidate is also consistent with a nonblended system, but is a late Mdwarf, which will always require radial velocity followup for confirmation since Mdwarfs can have similar radii as Jupiter mass planets. One other system consists of a candidate Neptune around a Mdwarf according to Moutou et al. (2009). We have therefore shown that 85% of the planet candidates can be rejected for the IRa01 field from the lightcurves alone. We propose to use this method on the Kepler database to study the fraction of real planets and to potentially increase the efficiency of followup. For long period candidates, possible nonzero eccentricity will affect the cut in planettostar ratio versus host star density, effectively increasing the sample size. However a single highresolution spectrum would be sufficient to determine the real host star density and estimate the size of transiting objects.
References
 Aigrain, S., Collier Cameron, A., Ollivier, M., et al. 2008, A&A, 488, L43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aigrain, S., Pont, F., Fressin, F., et al. 2009, A&A, 506, 425 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baglin, A., Auvergne, M., Boisnard, L., et al. 2006, 36th COSPAR Scientific Assembly, 36, 3749 [Google Scholar]
 Barge, P., Baglin, A., Auvergne, M., et al. 2008, A&A, 482, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barnes, J. W. 2007, PASP, 119, 986 [NASA ADS] [CrossRef] [Google Scholar]
 Batalha, N. M., Rowe, J. F., Gilliland, R. L., et al. 2010, ApJ, 713, L103 [NASA ADS] [CrossRef] [Google Scholar]
 Batalha, N. M., Borucki, W. J., Bryson, S. T., et al. 2011, ApJ, 729, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Borucki, W. J., Koch, D. G., Lissauer, J. J., et al. 2003, Proc. SPIE, 4854, 129 [Google Scholar]
 Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 728, 117 [NASA ADS] [CrossRef] [Google Scholar]
 Butler, R. P., Wright, J. T., Marcy, G. W., et al. 2006, ApJ, 646, 505 [NASA ADS] [CrossRef] [Google Scholar]
 Catala, C., Appourchaux, T., & Plato Mission Consortium, T. 2011, J. Phys. Conf. Ser., 271, 012084 [NASA ADS] [CrossRef] [Google Scholar]
 Carpano, S., Cabrera, J., Alonso, R., et al. 2009, A&A, 506, 491 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cox, A. N., & Pilachowski, C. A. 2000, Phys. Today, 53, 100000 [Google Scholar]
 Deeg, H. J., Gillon, M., Shporer, A., et al. 2009, A&A, 506, 343 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Devor, J., Charbonneau, D., O’Donovan, F. T., Mandushev, G., & Torres, G. 2008, ApJ, 135, 850 [CrossRef] [Google Scholar]
 Hastings, W. K. 1970, Biometrika, 57, 97 [Google Scholar]
 Geman, A., & Rubin, D. B. 1992, Statistical Science, 7, 457 [Google Scholar]
 Kipping, D. M. 2010, MNRAS, 407, 301 [NASA ADS] [CrossRef] [Google Scholar]
 Kwee, K. K., & van Woerden, H. 1956, Bull. Astron. Inst. Netherlands, 12, 327 [Google Scholar]
 Léger, A., Rouan, D., Schneider, J., et al. 2009, A&A, 506, 287 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Mandel, K., & Agol, E. 2002, ApJ, 580, L171 [NASA ADS] [CrossRef] [Google Scholar]
 Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., & Teller, E. 1953, J. Chem. Phys., 21, 1087 [NASA ADS] [CrossRef] [Google Scholar]
 Moutou, C., Pont, F., Bouchy, F., et al. 2009, A&A, 506, 321 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pont, F., Zucker, S., & Queloz, D. 2006, MNRAS, 373, 231 [NASA ADS] [CrossRef] [Google Scholar]
 Popper, D. M., & Etzel, P. B. 1981, AJ, 86, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Rauer, H., Queloz, D., Csizmadia, Sz., et al. 2009, A&A, 506, 281 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Seager, S., & MallénOrnelas, G. 2003, ApJ, 585, 1038 [NASA ADS] [CrossRef] [Google Scholar]
 Serra, P., Amblard, A., Temi, P., et al. 2011, ApJ, 740, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Sing, D. K. 2010, A&A, 510, A21 [Google Scholar]
 Smalley, B., Anderson, D. R., Collier Cameron, A., et al. 2011, A&A, 526, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Snellen, I. A. G., Koppenhoefer, J., van der Burg, R. F. J., et al. 2009, A&A, 497, 545 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tingley, B., Bonomo, A. S., & Deeg, H. J. 2011, ApJ, 726, 112 [NASA ADS] [CrossRef] [Google Scholar]
 Torres, G., Fressin, F., Batalha, N. M., et al. 2011, ApJ, 727, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Udalski, A., Szewczyk, O., Zebrun, K., et al. 2002, Acta Astron., 52, 317 [Google Scholar]
Online material
Appendix A: Comments on individual sources
In this Appendix we discuss in detail the sample of ten remaining CoRoT candidates, that were selected using the cut in impact parameter and were presented in Sect. 4 and Table 1. In Figs. A.1, we show for each candidate the blended light fraction k versus reduced χ^{2}, and the best fitting blended light models for k = 0.2, 0.5, 0.9, and 0.95. In Table A.1, we show bestmatching system parameters for the full CoRoT IRa01, assuming no blended light.
A.1. Selected planet candidates from the lightcurves alone
E211260102890318
We find a 2σ upper limit for blended light contribution of k < 0.14, therefore a blend scenario can be excluded for this source at high confidence using the lightcurve alone. In addition, by assuming that the host star is on the main sequence, its mean density points to a ~ 1.5 R_{Jup} radius, well in the range of known hot Jupiters. Of course, this source is exoplanet CoRoT1b (Barge et al. 2008).
E103300102912369
We find a 2σ upper limit for blended light contribution of k < 0.31 from its lightcurve, meaning that only a small contribution of blended light is tollerated. Assuming the host star is on the main sequence, its mean density points to a ~ 1.2 R_{Jup} radius for the secondary. This object is identified as exoplanet CoRoT4b (Aigrain et al. 2008). Eventhough the CoRoT4b host star is of similar brightness as CoRoT1b, the significantly longer orbital period, the residual variability (caused by a spotted rotating stellar photosphere) and the 1.8 times smaller transit depth are the causes of the lower confidence on blended light.
E202030102825481
The 2σ upper limit for blended light is k < 0.3 from its lightcurve. Radial velocity followup by the CoRoT team showed this to be an eclipsing binary of a lowmass M dwarf and a Gtype primary (Morales et al., in prep). Assuming the host star is on the main sequence, its mean density points to a ~1.7 R_{Jup} radius. Although not a planet, it is consistent with a nonblended system as found from our lightcurve fitting. Such systems always require RV followup since late M dwarfs and Jupitermass planets can have similar radii.
E217120102826302
We find a 2σ upper limit for blended light contribution of k < 0.93. We can therefore only exclude a high contribution of blended light for this shallow (2.4 mmag) transit. This means that at 2σ confidence the true eclipse depth is less than 2.4% in the presence of blended light. The fitted host star mean density points to an early type or evolved system. HARPS radial velocity followup has confirmed that the host star is an evolved fast rotator and Moutou et al. (2009) conclude that a triple system is the most probable scenario.
E141080102779966
Because of the poor signaltonoise of this transit and the relatively high impact parameter b = 0.8, the 2σ upper limit for blended light is k < 0.95, therefore only a very high contribution of blended light can be excluded for this candidate. Assuming the host star is on the main sequence, its density is slightly lower compared to the solar value, indicating a stellar radius of R_{1} ~ 1.2 R_{⊙}. However, spectroscopic followup with HARPS suggested that the host is a low mass (~0.8 M_{⊙}) star. No additional followup has thusfar been obtained by the CoRoT team.
A.2. Candidates rejected due to their large size
E146170102753331
The 2σ upper limit for blended light is k < 0.20, therefore a blend scenario can be excluded at high confidence for this source. Assuming the host star is on the main sequence, its very low density points to an early Btype primary with a K dwarf secondary. The planet hypothesis is rejected and no additional followup is therefore required judging from the lightcurve alone. Note that an orbit with an eccentricity of e = 0.5, orientated in the right way, could increase the estimated stellar density to that of the Sun, and decrease R_{2} to 2 R_{Jup}. This ambiguity can be easily removed by taking a single spectrum of the star, resolving its spectral type.
E224300102815260
We find a 2σ upper limit for blended light contribution of k < 0.44. Again, only a small contribution of blended light is tollerated. Assuming the host star is on the main sequence, its mean density, consistent with an A type or evolved star, points to a radius R_{2} > 2.5 R_{Jup}. Radial velocity followup by the CoRoT team showed this to be a single lined eclipsing binary of a fast rotating host star and an early type M dwarf (Moutou et al. 2009).
E240730102863810
For this source, we find a 2σ upper limit for blended light of k < 0.67. This object shows ~ 4% deep eclipses around a host star that is ~ 20% less dense than the sun. This candidate was introduced in the original list of Carpano et al. (2009), but is not mentioned in the followup paper of Moutou et al. (2009). With an anticipated secondary radius of ~ 2.1 R_{Jup} this object could still belong to the rare group of low mass stars or brown dwarfs. In the case of a stellar M5 secondary, the secondary eclipse would be detectable at ~ 3.5 mmag in depth.
E217360102855534
The 2σ upper limit for blended light is k < 0.62. The low mean density of the host star, consistent with a very early main sequence or evolved star, points to a > 2.0 R_{Jup} radius. Analysis of the lightcurve reveals a secondary eclipse at the 9σ level, which indicates the secondary is in fact a low mass star. CoRoT radial velocity followup has confirmed that the host star is a fast rotating early type star and the system is a single lined eclipsing binary.
E237240102759638
For this source, we find a 2σ upper limit for blended light of k < 0.78. Assuming the host star is on the main sequence, its very low density points to an A type primary, therefore R_{2} > 2.0 R_{Jup}. This object is listed both as a planet candidate and a binary by Carpano et al. (2009).
Fig. A.1
For each CoRoT IRa01 candidate: the blended light fraction k versus reduced χ^{2} (left panels) and the best fitting blended light models for k = 0.2, 0.5, 0.9, and 0.95. 
Fig. A.1
continued. 
Fig. A.1
continued. 
Fitting parameters for our blend models when applied to the CoRoT IRa01 sample, assuming k = 0.
All Tables
Fitting parameters for our blend models when applied to the CoRoT IRa01 sample, assuming k = 0.
All Figures
Fig. 1
Simulated lightcurve for a transiting exoplanet system consisting of a hot Jupiter in a 2.5 day orbit around a solar type star with impact parameter b = 0.2 (black dots). The solid curves show diluted binary models with bestfit parameters determined by MCMC, for blended light fraction k = [0.2, 0.5, 0.8, 0.95]. 

In the text 
Fig. 2
The photometric precision per 5 min required to exclude a blended light fraction k at 3σ for a hot Jupiter around a solar type host star (R_{2}/R_{1} = 0.1 and ρ_{ ∗ } = ρ_{⊙}), as a function of the system parameters b and k. The four solid curves are for impact parameters b = 0.2, 0.5, 0.8, and 0.95. The upper and lower horizontal dotted lines indicate the range in precision for objects in the IRa01 field, determined by Aigrain et al. (2009). 

In the text 
Fig. 3
As for Fig. 2, but for a 2.5 R_{Earth} SuperEarth planet around a solar type host in a 2.5 day orbital. We can exclude 80% blended light at the 3σ level at a moderate impact parameter of b = 0.5. The horizontal dashed line refers to the precision reached in the discovery lightcurve of CoRoT 7b, the first rocky SuperEarth planet (Leger et al. 2009). 

In the text 
Fig. 4
Fourier diagram of an example noise spectrum prior to lightcurve cleaning. Amplitude of the bestfitting sine curve on the vertical axis is plotted against frequency. Peaks around frequencies of 1.0 and ~14 are due to remaining systematics related to the satellite orbit and Earth’s rotational period. 

In the text 
Fig. 5
The MCMC solution for all IRa01 candidates in (b,(R_{p}/R_{ ∗ })) space. Note the strong parameter degeneracy at high b. Yellow CoRoT WinIDs are the two confirmed planets CoRoT1b and CoRoT4b, blue objects are confirmed blends from the followup work presented in Moutou et al. (2009), and red sources are either unsettled cases or confirmed genuine binaries with nonplanetary secondary masses from the radial velocity variations. 

In the text 
Fig. 6
Distribution of fitted impact parameter of the CoRoT IRa01 candidates. The distribution is strongly peaked around b = 1.0, indicating a significant population of (blended) EB contaminants. For a genuine planet distribution we would expect a flat histogram that falls off at high impact parameter. The dotted vertical line indicates the b = 0.85 cutoff we have proposed in this paper. 

In the text 
Fig. 7
Panel a): the reduced χ^{2} as function of blended light fraction k. The horizontal dashed line shows the 2,3 and 5σ rejection criteria. Our lightcurve model directly indicates an early type main sequence stellar host, with a 2σ upper limit for blended light of k ~ 20%. The low stellar density implies a large secondary radius, rejecting the planet hypothesis. Panel b): best fitting EB models with blended light fraction k = [0.2,0.5,0.90,0.95], clearly showing that solutions with low k are favoured. Note that an orbit with an eccentricity of e = 0.5, orientated in the right way, could increase the estimated stellar density to that of the Sun, and decrease R_{2} to 2 R_{Jup}. This ambiguity can be easily removed by taking a single spectrum of the star, resolving its spectral type. 

In the text 
Fig. 8
The R_{2}/R_{1} size ratio versus the log of the stellar density for the CoRoT candidates in the IRa01 field, assuming k = 0. The dotted lines mark the densities of A to M type mainsequence hosts. The five filled dots are the candidates that survive both our cuts in impact parameter and secondary size. The two confirmed transiting hot Jupiters CoRoT1b and CoRoT4b are shown as yellow filled dots. Open circles are the candidates we have excluded using our cuts. Blue circles indicate sources which have been identified as blended EBs by the CoRoT team followup, and red circles are either unsettled cases in the CoRoT followup or systems identified as genuine EBs through their radial velocities. The four solid curves indicate R_{2} = [0.10,0.50,1.0,2.0] R_{Jup}, assuming the main sequence massradius relation of Cox (2000) for the primary. 

In the text 
Fig. A.1
For each CoRoT IRa01 candidate: the blended light fraction k versus reduced χ^{2} (left panels) and the best fitting blended light models for k = 0.2, 0.5, 0.9, and 0.95. 

In the text 
Fig. A.1
continued. 

In the text 
Fig. A.1
continued. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.