Astronomy & Astrophysics manuscript no. w4 6 ACCEPTED c ESO 28 August, 28 WASP-4b and WASP-6b: two transiting hot Jupiters in.75-day and 9.3-day orbits A. M. S. Smith,2, D. R. Anderson 2, D. J. Armstrong 3, S. C. C. Barros 4, A. S. Bonomo 5, F. Bouchy 4, D. J. A. Brown 3, A. Collier Cameron 6, L. Delrez 7, F. Faedi 3, M. Gillon 7, Y. Gómez Maqueo Chew 3, G. Hébrard 8,9, E. Jehin 7, M. Lendl 7,, T. M. Louden 3, P. F. L. Maxted 2, G. Montagnier 8,9, M. Neveu-VanMalle,, H. P. Osborn 3, F. Pepe, D. Pollacco 3, D. Queloz,, J. W. Rostron 3, D. Segransan, B. Smalley 2, A. H. M. J. Triaud,2, O. D. Turner 2, S. Udry, S. R. Walker 3, R. G. West 3, and P. J. Wheatley 3 arxiv:48.887v3 [astro-ph.ep] Sep 24 N. Copernicus Astronomical Centre, Polish Academy of Sciences, Bartycka 8, -76 Warsaw, Poland. email:amss@camk.edu.pl 2 Astrophysics Group, Lennard-Jones Laboratories, Keele University, Keele, Staffordshire, ST5 5BG, UK 3 Department of Physics, University of Warwick, Coventry CV4 7AL, UK 4 Aix Marseille Université, CNRS, LAM (Laboratoire d Astrophysique de Marseille) UMR 7326, 3388, Marseille, France 5 INAF - Osservatorio Astrofisico di Torino, via Osservatorio 2, I-25 Pino Torinese, Italy 6 SUPA, School of Physics and Astronomy, University of St. Andrews, North Haugh, St. Andrews, Fife, KY6 9SS, UK 7 Institut d Astrophysique et de Géophysique, Université de Liège, Allée du 6 Août 7, Sart Tilman, Liège, Belgium 8 Institut d Astrophysique de Paris, UMR795 CNRS, Université Pierre & Marie Curie, 98bis boulevard Arago, 754, Paris, France 9 Observatoire de Haute-Provence, CNRS/OAMP, 487 Saint-Michel-l Observatoire, France Observatoire de Genève, Université de Genève, 5 Chemin des Maillettes, 29 Sauverny, Switzerland Cavendish Laboratory, J J Thomson Avenue, Cambridge CB3 HE, UK 2 Kavli Institute for Astrophysics & Space Research, Massachusetts Institute of Technology, Cambridge, MA 239, USA ABSTRACT We report the discovery from the WASP survey of two exoplanetary systems, each consisting of a Jupiter-sized planet transiting an th magnitude (V) main-sequence star. WASP-4b orbits its star in.75 d, whereas WASP-6b has the fourth-longest orbital period of any planet discovered by means of transits observed from the ground, orbiting every 9.29 d. Each planet is more massive than Jupiter (WASP-4b has a mass of.27±.5 M Jup, while WASP-6b has a mass of.93±.8 M Jup ). Both planets are just slightly larger than Jupiter, with radii of.4±.4 and.9±.4 R Jup for WASP-4 and WASP-6 respectively. No significant orbital eccentricity is detected in either system, and while this is not surprising in the case of the short-period WASP-4b, it is interesting in the case of WASP-6b, because many otherwise similar planets are known to have eccentric orbits. Key words. planetary systems planets and satellites: detection planets and satellites: fundamental parameters stars: individual: WASP-4 stars: individual: WASP-6. Introduction Transiting planets are vital for our understanding of planetary systems, as planetary radii and absolute masses can be measured. Recently, the large number of planets and planetary candidates discovered by the Kepler satellite has extended the parameter space of transiting planet discovery and led to major advances in our understanding of the statistics of the Galactic planetary population. Planets transiting bright stars are required to increase our knowledge of the range of properties exhibited by the nearby planetary population, as well as to make advances in our understanding of planetary formation and evolution. Examples of planet characterisation that require bright target stars include high-precision radial velocity measurements to measure the orbital obliquity, and measurement of planetary transmission and emission spectra to infer atmospheric properties (e.g. Winn 2). The Wide Angle Search for Planets (WASP; Pollacco et al. 26) and other wide-field ground-based surveys such as HATnet (Bakos et al. 22) are the leading discoverers of bright (8.5 V 2.5) transiting systems, with around 5 planets between them. Only a handful of the nearly one hundred planets discovered by WASP have orbital periods, P, greater than 5 d. This is because there is an intrinsic pile-up of planets at 3 5 days, and hence a relative dearth of planets at longer periods. This is further exacerbated by the fact that at long periods the probability that a planetary system is aligned such that transits are visible from Earth is significantly reduced. Furthermore, because a large number of transits must be observed to overcome the correlated noise present in wide-field survey photometry (e.g. Smith et al. 26), very long observational baselines are usually required to find longer-period systems. Recently Lendl et al. (24) reported the discovery of WASP-7b, the first WASP discovery with P > 9 d. Only two other systems with P > 9 d have been discovered from the ground by means of transits: HAT-P-5b (Kovács et al. 2) and HAT-P-7b (Howard et al. 22). Here we report the discovery of the fourth such planet, WASP-6b (P=9.29 d), along with a much shorter-period planet, WASP-4b (P =.75 d). Both planets orbit V= stars lying in the constellation Leo.
A. M. S. Smith et al.: WASP-4b and WASP-6b 2. Observations 2.. WASP photometry Both WASP-4 and WASP-6 were observed for four seasons by SuperWASP-N at the Observatorio del Roque de los Muchachos on La Palma, Spain, from 28 February 5 to 2 March 29. Each of the targets was also observed for two seasons by WASP-South at the South African Astronomical Observatory near Sutherland, South Africa, from 29 January 6 to 2 May 3. These observations resulted in 24632 data points on WASP-4 (3993 from SuperWASP-N and 639 from WASP- South) and 23624 measurements of WASP-6 (3 from SuperWASP-N and 523 from WASP-South). The WASP instruments consist of eight Canon 2mm f/.8 lenses, each equipped with an Andor 248 248 e2v CCD camera, on a single robotic mount. Further details of the instrument, survey and data reduction procedures are described in Pollacco et al. (26) and details of the candidate selection procedure can be found in Collier Cameron et al. (27) and Pollacco et al. (28). The data revealed the presence of transit-like signals with a period of around.75 days in the case of WASP-4, and around 9.29 days for WASP-6. The combined WASP light curves are shown binned and folded on the best-fitting orbital periods in the upper panels of Figures and 2. The light curves were also searched for additional transiting components with periods up to 2 d, using the method described in Smith et al. (29), but no such signals were found. 2.2. Spectroscopic follow-up Spectroscopic observations of both targets were conducted using the spectrograph mounted on the.93-m telescope of the Observatoire de Haute-Provence, France, and with on the.2-m Euler-Swiss telescope at La Silla, Chile. A total of 2 observations ( with, with ) were made of WASP-4 between 23 January 8 and 23 June 25. WASP-6 was observed 29 times (2 times with, 9 with ) between 23 January 6 and 24 February 8. was used in high-efficiency (R = 4) mode, with the CCD in slow-readout mode. A second fibre was placed on the sky to check for sky background contamination in the spectra of the target. Corrections of up to a few tens of m s were applied to some exposures contaminated by moonlight. More information about the instrument can be found in Bouchy et al. (29). The observations were conducted only in dark time to avoid moon light entering the fibre. Observations of thorium-argon emission line lamps were used to calibrate both the and spectra. The data were processed using the standard and data-reduction pipelines. The resulting radial velocity data are listed in Table and plotted in the lower panels of Figures and 2. In order to rule out the possibility that either system is a blended eclipsing binary, or that the RV variation is a result of stellar activity, we examined the bisector spans (e.g. Queloz et al. 2). For both stars, the bisector spans exhibit no significant correlation with radial velocity (Figure 3), as expected for true planetary systems. Table. Radial velocities of WASP-4 and WASP-6 Star BJD(UTC) RV σ RV BS Inst. 245 km s km s km s WASP-4 6364.666 28.386.4 -.82 COR. WASP-4 6365.6744 28.566.2 -.53 COR. WASP-4 6368.745 28.752.3 -.3 COR. WASP-4 637.672 28.45.3 -.27 COR. WASP-4 64.4762 28.523.9 -.57 COR. WASP-4 6438.4853 28.487.9 -.62 COR. WASP-4 6442.4977 28.76.8 -.4 COR. WASP-4 6444.5289 28.722. -.55 COR. WASP-4 645.574 28.336. -.69 COR. WASP-4 6469.4597 28.443.2 -.48 COR. WASP-4 63.5794 28.968.6 -.28 SOPH. WASP-4 6329.672 28.752. -.57 SOPH. WASP-4 633.5262 28.96. -.9 SOPH. WASP-4 633.582 28.898. -.35 SOPH. WASP-4 636.484 28.98.6 -.23 SOPH. WASP-4 6363.525 29.8.6 -.59 SOPH. WASP-4 6377.364 28.964.4 -.36 SOPH. WASP-4 6378.372 28.647.3 -.2 SOPH. WASP-4 64.442 28.646.7 -.45 SOPH. WASP-4 643.3674 28.737. -.35 SOPH. WASP-4 6423.3975 28.989.3 -.8 SOPH. WASP-6 6364.6843 7.27.24.65 COR. WASP-6 6365.7257 7.26.2 -.45 COR. WASP-6 6366.6264 7.353.46.5 COR. WASP-6 6368.7463 7.386.2.6 COR. WASP-6 6369.7726 7.363.2.4 COR. WASP-6 637.7735 7.268.22.94 COR. WASP-6 637.5952 7.54.9.5 COR. WASP-6 637.769 7.48.2.5 COR. WASP-6 6372.7689 7.75.8 -.22 COR. WASP-6 6373.784 7.25.2.58 COR. WASP-6 64.55 7.4.38.45 COR. WASP-6 6423.5873 7.39.4 -.28 COR. WASP-6 6424.5533 7.45.4.5 COR. WASP-6 6438.533 7.84.5 -.4 COR. WASP-6 644.4779 7.345.23 -.53 COR. WASP-6 6443.526 7.393.2.39 COR. WASP-6 6449.5284 7.234.6.7 COR. WASP-6 6469.4842 7.338.8.6 COR. WASP-6 648.4663 7.355.24 -.4 COR. WASP-6 6696.769 7.75.6 -.2 COR. WASP-6 6298.699 7.33.4.27 SOPH. WASP-6 6329.6333 7.242.2 -.3 SOPH. WASP-6 633.5696 7.347.7 -.35 SOPH. WASP-6 633.597 7.354.7.39 SOPH. WASP-6 636.495 7.34.3. SOPH. WASP-6 6363.5345 7.8.3 -. SOPH. WASP-6 6377.3763 7.38.67 -.99 SOPH. WASP-6 6378.3889 7.356.24.4 SOPH. WASP-6 643.3926 7.8.8.53 SOPH. Spectrograph used. COR.= on the Euler-Swiss telescope; SOPH.= on the.93-m telescope at the OHP. 2.3. Photometric follow-up High-precision photometric follow-up observations were made from ESO s La Silla Observatory in Chile and from the Observatorio del Roque de los Muchachos on La Palma, Spain. The telescopes used were the robotic Transiting Planets and Planetesimals Small Telescope (TRAPPIST; Jehin et al. 2) and the Euler-Swiss telescope with EulerCam (Lendl et al. 22) 2
A. M. S. Smith et al.: WASP-4b and WASP-6b Table 2. Observing log for follow-up transit photometry No. Date/UT Telescope Band full? WASP-4: (i) 23 Mar 29/3 TRAPPIST I + z Y (ii) 23 Apr 5/6 TRAPPIST I + z Y (iii) 23 Apr 2/3 TRAPPIST I + z Y (iv) 23 Jun 2/3 TRAPPIST I + z N (v) 23 Apr 2 Euler I Y (vi) 23 Apr 27 Euler I Y WASP-6: (vii) 23 May 3/4 Euler Gunn r N (viii) 23 May 3/4 TRAPPIST I + z N (ix) 24 Feb 25/26 TRAPPIST I + z N (x) 24 Feb 25/26 Liverpool V + R Y (xi) 24 Feb 25/26 Isaac Newton i N Indicates whether or not the complete transit was observed at La Silla; and the Liverpool Telescope (LT) with the Rapid Imaging Search for Exoplanets camera (RISE; Steele et al. 28; Gibson et al. 28) and the Isaac Newton Telescope (INT) with the Wide Field Camera (WFC) on La Palma. A summary of these observations is given in Table 2. Due to the relatively long transit duration of over five hours, it proved difficult to observe a full transit of WASP-6b. This was, however, achieved using the LT on the night of 24 February 25/26. Simultaneous observations from the nearby INT were stymied by dome-shutter problems during ingress, preventing observation of the full transit, although a partial transit was observed. 3. Determination of system parameters 3.. Stellar parameters The individual spectra were co-added to produce a single spectrum of each star with an average S/N of around 6: for WASP-4 and : for WASP-6. The analyses were performed using standard pipeline reduction products and the procedures given in Doyle et al. (23). The projected stellar rotation velocities, v sin i, were determined by fitting the profiles of several unblended Fe i lines. We assumed macroturbulent velocities (v mac ) of 2.5±.6 km s for WASP-4 and 4.7±.6 km s for WASP-6 using the asteroseismic-based calibration of Doyle et al. (24) and an instrumental resolution of 55. The resultant stellar parameters are given in Table 3. The lithium abundance of WASP-6 suggests an age of around 2 Gyr. The v sin i gives an upper limit to the age of.4 +.7.4 Gyr from the gyrochronological relation of Barnes (27). Our isochrone analysis (Fig. 6) gives an age estimate of 7±2 Gyr. Such discrepancies in the ages estimated from these two methods are not unusual, particularly for estimates derived from isochrone fitting (e.g. Brown 24). 3..3. Stellar distance We calculated the distances to WASP-4 and WASP-6 (see Table 3) using Tycho apparent V-band magnitudes of.2 ±. and.2 ±. respectively. The stellar luminosities were calculated from our best-fitting stellar effective temperatures and radii (Table 4), and we adopted bolometric corrections (from Flower 996) of.5 and.5, for WASP-4 and WASP- 6 respectively. We use the interstellar Na D lines and the calibration of Munari & Zwitter (997) to adopt colour excesses, E(B V) =.4±. and E(B V) =.2±. for WASP-4 and WASP-6 respectively, to correct for interstellar reddening. Table 3. Stellar parameters Parameter WASP-4 WASP-6 RA (J2.) h42m24.6s h5m43.3s Dec (J2.) +7 26 6.3 5 4 45.9 Identifier USNO-B. 974-234922 TYC 4927-63- T,eff / K 545±3 6±5 log g (cgs) 4.5±.2 4.±.2 ξ t / km s.9±..2±. v mac / km s 2.5±.6 4.7±.6 v sin i / km s.4±.7 6.3±.7 log A(Li) <.7 2.45 ±.2 Sp. Type G8 F9 Distance/pc 43± 283±2 Mass/M.2±.9.27±.5 Radius/R.93±.23.85±.52 [Fe/H] +.32±.9.9±.9 V-mag (Tycho).2 ±..2 ±. Note: The spectral type was estimated from T,eff using the table of (Gray 28, p. 57). The mass and radius were estimated using the Torres et al. (2) calibration. 3... Stellar activity We searched the WASP photometry for rotational modulations, using the sine-wave fitting algorithm described by Maxted et al. (2). We estimated the significance of periodicities by subtracting the fitted transit light curve and then repeatedly and randomly permuting the nights of observation. No significant periodic modulation was detected, to a 95 per cent confidence upper limit of 4 mmag in the case of WASP-4, and mmag for WASP-6. 3..2. Stellar age The low lithium abundance for WASP-4 suggests that the star is old (several Gyr, perhaps more than 5 Gyr; Sestito & Randich 25), and the low v sin i is consistent with this. From an analysis (Fig. 5) performed using the isochrones of Bressan et al. (22), we infer an age of around 3±2 Gyr. 3.2. Planetary system parameters For each system the photometric data were combined with the radial velocities and analysed simultaneously using the Markov Chain Monte Carlo (MCMC) method. We use the current version of the MCMC code described in Collier Cameron et al. (27) and Pollacco et al. (28), which uses the relation between stellar mass, effective temperature, density and metallicity appropriate for stars with masses between.2 and 3. M from Southworth (2). Briefly, the radial velocity data are modelled with a Keplerian orbit and the photometric transits are fitted using the model of Mandel & Agol (22) and limb-darkening was accounted for using a fourcoefficient, non-linear model, employing coefficients appropriate to the passbands from the tabulations of Claret (2, 24). The coefficients were determined using an initial interpolation in log g and [Fe/H] (values from Table 3), and an interpolation in 3
A. M. S. Smith et al.: WASP-4b and WASP-6b Relative flux.5.95.4.6.8.2.4.6 Relative flux.5.95.4.6.8.2.4.6. (i) TRAPPIST.99 (vii) EulerCam.98 (ii) TRAPPIST.98.96.97 (viii) TRAPPIST Relative flux.94 (iii) TRAPPIST Relative flux.96.95 (ix) TRAPPIST (iv) TRAPPIST.94.92 (v) EulerCam.93 (x) LT/RISE.9.92 (vi) EulerCam.9 (xi) INT/WFC.88 25.95.5.9 25.98.99..2 Relative radial velocity / m s - 2 5 5-5 - -5 Relative radial velocity / m s - 2 5 5-5 - -5-2 -2-25.4.6.8.2.4.6 Orbital phase Fig.. WASP-4 photometry and radial velocities. Upper panel: Phase-folded photometry from SuperWASP-N and WASP-South. Data are binned in phase with a bin-width equal to 2 s. Middle panel: High-precision follow-up photometry. Each light curve is offset in flux for clarity, and binned using a 2 s bin width. The light curves are identified with the same numbers as those in Table 2. Our best-fitting model is plotted with a solid line and TRAPPIST s meridian flips are marked with dashed vertical lines. Lower panel: Radial-velocity measurements from and, along with our bestfitting MCMC solution (solid line). The jitter added to the radial velocity uncertainties is not shown. The centre-of-mass velocity,γ=28.54842 km s, has been subtracted, as has the fitted offsetγ COR SOPH = 28.23 m s between the two datasets. 4-25.4.6.8.2.4.6 Orbital phase Fig. 2. WASP-6 photometry and radial velocities. Upper panel: Phase-folded photometry from SuperWASP-N and WASP-South. Data are binned in phase with a bin-width equal to 2 s. Middle panel: High-precision follow-up photometry. Each light curve is offset in flux for clarity, and binned using a 2 s bin width. The light curves are identified with the same numbers as those in Table 2. Our best-fitting model is plotted with a solid line and TRAPPIST s meridian flips are marked with dashed vertical lines. Lower panel: Radial-velocity measurements from and, along with our best-fitting MCMC solution (solid line). The point taken during transit and excluded from our MCMC analysis is shown in grey. The centreof-mass velocity,γ=7.24744 km s, has been subtracted, as has the fitted offsetγ COR SOPH = 59.57 m s between the two datasets.
A. M. S. Smith et al.: WASP-4b and WASP-6b T,eff at each MCMC step. The coefficient values corresponding to the best-fitting value of T,eff are given in Table 5. The MCMC proposal parameters we use are: the epoch of mid-transit, T c ; the orbital period, P; the transit duration, T 4 ; the fractional flux deficit that would be observed during transit in the absence of stellar limb-darkening, F; the transit impact parameter for a circular orbit, b=a cos(i P )/R ; the stellar reflex velocity semi-amplitude, K ; the stellar effective temperature, T,eff, and the stellar metallicity, [Fe/H]. When fitting for a non-zero orbital eccentricity, there are an additional two proposal parameters: e cosω and e sinω, where e is the orbital eccentricity, and ω is the argument of periastron (Anderson et al. 2). The best-fitting system parameters are taken to be the median values of the posterior probability distribution, and are recorded in Table 4. 3.2.. WASP-4 There is a star of approximately equal brightness to WASP- 4 (USNO-B. R =.82, c.f..28 for WASP-4), located 29 to the North. We found from our observations with the Euler Telescope that this neighbouring object (TYC 26-73-) is itself a visual binary, with the components separated from each other by around. Further, we found this binary to be unrelated to WASP-4, which has a systemic radial velocity of 28 km s compared to -3 km s for the binary. This conclusion is supported by proper motions from the UCAC4 catalogue (Zacharias et al. 22), which indicate that it is unlikely that WASP-4 (µ α = 3.3±2.7 mas yr ; µ δ = 2.3±.6 mas yr ) and the neighbouring object (µ α = 4.4± 6.6 mas yr ;µ δ = 4.8±2.2 mas yr ) are co-moving. Although this object is easily excluded from the photometric apertures used in our follow-up observations, it falls within the 3.5-pixel (= 48 ) WASP aperture, thus rendering the WASP photometry unreliable for measuring the transit depth. After determining the best-fitting orbital period from an MCMC analysis of all the available data, we exclude the WASP photometry from subsequent analysis. Instead, we fix the orbital period to the value previously determined and use MCMC to analyse only the follow-up photometry and the radial velocities. An initial fit for an eccentric orbit gave e =.92 +.23.66 (withω=5 +6 25 degrees). This is consistent with a circular orbit; according to the F-test of Lucy & Sweeney (97) there is a > 99 per cent chance that the apparent eccentricity could have arisen if the underlying orbit were circular. We also tried fitting for a drift in the radial velocities by fitting for a linear function in time. Such a trend may be indicative of a third body in the system, but fitting for radial acceleration, dγ/dt, is clearly not justified in this case, because we obtain dγ/dt= 24±28 m s yr (Figure 4). Linear functions of time were fitted to each light curve at each step of the MCMC, to remove systematic trends. We added jitter terms of 5 m s for and 6 m s for to the formal RV uncertainties in quadrature so as to obtain a spectroscopic reducedχ 2 of unity. 3.2.2. WASP-6 One of the radial velocity measurements was made during transit, and we exclude this point from our MCMC analysis because we do not model the Rossiter-McLaughlin effect. The excluded point is, however, plotted (in grey) in the lower Bisector span / m s - 2 - -2 WASP-4 WASP-6-2 - 2-2 - 2 Relative radial velocity / m s - Fig. 3. Radial velocity bisector span vs. relative radial velocity for WASP-4 (left) and WASP-6 (right). The indicated uncertainties in the radial velocities do not include the added jitter, and the uncertainties in the bisector span are taken to be twice the uncertainty in the radial velocity. The point taken during transit and excluded from our MCMC analysis is shown in grey. Residual radial velocity / m s - 4-4 2 6-6 -2 63 635 64 645 HJD - 245 (UTC) WASP-4 WASP-6 63 64 65 66 67 Fig. 4. Residual radial velocity as a function of time for WASP- 4 (upper) and WASP-6 (lower). There is no evidence of radial acceleration in either system. The indicated uncertainties in the radial velocities do not include the added jitter. The point taken during transit and excluded from our MCMC analysis is shown in grey. panel of Figure 2. It was not necessary to add any jitter to the RV uncertainties. Fitting for an eccentric orbit gives e=.3 +.32.22 (withω= 3 +3 82 degrees). This is again consistent with a circular orbit (the Lucy & Sweeney (97) F-test probability is 2 per cent). The best-fitting radial acceleration is dγ/dt= +3 8 m s yr, which is also consistent with zero (Figure 4). We therefore adopt, as for the WASP-4 system, a circular orbit and a constant systemic radial velocity. Linear functions of time were fitted to each light curve at each step of the MCMC, to remove systematic trends. In the case of the LT light curve (x), the function fitted was quartic, rather than linear in time. A function of CCD position, airmass and sky background instead of time was fitted to the INT light curve (xi). In each case, the addition of extra terms was justified by an improved value of the Bayesian Information Criterion (BIC; Schwarz 978). 5
Table 4. System parameters for WASP-4 and WASP-6 A. M. S. Smith et al.: WASP-4b and WASP-6b Parameter Symbol Unit WASP-4 WASP-6 Model parameters: Orbital period P d.755437 +.8 9.28975±..36 Epoch of mid-transit T c HJD, UTC 245646.26±.2 2456649.54972±.33 Transit duration T 4 d.7342±.56.22346 +.4.9 Planet-to-star area ratio F= R 2 P /R2....474±.2.642±.8 Transit impact parameter b....724±.4.27 +.49.87 Stellar orbital velocity semi-amplitude K m s 22.7±4.3 65.3±4.3 System velocity γ km s 28.5484±. 7.24744±.9 Velocity offset between and γ COR SOPH m s 28.2±. 59.57±.3 Stellar effective temperature T,eff K 5475±27 655±36 Stellar metallicity [Fe/H] dex +.32 ±.9.9 ±.9 Derived parameters: Ingress/egress duration T 2 = T 34 d.535±.75.685 +.8.39 Orbital inclination angle i P 83.63±.25 89.49 +.35.64 Orbital eccentricity (adopted) e... Orbital eccentricity (3-σ upper-limit).......6.3 Stellar mass M M.76±.49.92±.54 Stellar radius R R.963±.27.393 +.48.28 log (stellar surface gravity) log g (cgs) 4.53±.6 4.226 +..22 Stellar density ρ ρ.27±.7.445 +..39 Planet mass M P M Jup.272±.47.925±.76 Planet radius R P R Jup.37±.37.85 +.46.28 log (planet surface gravity) log g P (cgs) 3.353±.23 3.57 +.2.3 Planet density ρ P ρ J.866±.7.5 +..6 Scaled orbital major semi-axis a/r... 6.52±.3 4.2 +..43 Orbital major semi-axis a au.298±.45.97±.4 Planet equilibrium temperature T P,A= K 56±39 4±29 (uniform heat redistribution) The orbital period of WASP-4b and its uncertainty was determined from a fit to all the data (including photometry from the WASP instruments). The period was fixed in subsequent analyses (which excluded the WASP photometry). See Sec 3.2.. The following constant values are used: au=.49598 m; R = 6.9599 8 m; M =.9892 3 kg; R Jup = 7.492 7 m; M Jup =.89896 27 kg;ρ J = 24.67 kg m 3. Table 5. Stellar limb-darkening coefficients Star Instruments Observation bands Claret band a a 2 a 3 a 4 WASP-4 TRAPPIST Cousins I+ Sloan z Sloan z.82.828.22.527 WASP-4 Euler Cousins I Cousins I.82.835.38.573 WASP-6 WASP/ Euler/ LT Broad ( 4 7 nm)/ Gunn r/ V + R Cousins R.497.273.34.96 WASP-6 TRAPPIST Cousins I+ Sloan z Sloan z.592.34.337.99 WASP-6 INT Sloan i Sloan i.559.33.223.62 4. Results and Discussion We find that WASP-4 and WASP-6 are both main-sequence stars. WASP-4 has a mass and radius similar to that of the Sun, while WASP-6 is around 2 per cent more massive than the Sun, and has a radius around 4 per cent larger. Both planets are more massive than Jupiter (WASP-4b has a mass of.27±.5 M Jup, while WASP-6b has a mass of.93±.8 M Jup ). WASP-4b and WASP-6b are not bloated like many hot Jupiters. Indeed with radii of.4±.4 and.9±.4 R Jup respectively, neither planet is much larger than Jupiter. 4.. The circular nature of the orbit of WASP-6b It is interesting to compare WASP-6 to WASP-7 (Lendl et al. 24), a recently discovered transiting system consisting of a hot Jupiter in a. d orbit that is both significantly eccentric (e=.32±.23) and mis-aligned (ψ=69.6 +4.7 4. deg; β = 44± deg). By contrast, WASP-6b orbits in a circular or near-circular orbit (Sec. 3.2.2). The stars WASP-6 and WASP-7 (M =.26±.29M, R =.7 +.67.59 R, T,eff = 638±88 K) are also very similar to each other. Such planets may offer important evidence as to which of two proposed planetary migration scenarios is dominant. Diskdriven migration should produce giant planets in circular orbits with low obliquities, whereas planet-planet scattering and migration through the Kozai mechanism should lead to highly eccentric orbits, with large obliquities (see Lendl et al. 24 and references therein). These eccentric orbits may be subsequently circularised through tidal interactions with the star. The vast majority of known hot Jupiters lie in circular or near-circular orbits, and distinguishing between established circular orbits and once-eccentric orbits that have been relatively quickly made circular by the huge tidal forces at play is non- 6
A. M. S. Smith et al.: WASP-4b and WASP-6b (ρ * /ρ O ) -/3..9 6. Gyr 55 T *,eff / K 2. Gyr. M O.9 M O 5 Fig. 5. Modified Hertzsprung-Russell diagram for WASP-4 (black circle). Isochrones (solid green lines) from the models of Bressan et al. (22) are plotted for ages between. and 2. Gyr, spaced at. Gyr intervals. Also shown are evolutionary tracks (dashed orange lines) for.9,. and. M from Bertelli et al. (28). (ρ * /ρ O ) -/3.6.5.4.3.2 65 4. Gyr. M O 6 T *,eff / K.9 M O 2. Gyr 55 Fig. 6. Modified Hertzsprung-Russell diagram for WASP-6 (black circle). Isochrones (solid green lines) from the models of Bressan et al. (22) are plotted for ages between 4. and 2. Gyr, spaced at. Gyr intervals. Also shown are evolutionary tracks (dashed orange lines) for.9,. and. M from Bertelli et al. (28). trivial. Hot Jupiters with periods longer than 5 d, where the observed pile-up lies, may offer a clue to this puzzle, because the tidal forces acting to circularise the orbit are lesser. Using Equation () of Jackson et al. (28), we calculate the circularisation time-scale,τ e = ( ) de, e dt for the orbit of WASP- 6b, in terms of the tidal dissipation parameters for the planet, Q P, and for the star, Q,.46 τ e = ( QP )+.23 ( Q ) 5.5 6.5 Gyr, () Adopting Q P = 5.5 and Q = 6.5 (the best-fitting values from the study of Jackson et al. 28), we findτ e = 2.7 Gyr. Even in the extreme case that Q P = Q = 5, the circularisation time-scale is still as long as 4.6 Gyr, comparable to the age of the system. Hence it seems probable that the orbit of WASP- 6b has not been circularised from a highly eccentric starting point, but rather that the orbit has remained close to circular for the lifetime of the system. A large fraction of similar planets, however, lie on significantly eccentric orbits. Using the Exoplanets Data Explorer (Wright et al. 2), we find that there are 34 known planets of at least half the mass of Jupiter with orbital periods in the range 9 P d (including WASP-6b). Around half of these (8) are known to have significantly eccentric (e.) orbits, while the remainder have near-circular orbits, or their eccentricity is unknown. It would be interesting to see if the obliquity of the WASP- 6 system is also low, as well as the eccentricity. Measuring the sky-projected angle between the stellar spin axis and the orbital axis through the Rossiter-McLaughlin effect would give an insight into this. Such a measurement should prove possible; we predict an amplitude of around 26 m s for the Rossiter- McLaughlin effect, and a typical RV precision of 4 m s per -minute exposure with HARPS. Acknowledgements. We acknowledge the significant contributions made to the construction and operation of WASP-South by Coel Hellier, who also identified WASP-4 and WASP-6 as candidate planetary systems. WASP-South is hosted by the SAAO and SuperWASP by the Isaac Newton Group and the Instituto de Astrofísica de Canarias; we gratefully acknowledge their ongoing support and assistance. Funding for WASP comes from consortium universities and from the UK s Science and Technology Facilities Council (STFC). TRAPPIST is a project funded by the Belgian Fund for Scientific Research (FNRS) with the participation of the Swiss National Science Foundation. MG and EJ are FNRS Research Associates. The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias with financial support from the UK STFC. The Isaac Newton Telescope is operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias; observations were made under programme I/24A/9. AHMJT is a Swiss National Science Foundation fellow under grant number P3P2-47773. References Anderson, D. R., Collier Cameron, A., Hellier, C., et al. 2, ApJ, 726, L9 Bakos, G. Á., Lázár, J., Papp, I., Sári, P., & Green, E. M. 22, PASP, 4, 974 Barnes, S. A. 27, ApJ, 669, 67 Bertelli, G., Girardi, L., Marigo, P., & Nasi, E. 28, A&A, 484, 85 Bouchy, F., Hébrard, G., Udry, S., et al. 29, A&A, 55, 853 Bressan, A., Marigo, P., Girardi, L., et al. 22, MNRAS, 427, 27 Brown, D. J. A. 24, MNRAS, 442, 844 Claret, A. 2, A&A, 363, 8 Claret, A. 24, A&A, 428, Collier Cameron, A., Wilson, D. M., West, R. G., et al. 27, MNRAS, 38, 23 Doyle, A. P., Davies, G. R., Smalley, B., Chaplin, W. J., & Elsworth, Y. 24, MNRAS, in press, arxiv:48.3988 Doyle, A. P., Smalley, B., Maxted, P. F. L., et al. 23, MNRAS, 428, 364 Flower, P. J. 996, ApJ, 469, 355 Gibson, N. P., Pollacco, D., Simpson, E. K., et al. 28, A&A, 492, 63 Gray, D. F. 28, The Observation and Analysis of Stellar Photospheres, 3rd Edition (Cambridge University Press) Howard, A. W., Bakos, G. Á., Hartman, J., et al. 22, ApJ, 749, 34 Jackson, B., Greenberg, R., & Barnes, R. 28, ApJ, 678, 396 Jehin, E., Gillon, M., Queloz, D., et al. 2, The Messenger, 45, 2 Kovács, G., Bakos, G. Á., Hartman, J. D., et al. 2, ApJ, 724, 866 Lendl, M., Anderson, D. R., Collier-Cameron, A., et al. 22, A&A, 544, A72 Lendl, M., Triaud, A. H. M. J., Anderson, D. R., et al. 24, ArXiv e-prints For comparison, the equivalent value for WASP-4 is just 7.6 Myr. 7
Lucy, L. B. & Sweeney, M. A. 97, AJ, 76, 544 Mandel, K. & Agol, E. 22, ApJ, 58, L7 Maxted, P. F. L., Anderson, D. R., Collier Cameron, A., et al. 2, PASP, 23, 547 Munari, U. & Zwitter, T. 997, A&A, 38, 269 Pollacco, D., Skillen, I., Collier Cameron, A., et al. 28, MNRAS, 385, 576 Pollacco, D. L., Skillen, I., Cameron, A. C., et al. 26, PASP, 8, 47 Queloz, D., Henry, G. W., Sivan, J. P., et al. 2, A&A, 379, 279 Schwarz, G. 978, Ann. Stat., 6, 46 Sestito, P. & Randich, S. 25, A&A, 442, 65 Smith, A. M. S., Collier Cameron, A., Christian, D. J., et al. 26, MNRAS, 373, 5 Smith, A. M. S., Hebb, L., Collier Cameron, A., et al. 29, MNRAS, 398, 827 Southworth, J. 2, MNRAS, 47, 266 Steele, I. A., Bates, S. D., Gibson, N., et al. 28, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 74, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series Torres, G., Andersen, J., & Giménez, A. 2, A&A Rev., 8, 67 Winn, J. N. 2, Exoplanet Transits and Occultations, ed. S. Seager, 55 77 Wright, J. T., Fakhouri, O., Marcy, G. W., et al. 2, PASP, 23, 42 Zacharias, N., Finch, C. T., Girard, T. M., et al. 22, VizieR Online Data Catalog, 322, A. M. S. Smith et al.: WASP-4b and WASP-6b 8