Astronomy & Astrophysics manuscript no. rbatch c ESO 28 July 23, 28 Three irradiated and bloated hot Jupiters: WASP-76b, WASP-82b & WASP-9b R.G. West, J.-M. Almenara 2, D.R. Anderson 3, F. Bouchy 4, D. J. A. Brown 5, A. Collier Cameron 5, M. Deleuil 2, L. Delrez 6, A. P. Doyle 3, F. Faedi, A. Fumel 6, M. Gillon 6, G. Hébrard 4, C. Hellier 3, E. Jehin 6, M. Lendl 7, P.F.L. Maxted 3, F. Pepe 7, D. Pollacco, D. Queloz 7,8, D. Ségransan 7, B. Smalley 3, A.M.S. Smith 3,9, A.H.M.J. Triaud, 7,, and S. Udry 7 arxiv:3.567v [astro-ph.ep] 2 Oct 23 Department of Physics, University of Warwick, Coventry CV4 7AL, UK 2 Aix Marseille Université, CNRS, LAM (Laboratoire d Astrophysique de Marseille) UMR 7326, 3388, Marseille, France 3 Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK 4 Institut d Astrophysique de Paris, UMR 795 CNRS, Université Pierre & Marie Curie, France; Observatoire de Haute-Provence, CNRS/OAMP, 487, St Michel l Observatoire, France 5 SUPA, School of Physics and Astronomy, University of St. Andrews, North Haugh, Fife, KY6 9SS, UK 6 Institut d Astrophysique et de Géophysique, Université de Liège, Allée du 6 Août, 7, Bat. B5C, Liège, Belgium 7 Observatoire astronomique de l Université de Genève 5 ch. des Maillettes, 29 Sauverny, Switzerland 8 Cavendish Laboratory, J J Thomson Avenue, Cambridge, CB3 HE, UK 9 N. Copernicus Astronomical Centre, Polish Academy of Sciences, Bartycka 8, -76 Warsaw, Poland Department of Physics and Kavli Institute for Astrophysics & Space Research, Massachusetts Institute of Technology, Cambridge, MA 239, USA Preprint online version: July 23, 28 ABSTRACT We report three new transiting hot-jupiter planets discovered from the WASP surveys combined with radial velocities from OHP/SOPHIE and Euler/CORALIE and photometry from Euler and. All three planets are inflated, with radii.7.8 R Jup. All orbit hot stars, F5 F7, and all three stars have evolved, post-ms radii (.7 2.2 R ). Thus the three planets, with orbits of.8 3.9 d, are among the most irradiated planets known. This reinforces the correlation between inflated planets and stellar irradiation. Key words. stars: individual (WASP-76; BD+ 36) stars: individual (WASP-82) stars: individual (WASP-9) planetary systems. Introduction The naive expectation that a Jupiter-mass planet would have a one-jupiter radius has been replaced by the realisation that many of the hot Jupiters found by transit surveys have inflated radii. Planets as large as 2 R Jup have been found (e.g. WASP-7b, Anderson et al. 2; HAT-P-32b, Hartman et al. 2). It is also apparent that inflated planets are found preferentially around hot stars. For example Hartman et al. 22 reported three new HAT-discovered planets, with radii of.6.7 R Jup, all transiting F stars. Here we continue this theme by announcing three new hot Jupiters, again all inflated and all orbiting F stars. For a discussion of the radii of transiting exoplanets we refer the reader to the recent paper by Weiss et al. (23). It seems clear that stellar irradiation plays a large role in inflating hot Jupiters, since no inflated planets are known that receive less than 2 8 erg s cm 2 (Miller & Fortney 2; Demory & Seager 2). There is also an extensive literature discussing other mechanisms for inflating hot Jupiters, such as tidal dissipation (e.g. Leconte et al. 2, and references therein) and Ohmic dissipation (e.g. Batygin & Stevenson 2). Fellow of the Swiss National Science Foundation 2. Observations The three transiting-planet systems reported here are near the equator, and so have been observed by both the SuperWASP- North camera array on La Palma and by WASP-South at Sutherland in South Africa. Our methods all follow closely to those in previous WASP discovery papers. The WASP camera arrays are described in Pollacco et al. (26) while our planethunting methods are described in Collier-Cameron et al. (27a) and Pollacco et al. (27). Equatorial WASP candidates are followed up by obtaining radial velocities using the SOPHIE spectrograph on the.93- m telescope at OHP (as described in, e.g., Hébrard et al. 23) and the CORALIE spectrograph on the.2-m Euler telescope at La Silla (e.g., Triaud et al. 23). Higher-quality lightcurves of transits are obtained using on the.2-m telescope (e.g., Lendl et al. 23) and the robotic photometer at La Silla (e.g., Gillon et al. 23). The observations for our three new planets are listed in Table. 3. The host stars The stellar parameters for WASP-76, WASP-82 and WASP-9 were derived from the co-added RV spectra using the methods given in Doyle et al. (23). The excitation balance of the Fei lines was used to determine the effective temperature (T eff ). The
West et al.: WASP-76b, WASP-82b & WASP-9b Table. Observations Facility Date WASP-76: SuperWASP-North 28 Sep 2 Dec 2 8 points WASP-South 28 Jul 29 Dec 77 points OHP/SOPHIE 2 Sep 2 Dec 9 RVs Euler/CORALIE 22 Feb 22 Dec 8 RVs 2 Nov 6 I filter 22 Aug 25 I filter 22 Oct 3 Gunn r filter 22 Oct 3 I filter 22 Nov 2 I filter WASP-82: SuperWASP-North 28 Oct 2 Feb 5 points WASP-South 28 Oct 2 Jan 86 points OHP/SOPHIE 2 Dec 22 Feb 8 RVs Euler/CORALIE 22 Feb 23 Mar 2 RVs 22 Nov 2 Gunn r filter WASP-9: SuperWASP-North 24 May 2 Oct 4 8 points WASP-South 28 Jun 29 Oct 2 2 points Euler/CORALIE 2 Oct 22 Sep 5 RVs 22 Jun 3 I + z filter 22 Jul 28 Gunn r filter 22 Sep 3 I + z filter 23 Jun Gunn r filter surface gravity (log g) was determined from the ionisation balance of Feiand Feii. The Cailine at 6439Å and the NaiD lines were also used as log g diagnostics. Values of microturbulence (ξ t ) were obtained by requiring a null-dependence on abundance with equivalent width. The elemental abundances were determined from equivalent width measurements of several unblended lines. The quoted error estimates include that given by the uncertainties in T eff and log g, as well as the scatter due to measurement and atomic data uncertainties. The projected stellar rotation velocity (v sin I) was determined by fitting the profiles of several unblended Fe i lines. Macroturbulence was obtained from the calibration by Bruntt et al. (2). For WASP-76, the rotation rate (P=7.6±4. d) implied by the v sin I (assuming that the spin axis is perpendicular to us) gives a gyrochronological age of 5.3 +6. 2.9 Gyr using the Barnes (27) relation. The lithium age of several Gyr, estimated using results in Sestito & Randich (25), is consistent. For WASP-9, the rotation rate (P=.±.6 d) implied by the v sin I gives a gyrochronological age of 4.4 +8.4 2.4 Gyr. The T eff of this star is close to the lithium-gap (Böhm-Vitense, 24), and thus the lack of any detectable lithium in this star does not provide a usable age constraint. WASP-82 is too hot for reliable gyrochronological or lithium ages. We list in Tables 2 4 the proper motions of the three stars from the UCAC4 catalogue (Zacharias et al. 23). These are compatible with the stars being from the local thin-disc population. We also searched the WASP photometry of each star for rotational modulations by using a sine-wave fitting algorithm as described by Maxted et al. (2). For none of the three stars was a significant periodicity found (< mmag at 95%-confidence). 4. System parameters The radial-velocity and photometric data (Table ) were combined in a simultaneous Markov-chain Monte-Carlo (MCMC) analysis to find the system parameters (see Collier Cameron et al. 27b for an account of our methods). For limb-darkening we used the 4-parameter law from Claret (2), and list the resulting parameters in Tables 2 4. For WASP-76b and WASP-82b the radial-velocity data imply circular orbits with eccentricities less than.5 and.6 respectively. WASP-9 is a fainter star and WASP-9b is a lowermass planet, so, while the data are again compatible with a circular orbit, the limit on the eccentricity is weaker at.5. For all three we enforced a circular orbit in the MCMC analysis (see Anderson et al. 22 for the rationale for this). One of the WASP-82 RVs was taken during transit, and this point was given zero weight in the analysis. To translate transit and radialvelocity information (which give stellar density) into the star s mass and radius we need one additional mass radius constraint. Here we use the calibration presented by Southworth (2). The fitted parameters were thus T c, P, F, T 4, b, K, where T c is the epoch of mid-transit, P is the orbital period, F is the fractional flux-deficit that would be observed during transit in the absence of limb-darkening, T 4 is the total transit duration (from first to fourth contact), b is the impact parameter of the planet s path across the stellar disc, and K is the stellar reflex velocity semi-amplitude. The resulting fits are reported in Tables 2 to 4. 5. Discussion The three host stars, WASP-76, WASP-82 and WASP-9, are all F stars with temperatures of 625 65 K. Their metallicities ([Fe/H] =..2) and space velocities are compatible with the local thin-disk population. The stellar densities derived from the MCMC analysis, along with the temperatures from the spectral analysis, are shown on a modified H R diagram in Fig. 4. All three stars have inflated radii (R =.7 2.2 R ) and thus appear to have evolved significantly. The indicated ages of 2 Gyr are compatible with the estimates from gyrochronology (Section 2). The three planets also all have inflated radii (.7.8 R Jup ) and are thus bloated hot Jupiters. This is likely related to the high temperature and large radii of the host stars. Weiss et al. (23) have shown that the radii of hot-jupiter planets correlates well with their irradiation. Our three new planets fit their relationship, and are all at the extreme high-irradiation, high-radius end (their locations on Fig 4 of Weiss et al. being at fluxes 5., 5.2 and 2.9 9 erg cm 2, s, and radii 2., 8.4 and 9.6 R Earth for WASP-76b, WASP-82b and WASP-9b respectively). One caution, however, is that we have a bias against finding unbloated hot Jupiter around evolved, high-radius stars, since the transit depths would be low. Of particular note is that, at V= 9.5 and R P =.8 R Jup, WASP-76 is now the brightest known star transited by a planet larger than.5 R Jup. WASP-82 is not far behind at V=. and R P =.7 R Jup, comparable to WASP-79 (V=., R P =.7 R Jup ; Smalley et al. 22, and KOI-3 (V=., R P =.8 R Jup ; Santerne et al. 22). Thus the new discoveries will be useful for studying bloated hot Jupiters. For example, Triaud (2) suggests that the orbital inclinations of hot Jupiters are a function of system age. Given that radius changes of evolved systems give age constraints, WASP-76 and WASP-82 will be good systems for testing this idea. WASP-82 has a relatively small v sin I (Table 3) for its spectral type, which could indicate a mis-aligned orbit. 2
West et al.: WASP-76b, WASP-82b & WASP-9b Rel. mag Table 2. System parameters for WASP-76. -.6.6.6.8.2.4 WASP.99.98.97 Relative flux.96.95.94.93.92.9.9.92.94.96.98.2.4.6.8. SOPHIE CORALIE - Rel. RV (m s ) BD+ 36 SWASP J463.86+2422. 2MASS 46385+2429 RA = h 46m 3.86s, Dec = +2 42 2. (J2) V mag = 9.5 Rotational modulation < mmag (95%) pm (RA) 46.6 ±.7 (Dec), 39.9 ±.6 mas/yr Stellar parameters from spectroscopic analysis. Spectral type F7 625 ± log g 4.4 ±. ξt (km s ).4 ±. vmac (km s ) 4.±.3 v sin I (km s ) 3.3 ±.6 [Fe/H] +.23 ±. log A(Li) 2.28 ±. Distance 2 ± 2 pc Parameters from MCMC analysis. P (d).89886 ±. T c (HJD) (UTC) 245 67.8557 ±.34 T 4 (d).539 ±.8 +.8 T 2 = T 34 (d).54.3 2 2.89 ±.6 F = RP /R +. b.4.9 +.3 i ( ) 88..6 K (km s ).93 ±.8 γ (km s ).733 ±.2 e (adopted) (<.5 at 3σ) M (M ).46 ±.7 R (R ).73 ±.4 log g (cgs) 4.28 ±.5 +.8 ρ (ρ ).286.8 625 ± MP (MJup ).92 ±.3 +.6 RP (RJup ).83.4 log gp (cgs) 2.8 ±.2 ρp (ρj ).5 ±. a (AU).33 ±.5 T P,A= (K) 26 ± 4 Errors are σ; Limb-darkening coefficients were: z: a =.683, a2 =.349, a3 =.565, a4 =.286 rg : a =.593, a2 =.2, a3 =.327, a4 =.25 Bis. span -.6.8.2.4 Orbital phase Fig.. WASP-76b discovery data: (Top) The WASP data folded on the transit period. (Second panel) The binned WASP data with (offset) the follow-up transit lightcurves (ordered from the top as in Table ) together with the fitted MCMC model. (Third) The SOPHIE and CORALIE radial velocities with the fitted model. (Lowest) The bisector spans; the absence of any correlation with radial velocity is a check against transit mimics. 3
West et al.: WASP-76b, WASP-82b & WASP-9b Table 3. System parameters for WASP-82. Rel. mag -.6.6.6.8.2.4 Relative flux WASP.99.98.94.96.98.2.6 SOPHIE CORALIE - Rel. RV (m s ).4 Bis. span -.6.8.2 Orbital phase Fig. 2. WASP-82b discovery data (as for Fig. ) 4.4 SWASP J4538.56+5338. 2MASS 453856+5338 RA = 4h 5m 38.56s, Dec = + 53 38. (J2) V mag =. Rotational modulation <.6 mmag (95%) pm (RA) 7.5 ±.9 (Dec), 7.7 ±.7 mas/yr Stellar parameters from spectroscopic analysis. Spectral type F5 65 ± 8 log g 4.8 ±.9 ξt (km s ).5 ±. vmac (km s ) 5.±.3 v sin I (km s ) 2.6 ±.9 [Fe/H] +.2 ±. log A(Li) 3. ±.8 Distance 2 ± 3 pc Parameters from MCMC analysis. P (d) 2.75782 ±.3 T c (HJD) (UTC) 245 657.9898 ±.5 T 4 (d).277 ±.2 +.2 T 2 = T 34 (d).56.4.624 ±.2 F = R2P /R2 +.4 b.6. +.4 i ( ) 87.9.9 K (km s ).37 ±.9 γ (km s ) 23.62827 ±.7 e (adopted) (<.6 at 3σ) M (M ).63 ±.8 +.8 R (R ) 2.8.5 +.3 log g (cgs) 3.973.2 +.6 ρ (ρ ).58.4 649 ± MP (MJup ).24 ±.4 +.7 RP (RJup ).67.5 +.7 log gp (cgs) 3.7.32 +.7 ρp (ρj ).266.29 a (AU).447 ±.7 T P,A= (K) 29 ± 4 Errors are σ; Limb-darkening coefficients were: rg : a =.494, a2 =.424, a3 =.266, a4 =.436
West et al.: WASP-76b, WASP-82b & WASP-9b Table 4. System parameters for WASP-9. Rel. mag -.6.6.6.8.2.4 WASP.99 Relative flux.98.97.96.95.96.98.2.4 - Rel. RV (m s ) CORALIE - Bis. span SWASP J227.7+7323.7 2MASS 46385+2429 RA = 2h 2m 7.7s, Dec = +7 3 23.7 (J2) V mag =.7 Rotational modulation < mmag (95%) pm (RA).2 ±.4 (Dec), 8. ± 4.3 mas/yr Stellar parameters from spectroscopic analysis. Spectral type F6 644 ± 3 log g 4.32 ±.9 ξt (km s ).3 ±.2 vmac (km s ) 4.7±.3 v sin I (km s ) 6. ±.5 [Fe/H] +. ±.4 log A(Li) <.7 Distance 34 ± 6 pc Parameters from MCMC analysis. P (d) 3.96243 ±.3 T c (HJD) (UTC) 245 6235.5639 ±.5 T 4 (d).398 ±.22 T 2 = T 34 (d).33 ±.3 F = R2P /R2.7 ±.2 b.84 ±.3 i ( ) 82. ±.4 K (km s ).6 ±.6 γ (km s ) 4.36 ±.3 e (adopted) (<.5 at 3σ) M (M ).55 ±. R (R ).98 ±.9 log g (cgs) 4.33 ±.29 ρ (ρ ).2 ±.2 643 ± 3 MP (MJup ).63 ±.7 RP (RJup ).63 ±.9 log gp (cgs) 2.73 ±.6 ρp (ρj ).45 ±.27 a (AU).562 ±.2 T P,A= (K) 84 ± 5 Errors are σ; Limb-darkening coefficients were: I + z: a =.554, a2 =.4, a3 =.7, a4 =.86 rg : a =.476, a2 =.422, a3 =.226, a4 =.2 -.6.8.2.4 Orbital phase Fig. 3. WASP-9b discovery data (as for Fig. ) 5
West et al.: WASP-76b, WASP-82b & WASP-9b Santerne, A.,Moutou, C., Barros, S. C. C. et al. 22, A&A, 544, L2 Sestito P., Randich S., 25, A&A, 442, 65 Smalley, B., Anderson, D. R., Collier-Cameron, A. et al, 22, A&A, 547, A6 Smith A. M. S. et al. 22, AJ, 43, 8 Southworth J., 2, MNRAS, 47, 266 Triaud, A. H. M. J., 2, A&A, 534, L6 Triaud, A. H. M. J., Anderson, D. R., Collier Cameron, A. et al. 23, A&A, 55, A8 Weiss, L. M., Marcy, G. W.,Rowe, J. F. et al. 23, ApJ, 768, 4 Zacharias, N., Finch, C. T., Girard, T. M., Henden, A., Bartlett, J. L., Monet, D. G., Zacharias, M. I. 23, AJ, 45, 44 Fig. 4. Evolutionary tracks on a modified H R diagram (ρ versus T eff ). The green lines are for a metallicity of [Fe/H]=+.9; the dashed lines are isochrones for.7, 2. and 2.5 Gyr; the solid lines are mass tracks for.3 M and.4 M. The red lines are for a metallicity of [Fe/H]=+., with the same isochrones, and the mass track for.5 M. The models are from Girardi et al. (2). Acknowledgements. SuperWASP-North is hosted by the Issac Newton Group and the Instituto de Astrofísica de Canarias on La Palma while WASP-South is hosted by the South African Astronomical Observatory; we are grateful for their ongoing support and assistance. Funding for WASP comes from consortium universities and from the UK s Science and Technology Facilities Council. is funded by the Belgian Fund for Scientific Research (Fond National de la Recherche Scientifique, FNRS) under the grant FRFC 2.5.594.9.F, with the participation of the Swiss National Science Fundation (SNF). M. Gillon and E. Jehin are FNRS Research Associates. A.H.M.J. Triaud is a Swiss National Science Foundation Fellow under grant PBGEP2-45594 References Anderson, D. R., Collier Cameron, A., Gillon, M. et al. 22, MNRAS, 422, 988 Anderson, D. R., Hellier, C., Gillon, M. et al. 2, ApJ, 79, 59 Asplund, M., Grevesse, N., Sauval, A.J., & Scott, P. 29, ARA&A, 47, 48 Barnes, S.A. 27, ApJ, 669, 67 Batygin, K., Stevenson, D., J., 2, ApJ, 74, L238 Böhm-Vitense, E. 24, AJ, 28, 2435 Bruntt, H., Bedding, T.R., Quirion, P.-O., et al. 2b, MNRAS, 45, 97 Claret, A., 2, A&A, 363, 8 Collier Cameron, A., et al., 27a, MNRAS, 375, 95 Collier Cameron, A., et al., 27b, MNRAS, 38, 23 Demory, B.-O., Seager, S., 2, ApJS, 97, 2 Doyle, A.P., Smalley, B., Maxted, P.F.L., et al. 23, MNRAS, 428, 364 Gillon, M., Anderson, D. R., Collier-Cameron, A. et al. 23, A&A, 552, A82 Girardi, L., Bressan, A., Bertelli, G., Chiosi, C. 2, A&AS, 4, 37 Hartman, J. D., Bakos, G. Á., Béky, B. et al. 22, AJ, 44, 39 Hartman, J. D., Bakos, G. Á., Torres, G. et al. 2, ApJ, 742, 59 Hébrard, G., Collier Cameron, A., Brown, D. J. A. et al. 23, A&A, 549, A34 Leconte, J., Chabrier, G., Baraffe, I., Levrard, B. 2, A&A, 56, A64 Lendl, M., Anderson, D. R., Collier-Cameron, A. et al. 22, A&A, 544, A72 Maxted, P.F.L. et al. 2, PASP, 23, 547 Miller, N., Fortney, J. J. 2, ApJ, 736, L29 Pollacco, D., et al., 26, PASP, 8, 47 Pollacco, D., et al., 28, MNRAS, 385, 576 6