Fast Bayesian Inference for Exoplanet Discovery in Radial Velocity Data

Abstract

Inferring the number of planets $N$ in an exoplanetary system from radial velocity (RV) data is a challenging task. Recently, it has become clear that RV data can contain periodic signals due to stellar activity, which can be difficult to distinguish from planetary signals. However, even doing the inference under a given set of simplifying assumptions (e.g. no stellar activity) can be difficult. It is common for the posterior distribution for the planet parameters, such as orbital periods, to be multimodal and to have other awkward features. In addition, when $N$ is unknown, the marginal likelihood (or evidence) as a function of $N$ is required. Rather than doing separate runs with different trial values of $N$, we propose an alternative approach using a trans-dimensional Markov Chain Monte Carlo method within Nested Sampling. The posterior distribution for $N$ can be obtained with a single run. We apply the method to $\nu$ Oph and Gliese 581, finding moderate evidence for additional signals in $\nu$ Oph with periods of 36.11 $\pm$ 0.034 days, 75.58 $\pm$ 0.80 days, and 1709 $\pm$ 183 days; the posterior probability that at least one of these exists is 85%. The results also suggest Gliese 581 hosts many (7-15) "planets" (or other causes of other periodic signals), but only 4-6 have well determined periods. The analysis of both of these datasets shows phase transitions exist which are difficult to negotiate without Nested Sampling.

arXiv commentAccepted for publication in MNRAS. 9 pages, 12 figures. Code at http://www.github.com/eggplantbren/Exoplanet

References in corpus (24)

• The generalised Lomb-Scargle periodogram. A new formalism for the floating-mean and Keplerian periodograms
• No planet for HD 166435
• The Rossiter-McLaughlin effect and analytic radial velocity curves for transiting extrasolar planetary systems
• Parametrizing the exoplanet eccentricity distribution with the Beta distribution
• Stellar Activity Masquerading as Planets in the Habitable Zone of the M dwarf Gliese 581
• The HARPS search for southern extra-solar planets. X. A m sin i = 11 Mearth planet around the nearby spotted M dwarf GJ 674
• Exploring Exoplanet Populations with NASA's Kepler Mission
• A Sub-Earth-Mass Moon Orbiting a Gas Giant Primary or a High Velocity Planetary System in the Galactic Bulge
• Detecting extrasolar planets from stellar radial velocities using Bayesian evidence
• Bayesian Re-analysis of the Gliese 581 Exoplanet System
• Bayesian analysis of radial velocity data of GJ667C with correlated noise: evidence for only 2 planets
• Bayesian re-analysis of the radial velocities of Gliese 581. Evidence in favour of only four planetary companions
• Bayesian modeling of source confusion in LISA data
• Substellar Companions to Seven Evolved Intermediate-Mass Stars
• Recursive Pathways to Marginal Likelihood Estimation with Prior-Sensitivity Analysis
• Gaussian Process Modelling of Asteroseismic Data
• Probabilistic Catalogs for Crowded Stellar Fields
• A transdimensional Bayesian method to infer the star formation history of resolved stellar populations
• Bayesian Inference from Observations of Solar-Like Oscillations
• MOA-2013-BLG-220Lb: Massive Planetary Companion to Galactic-Disk Host
• Disentangling Overlapping Astronomical Sources using Spatial and Spectral Information
• Dynamical stability of the Gliese 581 exoplanetary system
• Vertical-likelihood Monte Carlo
• The Probabilities of Orbital-Companion Models for Stellar Radial Velocity Data

Cited by in corpus (24)

• Model selection and parameter inference in phylogenetics using Nested Sampling
• Nested sampling for physical scientists
• The Joker: A custom Monte Carlo sampler for binary-star and exoplanet radial velocity data
• Uncovering the planets and stellar activity of CoRoT-7 using only radial velocities
• Bayesian model selection without evidences: application to the dark energy equation-of-state
• Quantifying the Bayesian Evidence for a Planet in Radial Velocity Data
• kima: Exoplanet detection in radial velocities
• Comment on "Stellar activity masquerading as planets in the habitable zone of the M dwarf Gliese 581"
• Statistical methods for exoplanet detection with radial velocities
• Improving exoplanet detection capabilities with the false inclusion probability. Comparison with other detection criteria in the context of radial velocities
• BEBOP III. Observations and an independent mass measurement of Kepler-16 (AB) b -- the first circumbinary planet detected with radial velocities
• Improved Point Source Detection in Crowded Fields using Probabilistic Cataloging
• The stepping-stone sampling algorithm for calculating the evidence of gravitational wave models
• The HARPS search for southern extra-solar planets XL. Searching for Neptunes around metal-poor stars
• The need for a public forecast of stellar activity to optimise exoplanet radial velocity detections & transmission spectroscopy
• Scheduling strategies for the ESPRESSO follow-up of TESS targets
• DMPP-4: Candidate sub-Neptune mass planets orbiting a naked-eye star
• Computing Entropies With Nested Sampling
• DNest4: Diffusive Nested Sampling in C++ and Python
• PCAT-DE: Reconstructing point-like and diffuse signals in astronomical images using spatial and spectral information
• A continuous multiple hypothesis testing framework for optimal exoplanet detection
• Maximum Likelihood Estimation based on Random Subspace EDA: Application to Extrasolar Planet Detection
• Inferring planet occurrence rates from radial velocities
• BEBOP II: Sensitivity to sub-Saturn circumbinary planets using radial-velocities