RUN DMC: An efficient, parallel code for analyzing Radial Velocity Observations using N-body Integrations and Differential Evolution Markov chain Monte Carlo

Abstract

In the 20+ years of Doppler observations of stars, scientists have uncovered a diverse population of extrasolar multi-planet systems. A common technique for characterizing the orbital elements of these planets is Markov chain Monte Carlo (MCMC), using a Keplerian model with random walk proposals and paired with the Metropolis-Hastings algorithm. For approximately a couple of dozen planetary systems with Doppler observations, there are strong planet-planet interactions due to the system being in or near a mean-motion resonance (MMR). An N-body model is often required to accurately describe these systems. Further computational difficulties arise from exploring a high-dimensional parameter space ($\sim$7 x number of planets) that can have complex parameter correlations. To surmount these challenges, we introduce a differential evolution MCMC (DEMCMC) applied to radial velocity data while incorporating self-consistent N-body integrations. Our Radial velocity Using N-body DEMCMC (RUN DMC) algorithm improves upon the random walk proposal distribution of the traditional MCMC by using an ensemble of Markov chains to adaptively improve the proposal distribution. We describe the methodology behind the algorithm, along with results of tests for accuracy and performance. We find that most algorithm parameters have a modest effect on the rate of convergence. However, the size of the ensemble can have a strong effect on performance. We show that the optimal choice depends on the number of planets in a system, as well as the computer architecture used and the resulting extent of parallelization. While the exact choices of optimal algorithm parameters will inevitably vary due to the details of individual planetary systems, we offer recommendations for choosing the DEMCMC algorithm's algorithmic parameters that result in excellent performance for a wide variety of planetary systems.

14 pages, 7 figures, accepted for publication in ApJS

References in corpus (4)

• EXOFAST: A fast exoplanetary fitting suite in IDL
• Five Planets Orbiting 55 Cancri
• ExoFit: Orbital Parameters of Extra-solar Planets from Radial Velocities
• A Third Planet Orbiting HIP 14810

Cited by in corpus (52)

• First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole
• Data analysis recipes: Using Markov Chain Monte Carlo
• Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy
• Inferring probabilistic stellar rotation periods using Gaussian processes
• Secure TTV Mass Measurements: Ten Kepler Exoplanets between 3 and 8 Earth Masses with Diverse Densities and Incident Fluxes
• The mass of the Mars-sized exoplanet Kepler-138 b from transit timing
• Revisiting Proxima with ESPRESSO
• AT2017gfo: Bayesian inference and model selection of multi-component kilonovae and constraints on the neutron star equation of state
• Separating planetary reflex Doppler shifts from stellar variability in the wavelength domain
• The 55 Cancri Planetary System: Fully Self-Consistent N-body Constraints and a Dynamical Analysis
• Discovery of a Third Transiting Planet in the Kepler-47 Circumbinary System
• An Empirically Derived Three-Dimensional Laplace Resonance in the Gliese 876 Planetary System
• Detecting Ocean Glint on Exoplanets Using Multiphase Mapping
• Improving the Accuracy of Planet Occurrence Rates from Kepler using Approximate Bayesian Computation
• Diffuser-Assisted Infrared Transit Photometry for Four Dynamically Interacting \textit{Kepler} Systems
• The California Planet Survey IV: A Planet Orbiting the Giant Star HD 145934 and Updates to Seven Systems with Long-Period Planets
• Newly-Discovered Planets Orbiting HD~5319, HD~11506, HD~75784 and HD~10442 from the N2K Consortium
• Evidence for the volatile-rich composition of a 1.5-$R_{\oplus }$ planet
• New Constraints on Gliese 876 -- Exemplar of Mean-Motion Resonance
• $\mathtt{bajes}$: Bayesian inference of multimessenger astrophysical data, methods and application to gravitational-waves
• A new dynamical modeling of the WASP-47 system with CHEOPS observations
• Identifying inflated super-Earths and photo-evaporated cores
• Starspot mapping with adaptive parallel tempering I: Implementation of computational code
• The conjectured S-type retrograde planet in nu Octantis: more evidence including four years of iodine-cell radial velocities
• HD 20329b: An ultra-short-period planet around a solar-type star found by TESS
• TTV-determined Masses for Warm Jupiters and their Close Planetary Companions
• The central black hole in the dwarf spheroidal galaxy Leo I Not supermassive, at most an intermediate-mass candidate
• Dark Matter in Fractional Gravity I: Astrophysical Tests on Galactic Scales
• Allesfitter: Flexible Star and Exoplanet Inference From Photometry and Radial Velocity
• Orbital fitting of imaged planetary companions with high eccentricities and unbound orbits -- Application to Fomalhaut b and PZ Telescopii B
• Chasing extreme planetary architectures: I- HD196885Ab, a super-Jupiter dancing with two stars?
• Uniform Modelling of the Stellar Density of Thirteen Tidal Streams within the Galactic Halo
• Gravitationally Bound Gas Determines Star Formation in the Galaxy
• Validating Stellar Abundance Measurements from Multi-Resolution Spectroscopy
• Radii, masses, and transit-timing variations of the three-planet system orbiting the naked-eye star TOI-396
• UCLCHEMCMC: A MCMC Inference tool for Physical Parameters of Molecular Clouds
• The red giant branch tip in the SDSS, PS1, JWST, NGRST and Euclid photometric systems. Calibration in optical passbands using Gaia DR3 synthetic photometry
• The K2-24 planetary system revisited by CHEOPS
• Asteroseismology with PBjam 2.0: measuring dipole mode frequencies in coupling regimes from main sequence to low-luminosity red giant stars
• Angular Correlation Function from sample covariance with BOSS and eBOSS LRG
• Novel Cosmological Joint Constraints in Multidimensional Observables Space with Redshift-free Inferences
• Tidally Heated Sub-Neptunes, Refined Planetary Compositions, and Confirmation of a Third Planet in the TOI-1266 System
• jaxspec : a fast and robust Python library for X-ray spectral fitting
• Modified Hermite Integrators of Arbitrary Order
• Transit Timing Variations in HIP 41378: CHEOPS and TESS confirm a non-transiting sixth planet in the system
• Searching for rotation in X-COP galaxy clusters
• Mass Derivation of planets K2-21b and K2-21c from Transit Timing Variations
• The multi-planet system TOI-5624: Four transiting sub-Neptunes with an outer companion revealed by transit-timing variations
• Multi-component, axisymmetric dynamical models of dSphs based on distribution functions: inferences on dark matter and intermediate-mass black holes in Draco and Ursa Minor
• Pulsar electrodynamics inferred from frequency-dependent circular polarization diversity
• Optimal observational scheduling framework for binary and multiple stellar systems
• Bayesian inference in single-line spectroscopic binaries with a visual orbit