Issue 
A&A
Volume 573, January 2015



Article Number  A124  
Number of page(s)  14  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201424591  
Published online  09 January 2015 
SOPHIE velocimetry of Kepler transit candidates
XIV. A joint photometric, spectroscopic, and dynamical analysis of the Kepler117 system^{⋆,}^{⋆⋆}
^{1}
Aix Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de
Marseille), UMR 7326,
13388
Marseille,
France
email:
giovanni.bruno@lam.fr
^{2}
Centro de Astrofísica, Universidade do Porto,
Rua das Estrelas, 4150762
Porto,
Portugal
^{3}
Instituto de Astrofísica e Ciências do Espaço, Universidade do
Porto, CAUP, Rua das Estrelas,
PT 4150762
Porto,
Portugal
^{4}
Observatoire Astronomique de l’Université de Genève,
51 chemin des Maillettes,
1290
Versoix,
Switzerland
^{5}
INAF – Osservatorio Astrofisico di Torino, via Osservatorio
20, 10025
Pino Torinese,
Italy
^{6}
Observatoire de HauteProvence, 04670 Saint Michel l’Observatoire,
France
^{7}
Institut d’Astrophysique de Paris, UMR 7095 CNRS, Université
Pierre & Marie Curie, 98bis
boulevard Arago, 75014
Paris,
France
Received: 11 July 2014
Accepted: 8 November 2014
As part of our followup campaign of Kepler planets, we observed Kepler117 with the SOPHIE spectrograph at the Observatoire de HauteProvence. This F8type star hosts two transiting planets in nonresonant orbits. The planets, Kepler117 b and c, have orbital periods ≃ 18.8 and ≃ 50.8 days, and show transittiming variations (TTVs) of several minutes. We performed a combined Markov chain Monte Carlo (MCMC) fit on transits, radial velocities, and stellar parameters to constrain the characteristics of the system. We included the fit of the TTVs in the MCMC by modeling them with dynamical simulations. In this way, consistent posterior distributions were drawn for the system parameters. According to our analysis, planets b and c have notably different masses (0.094 ± 0.033 and 1.84 ± 0.18M_{J}) and low orbital eccentricities (0.0493 ± 0.0062 and 0.0323 ± 0.0033). The uncertainties on the derived parameters are strongly reduced if the fit of the TTVs is included in the combined MCMC. The TTVs allow measuring the mass of planet b although its radial velocity amplitude is poorly constrained. Finally, we checked that the best solution is dynamically stable.
Key words: planetary systems / stars: individual: Kepler117 / techniques: photometric / techniques: radial velocities / techniques: spectroscopic / methods: statistical
Appendix A is available in electronic form at http://www.aanda.org
Radial velocity tables are only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/573/A124
© ESO, 2015
1. Introduction
In the past few years, the number of known multiple planet systems detected by the Kepler space telescope has enormously increased. Multiple transiting planet systems have a low falsepositive probability: using conservative hypotheses, Lissauer et al. (2012) estimated a 1.12% probability of observing two false positives in the same system and a 2.25% probability for a system to host a planet and show the features of a false positive at the same time. Their estimation was based on the assumptions that false positives are randomly distributed among the Kepler targets and that there is no correlation between the probability of a target to host one or more detectable planets and display false positives.
At the time of writing, Kepler has detected 469 multiple planet systems^{1}. Among them, Kepler117 (also named KOI209) hosts the two transiting planets Kepler117 b and Kepler117 c. These planets were presented by Borucki et al. (2011) as candidates and validated by Rowe et al. (2014) with a confidence level of more than 99%, while radial velocity observations were still unavailable. After subjecting Kepler117 b and c to various falsepositive identification criteria, Rowe et al. used the statistical framework of Lissauer et al. (2012) (further refined in Lissauer et al. 2014) to promote them to bona fide exoplanets. Kepler117 b and c were found to have orbital periods ≃ 18.8 and ≃ 50.8 days and radii ≃ 0.72 and ≃ 1.04R_{J}.
Multiple planet systems offer insights on their dynamical history (e.g., Batygin & Morbidelli 2013) and can show transittiming variations (TTVs; Agol et al. 2005), especially if the planets are in mean motion resonance. TTVs are a powerful tool for detecting nontransiting planets, and for the determination of planetary masses (Holman et al. 2010; Nesvorný et al. 2013; Barros et al. 2014; Dawson et al. 2014) and can be a tracer of stellar activity, as well (Barros et al. 2013; Oshagh et al. 2013). Moreover, nondetected TTVs can cause an underestimation of the uncertainty on the stellar density derived from the photometry (Kipping 2014).
Using only the first quarter of the Kepler photometric data, Steffen et al. (2010) predicted TTVs to be observable for the Kepler117 system. At that time, the photometric time coverage was not sufficient to allow a verification. The TTVs were later confirmed by Mazeh et al. (2013). According to these authors, the ratio between the periodic modulation of the TTVs of the inner planet (b) and the orbital period of the outer one (c), P_{TTVs,b}/P_{c}, is ≃ 0.997. This ratio is the closest to 1 among the Kepler candidate twoplanets systems with TTVs presented in that paper. The similarity between the two periodicities is a strong indication that the two bodies are in the same system and thus is another argument for the validation of Kepler117 b and c.
In this paper, we included the information from the TTVs in the combined fit of the system parameters together with the photometry and radial velocities we acquired during our observation campaigns with the SOPHIE spectrograph at the Observatoire de HauteProvence. By fully exploiting the data, we obtained a precise measure of the masses and radii of the planets. In Sect. 2, the data acquisition and reduction is discussed. In Sect. 3, we describe the treatment of the stellar spectra, and in Sect. 4 we report on the stellar activity, the measurement of the TTVs, and the joint Bayesian fit of the system parameters. In Sect. 5, the results are discussed. The implications and conclusions are given in Sect. 6.
2. Observations and data reduction
2.1. Kepler photometric observations
Kepler117 was observed by the Kepler space telescope from quarter 1 to 17 between May 2009 and May 2013. The first three quarters were covered by a sampling of 29.4 min (longcadence data, LC), the following were sampled every 58.5 s (shortcadence data, SC). We chose to use the LC data only for quarters from 1 to 3, and relied on the SC data for the others. The light curves, already reduced by the Kepler pipeline (Jenkins et al. 2010), are publicly available on the Mikulski Archive for Space Telescopes (MAST)^{2}. We made use of the light curves corrected by the Presearch Data Conditioning (PDC) module, available in the light curve fits file.
For all the quarters and for both the LC and the SC data, the dispersion of the contamination of nearby stars, corrected for by the pipeline, is lower than 1%. The contamination value, then, was fixed in the following combined analysis (Sect. 4.2.2).
We isolated the photometric signal around every transit using a preliminary estimate of the ephemeris, following Rowe et al. (2014). The transits of the two planets sometimes superpose because of their different periods (≃ 18.8 days and ≃ 50.8 days). We discarded these overlapping transits because the software we used to fit the data sets (Sect. 4.2.2) does not yet include this modeling. No secondary eclipse was found, as expected from the relatively long periods.
We normalized the transits by fitting a secondorder polynomial to the flux outside of the transits and rejected the outliers through a 3σ clipping.
2.2. Spectroscopic observations
Kepler117 is part of our followup program of Kepler candidates (Bouchy et al. 2011; Santerne et al. 2012). We acquired 15 spectra of this star during two observing seasons, between July 2012 and November 2013, using the SOPHIE spectrograph at the 1.93 m telescope of the Observatoire de HauteProvence (Perruchot et al. 2008; Bouchy et al. 2013). The instrument was set in highefficiency mode, with a spectral resolution λ/ Δλ ~ 38 000. The exposures lasted from 1200 to 3600 s for a signaltonoise ratio (S/N) per pixel at 550 nm between 9 and 17. The spectra were reduced using the SOPHIE pipeline (Bouchy et al. 2009). The radial velocities (RVs) and their uncertainties were obtained through a Gaussian fit of the crosscorrelation function (CCF) with numerical masks corresponding to the F0, G2, and K5 spectral types. The final RVs were measured with the G2 mask because the spectral analysis showed Kepler117 to be close to a G star (i.e., Rowe et al. 2014, verified in Sect. 3). However, using different masks to compute the CCF did not result in a systematic difference between the RVs. The RV reference star HD185144 (Howard et al. 2010; Bouchy et al. 2013; Santerne et al. 2014) was used to correct the RVs by between ~5 and ~30 m s^{1}. Three spectra were affected by the moonlight: we corrected them for the RV of the Moon, as discussed in Baranne et al. (1996), Pollacco et al. (2008), and Hébrard et al. (2008). The charge transfer inefficiency effect was corrected for using the prescription of Santerne et al. (2012). The first three echelle orders at the blue edge of the spectrum were not used to calculate of the RVs because their low S/N degrades the precision of the measurements.
We rejected the point at BJD = 2 456 551.49295 because of its low S/N (9 at 550 nm, the lowest of all the set).
We checked for linear correlations between the bisector span of the CCF and the RVs, following Queloz et al. (2001) (Fig. 1). If linear correlations are observed, the planetary scenario is very likely to be rejected in favor of a blend. The Spearmanrankorder correlation coefficient between the bisector span and the RVs, excluding the points contaminated by the moonlight, is − 0.08 ± 0.32. The pvalue for this coefficient, with the null hypothesis of no correlation, is 0.98. Similarly, the Spearman correlation coefficient between the full width at half maximum of the CCFs of the spectra and their respective RV is 0.27 ± 0.30, with a pvalue of 0.39. The two diagnostics on the CCF are clearly compatible with the planetary scenario.
Log of SOPHIE radial velocity observations.
Fig. 1
Bisector span of the CCF plotted with respect to the radial velocity measurements (the mean RV has been subtracted). The red points indicate contamination by the Moon. The uncertainty on the bisector span of each point is twice the uncertainty on the RV for that point. 
3. Host star
To analyze the stellar atmosphere we used only spectra that were not affected by moonlight. The only spectrum at an S/N< 10 at 550 nm was discarded as well. The remaining twelve spectra were shifted according to their measured RVs, corrected for the cosmic rays neglected by the SOPHIE pipeline, coadded, and normalized. The final S/N in the continuum, at 550 nm, is ≃ 130 per resolution element.
We measured the stellar effective temperature T_{eff}, surface gravity log g, metallicity [Fe/H], and projected rotational velocity v sin i_{⋆} with the VWA software (Bruntt et al. 2010a,b, and references therein). This method is based on the fit of the metal spectral lines, especially the iron lines. The best parameters are those that minimize the correlation of the element abundances with the excitation potential and the equivalent width of the spectral lines. We obtained T_{eff} = 6260 ± 80 K, log g = 4.40 ± 0.11, and [ Fe / H ] = 0.10 ± 0.13, appropriate of an F8Vtype star.
The estimate of log g was confirmed with the pressuresensitive lines of CaI at 612.2 nm and the MgIb triplet. Finally, the couple v sin i_{⋆} and v_{macro} was jointly measured by fitting a rotational profile on a set of isolated spectral lines. This measure of v sin i_{⋆} (6 ± 2 km s^{1}) agrees with the one obtained with the fit of the CCF (Boisse et al. 2010): 6.8 ± 1.0 km s^{1}.
A first combined fit of the data sets (Sect. 4.2.2) showed a ~3σ difference between the spectroscopic log g (4.40 ± 0.11) and the one derived from the posteriors of the stellar parameters (4.102 ± 0.019). We compared our spectroscopic parameters with the other published ones. Everett et al. (2013) observed Kepler117 at the National Optical Astronomy Observatory (NOAO) Mayall 4m telescope on Kitt Peak with the RCSpec longslit spectrograph. They reported the result of two fits on the spectrum, obtaining log g = 4.26 ± 0.15 and 4.65 ± 0.15 (reaching, in this case, a parameter limit in their model). Rowe et al. (2014), instead, used the publicly available spectra recorded with the HIRES spectrograph at the Keck I telescope and found log g = 4.187 ± 0.150. This star has been observed with other telescopes as well, but the resulting parameters are published only in the Kepler Community Followup Observing Program (CFOP) online archive^{3} and not in the literature, so that we did not consider them. The log g is known to be a problematic parameter to measure accurately and is correlated with T_{eff} and [Fe/H]. In particular, a decrease in log g is usually reflected by a decrease in T_{eff} and in [Fe/H]. The complete set of the three parameters in the articles we referred to is reported in Table 2. The stellar densities derived from the SOPHIE spectrum, the HIRES spectrum (both calculated with the Dartmouth tracks), and those from the TTVs are shown in Fig. 2.
We were unable to identify a problem in the SOPHIE spectra or in our analysis method. We therefore chose to use the published combination of T_{eff}, log g and [Fe/H] whose log g is the closest to our posterior, which converges to a sharp distribution even in the tail of the large spectroscopic prior of SOPHIE. The final values we adopted are those of Rowe et al. (2014), which were used as priors in the Bayesian analysis.
Published spectroscopic parameters for Kepler117 compared with those of this work.
Fig. 2
Stellar densities derived from the spectroscopic parameters of the SOPHIE spectrum (dashed line) and the HIRES spectrum (Rowe et al. 2014) (continuous line). In red, the posterior distribution from the PASTIS analysis, shaded according to the 1, 2, and 3σ intervals. 
4. System analysis
4.1. Stellar activity
The light curve shows small periodic variations, arguably due to starspots. To identify the periodicities, we removed the transit features and computed the LombScargle periodogram (LSP: Press & Rybicki 1989) of the light curve (Fig. 3), finding a peak at 10.668 ± 0.028 days. The uncertainty is underestimated because it does not take into account the position of the spots on the stellar surface and the differential rotation. The ≃ 11 day periodicity was also isolated by the autocorrelation of the light curve (Fig. 3, bottom panel), precisely, by the main peak at 11.1 ± 1.4 days and the first two multiples at 21.8 ± 1.7 and 32.8 ± 1.9 days. The measurements agree with the stellar rotation period measured from the v sin i_{⋆} (Sect. 3) and the stellar radius R_{⋆} (Sect. 4.2.2), assuming that the rotation axis is perpendicular to the line of sight (11.6 ± 1.8 days). In conclusion, the peaks in the periodogram and the autocorrelation function can be considered as representative of the rotation of the host star, for which we conservatively adopted the photometric value with the largest uncertainty, that is, P_{⋆} = 11.1 ± 1.4 days. Comparing this with v sin i_{⋆} shows that the stellar inclination is compatible with 90°.
Fig. 3
Top panel: LombScargle periodogram of the light curve after removing the transits. Bottom panel: autocorrelation of the light curve. The green dotted lines indicate the Gaussian fit to the peaks. The red line corresponds to the maximum rotation period of the star deduced by the v sin i_{⋆} and the stellar radius. The red shadowed regions highlight the 1, 2, and 3σ confidence intervals for the rotation period. 
4.2. Combined fit of the system parameters
We performed a combined Bayesian analysis of photometry, radial velocities, stellar parameters, and TTVs. We begin by describing the measurement and modeling of the TTVs and then present the Bayesian analysis.
4.2.1. Transit time variations
After reducing the photometric data (Sect. 2.1), we fitted the transit times with a procedure similar to the one discussed in Barros et al. (2011). All the transits were fitted simultaneously to constrain the shape parameters, that is, the normalized separation of the planet a/R_{⋆}, the ratio of planettostar radius R_{p}/R_{⋆}, and the orbital inclination i. For each transit, the primary transit epoch T_{0} and three normalization parameters were also fitted to account for a quadratic trend with time. The derived transit times for each planet are given in Tables A.1 and A.2. After removing a linear ephemeris, the transit times of the two planets showed significant TTVs. The TTVs exhibited by planet b are ~4 times longer than those of planet c (≃ 28 min against ≃ 7 min). The TTV amplitude is proportional to the period of the perturbed planet and to the mass of the perturbing one (Agol et al. 2005; Holman & Murray 2005). Therefore, if the two planets had a similar mass, we would expect the outer one to show stronger TTVs. This indicates, then, that the outer planet is the heaviest.
The transit times were shifted to the mean ephemeris, and each transit was normalized using the derived normalization coefficients. To model the TTVs, we performed dynamical simulations with the mercury code, version 6.2 (Chambers 1999). The integrations were executed with a BulirschStoerwe algorithm. For each of them, we identified the transit times by interpolating the passage of the planets through the line of sight. We computed the TTVs by subtracting a linear fit to the transit times. As a compromise between execution time and accuracy of the TTVs with respect to the measured uncertainties, we set the simulations to cover the time span of the Kepler photometry, with a step of 0.4 days, that is, 1 / 47th of the lower orbital period.
Fig. 4
Top: phasefolded plot of the best transit model of planet b (left) and c (right), over the SC data. In black the model, in red the data binned every hundredth of orbital phase. Center: the same for the radial velocities. Bottom: the TTVs of planet b folded at the orbital period of planet c (left) and those of planet c folded at the first peak of its LombScargle periodogram (right, Sect. 5.1). For each plot, the lower panel shows the residuals as observed minus calculated (O−C) points. 
The simulations were included in the Bayesian fit that is described in the next paragraph.
4.2.2. Bayesian combined fit
The Kepler photometry, reduced as described in Sect. 2.1 and corrected for the TTVs as explained in Sect. 4.2.1, and the SOPHIE radial velocities were fitted together using the planet analysis and small transit investigation software (PASTIS) described in Díaz et al. (2014). This software has been primarily designed to calculate the Bayesian odds ratios between competing scenarios in planetary validation problems. PASTIS allows simultaneously modeling of several data sets and obtaining samples of the parameter posteriors with a Markov chain Monte Carlo (MCMC) algorithm.
An exploration phase was started at random points drawn from the priors listed in Table A.3. From the chains computed in this phase, we used the one with the highest likelihood for the starting values of the final MCMC set. The solution of the exploration phase with the highest likelihood has the lowest eccentricities.
To take into account the differences between the stellar models, we used four evolutionary tracks as input for the stellar parameters: Dartmouth (Dotter et al. 2008), PARSEC (Bressan et al. 2012), StarEvol (Palacios, priv. com.; Lagarde et al. 2012), and Geneva (Mowlavi et al. 2012). However, the intrinsic uncertainties in the models were not taken into account. We ran twentyfive chains of about 10^{5} steps for each of the stellar evolutionary tracks. At each step of the MCMC, the model light curves were oversampled and then binned by a factor 10, to correct for the distortions in the signal due to the finite integration time (Kipping 2010). We derived the stellar density ρ_{⋆} from the spectroscopic T_{eff}, log g, and [Fe/H] (Sect. 3) and set it, together with the spectroscopic T_{eff} and [Fe/H], as a jump parameter, with normal priors for all three of them. For each planet, we used Gaussian priors for the period P and the primary transit epoch T_{0} and noninformative priors for the argument of periastron ω, the inclination i, and the eccentricity e. We stress this last point: without imposing zero eccentricities, we consistently measured these key parameters by taking into account all the available information, TTVs included.
We used uniform priors for the coefficients of quadratic limbdarkening, for the planetarytostellar radius ratio R_{p}/R_{⋆}, and for the radial velocity amplitude K. For Kepler117 b, we fitted the longitude of the ascending node Ω, too, for which we imposed a uniform prior. The Ω of planet c was fixed at 180° because the symmetry of the problem allows freely choosing one of the two Ωs.
Planetary and stellar parameters with their 68.3% central confidence intervals.
We expressed the Kepler normalized flux offset, the systemic velocity, and the RV linear drift with uniform priors (separating LC and SC data in the photometry). Finally, we modeled the instrumental and astrophysical systematic sources of error with a jitter term for Kepler, two for the TTVs (one for each planet), and one for SOPHIE. A uniform prior was assigned to all the jitter terms.
After they were sampled, every posterior distribution was thinned according to its correlation length. A combined posterior distribution was derived by taking the same number of points from each stellar evolutionary track. This combined distribution gave the most probable values and the confidence intervals for the system parameters.
Finally, the derived log g and the posterior stellar radius R_{⋆}, T_{eff}, and [Fe/H], together with the magnitudes in Table A.4, were set as priors for another MCMC run to derive the distance of the system using the spectral energy distribution (SED). The magnitudes were fitted to sample the posterior distributions of the distance of the system, the interstellar extinction E(B − V), and the jitter of the SED. The model SED was interpolated from the PHOENIX/BTSettl synthetic spectral library (Allard et al. 2012), scaled with the distance, the stellar radius, and the reddening E(B − V), expressed through a Fitzpatrick (1999) extinction law. For both the distance and the reddening, noninformative priors were used.
Fig. 5
Model SED on the photometric bands. The residuals are shown in the lower panel. 
In Table 3 we present the mode and the 68.3% equaltailed confidence intervals of the system parameters. According to our analysis, Kepler117 A is a ≃ 5 Gyr old F8V star with two planets in loweccentricity orbits (0.0493 ± 0.0062 and 0.0323 ± 0.0033 for planet b and c), which differ widely in their mass, but less so in their radii: 0.094 ± 0.033M_{J}, 0.719 ± 0.024R_{J} for planet b and 1.84 ± 0.18M_{J}, 1.101 ± 0.035R_{J} for planet c. The planetary radii, in particular, agree with the estimate of Rowe et al. (2014): ≃ 0.72 ± 0.14R_{J} for planet b and ≃ 1.04 ± 0.20R_{J} for planet c.
Fig. 6
LombScargle periodogram of the measured (red) and modeled (black) TTVs for planet b (left) and c (right). The blue triangle in the plot on the left indicates the orbital period of planet c. 
Fig. 7
Upper left: probability density function (PDF) of the stellar density including or excluding the TTVs; to make the plot more readable, the PDF using the TTVs is divided by 2. The prior from spectroscopy (HIRES, Dartmouth) is shown for comparison. Upper right: the planetary eccentricities from different sets of data; the PDF using the TTVs is divided by 12. The color code is the same in the following plots. Central line: planetary masses and radii from the fit with and without TTVs. The PDF of the masses using the TTVs is divided by 5, that of the radii by 2. Bottom: orbital inclinations from different sets of data. 
We remark that the measured drift of the RVs is compatible with 0 km s^{1}; a nonzero drift would have been an indication of a possible third companion in the system that affected the amplitude of the TTVs.
5. Discussion
5.1. Modulation of the TTVs
We tested the robustness of our result by inspecting the periodic modulation of the measured and the modeled TTVs. To do this, we compared their LombScargle periodograms (Fig. 6). The main peaks coincide for both planets and also agree with the periodicities found by Mazeh et al. (2013; see also Ofir et al. 2014). The periodogram of the modeled TTVs of planet b reproduces that of the measured TTVs well. Some of the peaks of planet c, on the other hand, are due to noise. This can be explained by the different amplitude of the signal in the two cases.
5.2. Constraining the system parameters with the TTVs
The RV amplitude produced by Kepler117 b (6.5 ± 2.1 m s^{1}) is close to the sensitivity of SOPHIE for a ≃ 14.5 mag_{V} star. Indeed, the SOPHIE RVs alone do not have the precision required to measure this planet’s mass and can only provide an upper limit. Including the TTVs allowed us to accurately determine this parameter.
The combined fit of the TTVs reduces the uncertainties on other parameters, too. This is because the amplitude of the TTVs is mainly determined by orbital separations, periods, and eccentricities of the orbits and masses ratios of the planets (Agol et al. 2005). The strong constraint on the eccentricities, combined with the constraint on the stellar density (which is determined with a precision similar to the precisions achieved in asteroseismology), reduces the error bars on the planetary parameters. Once again, we stress that the derived uncertainties do not include the uncertainties on all the models of stellar atmospheres and evolutionary tracks we used.
To check the impact of the TTVs on the combined fit, we ran PASTIS without them. The different posteriors of the most affected parameters, with or without the TTVs, are compared in Figs. 7 and 8. The mass of planet b presents the most evident difference because of its poorly constrained RV amplitude: its value reaches from an upper limit (0.28 M_{J} at 68.3% confidence level) to a better constrained value (0.094 ± 0.033M_{J}). The difference is smaller for planet c because the amplitude of the RVs is larger and better fitted. However, its uncertainty is roughly reduced by 40%. This indicates that, if possible, including the TTVs in the combined fit is more effective than fitting them a posteriori, using a set of orbital parameters derived without considering them.
We remark that the mass of planet c found with the RVs alone and with the TTVs are fully compatible. Therefore, the observed TTVs are completely explained by the two planets, within the data error bars. This agrees with the absence of any RV drift (Sect. 4.2.2).
Fig. 8
Massradius diagram for the solutions without (cyan for planet b, magenta for planet c) and with (blue for planet b, red for planet c) TTVs. The blue and red solutions are those indicated in Table 3. The colors, from the center to the edge of the regions, correspond to the 39.3%, 86.5%, and 98.9% joint confidence intervals. Jupiter and Saturn (labeled J and S) are marked for comparison. The other planet parameters were taken from Wright et al. (2011). 
The transit signature is degenerate with respect to the stellar hemisphere the planet covers, while the TTVs are not. In a twoplanet scenario, this can lead to strong correlations between the two inclinations. While running PASTIS, we therefore constrained one of the transits in one of the hemispheres and left the other free to vary. As the inclination of planet b is lower than that of planet c, the inclination of planet b was limited to one hemisphere (50°<i< 90°) and that of planet c was left free to vary between both (89°<i< 91°). In spite of this, our fit allowed both hemispheres to be transited by planet c because the final inclination of its orbit is almost symmetric with respect to the stellar equator (Fig. 7, bottom line). The solutions with i> 90° are compatible with those without TTVs (all < 90°) at 1σ. In particular, for planet c, we found 89.64 ± 0.10° with TTVs and i = 89.75 ± 0.13° without them.
Using the stellar inclination (Sect. 4.1) and the system parameters, we calculated the expected amplitude of the RossiterMcLaughlin effects following Eq. (11) of Gaudi & Winn (2007). For planet b, we found 10.9 ± 3.0 m s^{1}, for planet c 79 ± 13 m s^{1}. Measuring the spinorbit misalignment would then be possible for planet c, but the transit duration (≃ 11 h) would require a joint effort from different locations to cover a whole transit.
The difference between the resulting longitude of the ascending node Ω for planet b (177.9 ± 5.6°) and that of planet c (fixed to 180°) is compatible with 0°. Combined with the similar inclinations, this implies two almost coplanar orbits. As most of the Kepler planetary systems (Fabrycky et al. 2014), Kepler117 clearly has a flat configuration of the orbits.
We ran an MCMC set without the RVs to determine the reliability of the fit. As expected, the posterior distributions are the same as with the RVs. Systems with lowmass planets presenting TTVs, which are challenging for the RV observations, would benefit from the approach used in this paper.
Finally, we verified that the configuration of the most probable solution is dynamically stable. We ran mercury over a time span of 10 Myr (Fig. 9). The semimajor axes, eccentricities, and orbital inclinations oscillate over a time scale of around 200 years, but all the parameters are stable in the long term.
5.3. A third nontransiting companion around Kepler117?
The possibility of a third nontransiting companion can be probed with the RVs and the TTVs. As already mentioned in Sect. 4.2.2, the absence of stellar drift in the RVs brings no evidence of the possible presence of a third nontransiting planet in the system. Moreover, the agreement between the mass of planet c, found with the RVs and with or without the TTVs (Sect. 5.2), shows that the TTVs are not affected by a nontransiting body.
A more precise constraint can be obtained by subtracting the modeled RVs of the two planets from the RV measurements. We folded the residuals for several periods and fitted them with a sinusoid. The amplitude of the sinusoid and the mass of the star (Table 3) allow extracting the maximum mass of the possible companion.
The result is plotted in Fig. 10 (filled regions) for the 68.27%, 95.45%, and 99.73% confidence intervals. The RVs allow the presence of a Jupitermass planet for some orbital periods. The TTVs, however, impose a stronger constraint, since including a third body with the combination of mass and period allowed by the RVs (with the simplifying assumptions of a circular orbit and 90° inclination) would not fit the TTVs. The black line in the plot represents a 3σ difference in the residuals between the fit of the TTVs with a third planet and the best solution with two.
Fig. 9
Evolution of semimajor axes (top), eccentricities (center), and orbital inclinations (bottom) over a 10 Myr simulation of the most probable solution. The respective mean has been subtracted from the two semimajor axes. On the left, a zoom on the first 300 yr; on the right, the variation intervals of the parameters. In blue planet b, in red planet c. The shaded regions in the left panel of the inclinations correspond to the values resulting in a transit (see Eq. (7) of Winn 2010). 
Therefore, under some simplifying assumptions, the presence of a nondetected third companion above ~0.1M_{J} on an orbit shorter than ~100 days, as well as that of a giant companion with an orbit shorter than ~250 days, is very unlikely.
6. Summary and conclusions
We presented the combined analysis of the Kepler photometry, the TTVs, the SOPHIE RVs, and the spectroscopic observations of Kepler117. This allowed us to measure the stellar, planetary, and orbital parameters of the system. According to our analysis, Kepler117 A is an F8type mainsequence star that is about 5 Gyr old. The system is composed of two planets. Kepler117 c, the outer one, has a period of ≃ 51 days and a ~2M_{J} mass. Kepler117 b, the inner one, has a period of ≃ 19 days and a ~30 M_{Earth} mass. The latter produces a RV semiamplitude of 6.5 ± 2.1 m s^{1}, close to the limit of sensitivity of SOPHIE for faint magnitudes. Therefore, its mass and eccentricity cannot be obtained from the RVs alone, even though the derived upper limits confirm that the transiting body belongs to the planet realm.
Our analysis shows that the inclusion of the TTVs in the combined fit allows tightly constraining the mass of the lighter planet. Taking into account the TTVs in the fit also allows a better determination of the other system parameters. The stellar density was tightly constrained despite a loose prior on the spectroscopic log g. The planetary radii were strongly constrained, as were the orbital eccentricities, even if small (around ~0.03−0.05 for both orbits). Measuring the eccentricity accurately is important for testing the dynamical models of young systems with giant planets. Simulations show that the complex evolution of systems with two planets can be the result of ejection or merging in systems with three planets and can lead to stable, resonant, and loweccentricity orbits (e.g., Lega et al. 2013).
Fig. 10
Maximum mass of a third possible companion as a function of the orbital period, based on the RV (orange filled areas) and TTV (blue line) observations. The color, from darker to lighter, corresponds to the 68.27%, 95.45%, and 99.73% confidence intervals from the RV data. The blue line represents a 3σ difference in the residuals on the TTVs from the best solution with two planets. 
While the RV and TTVmeasured mass agree for planet c, the same comparison for planet b would benefit from a spectrograph with higher sensitivity than SOPHIE. Only a small part of the planets known to date have their mass measured with both RVs and TTVs (e.g., Barros et al. 2014). As observed by Weiss & Marcy (2014), planets smaller that four Earth radii with a mass measured with TTVs are systematically lower in the massradius diagram than those discovered by RV surveys. This could be due to nondetected companions that might dampen the TTVs, causing a systematic underestimation of the masses, or to a lower density of the planets that show TTVs. Indeed, in the compact multiplanetary systems that are likely to produce observable TTVs, planets with lower masses for a given size are more likely to reach stable orbits (JontofHutter et al. 2014). In this context, it is remarkable that the RV and the TTVmeasured mass of Kepler117 c agree. In addition, under some simplifying assumptions, the TTVs almost exclude a nondetected ~0.1M_{J} companion with an orbit shorter than ~100 days, as well as a giant companion with an orbit shorter than ~250 days.
We cannot exclude that while in this particular case the conditions are fulfilled for the combined fit to be effective, this is not the case in general. In fact, TTVs with sufficiently high amplitude are necessary. This system exhibits significant TTVs even though the orbital period ratio of the planets is far from an exact loworder mean motion resonance, for which strong TTVs are expected (e.g., Lithwick & Wu 2012). The orbital period ratio between planet c and b (≃ 2.7) places this system on the wide side of the 5:2 mean motion resonance. The overabundance of systems with period ratios some percent higher than resonant values than those with a ratio slightly lower than these ones has been well established for the Kepler systems in the case of firstorder resonances, that is, 2:1 or 3:2 (e.g., see Lissauer et al. 2011). The explanation is given in terms of tidal dissipation related to diskplanet or starplanet interactions for closein orbits (Lithwick & Wu 2012; Batygin & Morbidelli 2013; Delisle et al. 2014), causing the orbital periods to diverge. To date, this pilingup appears only for some firstorder resonances in the Kepler systems. Instead, the rest of the periodratio distribution, including higher order resonances, remains flat (e.g., Batygin & Morbidelli 2013; Fabrycky et al. 2014).
We verified that a system with the most probable solution of our analysis is dynamically stable. However, we noted the eccentricities and the inclinations show small oscillations, that do not affect the stability of the system. We found the planetary orbits to be almost coplanar. This places Kepler117 in the most common population of the Kepler multiplanetary systems with a flat configuration, as highlighted by Fabrycky et al. (2014).
In conclusion, a deeper understanding of the dynamics of orbital resonances is needed to better reconstruct the history of Kepler117, which adds valuable information to our knowledge of multiplanetary systems.
Acknowledgments
This paper includes data collected by the Kepler mission. Funding for the Kepler mission is provided by the NASA Science Mission directorate. We made use of the Mikulski Archive for Space Telescopes (MAST). Support for MAST for nonHST data is provided by the NASA Office of Space Science via grant NNX09AF08G and by other grants and contracts. We thank the technical team at the Observatoire de HauteProvence for their support with the SOPHIE instrument and the 1.93 m telescope and in particular for the essential work of the night assistants. Financial support for the SOPHIE observations comes from the Programme National de Planetologie (PNP) of CNRS/INSU, France is gratefully acknowledged. We also acknowledge support from the French National Research Agency (ANR08 JCJC010201). The team at LAM acknowledges support by CNES grants 98761 (SCCB), 426808 (CD), and 251091 (JMA). A.S. acknowledge the support from the European Research Council/European Community under the FP7 through Starting Grant agreement number 239953. A.S. is supported by the European Union under a Marie Curie IntraEuropean Fellowship for Career Development with reference FP7PEOPLE2013IEF, number 627202. A.S.B. acknowledges funding from the European Union Seventh Framework Programme (FP7/2007–2013) under Grant agreement No. 313014 (ETAEARTH). We thank John Chambers for his explanations about the use of mercury and Rosemary Mardling for the fruitful discussions about the dynamic of threebody systems. This research was made possible through the use of data from different surveys: the AAVSO Photometric AllSky Survey (APASS), funded by the Robert Martin Ayers Sciences Fund; the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation; the Widefield Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. This research has made reference to the Exoplanet Orbit Database and the Exoplanet Data Explorer at exoplanets.org.
References
 Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567 [NASA ADS] [CrossRef] [Google Scholar]
 Allard, F., Homeier, D., & Freytag, B. 2012, Roy. Soc. London Philos. Trans. Ser. A, 370, 2765 [Google Scholar]
 Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barros, S. C. C., Pollacco, D. L., Gibson, N. P., et al. 2011, MNRAS, 416, 2593 [NASA ADS] [CrossRef] [Google Scholar]
 Barros, S. C. C., Boué, G., Gibson, N. P., et al. 2013, MNRAS, 430, 3032 [NASA ADS] [CrossRef] [Google Scholar]
 Barros, S. C. C., Díaz, R. F., Santerne, A., et al. 2014, A&A, 561, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Batygin, K., & Morbidelli, A. 2013, AJ, 145, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Boisse, I., Eggenberger, A., Santos, N. C., et al. 2010, A&A, 523, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 728, 117 [NASA ADS] [CrossRef] [Google Scholar]
 Bouchy, F., Hébrard, G., Udry, S., et al. 2009, A&A, 505, 853 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bouchy, F., Bonomo, A. S., Santerne, A., et al. 2011, A&A, 533, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bouchy, F., Díaz, R. F., Hébrard, G., et al. 2013, A&A, 549, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Bruntt, H., Bedding, T. R., Quirion, P.O., et al. 2010a, MNRAS, 405, 1907 [NASA ADS] [Google Scholar]
 Bruntt, H., Deleuil, M., Fridlund, M., et al. 2010b, A&A, 519, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chambers, J. E. 1999, MNRAS, 304, 793 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
 Dawson, R. I., Johnson, J. A., Fabrycky, D. C., et al. 2014, ApJ, 791, 89 [NASA ADS] [CrossRef] [Google Scholar]
 Delisle, J.B., Laskar, J., & Correia, A. C. M. 2014, A&A, 566, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Díaz, R. F., Almenara, J. M., Santerne, A., et al. 2014, MNRAS, 441, 983 [NASA ADS] [CrossRef] [Google Scholar]
 Dotter, A., Chaboyer, B., Jevremović, D., et al. 2008, ApJS, 178, 89 [NASA ADS] [CrossRef] [Google Scholar]
 Everett, M. E., Howell, S. B., Silva, D. R., & Szkody, P. 2013, ApJ, 771, 107 [NASA ADS] [CrossRef] [Google Scholar]
 Fabrycky, D. C., Lissauer, J. J., Ragozzine, D., et al. 2014, ApJ, 790, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Fitzpatrick, E. L. 1999, PASP, 111, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Gaudi, B. S., & Winn, J. N. 2007, ApJ, 655, 550 [NASA ADS] [CrossRef] [Google Scholar]
 Hébrard, G., Bouchy, F., Pont, F., et al. 2008, A&A, 488, 763 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Holman, M. J., Fabrycky, D. C., Ragozzine, D., et al. 2010, Science, 330, 51 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Howard, A. W., Marcy, G. W., Johnson, J. A., et al. 2010, Science, 330, 653 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Jenkins, J. M., Caldwell, D. A., Chandrasekaran, H., et al. 2010, ApJ, 713, L87 [NASA ADS] [CrossRef] [Google Scholar]
 JontofHutter, D., Lissauer, J. J., Rowe, J. F., & Fabrycky, D. C. 2014, ApJ, 785, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Kipping, D. M. 2010, MNRAS, 408, 1758 [NASA ADS] [CrossRef] [Google Scholar]
 Kipping, D. M. 2014, MNRAS, 440, 2164 [NASA ADS] [CrossRef] [Google Scholar]
 Lagarde, N., Decressin, T., Charbonnel, C., et al. 2012, A&A, 543, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lega, E., Morbidelli, A., & Nesvorný, D. 2013, MNRAS, 431, 3494 [NASA ADS] [CrossRef] [Google Scholar]
 Lissauer, J. J., Ragozzine, D., Fabrycky, D. C., et al. 2011, ApJS, 197, 8 [NASA ADS] [CrossRef] [Google Scholar]
 Lissauer, J. J., Marcy, G. W., Rowe, J. F., et al. 2012, ApJ, 750, 112 [NASA ADS] [CrossRef] [Google Scholar]
 Lissauer, J. J., Marcy, G. W., Bryson, S. T., et al. 2014, ApJ, 784, 44 [NASA ADS] [CrossRef] [Google Scholar]
 Lithwick, Y., & Wu, Y. 2012, ApJ, 756, L11 [NASA ADS] [CrossRef] [Google Scholar]
 Mazeh, T., Nachmani, G., Holczer, T., et al. 2013, ApJS, 208, 16 [NASA ADS] [CrossRef] [Google Scholar]
 Mowlavi, N., Eggenberger, P., Meynet, G., et al. 2012, A&A, 541, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nesvorný, D., Kipping, D., Terrell, D., et al. 2013, ApJ, 777, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Ofir, A., Dreizler, S., Von Essen, C., & Aharonson, O. 2014, CoRoT3KASC7 Symposium: The Space Photometry Revolution, eds. J. Ballot, & R. A. Garcia, EPJ Web Conf., in press [Google Scholar]
 Oshagh, M., Santos, N. C., Boisse, I., et al. 2013, A&A, 556, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Perruchot, S., Kohler, D., Bouchy, F., et al. 2008, in SPIE Conf. Ser., 7014 [Google Scholar]
 Pollacco, D., Skillen, I., Collier Cameron, A., et al. 2008, MNRAS, 385, 1576 [NASA ADS] [CrossRef] [Google Scholar]
 Press, W. H., & Rybicki, G. B. 1989, ApJ, 338, 277 [NASA ADS] [CrossRef] [Google Scholar]
 Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rowe, J. F., Bryson, S. T., Marcy, G. W., et al. 2014, ApJ, 784, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Santerne, A., Díaz, R. F., Moutou, C., et al. 2012, A&A, 545, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santerne, A., Hébrard, G., Deleuil, M., et al. 2014, A&A, 571, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
 Steffen, J. H., Batalha, N. M., Borucki, W. J., et al. 2010, ApJ, 725, 1226 [NASA ADS] [CrossRef] [Google Scholar]
 Weiss, L. M., & Marcy, G. W. 2014, ApJ, 783, L6 [NASA ADS] [CrossRef] [Google Scholar]
 Winn, J. N. 2010, Exoplanets (Tucson: University of Arizona Press) [arXiv:1001.2010] [Google Scholar]
 Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
 Wright, J. T., Fakhouri, O., Marcy, G. W., et al. 2011, PASP, 123, 412 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Appendix A: Additional figures and tables
