arxiv: v1 [astro-ph.ep] 7 Jan 2015

Transcription

1 Astronomy & Astrophysics manuscript no. koi c ESO 2018 April 9, 2018 SOPHIE velocimetry of Kepler transit candidates. XV. KOI-614b, KOI-206b, and KOI-680b: a massive warm Jupiter orbiting a G0 metallic dwarf and two highly inflated planets with a distant companion around evolved F-type stars. J.M. Almenara 1,2,3, C. Damiani 1, F. Bouchy 1, M. Havel 4, G. Bruno 1, G. Hébrard 5,6, R.F. Diaz 1,7, M. Deleuil 1, S.C.C. Barros 1, I. Boisse 1, A.S. Bonomo 8, G. Montagnier 5,6, and A. Santerne 9,10,1 arxiv: v1 [astro-ph.ep] 7 Jan Aix Marseille Université, CNRS, LAM (Laboratoire d Astrophysique de Marseille) UMR 7326, 13388, Marseille, France 2 Univ. Grenoble Alpes, IPAG, F Grenoble, France 3 CNRS, IPAG, F Grenoble, France 4 NASA Postdoctoral Program Fellow, Ames Research Center, P.O. Box 1, Moffett Field, CA Observatoire de Haute Provence, Saint Michel l Observatoire, France 6 Institut d Astrophysique de Paris, UMR7095 CNRS, Université Pierre & Marie Curie, 98bis boulevard Arago, Paris, France 7 Observatoire Astronomique de l Université de Genève, 51 chemin des Maillettes, 1290 Versoix, Switzerland 8 INAF - Osservatorio Astrofisico di Torino, via Osservatorio 20, Pino Torinese, Italy 9 Centro de Astrofísica, Universidade do Porto, Rua das Estrelas, Porto, Portugal 10 Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, PT Porto, Portugal ABSTRACT We report the validation and characterization of three new transiting exoplanets using SOPHIE radial velocities: KOI-614b, KOI-206b, and KOI-680b. KOI-614b has a mass of 2.86±0.35 M Jup and a radius of R Jup, and it orbits a G0, metallic ([Fe/H]=0.35±0.15) dwarf in 12.9 days. Its mass and radius are familiar and compatible with standard planetary evolution models, so it is one of the few known transiting planets in this mass range to have an orbital period over ten days. With an equilibrium temperature of T eq = 1000 ± 45 K, this places KOI-614b at the transition between what is usually referred to as hot and warm Jupiters. KOI-206b has a mass of 2.82 ± 0.52 M Jup and a radius of 1.45 ± 0.16 R Jup, and it orbits a slightly evolved F7-type star in a 5.3-day orbit. It is a massive inflated hot Jupiter that is particularly challenging for planetary models because it requires unusually large amounts of additional dissipated energy in the planet. On the other hand, KOI-680b has a much lower mass of 4 ± 0.15 M Jup and requires less extra-dissipation to explain its uncommonly large radius of 1.99 ± 0.18 R Jup. It is one of the biggest transiting planets characterized so far, and it orbits a subgiant F9-star well on its way to the red giant stage, with an orbital period of 8.6 days. With host stars of masses of 1.46 ± 0.17 M and 1.54 ± 9 M, respectively, KOI-206b, and KOI-680b are interesting objects for theories of formation and survival of short-period planets around stars more massive than the Sun. For those two targets, we also find signs of a possible distant additional companion in the system. Key words. stars: planetary systems techniques: photometry techniques: radial velocities techniques: spectroscopic 1. Introduction Since 2010, we have been conducting a radial-velocity followup of Kepler (Borucki et al. 2010) transiting planetary candidates with the SOPHIE spectrograph (Perruchot et al. 2008; Bouchy et al. 2013). Transits combined with radial velocity measurements enable us to infer the masses and radii of the star and planet up to a one-parameter degeneracy (Wright & Gaudi 2013). To derive the absolute planetary parameters, an independent stellar parameter is needed, which can be estimated with the aid of stellar models. With the mass and radius of the planet, we can determine the mean density that gives insight into the composition and formation of the objects. In addition, such a survey allows the rate of false-positive of the Kepler close-in giant exoplanet candidates to measured (Santerne et al. 2012). In this paper, we present the validation and characterization of three new transiting exoplanets: KOI-614b, KOI-206b, and Based on observations made with SOPHIE on the 1.93-m telescope at the Observatoire de Haute-Provence (CNRS), France. KOI-680b. Two of them, KOI-206b and KOI-680b, have already been announced by our group to be planetary companions in Santerne et al. (2012). The third one, KOI-614b, was not in the sample considered by Santerne et al. (2012) owing to its transit depth of 0.39%, slightly lower than the % selection cut-off. The structure of the paper is as follows. We present the Kepler and SOPHIE observations (Section 2), the derivation of the spectral parameters (Section 3), the stellar and planetary parameter estimation (Section 4), the planetary evolution models (Section 5), and finally the discussion (Section 6) and the conclusions (Section 7) of our work. 2. Observations 2.1. Photometric detection with Kepler The three targets analyzed in this paper were identified as KOI (Kepler Object of Interest) by Borucki et al. (2011). Their IDs, coordinates, and magnitudes are reported in Table 2. Two of them were found to be hosting a transiting companion of plane- 1

