Detecting extrasolar planets from stellar radial velocities using Bayesian evidence

Stellar radial velocity (RV) measurements have proven to be a very successful method for detecting extrasolar planets. Analysing RV data to determine the parameters of the extrasolar planets is a significant statistical challenge owing to the presence of multiple planets and various degeneracies between orbital parameters. Determining the number of planets favoured by the observed data is an even more difficult task. Bayesian model selection provides a mathematically rigorous solution to this problem by calculating marginal posterior probabilities of models with different number of planets, but the use of this method in extrasolar planetary searches has been hampered by the computational cost of the evaluating Bayesian evidence. None the less, Bayesian model selection has the potential to improve the interpretation of existing observational data and possibly detect yet undiscovered planets. We present a new and efficient Bayesian method for determining the number of extrasolar planets, as well as for inferring their orbital parameters, without having to calculate directly the Bayesian evidence for models containing a large number of planets. Instead, we work iteratively and at each iteration obtain a conservative lower limit on the odds ratio for the inclusion of an additional planet into the model. We apply this method to simulated data sets containing one and two planets and successfully recover the correct number of planets and reliable constraints on the orbital parameters. We also apply our method to RV measurements of HD 37124, 47 Ursae Majoris and HD 10180. For HD 37124, we confirm that the current data strongly favour a three-planet system. We find strong evidence for the presence of a fourth planet in 47 Ursae Majoris, but its orbital period is suspiciously close to 1 yr, casting doubt on its validity. For HD 10180 we find strong evidence for a six-planet system.

[1]  J. Valenti,et al.  Spectroscopic Properties of Cool Stars (SPOCS). I. 1040 F, G, and K Dwarfs from Keck, Lick, and AAT Planet Search Programs , 2005 .

[2]  Maciej Konacki,et al.  Orbital configurations and dynamical stability of multiplanet systems around sun-like stars HD 202206, 14 Herculis, HD 37124, and HD 108874 , 2006 .

[3]  David J. C. MacKay,et al.  Information Theory, Inference, and Learning Algorithms , 2004, IEEE Transactions on Information Theory.

[4]  R. Paul Butler,et al.  A Planet Orbiting 47 Ursae Majoris , 1996 .

[5]  R. Paul Butler,et al.  Seven New Keck Planets Orbiting G and K Dwarfs , 2003 .

[6]  C. G. Tinney,et al.  Catalog of nearby exoplanets , 2006 .

[7]  R. Paul Butler,et al.  A Second Planet Orbiting 47 Ursae Majoris , 2002 .

[8]  M. Hobson,et al.  A Bayesian approach to discrete object detection in astronomical data sets , 2002, astro-ph/0204457.

[9]  O. Lahav,et al.  exofit: orbital parameters of extrasolar planets from radial velocities , 2008, 0805.3532.

[10]  N. Lomb.  Least-squares frequency analysis of unequally spaced data , 1976 .

[11]  F. Feroz,et al.  Multimodal nested sampling: an efficient and robust alternative to Markov Chain Monte Carlo methods for astronomical data analyses , 2007, 0704.3704.

[12]  Jean-Luis Lizon,et al.  Setting New Standards with HARPS , 2003 .

[13]  F. Feroz,et al.  MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics , 2008, 0809.3437.

[14]  F. Feroz,et al.  Cluster detection in weak lensing surveys , 2008, 0810.0781.

[15]  F. Feroz,et al.  Bayesian modelling of clusters of galaxies from multifrequency‐pointed Sunyaev–Zel'dovich observations , 2008, 0811.1199.

[16]  U. von Toussaint,et al.  Bayesian inference and maximum entropy methods in science and engineering , 2004 .

[17]  Jeffrey D. Scargle,et al.  Statistical challenges in modern astronomy II , 1997 .

[18]  P. Gregory.  A Bayesian Analysis of Extrasolar Planet Data for HD 73526 , 2005 .

[19]  A. Liddle,et al.  Information criteria for astrophysical model selection , 2007, astro-ph/0701113.

[20]  J. Scargle.  Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data , 1982 .

[21]  P. C. Gregory,et al.  A Bayesian periodogram finds evidence for three planets in HD 11964 , 2007, 0709.0970.

[22]  F. Feroz,et al.  Use of the MULTINEST algorithm for gravitational wave data analysis , 2009, 0904.1544.

[23]  R. Paul Butler,et al.  Five New Multicomponent Planetary Systems , 2005 .

[24]  Debra A. Fischer,et al.  A Bayesian Periodogram Finds Evidence for Three Planets in 47 Ursae Majoris , 2010, 1003.5549.

[25]  Structure and Evolution of Nearby Stars with Planets. II. Physical Properties of ~1000 Cool Stars from the SPOCS Catalog , 2006, astro-ph/0607235.

[26]  M. Tuomi,et al.  Bayesian analysis of the radial velocities of HD 11506 reveals another planetary companion , 2009, 0902.2997.

[27]  P. C. Gregory.  A Bayesian Kepler periodogram detects a second planet in HD 208487 , 2006 .

[28]  R. Trotta.  Applications of Bayesian model selection to cosmological parameters , 2005, astro-ph/0504022.

[29]  Harold F. Levison,et al.  A SEARCH FOR MULTI-PLANET SYSTEMS USING THE HOBBY–EBERLY TELESCOPE , 2009, 0903.0652.