Transittiming variations for planet b.
Transittiming variations for planet c.
Prior distributions used in the combined fit with PASTIS.
Magnitudes of the Kepler117 system.
Fig. A.1
TTVs of planet b (left) and c (right) as a function of time. In red the data, in black the fit. The lower panels show the residuals. 
Fig. A.2
Pyramid of the combined posterior distributions. 
All Tables
Published spectroscopic parameters for Kepler117 compared with those of this work.
All Figures
Fig. 1
Bisector span of the CCF plotted with respect to the radial velocity measurements (the mean RV has been subtracted). The red points indicate contamination by the Moon. The uncertainty on the bisector span of each point is twice the uncertainty on the RV for that point. 

In the text 
Fig. 2
Stellar densities derived from the spectroscopic parameters of the SOPHIE spectrum (dashed line) and the HIRES spectrum (Rowe et al. 2014) (continuous line). In red, the posterior distribution from the PASTIS analysis, shaded according to the 1, 2, and 3σ intervals. 

In the text 
Fig. 3
Top panel: LombScargle periodogram of the light curve after removing the transits. Bottom panel: autocorrelation of the light curve. The green dotted lines indicate the Gaussian fit to the peaks. The red line corresponds to the maximum rotation period of the star deduced by the v sin i_{⋆} and the stellar radius. The red shadowed regions highlight the 1, 2, and 3σ confidence intervals for the rotation period. 