2 J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates tary nature with high probability, namely KOI-206 and KOI-680 (vetting flag of 2), whereas KOI-614 was considered a moderate probability candidate because of the V-shaped transit (vetting flag of 3). There is no significant spatial shift at the time of the transits in the Kepler photometric centroid for the three targets. No other transits with different periods were detected in any of the light curves, so there are no signs of multiple transiting systems. At the time of the beginning of the analysis, the Kepler light curves of the first 15 quarters are publicly available at the MAST archive 1. Long cadence (29.4 minutes sampling) data only is available for KOI-614, but KOI-206 and KOI-680 have short-cadence data (1-minute sampling) available in 3 of the 15 first quarters. We used the light curve of quarters Q1 to Q15 calibrated through the Kepler Science Pipeline and corrected for instrumental artifacts and other known systematics errors (Jenkins et al. 2010). The three light curves present transits with depths of about 0.5% and a typical uncertainty for individual points at the level of ppm for the long-cadence data. Besides the transits, the three light curves show signs of stellar variability with flux modulations of about 90, 75, and 45 ppm amplitudes for KOI-206, KOI-614, and KOI-680, respectively. As a result, we computed their Lomb-Scargle periodograms (Press & Rybicki 1989) using the long-cadence sampling after subtracting the transits (Fig. 1, top panels). They all present peaks compatible with the rotational period of the stars as expected from their measured projected rotational velocities (v sin i, Section 3) and radii (Section 4), which give upper limits at 10.5 ± 1.5, , and 27.3 ± 5.7 days for KOI-206, KOI- 614, and KOI-680, respectively. Given the amplitude of the variations and their period, we conclude that those peaks are likely due to spots covering the stellar surface. The evolution of the spots, which vary in number and position with time, may scramble the periodicity of the signal, so we divided the light curves in chunks of the same duration, chosen to span over about five expected rotation period, and we again computed the Lomb- Scargle in each one. For each light curve, the resulting mean periodograms were plotted in Fig. 1 (middle panels). We also give the autocorrelation functions computed on the whole light curve (Fig. 1, bottom panels). The main period and its harmonics are clearly detectable for the three targets. All methods are in good agreement, and we estimate stellar rotational periods of 7.95 ± 0.90, ± 0.10, and ± 3 days for KOI-206, KOI-614, and KOI-680, respectively. Finally, owing to the high level of noise at high frequency, we cannot detect any evidence of pulsations in the expected range for those stars (Belkacem et al. 2011) Radial velocities with SOPHIE The spectroscopic follow-up of the three targets was performed with the SOPHIE spectrograph (Perruchot et al. 2008) after its upgrade made in 2011 (Bouchy et al. 2013). The three Kepler targets were observed from July 2011 to August 2012 with a few additional measurements made in September 2013 for KOI-680. To increase the throughput for these faint targets, the high-efficiency mode of SOPHIE was used (resolving power R ), as well as the slow read-out mode of the detector. The exposure times ranged from 15 to 60 min depending on the targets and the weather conditions, in order to keep the signal-tonoise ratio as constant as possible for a given target. While the first fiber pointed at the target, the sky background was acquired 1 search/search.php Table 1. SOPHIE radial velocity measurements. BJD RV ±1 σ bisect. exp. S/N [km s 1 ] [km s 1 ] [km s 1 ] [s] KOI KOI KOI : bisector spans; associated error bars are twice those of RVs. : signal-to-noise ratio per pixel at 550 nm. : measurements corrected for moonlight pollution. using the second fiber aperture of SOPHIE pointed at 2 arcmin from the target. The spectra were extracted using the SOPHIE pipeline (Bouchy et al. 2009), which includes localization of the spectral orders on the 2D-images, optimal order extraction, cosmic-ray rejection, wavelength calibration, and corrections of flat field. They were corrected for the spectrograph drifts, measured from thorium calibration exposures every two to three hours during the night. Radial velocities were obtained from a Gaussian fit to the weighted cross-correlation function (CCF) with F0-type numerical masks for the three targets (Baranne et al. 1996; Pepe et al. 2002) and are reported in Table 1. Some spectra were contaminated by moonlight and were corrected following the method described in Pollacco et al. (2008) and Hébrard et al. (2008) using the background sky obtained with the second fiber. The corrected exposures are identified in Table 1. Radial velocities were also corrected using the RV reference star HD (Howard et al. 2010; Bouchy et al. 2013; Santerne et al. 2014) 2

3 J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates Period [days] Lomb-Scargle Period [days] Autocorrelation Lomb-Scargle Lomb-Scargle Lomb-Scargle 5 Autocorrelation Autocorrelation Lomb-Scargle Lomb-Scargle Period [days] Fig. 1. From left to right: periodograms of KOI-206, KOI-614, and KOI-680. Top panels: Lomb-Scargle periodogram of the complete light curve. Middle panels: mean Lomb-Scargle periodogram of the light curve divided in chunks. Bottom panel: light curve autocorrelation. The green dashed lines indicate the center of Gaussians fitted to the data to derive the main frequency and its harmonics. The red vertical line corresponds to the maximum rotation period deduced from the v sin i? and the star radius, and the shaded areas are the associated 1, 2, and 3-σ upper limit regions from the darkest to the lighter ones, respectively Bisector span [km s 1 ] Bisector span [km s 1 ] Bisector span [km s 1 ] Radial velocity [km s 1 ] Radial velocity [km s 1 ] Radial velocity [km s 1 ] Fig. 2. Bisector span as a function of the detrended radial velocities with 1-σ error bars for KOI-206, KOI-614, and KOI-680 (from left to right). Points affected and corrected by Moon contamination are indicated in blue. The ranges here have the same extents in the x- and y-axes. which was systematically observed the same nights of the three Kepler targets. The RV dispersion over the two years of measurements on this constant star is about 10 m/s but with some variations at the level of m/s over a few days that seem to be correlated with the dome temperature. The correction by this constant star slightly improve the residual of the RVs, but without any significant change of the parameters. The error bars on the radial velocities were computed from the CCF using the method presented by Boisse et al. (2010). For faint targets such as these, some spectral domains are too noisy and degrade the accuracy of the radial-velocity measurement, therefore the number of spectral orders used in the CCF was adjusted to decrease the dispersion of the measurements. The radial velocity measurements are displayed in the lower plots of Fig. 3, together with their best Keplerian fits and the residuals. They show variations in phase with the transit ephemeris derived from Kepler and with semi-amplitudes implying a companion in the planet-mass range. For KOI-614b, no significant linear drift, mask effect, or correlation between the bissector span (Fig. 2, middle panel) were detected, reinforcing the conclusion that the radial-velocity vari- ations are not caused by spectral-line profile changes attributable to blends or stellar activity. On the other hand, the radial velocities of KOI-206b and KOI-680b show a significant linear drift of 500 ± 200 and 43.9 ± 9.2 m s 1 yr 1, suggesting the presence of a distant companion in those two systems of minimum mass of about 4.3 ± 1.6 and 2.9 ± MJup respectively. The nature of these companions is difficult to assess with the available data. Using different stellar masks to compute the CCF produced no significant effect in both cases. After removing for each target the two measurements affected by the moonlight, the hypothesis of a null Spearman s rank correlation coefficient between the bisector span and radial velocities cannot be rejected at the 5% level for both targets, supporting the hypothesis that the companion has no detectable effect on the CCF. We find the same result when the long-term drift is subtracted from the radial velocities as can be seen in Fig. 2 (left and right panels), where the CCF bisector spans show neither variations nor trends as a function of radial velocity. Neither is there a significative correlation between the radial-velocity long-term drift and the CCF bisector spans (Spearman s rank correlation p-value of 9 for KOI-206 and 0.32 for KOI-680). 3

