The Efficiency of Geometric Samplers for Exoplanet Transit Timing Variation Models

Abstract

Transit timing variations (TTVs) are a valuable tool to determine the masses and orbits of transiting planets in multi-planet systems. TTVs can be readily modeled given knowledge of the interacting planets' orbital configurations and planet-star mass ratios, but such models are highly nonlinear and difficult to invert. Markov chain Monte Carlo (MCMC) methods are often used to explore the posterior distribution for model parameters, but, due to the high correlations between parameters, nonlinearity, and potential multi-modality in the posterior, many samplers perform very inefficiently. Therefore, we assess the performance of several MCMC samplers that use varying degrees of geometric information about the target distribution. We generate synthetic datasets from multiple models, including the TTVFaster model and a simple sinusoidal model, and test the efficiencies of various MCMC samplers. We find that sampling efficiency can be greatly improved for all models by sampling from a parameter space transformed using an estimate of the covariance and means of the target distribution. No one sampler performs the best for all datasets, but several samplers, such as Differential Evolution Monte Carlo and Geometric adaptive Monte Carlo, have consistently efficient performance. For datasets with near Gaussian posteriors, Hamiltonian Monte Carlo samplers with 2 or 3 leapfrog steps obtained the highest efficiencies. Based on differences in effective sample sizes per time, we show that the right choice of sampler can improve sampling efficiencies by several orders of magnitude.

References in corpus (11)

• emcee: The MCMC Hammer
• Architecture and Dynamics of Kepler's Candidate Multiple Transiting Planet Systems
• Validation of Kepler's Multiple Planet Candidates. III: Light Curve Analysis & Announcement of Hundreds of New Multi-planet Systems
• Quantifying the Uncertainty in the Orbits of Extrasolar Planets
• Kepler-36: A Pair of Planets with Neighboring Orbits and Dissimilar Densities
• Extracting Planet Mass and Eccentricity From TTV data
• Transit Timing Observations from Kepler. IX. Catalog of the Full Long-Cadence Data Set
• TTVFast: An efficient and accurate code for transit timing inversion problems
• Secure TTV Mass Measurements: Ten Kepler Exoplanets between 3 and 8 Earth Masses with Diverse Densities and Incident Fluxes
• Forward-Mode Automatic Differentiation in Julia
• Transit timing to first order in eccentricity

Cited by in corpus (4)

• A Simplified Photodynamical Model for Planetary Mass Determination in Low-Eccentricity Multi-Transiting Systems
• The Exoplanet Edge: Planets Don't Induce Observable TTVs with a Dominant TTV Period Faster than Half their Orbital Period
• A differentiable N-body code for transit timing and dynamical modelling -- II. Photodynamics
• Nauyaca: a new tool to determine planetary masses and orbital elements through transit timing analysis