In the text 
Fig. 4
Top: phasefolded plot of the best transit model of planet b (left) and c (right), over the SC data. In black the model, in red the data binned every hundredth of orbital phase. Center: the same for the radial velocities. Bottom: the TTVs of planet b folded at the orbital period of planet c (left) and those of planet c folded at the first peak of its LombScargle periodogram (right, Sect. 5.1). For each plot, the lower panel shows the residuals as observed minus calculated (O−C) points. 

In the text 
Fig. 5
Model SED on the photometric bands. The residuals are shown in the lower panel. 

In the text 
Fig. 6
LombScargle periodogram of the measured (red) and modeled (black) TTVs for planet b (left) and c (right). The blue triangle in the plot on the left indicates the orbital period of planet c. 

In the text 
Fig. 7
Upper left: probability density function (PDF) of the stellar density including or excluding the TTVs; to make the plot more readable, the PDF using the TTVs is divided by 2. The prior from spectroscopy (HIRES, Dartmouth) is shown for comparison. Upper right: the planetary eccentricities from different sets of data; the PDF using the TTVs is divided by 12. The color code is the same in the following plots. Central line: planetary masses and radii from the fit with and without TTVs. The PDF of the masses using the TTVs is divided by 5, that of the radii by 2. Bottom: orbital inclinations from different sets of data. 