4 J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates The flux of identified (separations 2 ) contaminant stars in the Kepler photometric mask is estimated for each quarter, and the light curves are corrected. There is still the possibility that the light curve and the CCF are contaminated by stars closer than 2. This contaminant star can be physically related to the target or by chance close to the line of sight. We estimated the probability of the latter being 1.9, 1.8, and 2.0% for KOI-206, KOI-614, and KOI-680, respectively, using the Besançon model (Robin et al. 2003). The magnitude of the stars in the simulation is limited according to the observed transit depth, and its distance to the target to less than 2. We try to constrain the contaminant adding synthetic spectra to the co-added SOPHIE spectrum. The synthetic secondary CCF peak is detectable for early than K0V-K3V stars at the same distance of the target (or equivalent flux ratio for an unrelated contaminating star), if this CCF peak is well detached from the main one in radial velocity. That translated into an upper limit for the possible flux contamination in the light curve and CCF of 5% for KOI-206, 10% for KOI-614, and 3% for KOI-680. Assuming the RV drift observed for KOI-206 and KOI-680 is caused by a companion in a circular orbit, we estimated the allowed masses at each period fitting the RV residuals after having subtracted the Keplerian (Section 4). The result is presented in Fig. 11. We estimated the probability that the companion allowed by the drift is stellar in 34% for KOI-206 and 31% for KOI-680 using the statistics of Raghavan et al. (2010). As a conclusion, for the three targets, the radial-velocity variations detected at the Kepler period are compatible with a companion of a planetary nature, and in the case of for KOI-206b and KOI- 680b, there is no evidence that their estimated parameters are significantly affected by the possible additional companions. 3. Spectroscopic parameters For KOI-206 and KOI-680 we used the SOPHIE spectra described in Section 2.2. For KOI-614 we used spectra obtained with the ESPaDOnS spectrograph. The target was observed in object+sky mode, with a spectral resolution λ/ λ (Program 12AF95). The six observations ran from June 30 to July 5, 2012, with exposure times ranging from 3515 s for the first one to 1715 s for the five others. Each spectrum was corrected for its radial velocity shift, and then the respective spectra were co-added. The final spectrum has a signal-to-noise ratio per element of resolution at 5500Å of 50, 90, and 250 for KOI-614, KOI- 206, and KOI-680, respectively. We derived the effective temperature T eff, surface gravity log g, metallicity [Fe/H], and microturbulence v micro of the stars with the VWA package (Bruntt et al. 2004, 2008, 2010b). We used a grid of MARCS models (Gustafsson et al. 2008) to reach the ionization equilibrium over a large selection of spectral lines. Then, we fitted a rotational profile on a set of isolated metallic lines, solving for the projected rotational velocity v sin i and the macroturbulence v macro. For KOI-206 and KOI-680, instead, the v sin i was measured through the Gaussian width of the cross-correlation function (CCF) (Boisse et al. 2010, and references therein), because of the paucity of isolated lines to be fitted with VWA. The spectrum of KOI-614 did not present any particular problem, and the relative atmospheric parameters are presented in Table 2. For KOI-206, we found a considerably high value for v micro (2.8 ± km s 1 ). The large uncertainty can be influenced by its degeneracy with v macro, which is problematic for constraining a spectrum with a low signal-to-noise ratio. However, the measurement favors a value higher than 2 km s 1, indicating an evolved object (e.g., Gray 2005). We imposed different values for v macro (between 5.0 and 7.5 km s 1 ) and fitted all the other atmospheric parameters. In this way, we verified that the latter are not significantly affected. Eventually, we measured v sin i = 11 ± 1 km s 1 with the CCF. In conclusion, the ionization balance required T eff = 6340 ± 140 K, log g = 4.0 ± 0.3 dex, and [Fe/H]= 6 ± 0.19 dex. For KOI-680, the analysis yielded T eff = 6090 ± 110 K, log g = 3.5 ± 0.1 dex, [Fe/H] = 0.17 ± 0.10 dex, and v micro = 1.3 ± km s 1. Because of the loose constraint we obtained on v macro with VWA, we chose to fix it through the T eff, according to the relation in Bruntt et al. (2010a), to 3.6 km s 1. Finally, from the CCF we measured v sin i = 6 ± 1 km s 1. We used the SME package from Valenti & Piskunov (1996) and Valenti & Fischer (2005) to check the reliability of the results for KOI-206 and KOI-680. We verified T eff on H α, the log g on the pressure-sensitive lines Caλ6122, Caλ6162, and Caλ6439, and on the MgIb triplet. Fixing the T eff, log g, and [Fe/H] on KOI- 206, we fitted v macro 7 km s 1 and v micro 2.8 km s Modeling with PASTIS We use the Planet Analysis and Small Transit Investigation Software (PASTIS, Díaz et al. 2014), that allows the light curve, the radial velocity, and the spectral energy distribution (SED) to be fit with a Markov Chain Monte Carlo (MCMC) algorithm. The photometric magnitudes used for the SED fit are listed in Table 2. Only the sections of the light curve around the transits are used, after normalization to a linear trend fitted outside the transit. A 3-σ clipping is performed to eliminate outliers, and the measured flux is corrected for the contamination factors provided for each quarter at MAST archive. PASTIS uses the EBOP code (Nelson & Davis 1972; Etzel 1981; Popper & Etzel 1981) extracted from the JKTEBOP package (Southworth 2011) to model the photometric transits. Both long- and short-cadence (when available) are used, with an oversampling factor of 10 for the long-cadence data to account for the long integration time when comparing with the model. Radial velocity curves were simultaneously fitted with eccentric Keplerian orbit and a linear drift. We also account for additional sources of Gaussian noise in the light curves, radial velocities, and SED by fitting a jitter value to each data set. We used the stellar evolutionary tracks from STAREVOL (Palacios, priv. com.), Dartmouth (Dotter et al. 2008), and PARSEC (Bressan et al. 2012a) to estimate the stellar parameters, and the PHOENIX/BT-Settl synthetic spectral library (Allard et al. 2012) to model the SED. The interpolated spectrum is scaled to a given distance and corrected from interstellar extinction. The distance d and the color excess E (B V) were included as free parameters in our model. The model has therefore 20 jump parameters for KOI-614b, and two more for KOI-206b and KOI-680b that have short-cadence data, indicated in Table 2. We used non-informative priors (uniform or Jeffreys distributions) except for T eff, [Fe/H], and ρ, for which we used spectroscopic estimations as normal and asymmetric normal distribution priors, and P and T c, for which we used the determination by Borucki et al. (2011) as normal distribution priors, but the width of the distribution was increased by an order of magnitude to ensure that it would not bias the results. It is a well known problem that the uncertainty in the stellar parameters obtained from evolutionary models may be underestimated. To account for this, we use different evolutionary models as input for the stellar parameters, and the combined outcome will naturally enlarge the errors in the stellar param- 4

