Source author record

Benjamin E. Nelson

Benjamin E. Nelson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.IM

Catalog footprint

What is connected

4works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2019arXiv

Quantifying the Bayesian Evidence for a Planet in Radial Velocity Data

We present results from a data challenge posed to the radial velocity (RV) community: namely, to quantify the Bayesian "evidence" for n={0,1,2,3} planets in a set of synthetically generated RV datasets containing a range of planet signals. Participating teams were provided the same likelihood function and set of priors to use in their analysis. They applied a variety of methods to estimate Z, the marginal likelihood for each n-planet model, including cross-validation, the Laplace approximation, importance sampling, and nested sampling. We found the dispersion in Z across different methods grew with increasing n-planet models: ~3 for 0-planets, ~10 for 1-planet, ~100-1000 for 2-planets, and >10,000 for 3-planets. Most internal estimates of uncertainty in Z for individual methods significantly underestimated the observed dispersion across all methods. Methods that adopted a Monte Carlo approach by comparing estimates from multiple runs yielded plausible uncertainties. Finally, two classes of numerical algorithms (those based on importance and nested samplers) arrived at similar conclusions regarding the ratio of Zs for n and (n+1)-planet models. One analytic method (the Laplace approximation) demonstrated comparable performance. We express both optimism and caution: we demonstrate that it is practical to perform rigorous Bayesian model comparison for <=3-planet models, yet robust planet discoveries require researchers to better understand the uncertainty in Z and its connections to model selection.

preprint2015arXiv

An Empirically Derived Three-Dimensional Laplace Resonance in the Gliese 876 Planetary System

We report constraints on the three-dimensional orbital architecture for all four planets known to orbit the nearby M dwarf Gliese 876 based solely on Doppler measurements and demanding long-term orbital stability. Our dataset incorporates publicly available radial velocities taken with the ELODIE and CORALIE spectrographs, HARPS, and Keck HIRES as well as previously unpublished HIRES velocities. We first quantitatively assess the validity of the planets thought to orbit GJ 876 by computing the Bayes factors for a variety of different coplanar models using an importance sampling algorithm. We find that a four-planet model is preferred over a three-planet model. Next, we apply a Newtonian MCMC algorithm to perform a Bayesian analysis of the planet masses and orbits using an n-body model in three-dimensional space. Based on the radial velocities alone, we find that a 99% credible interval provides upper limits on the mutual inclinations for the three resonant planets ($Φ_{cb}<6.20^\circ$ for the "c" and "b" pair and $Φ_{be}<28.5^\circ$ for the "b" and "e" pair). Subsequent dynamical integrations of our posterior sample find that the GJ 876 planets must be roughly coplanar ($Φ_{cb}<2.60^\circ$ and $Φ_{be}<7.87^\circ$), suggesting the amount of planet-planet scattering in the system has been low. We investigate the distribution of the respective resonant arguments of each planet pair and find that at least one argument for each planet pair and the Laplace argument librate. The libration amplitudes in our three-dimensional orbital model supports the idea of the outer-three planets having undergone significant past disk migration.

preprint2014arXiv

The 55 Cancri Planetary System: Fully Self-Consistent N-body Constraints and a Dynamical Analysis

We present an updated study of the planets known to orbit 55 Cancri A using 1,418 high-precision radial velocity observations from four observatories (Lick, Keck, Hobby-Eberly Telescope, Harlan J. Smith Telescope) and transit time/durations for the inner-most planet, 55 Cancri "e" (Winn et al. 2011). We provide the first posterior sample for the masses and orbital parameters based on self-consistent n-body orbital solutions for the 55 Cancri planets, all of which are dynamically stable (for at least $10^8$ years). We apply a GPU version of Radial velocity Using N-body Differential evolution Markov Chain Monte Carlo (RUN DMC; B. Nelson et al. 2014) to perform a Bayesian analysis of the radial velocity and transit observations. Each of the planets in this remarkable system has unique characteristics. Our investigation of high-cadence radial velocities and priors based on space-based photometry yields an updated mass estimate for planet "e" ($8.09\pm0.26$ M$_\oplus$), which affects its density ($5.51\pm^{1.32}_{1.00}$ g cm$^{-3}$) and inferred bulk composition. Dynamical stability dictates that the orbital plane of planet "e" must be aligned to within $60^o$ of the orbital plane of the outer planets (which we assume to be coplanar). The mutual interactions between the planets "b" and "c" may develop an apsidal lock about $180^o$. We find 36-45% of all our model systems librate about the anti-aligned configuration with an amplitude of $51^o\pm^{6^o}_{10^o}$. Other cases showed short-term perturbations in the libration of $\varpi_b-\varpi_c$, circulation, and nodding, but we find the planets are not in a 3:1 mean-motion resonance. A revised orbital period and eccentricity for planet "d" pushes it further toward the closest known Jupiter analog in the exoplanet population.

preprint2013arXiv

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

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.