Orbits for the Impatient: A Bayesian Rejection-sampling Method for Quickly Fitting the Orbits of Long-period Exoplanets

We describe a Bayesian rejection sampling algorithm designed to efficiently compute posterior distributions of orbital elements for data covering short fractions of long-period exoplanet orbits. Our implementation of this method, Orbits for the Impatient (OFTI), converges up to several orders of magnitude faster than two implementations of MCMC in this regime. We illustrate the efficiency of our approach by showing that OFTI calculates accurate posteriors for all existing astrometry of the exoplanet 51 Eri b up to 100 times faster than a Metropolis-Hastings MCMC. We demonstrate the accuracy of OFTI by comparing our results for several orbiting systems with those of various MCMC implementations, finding the output posteriors to be identical within shot noise. We also describe how our algorithm was used to successfully predict the location of 51 Eri b six months in the future based on less than three months of astrometry. Finally, we apply OFTI to ten long-period exoplanets and brown dwarfs, all but one of which have been monitored over less than 3% of their orbits, producing fits to their orbits from astrometric records in the literature.

[1]  Daniel C. Fabrycky,et al.  Submitted to ApJ , 1999 .

[2]  A.-M. Lagrange,et al.  Constraints on planets around β Pic with Harps radial velocity data , 2012, 1202.2579.

[3]  Bertrand Mennesson,et al.  DISCOVERY OF A LOW-MASS COMPANION AROUND HR 3549 , 2015, 1509.05689.

[4]  Howard Isaacson,et al.  ORBITAL ARCHITECTURES OF PLANET-HOSTING BINARIES. I. FORMING FIVE SMALL PLANETS IN THE TRUNCATED DISK OF KEPLER-444A , 2015, 1512.03428.

[5]  Gautam Vasisht,et al.  THE κ ANDROMEDAE SYSTEM: NEW CONSTRAINTS ON THE COMPANION MASS, SYSTEM AGE, AND FURTHER MULTIPLICITY , 2013, 1309.3372.

[6]  K. Y. L. Su,et al.  Steady State Evolution of Debris Disks around A Stars , 2007 .

[7]  Daniel Foreman-Mackey,et al.  The Joker: A Custom Monte Carlo Sampler for Binary-star and Exoplanet Radial Velocity Data , 2016, 1610.07602.

[8]  Gautam Vasisht,et al.  THE ORBIT AND TRANSIT PROSPECTS FOR β PICTORIS b CONSTRAINED WITH ONE MILLIARCSECOND ASTROMETRY , 2016, 1607.05272.

[9]  Jungyeon Cho,et al.  IMBALANCED RELATIVISTIC FORCE-FREE MAGNETOHYDRODYNAMIC TURBULENCE , 2013, 1312.6128.

[10]  Virginie Faramaz,et al.  An independent determination of Fomalhaut b’s orbit and the dynamical effects on the outer dust belt , 2013, 1311.5035.

[11]  Andrew W. Serio,et al.  First light of the Gemini Planet Imager , 2014, Proceedings of the National Academy of Sciences.

[12]  Eric Mamajek,et al.  The Planetary Mass Companion 2MASS 1207–3932B: Temperature, Mass, and Evidence for an Edge-on Disk , 2006, astro-ph/0610550.

[13]  D. Mouillet,et al.  A giant planet candidate near a young brown dwarf - Direct VLT/NACO observations using IR wavefront sensing , 2004 .

[14]  Luca Casagrande,et al.  High precision abundances of the old solar twin HIP 102152: Insights on Li depletion from the oldest sun , 2013 .

[15]  Jason J. Wang,et al.  Discovery and spectroscopy of the young jovian planet 51 Eri b with the Gemini Planet Imager , 2015, Science.

[16]  Laird M. Close,et al.  THE GEMINI NICI PLANET-FINDING CAMPAIGN: DISCOVERY OF A SUBSTELLAR L DWARF COMPANION TO THE NEARBY YOUNG M DWARF CD−35 2722 , 2011, 1101.2893.

[17]  F. Fekel,et al.  The properties of close multiple stars , 1981 .

[18]  Inseok Song,et al.  HST NICMOS Imaging of the Planetary-mass Companion to the Young Brown Dwarf 2MASSW J1207334-393254 , 2006 .

[19]  David A. Golimowski,et al.  THE MASS–LUMINOSITY RELATION IN THE L/T TRANSITION: INDIVIDUAL DYNAMICAL MASSES FOR THE NEW J-BAND FLUX REVERSAL BINARY SDSS J105213.51+442255.7AB , 2015, 1503.06212.

[20]  Eric B. Ford,et al.  Improving the Efficiency of Markov Chain Monte Carlo for Analyzing the Orbits of Extrasolar Planets , 2005, astro-ph/0512634.

[21]  Gautam Vasisht,et al.  ASTROMETRIC CONFIRMATION AND PRELIMINARY ORBITAL PARAMETERS OF THE YOUNG EXOPLANET 51 ERIDANI b WITH THE GEMINI PLANET IMAGER , 2015, 1509.07514.

[22]  Dmitry Savransky,et al.  WFIRST-AFTA coronagraph science yield modeling with EXOSIMS , 2015, 1511.02869.

[23]  Serge Correia,et al.  DIRECT IMAGING AND SPECTROSCOPY OF A CANDIDATE COMPANION BELOW/NEAR THE DEUTERIUM-BURNING LIMIT IN THE YOUNG BINARY STAR SYSTEM, ROXs 42B , 2013, 1310.4825.

[24]  D. Mouillet,et al.  Giant Planet Companion to 2MASSW J1207334-393254 , 2005 .