5 Table 2. KOI-206, KOI-614, and KOI-680 system parameters. J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates Object Kepler-433 Kepler-434 Kepler-435 Kepler object of interest KOI-206 KOI-614 KOI-680 Kepler input catalog KIC KIC KIC USNO-A2 ID MASS ID Coordinates RA (J2000) [hh:mm:ss.sss] 19:50: :34: :29: Dec (J2000) [dd:mm:ss.ss] +40:58: :55: :11:51 Magnitudes Kepmag (AB) d51mag (AB) g Sloan-Gunn (AB) r Sloan-Gunn (AB) i Sloan-Gunn (AB) z Sloan-Gunn (AB) J 2MASS (Vega) ± ± ± 20 H 2MASS (Vega) ± ± ± 19 Ks 2MASS (Vega) ± ± ± 22 WISE W1 (Vega) ± ± 23 WISE W2 (Vega) ± ± 24 WISE W3 (Vega) ± Spectroscopic parameters Effective temperature, T eff [K] 6340 ± ± ± 110 Surface gravity, log g [cgs] 4.0 ± ± ± 0.1 Metallicity, [Fe/H] [dex] 6 ± ± ± 0.10 Stellar rotational velocity, v sin i [ km s 1 ] 11 ± ± 6 ± 1 Microturbulent velocity, v micro [ km s 1 ] 2.8 ± ± 1.3 ± Macroturbulent velocity, v macro [ km s 1 ] 5 ± ± 3.6 (fixed) Result from light curve, radial velocity, and SED combined analysis Planet orbital period, P [days] ± ± ± Mid-transit time, T c [BJD] ± ± ± cov(p, T c ) [days 2 ] Orbital eccentricity, e ± 79, < ± 72, < ± 77, < 0.30 Argument of periastron, ω [deg] ± ± 36 Orbit inclination, i [deg] ± 0.52 Orbital semi-major axis, a [AU] 679 ± ± ± 018 semi-major axis / radius of the star, a/r 6.44 ± ± ± 0.51 Radius ratio, k = R p /R 6590 ± ± 0020 Linear limb darkening coefficient, u a ± 22 0 ± ± 24 Quadratic limb darkening coefficient, u b 20 ± 45 1 ± ± 42 Transit duration [h] ± ± Impact parameter, b 82 ± ± ± 24 Systemic velocity, V r [ km s 1 ] ± ± ± 079 Radial velocity linear drift [ m s 1 yr 1 ] -500 ± ± ± 9.2 Radial velocity semi-amplitude, K [ m s 1 ] 259 ± ± ± 12 Jitter Kepler long cadence, σ Kepler,LC [ppm] 65 ± ± Jitter Kepler short cadence, σ Kepler,SC [ppm] 88 ± ± 32 Jitter radial velocity, σ RV [ m s 1 ] 34 ± ± 10 Jitter SED, σ SED [mag] 56 ± ± 16 Effective temperature, T eff [K] 6360 ± ± ± 94 Metallicity, [Fe/H] [dex] -1 ± 0 5 ± ± 0.11 Stellar density, ρ [ρ ] ± 33 6 ± ± 11 Star mass, M [M ] 1.46 ± ± ± 88 Star radius, R [R ] 2.26 ± ± ± 0.30 Deduced stellar surface gravity, log g [cgs] ± ± Age of the star [Gyr] 2.67 ± ± ± 2 Luminosity of the star, log(l/l ) [dex] 8 ± ± ± 89 Distance of the system [pc] 1870 ± ± ± 200 Color excess, E (B V) [mag] 10 ± ± ± 25 Planet mass, M p [M Jup ] 2.82 ± ± ± 0.15 Planet radius, R p [R Jup ] 1.45 ± ± 0.18 Planet mean density, ρ p [g cm 3 ] 1.13 ± ± ± 37 Planet surface gravity, log g p [cgs] ± 0.11 Planet equilibrium temperature, T eq [K] 1776 ± ± ± 70 Not used in the SED fitting. MCMC jump parameter. upper limit, 99% confidence. reference time BJD T eq = T eff (1 A) 1/4 R, using an albedo A = 0. 2a M = kg, R = m, M Jup = kg, R Jup = m 5