In the text 
Fig. 8
Massradius diagram for the solutions without (cyan for planet b, magenta for planet c) and with (blue for planet b, red for planet c) TTVs. The blue and red solutions are those indicated in Table 3. The colors, from the center to the edge of the regions, correspond to the 39.3%, 86.5%, and 98.9% joint confidence intervals. Jupiter and Saturn (labeled J and S) are marked for comparison. The other planet parameters were taken from Wright et al. (2011). 

In the text 
Fig. 9
Evolution of semimajor axes (top), eccentricities (center), and orbital inclinations (bottom) over a 10 Myr simulation of the most probable solution. The respective mean has been subtracted from the two semimajor axes. On the left, a zoom on the first 300 yr; on the right, the variation intervals of the parameters. In blue planet b, in red planet c. The shaded regions in the left panel of the inclinations correspond to the values resulting in a transit (see Eq. (7) of Winn 2010). 

In the text 
Fig. 10
Maximum mass of a third possible companion as a function of the orbital period, based on the RV (orange filled areas) and TTV (blue line) observations. The color, from darker to lighter, corresponds to the 68.27%, 95.45%, and 99.73% confidence intervals from the RV data. The blue line represents a 3σ difference in the residuals on the TTVs from the best solution with two planets. 

In the text 
Fig. A.1
TTVs of planet b (left) and c (right) as a function of time. In red the data, in black the fit. The lower panels show the residuals. 

In the text 
Fig. A.2
Pyramid of the combined posterior distributions. 

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.