SF2A 2016 L. Cambr´esy, M. Deleuil, E. P´econtal, C. Reyl´e, J. Richard, L. Tresse and I. Vauglin (eds) FOLLOW-UP AND CHARACTERIZATION OF THE TESS EXOPLANETS WITH SOPHIE, SPIROU, AND JWST N. Crouzet1, X. Bonfils2, X. Delfosse2, I. Boisse3, G. H´ebrard4, T. Forveille2, J.-F. Donati5, F. Bouchy3,6, C. Moutou3,7, R. Doyon8,9, E. Artigau8,9, L. Albert8,9, L. Malo7,8,9, 7 A. Lecavelier des Etangs4 and A. Santerne3 1 0 2 Abstract. n The NASA TESS mission will deliver hundreds of transiting exoplanet candidates orbiting bright stars. a The spectrometers SOPHIE at OHP and SPIRou at CFHT will be ideal to obtain radial velocities of J these candidates, confirm their nature, and derive the planets’ masses. These measurements will be crucial 3 to deliver the best targets for atmospheric characterization with JWST. Here, we calculate the required 1 observingtimewithSOPHIE,SPIRou,andJWST foreachoftheTESS targetsinordertopreparefollow- up observations. To infer their potential for JWST, we restrict the calculations to the case of transmission ] P spectroscopywithNIRISS.Theradialvelocityfollow-upofthegiantplanets(Rp >4RE)couldbeachieved E withSOPHIE,withamedianobservingtimeof3.47hourspertarget,andatotalobservingtimeof305hours . thatincludesthe80%mostfavorablecases. Severalsmallplanets(Rp <4RE)couldalsobeconfirmed,but h most of them would require an unrealistic time investment. On the other hand, SPIRou is ideally suited p to the follow-up of the small planets, with a median observing time of 2.65 hours per target, and a median - o observing time of 4.70 hours for the terrestrial planets in the habitable zone (Rp <2RE, S <2SE). With r JWST, the 10% most favorable small planets (184 planets) have a median observing time of 16.2 hours, t s and the 10% most favorable habitable zone terrestrial planets (7 planets) have a median observing time of a 59.7 hours. Overall, this study will help define a follow-up strategy and prepare observation programs with [ SOPHIE and SPIRou before the first planet candidates are delivered by TESS. 1 v Keywords: Surveys, Planets and satellites: detection, Planets and satellites: atmospheres, Methods: 9 observational, Techniques: radial velocities, Techniques: spectroscopic 3 5 3 1 Introduction 0 . 1 Exoplanets transiting in front of bright stars offer the best prospects for their physical characterization: by 0 deliveringmorephotons, brightstarsincreasetheefficiencyoffollow-upobservationsandmoreinformationcan 7 be obtained on the systems compared to fainter stars. TESS (Transiting Exoplanet Survey Satellite, launch: 1 2017, PI: G. Ricker, Ricker et al. 2014) is a NASA mission to detect transiting exoplanets around bright : v stars (4<I <13) by high precision wide field photometry. One year later, JWST (James Webb Space mag i X r a 1 Dunlap Institute for Astronomy and Astrophysics, University of Toronto, 50 St Georges Street, Toronto, ON M4Y 2K2, Canada 2 Universit´eGrenobleAlpes,CNRS,IPAG,38000Grenoble,France 3 Aix-MarseilleUniversit´e,CNRS,LAM(Laboratoired’AstrophysiquedeMarseille)UMR7326,13388Marseille,France 4 Institut d’Astrophysique de Paris, UMR 7095 CNRS, Universit´e Pierre & Marie Curie, 98bis boulevard Arago, 75014 Paris, France 5Universit´edeToulouse,UPS-OMP,IRAP,14avenueE.Belin,Toulouse,F-31400France&CNRS,IRAP/UMR5277,Toulouse, 14avenueE.Belin,F-31400France 6 ObservatoiredeGen`eve,Universit´edeGen`eve,51ch. desMaillettes,1290Sauverny,Switzerland 7 Canada-France-HawaiiTelescopeCorporation,CNRS,65-1238MamalahoaHwy,Kamuela,HI96743,USA 8 D´epartementdephysiqueandObservatoireduMont-M´egantic,Universit´edeMontr´eal,Montr´eal,QCH3C3J7,Canada 9 InstitutdeRecherchesurlesExoplan`etes(iREx),Universit´edeMontr´eal,D´epartementdePhysique,C.P.6128Succ. Centre- ville,Montr´eal,QCH3C3J7,Canada (cid:13)c Soci´et´eFrancaised’Astronomieetd’Astrophysique(SF2A)2016 238 SF2A 2016 Telescope,launch: 2018)willprovideauniqueopportunitytoprobetheatmospheresoftransitingexoplanetsby spectroscopyatinfraredwavelengths. MostoftheexoplanettargetsforJWST areexpectedtocomefromTESS. With such a tight schedule, efficient follow-up of the TESS planet candidates is crucial. The spectrometers SOPHIE (Bouchy et al. 2009) at OHP (Observatoire de Haute-Provence) and SPIRou (Delfosse et al. 2013) at CFHT (Canada-France-Hawaii Telescope) will be ideal to confirm the planetary nature of these candidates and tomeasuretheirmassesbyradialvelocities. SOPHIEhasbeensuccessfullyusedforthefollow-upofCoRoT and Kepler candidates,andtheTESS targetswillbesignificantlybrighter. SOPHIEoperatesatvisiblewavelengths (3872-6943 ˚A). SPIRou is a near-infrared spectro-polarimeter that will enable radial velocity measurements at high precision (< 1 ms−1) and high resolving power (> 70,000) in the 0.98-2.35 µm bandpass in one shot. One of its main goals is the detection of low-mass exoplanets around low-mass stars, and the first light is planned for the end of 2017, early 2018. The small size of M dwarfs also favours the detection and characterization of small planets through their transits, and TESS will detect many more planets around M dwarfs than Kepler andK2 thankstoitsalmostfullskycoverage. Inthiswork,wepresentestimatesoftheobservingtimerequired with SOPHIE and SPIRou to follow-up the TESS simulated planets in radial velocity, and the time required with JWST to characterize their atmospheres. This study will be useful to define follow-up strategies and to prepare observation campaigns before the first TESS candidates are delivered. The observation chain and the simulations are described in Section 2. Section 3 gives estimates of the required observing time per target. Section 4 presents the cumulative observing time with SOPHIE and SPIRou for each planet category in order to prepare follow-up programs. 2 Observation chain and description of the simulations 2.1 Transiting planet candidates with TESS The main goal of the NASA TESS mission is to detect transiting exoplanets orbiting around bright stars. We use the catalog of simulated TESS detections presented in Sullivan et al. (2015) as an input of our simulations. We consider the expected yield of planets around the 2×105 target stars: about 70 Earths, 486 Super-Earths, 1111 Sub-Neptunes, and 67 giant planets; these numbers are only indicative as they differ for each simulation. Wedonotconsiderdetectionsfromthefullframeimages(seeFigure18ofSullivanetal.2015,formoredetails). Thiscatalogprovidesseveralphysicalparametersofthesimulatedsystems: planetaryradius R , orbitalperiod p P, insolation S, host star radius R , effective temperature T , magnitude in several bandpasses, and radial (cid:63) eff velocity semi-amplitude K. We derive other quantities: the planet mass M using the empirical mass-radius p relation provided by Weiss et al. (2013), the host star mass M , the semi-major axis a, the transit depth δ and s duration τ, the planet equilibrium temperature T , density ρ , surface gravity logg , and atmospheric scale eq p p height H. We assume equatorial transits (inclination i=0) and circular orbits (eccentricity e=0). We do not consider planet multiplicity in this study. 2.2 Radial velocity follow-up with SOPHIE and SPIRou Radial velocity follow-up of the TESS candidates will be necessary to confirm their planetary nature and measure their masses. The spectrometers SOPHIE in the visible and SPIRou in the near-infrared are ideal in that purpose. For each TESS target, we calculate the observing time t that would be necessary to measure RV K with a given signal to noise ratio S/N with these two instruments. We use a reference radial velocity precision K of 6 ms−1 at a magnitude V of 12 in an exposure time t of 15 minutes for SOPHIE in the ref ref ref high resolution mode (G. H´ebrard, private communication), and an expected precision K of 1 ms−1 at a ref magnitude J of 9.5 in an exposure time t of 15 minutes for SPIRou (Delfosse et al. 2013). We scale these ref ref reference precisions to the magnitudes of the TESS targets according to Poisson noise. We impose a minimum exposuretimeof10minutestoaveragestellarpulsations,andafloorprecisionperpointof2ms−1 forSOPHIE and 1 ms−1 for SPIRou. We assume that the measurements are made at quadratures of the radial velocity curves, and we consider that a planet is detected if K is measured with a S/N of 3. In this study, we do not consider stellar rotation and stellar activity which may limit the radial velocity precision. 2.3 Atmospheric characterization with JWST Theatmosphereoftransitingexoplanetscanbecharacterizedbytransmissionandemissionspectroscopyduring their transits and eclipses, respectively. Here, we consider the case of transit spectroscopy of the TESS planets Characterization of the TESS exoplanets 239 with JWST, and we calculate the observing time t that is necessary to detect molecules in their atmo- JWST spheres. We use the JWST/NIRISS instrument in the SOSS mode (Single-Object Slitless Spectroscopy) which is dedicated to exoplanet spectroscopy, through the online NIRISS SOSS 1-D Simulator∗. We run a reference simulation using a set of parameters corresponding to the mean of the TESS small planets (R <4R , where p E R is the Earth’s radius). Then, we scale the results to each TESS target: we calculate the target’s flux in the E NIRISS bandpass from its T and J-band magnitude using a blackbody spectral energy distribution, and we eff assume that the noise scales as the Poisson noise. We assume the same observing time in and out of transit, andweusethe1storderspectrumonly. Theobservationefficiency,definedasintegrationtime/totalobserving time, is 33% in the reference simulation. We keep this value for all the targets for simplicity, because it is opti- mized during the simulation and the simulator can only be ran on individual targets. A future improvement, if available, would be to run the simulator on all the TESS targets. For the atmospheric signal, we calculate the fractional loss of light δ due to molecular absorption from an atmospheric annulus of thickness H: H dδ 2HR δ = ×H = p (2.1) H dR R2 p (cid:63) Thismetricallowsustoestimatethecharacteristicamplitudeofamolecularabsorptionfeaturewithoutrelying on a theoretical transmission spectrum. Instead, we consider a single spectral feature with an amplitude δ , H a width ∆λ = 100 nm, and a central wavelength λ = 1.8 µm corresponding to the center of the NIRISS c SOSS order 1 spectrum. Then, we calculate the observing time that is necessary to detect this 1-scale height amplitude feature with a S/N = 1. In practice, a molecular feature is expected to span several scale heights and a spectrum may contain several features, which would increase the S/N (for example, an amplitude of 3 scale heights would yield S/N =3). These calculations are made for each wavelength tabulated in the NIRISS simulation output and we take the median of the observing times as an estimate for t . We also set a JWST lower limit for t as twice the transit duration, because the observations will span at least one full transit JWST including an out-of-transit baseline. 3 Observing time per target 3.1 SOPHIE and SPIRou We split the target star population into those harbouring giant planets (R > 4R ) and those harbouring p E small planets (R < 4R ) including the terrestrial planets in the habitable zone (R < 2R , S < 2S ), p E p E E where R and S are the Earth radius and insolation respectively. In this definition, the habitable zone has no E E outward limit, but most of the TESS targets will be observed for 30 days only which will reduce the number of long period, cold planets. We consider only the targets that can be observed with an airmass lower than 2 (−16.07◦ <DEC <+90◦ for SOPHIE, −40.17◦ <DEC <+79.83◦ for SPIRou). Figure 1 shows the required observing time t as a function of stellar magnitude. The follow-up of the giant planets could be achieved RV withSOPHIE,withamediant of3.47hpertarget,andatotaltimeof305hafterexcluding17planetswith RV t > 25 h (the total time varies significantly depending on the limit adopted for t ). Some small planets RV RV could also be confirmed by SOPHIE, with a median of 5.92 h and a maximum of 10.71 h for the 10% most favorable cases (103 planets). However, most of the small planets would require an unrealistic time investment, with a median of 116 h per target. On the other hand, SPIRou is ideally suited to the follow-up of the small planets, with a median t of 2.65 h, and a median of only 4.70 h for the terrestrial planets in the habitable RV zone. As limitations, these calculations do not take into account the observational overheads, and constraining the eccentricity will require taking some data points out of quadratures. Thus, the actual observing times will be slightly larger than the estimates given here. 3.2 JWST To calculate H and δ , we consider a H/He-dominated atmosphere for the giant planets (µ=2.32 u, where µ H is the atmospheric mean molecular mass), a H O-dominated atmosphere for the hot and warm small planets 2 (µ = 18 u), and an Earth-like atmosphere for the cold small planets (µ = 29 u). We define hot, warm, and cold according to the planets’ insolation: S > 10S , 2S < S < 10S , and S < 2S respectively (with this E E E E ∗http://maestria.astro.umontreal.ca/niriss/simu1D/simu1D.php 240 SF2A 2016 Fig. 1. Observing time required with SOPHIE (left) and SPIRou (right) for each TESS target as a function of the V and K magnitude, respectively. Giant planets are indicated as open circles, small planets are indicated as filled black s circles, and terrestrial planets in the habitable zone are indicated as filled red circles. See text for definitions. definition, “cold” is equivalent to “habitable zone”). We consider all the targets, with no restriction on DEC. Figure 2 shows the observing time t required to detect the molecular absorption feature. Interestingly, JWST we find that t increases only weakly with the K magnitude for the small planets. This is because R JWST s (cid:63) decreasesasafunctionofK (Mdwarfsarefainterthanmoremassive,largerstars)andasaresultδ increases s H with K (the relative transit depth is larger for small stars). We note that R and R appear uncorrelated for s (cid:63) p the TESS simulated small planets. Whether this only weak correlation between the time required for transit spectroscopyandthestellarmagnitudeisgenerallytrueforplanetsaroundM-dwarfsremainstobeinvestigated. The median observing time t is 13.5 h for the giant planets and is set by twice the transit duration. The JWST small planets have a median of 201 h, and the 10% most favorable ones (184 planets) have a median of 16.2 h andamaximumof28.6h. The10%mostfavorablehabitablezoneterrestrialplanets(7planets)haveamedian of 59.7 h and a maximum of 80.3 h. These planets are highly valuable but will require a significant JWST time investment, thus identifying the best ones with SPIRou is crucial. If we assume H/He-dominated atmospheres for all the planets, the median observing time drops to 8 h for the small planets and is limited by twice the transit duration, and to 13.1 h for the habitable zone terrestrial planets. Fig.2.ObservingtimerequiredwithJWST NIRISSforeachTESS planetasafunctionoftheK magnitude. Thecolor s codeisthesameasinFigure1. Wesetalowerlimitont astwicethetransitduration,whichisthelimitingfactor JWST for the giant planets. We show the full TESS sample, although targets requiring more than ∼100 hours will probably not be observed by JWST. Characterization of the TESS exoplanets 241 4 Preparing radial velocity follow-up programs In this section, we estimate the cumulative observing time that would be necessary with SOPHIE and SPIRou for the radial velocity follow-up of the TESS planet candidates, in order to prepare and optimize observation campaigns. We split the planets into giants (R >4R ), sub-Neptunes and Super-Earths (2R <R <4R ), p E E p E and terrestrial planets (R < 2R ), and into hot, warm, and cold planets as defined in Section 3.2. We sort p E the planets by increasing t and compute its cumulative distribution. The results are shown in Figure 3. We RV exclude planets with t >25 hours, which would extend the cumulative distributions on the right hand side. RV These diagrams can be used to define a follow-up strategy: they give directly the number and type of planets thatcouldbefollowed-upinagivenamountoftelescopetime. PlacingtheTESS candidatesonthesediagrams as they are being discovered (from an initial guess) will also help decide on their follow-up. Finally, they are useful to define a combined strategy between SOPHIE and SPIRou; an obvious split would be to observe the giant planets with SOPHIE and the small planets around M-dwarfs with SPIRou. These estimates are based on the TESS planet simulations: the uncertainties in the number of planets could be as large as 50% (Sullivan et al. 2015) and we do not account for false positive discrimination. As a future improvement, the targets’ potential for observations with JWST could be included quantitatively in defining the radial velocity follow-up samples, for example by means of a unique merit function taking into account t and t . RV JWST Fig. 3. Cumulative number of TESS planets that can be followed-up in radial velocity as a function of the cumulative observing time for SOPHIE (top) and SPIRou (bottom), for the giant planets (left), Sub-Neptunes and Super-Earths (middle), and terrestrial planets (right). We show the full sample (grey), as well as the hot (red), warm (yellow), and cold (blue) samples. The planets are sorted by increasing t . See text for definitions. RV 5 Conclusion The observing time estimates provided in this study will help define the TESS planet candidate samples to observe in radial velocity with SOPHIE and SPIRou as well as their potential for atmospheric characterization withJWST.Practically,thisapproachwillhelpusdecidewhichcandidatestoobservewithSOPHIEandSPIRou while they are being identified by TESS and deliver the best targets for JWST during the process. This could be particularly important in the context of a large SPIRou program at CFHT, the SPIRou Legacy Survey, which would aim at ∼500 nights in the first 5 years of operation. This program would be shared into several science topics and could include the follow-up of the TESS low-mass exoplanet candidates around M-dwarfs. Finally, in the coming years, SOPHIE may be upgraded with a new detector in order to improve its sensitivity and extend its wavelength coverage. This would result in a gain equivalent to one magnitude for mid- and early-M-dwarfs. 242 SF2A 2016 TheDunlapInstituteisfundedthroughanendowmentestablishedbytheDavidDunlapfamilyandtheUniversityofToronto. References Bouchy, F., H´ebrard, G., Udry, S., et al. 2009, A&A, 505, 853 Delfosse, X., Donati, J.-F., Kouach, D., et al. 2013, in SF2A-2013: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, ed. L. Cambresy, F. Martins, E. Nuss, & A. Palacios, 497–508 Ricker,G.R.,Winn,J.N.,Vanderspek,R.,etal.2014,inProc.SPIE,Vol.9143,SpaceTelescopesandInstrumentation 2014: Optical, Infrared, and Millimeter Wave, 914320 Sullivan, P. W., Winn, J. N., Berta-Thompson, Z. K., et al. 2015, ApJ, 809, 77 Weiss, L. M., Marcy, G. W., Rowe, J. F., et al. 2013, ApJ, 768, 14