6 J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates eters from the differences between them. For each of the targets, we therefore fit the data three times using a different evolutionary model: STAREVOL, Dartmouth, and PARSEC. To ensure a broad exploration of the parameter space, we compute 30 chains of 10 6 steps, starting at random points drawn from the joint prior. After removing the burn-in interval of each chain, we obtain a merged chain by thinning each chain. The thinning factor is determined by the maximum correlation length found among all parameters in each chain (e.g. Tegmark et al. 2004). This is done to obtain independent points to form the posterior distributions. Then we take the same number of samples of each merged chain from each stellar model, obtaining the combined merged chain. We obtain 3168, 729, and 4564 independent points in the combined merged chain for KOI-206, KOI-614, and KOI-680, respectively, which we use to obtain the estimated value and the 68.3% central confidence interval, both for jump and derived parameters (listed in Table 2). The model of maximum likelihood is plotted together with the data in Fig. 3. The correlation distributions and histograms are shown in Figs. 12, 13, and 14 for the MCMC jump parameters. We found that the distributions of the posteriors obtained from the different stellar models are compatible within 1-σ. The stellar evolution tracks from STAREVOL in a Hertzsprung-Russell diagram are given in Fig. 4, along with the luminosity-t eff joint posterior distributions to show the evolutionary stage of the host stars. 5. Planetary evolution models KOI-614b, KOI-206b, and KOI-680b are all different from each other: with radii of 1.13, 1.45, and up to 1.99 R Jup, and masses of 2.86, 2.82, and down to 4 M Jup, they occupy distinct regions of the mass-radius diagram as illustrated in Fig. 5. To estimate the structure and composition of the three planets, we computed the joint evolution of both planet and host star for each system, combining the CEPAM (Guillot & Morel 1995; Guillot 2010) and PARSEC (Bressan et al. 2012b) evolutionary models using SET (see Guillot & Havel 2011; Havel et al. 2011; see also Almenara et al. 2013). Since stellar and planetary absolute parameters are ultimately model dependent, the values we Fig. 5. Radius versus mass of all known transiting planets as of February 2014 (see Wright et al. 2011). The three planets presented in this paper are highlighted in red and labeled. obtained are solely based on observations. Through the use of SET s MCMC algorithm, we thus obtained the posterior probability distributions of the bulk composition of the planets (i.e., their core mass), as well as an independent estimation of planetary and stellar parameters. The latter are entirely consistent with those presented in Section 4. The results for KOI-206b, KOI-614b, and KOI-680b are presented in terms of planetary radii as a function of age in Fig. 6: the regions show the 1, 2, and 3-σ constraints from modeling the star and transit, while the lines show a subset of planetary models for the nominal mass and equilibrium temperature of the planet and different compositions. All our planetary evolution models assume the following: the planet is made of a central rocky 2 core and a solar-composition 3 envelope. Of course we do not know whether all heavy elements are in the core or mixed in the envelope or know their exact composition. However, Baraffe et al. (2008) show that having the heavy elements mixed in the envelope generally leads to smaller radii. Therefore our models should provide an upper limit for Z, the fraction of heavy elements in the planet. In addition, we do not expect the envelope to have a solar composition (it is not even the case for Jupiter or Saturn; see, e.g., Guillot & Gautier 2014), but the uncertainties in the models themselves (EOS, opacities, atmospheric model) play a larger role than the small adjustments 4 in the composition of the envelope due to the different stellar metallicities (Z 13 for KOI-680 and Z 42 for KOI-614). For consistency with previous modeling of exoplanets, planetary evolution models are calculated in two cases: (i) a standard case where only the irradiation from the star is accounted for; (ii) a dissipated-energy model in which, in addition to the standard case, a fraction (5%) of the incoming stellar flux is assumed to be converted into kinetic energy 5 and then dissipated at the center of the planet (see, e.g., Guillot & Showman 2002 and Spiegel & Burrows 2013 for detailed discussions). With a radius of 1.13 R Jup and a mass of 2.86 M Jup, both standard and dissipated models provide solutions for the planetary radius of KOI-614b that match the available constraints, resulting in a core mass of M. As seen in Fig. 6, models with no core (i.e., pure H/He planet) land in the middle of the 1-σ region, but the high metallicity of the star (+0.35) tends to indicate that the planet should probably have a significant number of heavy elements (Guillot et al. 2006). On the other hand, KOI-206b has a similar mass of 2.82 M Jup, but a radius of 1.45 R Jup, and an inferred density of 1.13 g cm 3. For such a high radius and nearly three times the mass of Jupiter, KOI-206b does not have many similar siblings: the three closest in the mass-radius diagram (Fig. 5) are CoRoT- 2b (Alonso et al. 2008), CoRoT-18b (Hébrard et al. 2011), and Kepler-17b (Désert et al. 2011; Bonomo et al. 2012). The former two, for instance, have been modeled and are still challenging the current planetary models (see, e.g., Guillot & Havel 2011; Moutou et al. 2013). As Fig. 6 shows, this is also the case of KOI-206b: standard models barely succeed at matching the 2- σ age-radius constraint, while dissipated-energy models are just reaching the bottom of the 1-σ region. Although an indepen- 2 38% SiO 2, 25% MgO, 25% FeS, 12% FeO 3 70% H, 28% He, 2% others (Z) 4 as compared to the expected, much higher Z values in planets than in stars (Guillot et al. 2006) 5 Note the use of ke labels in figures: 1% ke = 5% of the incident flux at substellar point or 1% of the irradiating flux averaged over the planet s area converted into kinetic energy (ke). 6

