We present high-precision photometry of five consecutive transits of WASP-18, an extrasolar planetary system with one of the shortest orbital periods known. Through the use of telescope defocussing we achieve a photometric precision of 0.47–0.83 mmag per observation over complete transit events. The data are analysed using the jktebop code and three different sets of stellar evolutionary models. We find the mass and radius of the planet to be and (statistical and systematic errors) respectively. The systematic errors in the orbital separation and the stellar and planetary masses, arising from the use of theoretical predictions, are of a similar size to the statistical errors and set a limit on our understanding of the WASP-18 system. We point out that seven of the nine known massive transiting planets () have eccentric orbits, whereas significant orbital eccentricity has been detected for only four of the 46 less massive planets. This may indicate that there are two different populations of transiting planets, but could also be explained by observational biases. Further radial velocity observations of low-mass planets will make it possible to choose between these two scenarios.
Physical properties of WASP-18]Physical properties of the 0.94-day period transiting planetary system WASP-18
Astrophysics Group, Keele University, Newcastle-under-Lyme, ST5 5BG, UK
Armagh Observatory, College Hill, Armagh, BT61 9DG, Northern Ireland, UK
Niels Bohr Institute and Centre for Star and Planet Formation, University of Copenhagen, Juliane Maries vej 30, 2100 Copenhagen Ø, Denmark
SUPA, University of St Andrews, School of Physics & Astronomy, North Haugh, St Andrews, KY16 9SS, UK
Department of Physics & Astronomy, Aarhus University, Ny Munkegade, 8000 Aarhus C, Denmark
Astronomisches Rechen-Institut, Zentrum für Astronomie, Universität Heidelberg, Mönchhofstraße 12-14, 69120 Heidelberg, Germany
Dipartimento di Fisica “E. R. Caianiello”, Università di Salerno, Baronissi, Italy
Instituto Nazionale di Fisica Nucleare, Sezione di Napoli, Italy
Deutsches SOFIA Institut, NASA Ames Research Center, Mail Stop 211-3, Moffett Field, CA 94035, USA
Institut für Astrophysik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
Institut d’Astrophysique et de Géophysique, Université de Liège, 4000 Liège, Belgium
Department of Physics, Sharif University of Technology, Tehran, Iran
European Southern Observatory, Casilla 19001, Santiago 19, Chile
The recent discovery of the transiting extrasolar planetary system WASP-18 (Hellier et al., 2009, hereafter H09) lights the way towards understanding the tidal interactions between giant planets and their parent stars. WASP-18 b is one of the shortest-period ( d) and most massive () extrasolar planets known. These properties make it an unparallelled indicator of the tidal dissipation parameters (Goldreich & Soter, 1966) for the star and the planet (Jackson et al., 2009). A value similar to that observed for Solar system bodies (–; Peale 1999) would cause the orbital period of WASP-18 to decrease at a sufficient rate for the effect to be observable within ten years (H09).
In this work we present high-precision follow-up photometry of WASP-18, obtained using telescope-defocussing techniques (Southworth et al., 2009a, b) which give a scatter of only 0.47 to 0.83 mmag per observation. These are analysed to yield improved physical properties of the WASP-18 system, with careful attention paid to statistical and systematic errors. The quality of the light curve is a critical factor in measurement of the physical properties of transiting planets (Southworth, 2009). In the case of WASP-18 the systematic errors arising from the use of theoretical stellar models are also important, and are a limiting factor in the understanding of this system.
2. Observations and data reduction
|Date||Start time (UT)||End time (UT)||Exposure time (s)||Filter||Airmass||Moon||Distance ()||Scatter (mmag)|
|2009 09 07||06:41||10:05||93||80.0||1.05 1.27||0.860||61.9||0.83|
|2009 09 08||05:00||10:04||140||80.0||1.15 1.04 1.28||0.784||67.1||0.68|
|2009 09 09||03:55||10:00||169||80.0||1.31 1.04 1.27||0.696||73.3||0.51|
|2009 09 10||03:25||08:29||141||80.0||1.41 1.04 1.10||0.597||80.2||0.56|
|2009 09 11||04:07||07:12||86||80.0||1.25 1.04||0.492||87.6||0.47|
We observed five consecutive transits of WASP-18 on the nights of 2009 September 7–11, using the 1.54 m Danish Telescope at ESO La Silla with the focal-reducing imager DFOSC. The plate scale of this setup is 0.39 pixel. The CCD was windowed down to an area of 128 in order to decrease the readout time to 51 s. This window was chosen to include WASP-18 (HD 10069, spectral type F6 V, , ) and a good comparison star (HD 10179, spectral type F3 V, , ).
A Johnson filter and exposure times of 80 s were used (instead of our usual Cousins and 120 s), to obtain lower count rates from the target and comparison stars. We defocussed the telescope to a point spread function (PSF) diameter of 90 pixels (35), to limit the peak counts to 45 000 per pixel for WASP-18 and also lower the flat-fielding noise. This focus setting was used for all observations and the pointing of the telescope was maintained using autoguiding. An observing log is given in Table 1.
We took several images with the telescope properly focussed in order to check that there were no nearby stars contaminating the PSF of WASP-18. There are five detectable objects nearby, but all are at least 100 pixels distant. The brightest is 6.75 mag fainter than WASP-18 and the other four are more than 9.5 mag fainter. We conclude that the PSF of WASP-18 is not contaminated by any star detectable with our equipment.
Data reduction was performed in the same way as in Southworth et al. (2009a, b). In short, we run a reduction pipeline written in idl
We calculated differential magnitudes for WASP-18 using an ensemble of four comparison stars, of which HD 10179 is by far the brightest. The light curve for each night was normalised to zero differential magnitude by fitting a straight line to the observations taken outside transit, whilst simultaneously optimising the weights of the four comparison stars. We have applied bias and flat-field corrections to the images, but find that this does not have a significant effect on the photometry. The individual light curves are shown in Fig. 1, and the full 629 datapoints are given in Table 2. The scatter in the final light curves varies from 0.47 to 0.83 mmag per point, and is higher for data taken when the moon was bright (Table 1).
|Midpoint of||Relative||Error in|
3. Light curve analysis
The light curves of WASP-18 were analysed with the jktebop
WASP-18 b has a slightly eccentric orbit, which should be taken into account as it has a small effect on the shape of the transit. However, for transit light curves the effect is in general too subtle to include it as a fitted parameter (Kipping, 2008). We could fix eccentricity, , and periastron longitude, , to the values obtained from the velocity variation of the parent star (H09), but this would neglect their measurement uncertainties. We have therefore modified jktebop to allow the inclusion of and as fitted parameters constrained by the known values and uncertainties. In practise we use and , as these two parameters are only weakly correlated with each other. In the current case the uncertainties in the and values have a very minor effect on our results.
We have incorporated the H09 photometry of WASP-18 in order to obtain the most precise ephemeris. This was done by including the reference transit epoch from H09 as an observed quantity, and then fitting for transit epoch and orbital period directly (see Southworth et al., 2007). We chose a new transit epoch which is close to the midpoint of our own observations, so is essentially based on just the data presented in this work. The final eclipse ephemeris is
where is the transit midpoint, is the number of orbital cycles after the reference epoch, and quantities in parentheses denote the uncertainty in the final digit of the preceding number.
|Linear LD law||Quadratic LD law||Square-root LD law||Logarithmic LD law||Cubic LD law|
|All LD coefficients fixed|
|0.60 fixed||0.40 fixed||0.20 fixed||0.70 fixed||0.40 fixed|
|0.30 fixed||0.60 fixed||0.25 fixed||0.15 fixed|
|Fitting for the linear LD coefficient and fixing the nonlinear LD coefficient|
|0.30 fixed||0.60 fixed||0.25 fixed||0.15 fixed|
|Fitting for the linear LD coefficient and perturbing the nonlinear LD coefficient|
|0.30 perturbed||0.60 perturbed||0.25 perturbed||0.15 perturbed|
|Fitting for both LD coefficients|
The limb darkening (LD) of WASP 18 A was accounted for using five different parametric laws (see Southworth, 2008). Theoretical LD coefficients were obtained by bilinear interpolation, to the known effective temperature () and surface gravity of the parent star, in the tables of Van Hamme (1993), Claret (2000, 2004a) and Claret & Hauschildt (2003). We obtained solutions with the LD coefficients fixed at the theoretical values, with the linear coefficient fitted for and the nonlinear coefficient fixed (and optionally perturbed by 0.05 on a flat distribution in the error analyses), and with both LD coefficients included as fitted parameters. The full set of solutions is given in Table 3.
The uncertainties of the light curve parameters were assessed using Monte Carlo simulations (Southworth et al., 2004c, 2005b). The importance of red noise was checked using a residual-permutation approach, and found to be minor. The solutions with LD coefficients fixed to theoretical values are poorer than those where one or both LD coefficients are fitted parameters. We adopt the mean of the solutions with non-linear LD and both LD coefficients fitted, as these are the most internally consistent. The final uncertainties come from Monte Carlo solutions but include contributions from the residual-permutation analyses and the (minor) variation between the solutions with different LD laws. We find the fractional radius
4. The physical properties of WASP-18
|Cambridge models||Y models||Claret models||Final result (this work)||Hellier et al. (2009)|
|Planet surface gravity||( m s)||191||17||191||17||191||17||191||17|
|Stellar age||(Gyr)||0.0 – 0.6||0.0 – 2.1||0.0 – 2.0||0.5 – 1.5|
The physical properties of a transiting planetary system cannot in general be calculated purely from observed quantities. The most common way to overcome this difficulty is to impose predictions from theoretical stellar evolutionary models onto the parent star. We have used tabulated predictions from three sources: Claret (Claret, 2004b, 2005, 2006, 2007), Y (Demarque et al., 2004) and Cambridge (Pols et al., 1998; Eldridge & Tout, 2004). This allows the assessment of the systematic errors caused by using stellar theory.
We began with the parameters measured from the light curve and the observed velocity amplitude of the parent star, m s (H09). These were augmented by an estimate of the velocity amplitude of the planet, , to calculate preliminary physical properties of the system. We then interpolated within one of the grids of theoretical predictions to find the expected radius and of the star for the preliminary mass and the measured metal abundance (; H09). was then iteratively refined to minimise the difference between the model-predicted radius and , and the calculated radius and measured ( K; H09). This was done for a range of ages for the star, and the best overall fit retained as the optimal solution. Finally, the above process was repeated whilst varying every input parameter by its uncertainty to build up a complete error budget for each output parameter (Southworth et al., 2005a). A detailed description of this process can be found in Southworth (2009).
Table 4 shows the results of these analyses. The physical properties calculated using the Claret and Y sets of stellar models are in excellent agreement, but those from the Cambridge models are slightly discrepant. This causes systematic errors of 4% in the stellar mass and 2% in the planetary mass, both of similar size to the corresponding statistical errors. The quality of our results is therefore limited by our theoretical understanding of the parent star. Our final results are in good agreement with those of H09 (Table 4), but incorporate a more comprehensive set of uncertainties.
These results give the equilibrium temperature of the planet to be one of the highest for the known planets:
where is the Bond albedo and is the heat redistribution factor. This equilibrium temperature, and the closeness to its parent star, make WASP-18 b a good target for the detection of thermal emission and reflected light.
We have presented high-quality observations of five consecutive transits by the newly-discovered planet WASP-18 b, which has one of the shortest orbital periods of all known transiting extrasolar planetary systems (TEPs). Our defocussed-photometry approach yielded scatters of between 0.47 and 0.83 mmag per point in the final light curves. These data were analysed using the jktebop code, which was modified to include the spectroscopically derived orbital eccentricity in a statistically correct way. The light curve parameters were then combined with the predictions of theoretical stellar evolutionary models to determine the physical properties of the planet and its host star.
A significant source of uncertainty in our results stems from the use of theoretical models to constrain the physical properties of the star. Further uncertainty comes from observed and , for which improved values are warranted. However, the systematic error from the use of stellar theory is an important uncertainty in the masses of the star and planet. This is due to our fundamentally incomplete understanding of the structure and evolution of low-mass stars. As with many other transiting systems (e.g. WASP-4; Southworth et al. 2009b), our understanding of the planet is limited by our lack of understanding of the parent star.
We confirm and refine the physical properties of WASP 18 found by H09. WASP-18 b is a very massive planet in an extremely short-period and eccentric orbit, which is a clear indicator that the tidal effects in planetary systems are weaker than expected (see H09). Long-term follow-up studies of WASP-18 will add progressively stricter constraints on the orbital decay of the planet and thus the strength of these tidal effects.
We now split the full sample of known (i.e. published) TEPS into two classes according to planetary mass. The mass distribution of transiting planets shows a dearth of objects with masses in the interval 2.0–3.1. There are nine planets more massive than this and 46 less massive. Seven of the nine high-mass TEPs have eccentric orbits (HAT-P-2, Bakos et al. 2007; HD 17156, Barbieri et al. 2007; HD 80606, Laughlin et al. 2009; WASP-10, Christian et al. 2009; WASP-14, Joshi et al. 2009; WASP-18; XO-3, Johns-Krull et al. 2008), and the existing radial velocity observations of the remaining two cannot rule out an eccentricity of or lower (CoRoT-Exo-2, Alonso et al. 2008; OGLE-TR-L9, Snellen et al. 2009). By comparison, only four of the 46 low-mass TEPs have a significant (Lucy & Sweeney, 1971) orbital eccentricity measurement.
These numbers imply that the more massive TEPs are a different population to the less massive ones; Fisher’s exact test (Fisher, 1922) returns a probability lower than of the null hypothesis (although this does not account for our freedom to choose the dividing line between the two classes). This indicates that the two types of TEPs have a different internal structure, formation mechanism, or evolution, a suggestion which is supported by observations of misalignment between the spin and orbital axes of TEPs (Johnson et al., 2009).
There is, however, a bias at work here. The more massive TEPs cause a larger radial velocity signal in their parent star (), so a given set of radial velocity measurements can detect smaller eccentricities (see also Shen & Turner, 2008). The eccentricity of the WASP-18 system is in fact below the detection limit of existing observations of most TEPs. We therefore advocate the acquisition of additional velocity data for the known low-mass TEPs, in order to equalise the eccentricity detection limits between the two classes of TEPs. These observations would allow acceptance or rejection of the hypothesis that more massive TEPs represent a fundamentally different planet population to their lower-mass brethren.
- Royal Society University Research Fellow
- Based on data collected by MiNDSTEp with the Danish 1.54 m telescope at the ESO La Silla Observatory
- The acronym idl stands for Interactive Data Language and is a trademark of ITT Visual Information Solutions. For further details see http://www.ittvis.com/ProductServices/IDL.aspx.
- jktebop is written in fortran77 and the source code is available at http://www.astro.keele.ac.uk/jkt/
- Fractional radius is the radius of a component of a binary system expressed as a fraction of the orbital semimajor axis. The utility of this quantity is that it is measureable from light curve data alone.
- Alonso, R., et al., 2008, A&A, 482, L21
- Bakos, G. Á., et al., 2007, ApJ, 670, 826
- Barbieri, M., et al., 2007, A&A, 476, L13
- Christian, D. J., et al., 2009, MNRAS, 392, 1585
- Claret, A., 2000, A&A, 363, 1081
- Claret, A., 2004a, A&A, 428, 1001
- Claret, A., 2004b, A&A, 424, 919
- Claret, A., 2005, A&A, 440, 647
- Claret, A., 2006, A&A, 453, 769
- Claret, A., 2007, A&A, 467, 1389
- Claret, A., Hauschildt, P. H., 2003, A&A, 412, 241
- Demarque, P., Woo, J.-H., Kim, Y.-C., Yi, S. K., 2004, ApJS, 155, 667
- Eldridge, J. J., Tout, C. A., 2004, MNRAS, 353, 87
- Fisher, R. A., 1922, Journal of the Royal Statistical Society, 85, 87
- Goldreich, P., Soter, S., 1966, Icarus, 5, 375
- Hellier, C., et al., 2009, Nature, 460, 1098
- Jackson, B., Barnes, R., Greenberg, R., 2009, ApJ, 698, 1357
- Johns-Krull, C. M., et al., 2008, ApJ, 677, 657
- Johnson, J. A., Winn, J. N., Albrecht, S., Howard, A. W., Marcy, G. W., Gazak, J. Z., 2009, PASP, 121, 1104
- Joshi, Y. C., et al., 2009, MNRAS, 392, 1532
- Kipping, D. M., 2008, MNRAS, 389, 1383
- Laughlin, G., Deming, D., Langton, J., Kasen, D., Vogt, S., Butler, P., Rivera, E., Meschiari, S., 2009, Nature, 457, 562
- Lucy, L. B., Sweeney, M. A., 1971, AJ, 76, 544
- Peale, S. J., 1999, ARA&A, 37, 533
- Pols, O. R., Schroder, K.-P., Hurley, J. R., Tout, C. A., Eggleton, P. P., 1998, MNRAS, 298, 525
- Shen, Y., Turner, E. L., 2008, ApJ, 685, 553
- Snellen, I. A. G., et al., 2009, A&A, 497, 545
- Southworth, J., 2008, MNRAS, 386, 1644
- Southworth, J., 2009, MNRAS, 394, 272
- Southworth, J., Maxted, P. F. L., Smalley, B., 2004a, MNRAS, 349, 547
- Southworth, J., Maxted, P. F. L., Smalley, B., 2004b, MNRAS, 351, 1277
- Southworth, J., Zucker, S., Maxted, P. F. L., Smalley, B., 2004c, MNRAS, 355, 986
- Southworth, J., Maxted, P. F. L., Smalley, B., 2005a, A&A, 429, 645
- Southworth, J., Smalley, B., Maxted, P. F. L., Claret, A., Etzel, P. B., 2005b, MNRAS, 363, 529
- Southworth, J., Bruntt, H., Buzasi, D. L., 2007, A&A, 467, 1215
- Southworth, J., et al., 2009a, MNRAS, 396, 1023
- Southworth, J., et al., 2009b, MNRAS, in press (preprint arXiv:0907:3356)
- Stetson, P. B., 1987, PASP, 99, 191
- Van Hamme, W., 1993, AJ, 106, 2096