An Affine-Invariant Sampler for Exoplanet Fitting and Discovery in Radial Velocity Data

Abstract

Markov Chain Monte Carlo (MCMC) proves to be powerful for Bayesian inference and in particular for exoplanet radial velocity fitting because MCMC provides more statistical information and makes better use of data than common approaches like chi-square fitting. However, the non-linear density functions encountered in these problems can make MCMC time-consuming. In this paper, we apply an ensemble sampler respecting affine invariance to orbital parameter extraction from radial velocity data. This new sampler has only one free parameter, and it does not require much tuning for good performance, which is important for automatization. The autocorrelation time of this sampler is approximately the same for all parameters and far smaller than Metropolis-Hastings, which means it requires many fewer function calls to produce the same number of independent samples. The affine-invariant sampler speeds up MCMC by hundreds of times compared with Metropolis-Hastings in the same computing situation. This novel sampler would be ideal for projects involving large datasets such as statistical investigations of planet distribution. The biggest obstacle to ensemble samplers is the existence of multiple local optima; we present a clustering technique to deal with local optima by clustering based on the likelihood of the walkers in the ensemble. We demonstrate the effectiveness of the sampler on real radial velocity data.

24 pages, 7 figures, accepted to ApJ

References in corpus (10)

• MCMC using Hamiltonian dynamics
• The Occurrence and Mass Distribution of Close-in Super-Earths, Neptunes, and Jupiters
• Inferring the eccentricity distribution
• On the Eccentricity Distribution of Exoplanets from Radial Velocity Surveys
• Observational biases in determining extrasolar planet eccentricities in single-planet systems
• Precise radial velocities of giant stars: I. Stable stars
• Precise radial velocities of giant stars. IV. A correlation between surface gravity and radial velocity variation and a statistical investigation of companion properties
• The Posterior Distribution of sin(i) for Exoplanets with M_T sin(i) Determined from Radial Velocity Data
• The Lyot Project Direct Imaging Survey of Substellar Companions: Statistical Analysis and Information from Nondetections
• Planets around Giant Stars

Cited by in corpus (31)

• emcee: The MCMC Hammer
• PS1-10jh: The Disruption of a Main-Sequence Star of Near-Solar Composition
• Bayesian parameter estimation for effective field theories
• Reconnaissance of the HR 8799 Exosolar System II: Astrometry and Orbital Motion
• Pyaneti: a fast and powerful software suite for multi-planet radial velocity and transit fitting
• RUN DMC: An efficient, parallel code for analyzing Radial Velocity Observations using N-body Integrations and Differential Evolution Markov chain Monte Carlo
• The surprising accuracy of isothermal Jeans modelling of self-interacting dark matter density profiles
• Asteroseismic inference on the spin-orbit misalignment and stellar parameters of HAT-P-7
• The Black Hole Mass Function Derived from Local Spiral Galaxies
• The discovery of WASP-151b, WASP-153b, WASP-156b: Insights on giant planet migration and the upper boundary of the Neptunian desert
• Bayesian parameter constraints for neutron star masses and radii using X-ray timing observations of accretion-powered millisecond pulsars
• Probing Final Stages of Stellar Evolution with X-Ray Observations of SN 2013ej
• The density of mid-sized Kuiper belt objects from ALMA thermal observations
• Constraining the Orbit and Mass of epsilon Eridani b with Radial Velocities, Hipparcos IAD-Gaia DR2 Astrometry, and Multi-epoch Vortex Coronagraphy Upper Limits
• Stellar Disc Truncations and Extended Haloes in Face-on Spiral Galaxies
• Reanalysis of radial velocity data from the resonant planetary system HD128311
• WD 1202-024: The Shortest-Period Pre-Cataclysmic Variable
• The HARPS search for southern extra-solar planets. XXXIX. HD175607 b, the most metal-poor G dwarf with an orbiting sub-Neptune
• Inferring heat recirculation and albedo for exoplanetary atmospheres: Comparing optical phase curves and secondary eclipse data
• Searching for comets on the World Wide Web: The orbit of 17P/Holmes from the behavior of photographers
• Modeling the Formation of the Family of the Dwarf Planet Haumea
• The Size and Shape of the Milky Way Disk and Halo from M-type Brown Dwarfs in the BoRG Survey
• Beyond Point Masses. II. Non-Keplerian Shape Effects are Detectable in Several TNO Binaries
• Marginal likelihoods of distances and extinctions to stars: computation and compact representation
• Chasing extreme planetary architectures: I- HD196885Ab, a super-Jupiter dancing with two stars?
• Reconstructing the Cosmic Expansion History with a Monotonicity Prior
• SNITCH: Seeking a simple, informative star formation history inference tool
• Closing the Evidence Gap: reddemcee, a Fast Adaptive Parallel Tempering Sampler
• Bayesian inference in single-line spectroscopic binaries with a visual orbit
• Optimal observational scheduling framework for binary and multiple stellar systems
• Dynamical Interactions in the Planetary System GJ4276