7 J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates KOI-206 KOI-614 KOI log10 (L/L ) 3 log10 (L/L ) log10 (L/L ) Fig. 3. Model of the maximum likelihood step of the merged chain of KOI-206, KOI-614, and KOI-680 (from left to right columns) plotted with the data. From top to bottom: SED, light curve, radial velocity in time, and radial velocity in phase. In the bottom panels: residuals after subtracting the model to the observed data. SED plot: the solid line plots the PHOENIX/BT-Settl interpolated synthetic spectrum, in red circles the absolute photometric observations listed in Table 2, in blue open circles the result of integrated the synthetic spectrum in the observed bandpasses. Light curve plot: folded transit observed by Kepler in long-the black solid line plots the Keplerian model a for the radial velocities and red circles the SOPHIE observations Teff [K] Teff [K] Teff [K] Fig. 4. Stellar evolution tracks from STAREVOL, from the zero-age main-sequence stage to 25 Gyr age, in the luminosity-effective temperature plane for [Fe/H] = 0, 5, and -0 (from left to right, solid line) and [Fe/H] = 0.10, 0.35, and (from left to right, dashed line) with the joint posterior distributions of the MCMC (98.9% joint confidence region is denoted in red) for KOI-206, KOI-614, and KOI-680, respectively. The mass in solar masses is annotated at the beginning of the main sequence of each track. 7

8 J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates R p [R Jup ] M p =2.82 M jup T eq =1776 K ɛ =1% ke = erg/s KOI-206b 50% ke 0 M 10% ke 0 M age [Ga] 0 M 100 M 199 M 299 M R p [R Jup ] ɛ =1% ke = erg/s M p =2.86 M jup T eq =1000 K KOI-614b 397 M age [Ga] 0 M 101 M 202 M 303 M 404 M R p [R Jup ] ɛ =1% ke = erg/s M p =4 M jup T eq =1729 K KOI-680b 32% ke 0 M 29% ke 0 M 23% ke 0 M 16% ke 0 M 10% ke 0 M age [Ga] 0 M 15 M 30 M 44 M 59 M 74 M Fig. 6. Evolution of the planet radius as a function of the age. KOI-206b with mass 2.82 M Jup and effective temperature 1776 K is plotted to the left, KOI-614b (2.86 M Jup and 1000 K) in the center plot, and KOI-680b (4 M Jup and 1729 K) in the right plot. The 68.3%, 95.5%, and 99.7% joint confidence regions are denoted by black, dark gray, and light gray areas, respectively. The curves represent the thermal evolution of a respective planet mass and equilibrium temperature. Text labels indicate the amount of heavy elements in the planet (its core mass in Earth masses). Dashed lines represent planetary evolution models for which 5% of the incoming stellar flux ( ɛ) is dissipated into the core of the planet, whereas plain lines do not account for this dissipation (standard models). Dashed-dotted lines (gray shades) are coreless models with higher amounts of dissipated energy, as labeled. dent MCMC 1-D distribution gives a core mass of 0 M with a 68.3% confidence interval [0, 74] M, thus promoting solutions for which the planet has a much lower radius. Extra models with incredibly large amounts of energy dissipated in the planet have to be created to explain the observed radius at 1- σ. Up to 75% of the total light reaching the planet would need to be dissipated for these models to work out a solution. There is, to date, no known physical process that would allow such efficiency. However, the recent work of Spiegel & Burrows (2013) suggests that dissipating the energy in the upper atmosphere instead of the deep interior maintains the planet s hot state (i.e., large radius) longer for ages below 2-3 Gyr. Last but not least, the use of a newer H/He equation-of-state (EOS) may contribute significantly to solving the problem: while our models are using SCVH (Saumon-Chabrier-Van Horn, Saumon et al. 1995), the EOS from Militzer & Hubbard (2013) may produce planets with bigger radius, up to about R Jup larger than ours. Last of all, with a radius of 1.99 R Jup and a mass of 4 M Jup, KOI-680b is among the largest planets ever known. Only two planets have a higher radius, HAT-P-32b (Hartman et al. 2011) and WASP-17 (Anderson et al. 2011) (see Fig. 5), and as few as eleven have a radius larger than 1.7 R Jup. WASP-78b (Smalley et al. 2012) is the only one of similar mass and radius around an evolved star with low metallicity: unfortunately, it has not been modeled yet. As expected, our standard planetary models cannot account for such a high radius. Even with the dissipation of 1% of the available irradiation, we hardly reach the 3-σ constraint region. Once again, we had to produce some more extreme models with unphysical amounts of energy dissipated into the deep planet to achieve the 2 R Jup limit. Interestingly, we note that dissipating just 10% to 40% of the incident stellar flux is enough to match the 1-σ region, possibly less than what is required for KOI-206b. Previously suggested improvements over our models for KOI-206b should apply for this planet. On the other hand, KOI-680b is below a Jupiter mass and should therefore be easier to re-inflate. 6. Discussion 6.1. KOI-614b KOI-614b was a moderate probability candidate, but the radial velocity measurements presented in this paper establish the Fig. 7. Planetary mass as a function of the orbital period for known transiting giant exoplanets (Wright et al. 2011). The three planets presented in this paper are highlighted in red and labeled. The dotted line represents the Roche limit, and the dashed line corresponds to twice this value, which would be the final orbital period of an initial eccentric long-period planet that would have completed its orbit. planetary nature of the companion, its high-impact parameter (b 0.98) explaining the V-shaped transits. With an equilibrium temperature of T eq = 1000 ± 45 K, KOI-614b is at the transition between what is usually referred to as hot and warm Jupiters and is a good candidate for testing migration theories. KOI-614b is one of the few transiting giant planets with an orbital period longer than ten days, as can be seen in Fig. 7. Unlike hot Jupiters, this kind of planets has a pericenter distance that is too large to allow efficient tidal dissipation to induce migration through the high-eccentricity migration mechanism (Dong et al. 2014). In order to access the close pericenter required for migration and circularization by tidal dissipation of the eccentricity, they need to be accompanied by a strong enough perturber to overcome the precession caused by General Relativity. KOI- 614b has a low eccentricity (e 0.32), but there is no trace of an additional massive companion in the system. Using the radial 8