[25]  A. Vigan,et al.  Spectroscopy across the brown dwarf/planetary mass boundary I. Near-infrared JHK spectra , , 2012, 1201.3921.

[26]  Jessica R. Lu,et al.  Measuring Distance and Properties of the Milky Way’s Central Supermassive Black Hole with Stellar Orbits , 2008, 0808.2870.

[27]  Laird M. Close,et al.  NEW SPATIALLY RESOLVED OBSERVATIONS OF THE T Cha TRANSITION DISK AND CONSTRAINTS ON THE PREVIOUSLY CLAIMED SUBSTELLAR COMPANION , 2015, 1501.01964.

[28]  C. A. Grady,et al.  DIRECT IMAGING OF A COLD JOVIAN EXOPLANET IN ORBIT AROUND THE SUN-LIKE STAR GJ 504 , 2013, 1307.2886.

[29]  Timothy D. Brandt,et al.  Characterization of the gaseous companion κ Andromedae b - New Keck and LBTI high-contrast observations , 2013, 1308.3859.

[30]  W. Press,et al.  The time delay of gravitational lens 0957+561. II: Analysis of radio data and combined optical-radio analysis , 1992 .

[31]  R. J. Terrile,et al.  A Candidate Substellar Companion to HR 7329 , 2000 .

[32]  G. Carraro,et al.  A brown dwarf companion to the intermediate-mass star HR 6037 , 2010, 1009.4650.

[33]  Laird M. Close,et al.  A UNIFORM ANALYSIS OF 118 STARS WITH HIGH-CONTRAST IMAGING: LONG-PERIOD EXTRASOLAR GIANT PLANETS ARE RARE AROUND SUN-LIKE STARS , 2009, 0909.4531.

[34]  Marc J. Kuchner,et al.  The Gemini NICI Planet-Finding Campaign: The Orbit of the Young Exoplanet β Pictoris b , 2014 .

[35]  O. Ilbert,et al.  REST-UV ABSORPTION LINES AS METALLICITY ESTIMATOR: THE METAL CONTENT OF STAR-FORMING GALAXIES AT z ∼ 5 , 2015, 1512.00018.

[36]  G. Chauvin,et al.  Mid-infrared imaging of brown dwarfs in binary systems , 2007, 0712.1887.

[37]  T. Fusco,et al.  First light of the VLT planet finder SPHERE II. The physical properties and the architecture of the young systems PZ Telescopii and HD 1160 revisited , 2015, 1511.04072.

[38]  Wolfgang Brandner,et al.  Infrared spectrum and proper motion of the brown dwarf companion of HR 7329 in Tucanae , 2001 .

[39]  Daniel Foreman-Mackey,et al.  emcee: The MCMC Hammer , 2012, 1202.3665.

[40]  Olivier Guyon,et al.  Science yield estimate with the Wide-Field Infrared Survey Telescope coronagraph , 2016 .

[41]  Jason J. Wang,et al.  CONSTRAINTS ON THE ARCHITECTURE OF THE HD 95086 PLANETARY SYSTEM WITH THE GEMINI PLANET IMAGER , 2016, 1604.05139.

[42]  Adam L. Kraus,et al.  SEARCHING FOR SCATTERERS: HIGH-CONTRAST IMAGING OF YOUNG STARS HOSTING WIDE-SEPARATION PLANETARY-MASS COMPANIONS , 2016, 1606.06744.

[43]  Julien H. Girard,et al.  Direct-imaging discovery of a 12–14 Jupiter-mass object orbiting a young binary system of very low-mass stars , 2013, 1303.4525.

[44]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[45]  Laird M. Close,et al.  THE GEMINI NICI PLANET-FINDING CAMPAIGN: THE FREQUENCY OF GIANT PLANETS AROUND YOUNG B AND A STARS , 2013, Proceedings of the International Astronomical Union.

[46]  Alycia J. Weinberger,et al.  IDENTIFYING THE YOUNG LOW-MASS STARS WITHIN 25 pc. II. DISTANCES, KINEMATICS, AND GROUP MEMBERSHIP , 2012, 1207.5074.

[47]  C. A. Grady,et al.  NEAR-INFRARED MULTI-BAND PHOTOMETRY OF THE SUBSTELLAR COMPANION GJ 758 B , 2010, 1011.5505.

[48]  Laird M. Close,et al.  THE GEMINI NICI PLANET-FINDING CAMPAIGN: DISCOVERY OF A MULTIPLE SYSTEM ORBITING THE YOUNG A STAR HD 1160 , 2012, 1202.2854.

[49]  Jena,et al.  Further deep imaging of HR 7329 A (η Tel A) and its brown dwarf companion B , 2011, 1106.1388.

[50]  Brendan P. Bowler,et al.  Imaging Extrasolar Giant Planets , 2016, 1605.02731.

[51]  Jason J. Wang,et al.  β PICTORIS’ INNER DISK IN POLARIZED LIGHT AND NEW ORBITAL PARAMETERS FOR β PICTORIS b , 2015, 1508.04787.

[52]  Eric B. Ford.  Quantifying the Uncertainty in the Orbits of Extrasolar Planets with Markov Chain Monte Carlo , 2004 .

[53]  D. Fantinel,et al.  Characterizing HR 3549 B using SPHERE , 2016, 1606.06654.

[54]  Adam L. Kraus,et al.  THREE WIDE PLANETARY-MASS COMPANIONS TO FW TAU, ROXs 12, AND ROXs 42B , 2013, 1311.7664.

[55]  Julien Lozi,et al.  SCExAO AND GPI Y JH BAND PHOTOMETRY AND INTEGRAL FIELD SPECTROSCOPY OF THE YOUNG BROWN DWARF COMPANION TO HD 1160 , 2016, 1610.05786.