9 J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates greater heavy element content can be obtained with particular disk conditions (the case of monarchical growth, Mordasini et al. 2009; Weidenschilling 2005) or by the collision with another planet(s) or planetesimals after substantial accretion of the envelope (Ikoma et al. 2006). In the latter case, the disk would have fostered the formation of multiple embryos, which would provide a natural way to stop the migration of the planet at its current semi-major axis (a AU). Both migration scenarios require the presence of additional companion(s) in the system. High-eccentricity migration mechanisms require a massive, not very distant perturber to bring the warm Jupiter into a Kozai-Lidov cycle, type II migration requires multiple planets to clear a wide inner hole and prematurely stop the migration of the planet. KOI-614b is thus an interesting candidate for future observations with the aim of finding additional low-mass, short-period planetary companions or a more distant and massive perturber. Fig. 8. Planetary density as a function of mass for known transiting giant exoplanets (Wright et al. 2011). The three planets presented in this paper are highlighted in red and labeled. velocity residuals of our fit, we can exclude the presence of a companion of more than 2.0 M Jup on orbital periods less than 200 days, at the 3-σ level. Given the length of our observation period, we cannot make any conclusion about more distant companions. In addition, Mazeh et al. (2013) report no significant transit timing variations for this target. The high-eccentricity migration scenario cannot be favored yet. Considering the mass of 3 M Jup of KOI-614b, its formation would have quickly opened a gap in the protoplanetary disk, and another migration scenario would involve type II migration. However, as shown by Burkert & Ida (2007), whatever the disk lifetime, type II migration of a single Jovian planet in a truncated disk naturally leads to the creation of period gap between 8 a AU in the distribution semi-major axis, which is especially pronounced for systems with stellar masses M 1.2 M. As this is precisely the case of this system, the migration of the planet within the disk could be doubted. Nevertheless, while single planets open only narrow gaps, multiple planets can clear wide inner holes (Dodson-Robinson & Salyk 2011), thus preventing the migration of the Jovian planet down to a short-period orbit. As previously said, our observations exclude the presence of Jovian companions over orbital periods < 100 days, but we cannot exclude the presence of less massive, terrestrial planets in the system. In this regard, the composition of the planet can provide additional clues. The estimated radius of 1.13 R Jup makes KOI-614b a high-density giant planet. As can be seen in Fig.8, this mass and density are common in the population of the transiting exoplanets discovered so far 6. Within our estimated uncertainties, we have established that the dimensions of this object can encompass internal compositions ranging from pure hydrogen to containing a solid core of heavy elements up to 150 M. Even if the pure hydrogen scenario cannot be discarded, some authors (see, e.g., Guillot et al. 2006) have noted that there is observational evidence of a correlation between star metallicity and heavy elements in the planet. Considering that KOI-614b is orbiting a main-sequence, metallic ([Fe/H] = 0.35 dex), G0V star, it is reasonable to favor a high heavy-element content. In the core-accretion scenario, the canonical core mass is about 10 M (Mizuno 1980). A 6 as of Februray 2014 at KOI-206 Like KOI-614b, KOI-206b is a relatively massive planet of about 2.8 M Jup, but its radius of 1.45 ± 0.16 R Jup implies a different inner structure with a Jupiter-like density. Although the much tighter orbit of KOI-206b (a 679 AU) puts it in the range of hot Jupiters, and although the higher effective temperature of its host star does imply a stronger irradiation than for KOI-614b, its large radius is particularly challenging for planetary evolution models. The radius of a transiting planet is one of the most robust measurements because it only relies on the knowledge of the transit depth and of the stellar radius (with a second-order dependence on limb darkening, and perhaps light curve contamination). Here, the stellar radius is determined using evolutionary models, which are known to be inaccurate, and yet this is partially taken into account by our fitting method by combining several sets of evolutionary tracks in the final result. Concerning the transit depth, it is precisely measured at the 210-σ level. Nonetheless, the light curve of KOI-206b exhibits rotational modulations due to spot coverage (cf. Sect. 2.1). There is no sign of spot occultation during the transits of the planet, but high-latitude spot could still affect the estimation of the radius. In particular, the presence of such a group of spots could be indistinguishable in the light curve and would induce an overestimation of the radius of the planet. As shown in Fig. 6, typical models of bloated hot Jupiters, where 5% of the incoming stellar flux is dissipated in the core, would agree better with a planetary radius in the bottom of the 1-σ confidence region. Still, to reduce the estimated planetary radius to this value, i.e. to about 1.3 R Jup, we would need to assume that the value of the flux continuum is actually a factor 1.26 more than the value we are using. It would require a polar spot of angular radius of 25, in the very conservative assumption that the contrast of the spot is 0 and the rotational axis of the star is aligned on the line of sight. More realistically, if we assume a surface-to-spot contrast of about 7 (Lanza et al. 2003) and the limb darkening coefficients of our best model, this would require a polar spot with a radius of at least 59 (estimated using Kipping 2012). Such a spot could be betrayed by the activity level estimated with the emission cores of the CaII lines. Unfortunately, the signalto-noise ratio of the SOPHIE spectra in the blue part of the spectrum is too low to estimate the activity with CaII H and K lines. Either way, the dimensions of such a spot is rarely observed (Lanza et al. 2003), and it is very unlikely that our plan- 9

10 J.M. Almenara et al.: SOPHIE velocimetry of Kepler transit candidates Fig. 9. Planetary radius as a function of stellar radius for known transiting exoplanets (Wright et al. 2011). The dotted lines give the loci of some typical transit depth given in labels. The three planets presented in this paper are highlighted in red and labeled. etary radius estimation would be overestimated because of such an effect. Besides, the stellar density inferred from the transit (Table. 2) is in good agreement with the stellar parameter determined by the spectral analysis, and there is no evidence of a false positive. At the time of writing, KOI-206b is the transiting planet detected around the host with the second largest radius so far 7 just after KOI-680b, as shown in Fig. 9. There are only a few other known bloated planets in this mass range (see Fig. 5), and KOI- 206b is the one with the lowest density in this group (see Fig. 8). For such a mass, the amount of extra dissipated energy to produce that kind of density is preposterous. KOI-206b orbits an evolved F7-type star that has recently left the main sequence (see Fig. 4), and the dynamic aspects related to the rapid and recent evolution of the host might have effects that are unaccounted for. As noted by Burrows et al. (2003), there is a difference between the observed transit radius (at mbar levels) and the modeled radius, R. Since the evolution of a giant planet is largely controlled by its atmosphere, the planetary evolution code we used, CEPAM, has a dedicated module that computes the exact P-T profile of the atmosphere up to mbar levels (see Parmentier et al for more details on the atmospheric model and, e.g., Guillot 2008 for an estimation of R). Therefore, our results for modeling the interior structure of the three giant planets can be directly compared to the observations: R cannot be the explanation for the discrepancy between theoretical models and observed radius KOI-680b The host star of KOI-680b shares many properties with the one of KOI-206b, but with a mass of only 4 ± 0.15 M Jup and a radius of 1.99 ± 0.18 R Jup, KOI-680b is one of the biggest (Fig. 5) and least dense (Fig. 8) planets detected to date. Such properties exclude a high concentration of heavy elements in its interior (Fig. 6), which correlates well with the low metallicity of the host star ([Fe/H] = 0.17 ± 0.10 dex). Although standard planetary evolution models fail to reproduce its observed characteristics, the low mass of KOI-680b also means that it is easier to 7 see note 6 Fig. 10. Kiel diagram for the stars of known transiting (in blue) and radial velocity (in orange) exoplanets (Wright et al. 2011). KOI-614, KOI-206, and KOI-680 are drawn in black and labeled. We also plot the STAREVOL evolutionary tracks for stars of 1, 1.5, and 2 M and solar metallicity. inflate, and improvements in the equation of state might easily solve the problem (Militzer & Hubbard 2013). Not only is the planet one of the biggest ever observed, but it also orbits the bigger star observed with the transit technique (Fig. 9). Indeed the host star of KOI-680b is an F9 subgiant, well on its way toward the beginning of the red giant branch (RGB). Although it is slightly younger than the host of KOI-206b, it is also more massive (M = ± 88 M ) and thus more evolved. As can be seen in Fig. 10, the host star of KOI-206b is observed in a phase of rapid stellar evolution, where no other planets have been detected before. With an orbital semi-major axis of 948±018 AU, KOI- 680b is doomed to be engulfed by the star, even neglecting tidal interaction, the stellar radius will catch up with the current orbit of the planet within the next 260 Myr or so. We note that we detect a companion of minimum mass 2.9± M Jup with minimum orbital period of about 790 days (i.e., about 1.9 AU). Further radial velocity follow-up is required to determine the period and the planetary nature of the companion. 7. Conclusion We have presented the detection and characterization of KOI- 614b, KOI-206b, and KOI-680b, which are three new transiting, giant extrasolar planets. They were first detected as promising candidates by the Kepler team from the Kepler light curve. We established their planetary nature with radial velocity followup using the spectrograph SOPHIE and characterized the systems through a combined MCMC Bayesian fit of photometric and spectroscopic datasets. We find that KOI-614b is a dense, warm Jupiter, with a mass of 2.86 ± 0.35 M Jup and a radius of R Jup, and it is orbiting a G0, metallic ([Fe/H]=0.35 ± 0.15) dwarf in 12.9 days. KOI-206b has a mass of 2.82 ± 0.52 M Jup and a radius of 1.45 ± 0.16 R Jup, and it orbits a slightly evolved F7-type star in a 5.3-day orbit. Orbiting a F9 subgiant in 8.6 days, KOI-680b has a much lower mass of 4 ± 0.5 M Jup and an uncommonly large radius of 1.99 ± 0.18 R Jup. With host stars of masses of 1.46 ± 0.17 M and 1.54 ± 9 M and 10