orvara: An Efficient Code to Fit Orbits using Radial Velocity, Absolute, and/or Relative Astrometry

Abstract

We present an open-source Python package, Orbits from Radial Velocity, Absolute, and/or Relative Astrometry (orvara), to fit Keplerian orbits to any combination of radial velocity, relative astrometry, and absolute astrometry data from the Hipparcos-Gaia Catalog of Accelerations. By combining these three data types, one can measure precise masses and sometimes orbital parameters even when the observations cover a small fraction of an orbit. orvara achieves its computational performance with an eccentric anomaly solver five to ten times faster than commonly used approaches, low-level memory management to avoid python overheads, and by analytically marginalizing out parallax, barycenter proper motion, and the instrument-specific radial velocity zero points. Through its integration with the Hipparcos and Gaia intermediate astrometry package htof, orvara can properly account for the epoch astrometry measurements of Hipparcos and the measurement times and scan angles of individual Gaia epochs. We configure orvara with modifiable .ini configuration files tailored to any specific stellar or planetary system. We demonstrate orvara with a case study application to a recently discovered white dwarf/main sequence (WD/MS) system, HD 159062. By adding absolute astrometry to literature RV and relative astrometry data, our comprehensive MCMC analysis improves the precision of HD 159062B's mass by more than an order of magnitude to $0.6083_{-0.0073}^{+0.0083}{M}_{\odot }$. We also derive a low eccentricity and large semimajor axis, establishing HD 159062AB as a system that did not experience Roche lobe overflow.

arXiv comment24 pages, 5 figures, 5 tables. AJ accepted with minor changes. orvara is available at https://github.com/t-brandt/orvara

References in corpus (83)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Astropy: A Community Python Package for Astronomy
• emcee: The MCMC Hammer
• The NumPy array: a structure for efficient numerical computation
• Gaia Data Release 2. Summary of the contents and survey properties
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• The Gaia mission
• Validation of the new Hipparcos reduction
• Gaia Data Release 2: The astrometric solution
• A Survey of Stellar Families: Multiplicity of Solar-Type Stars
• Direct Imaging of Multiple Planets Orbiting the Star HR 8799
• Evolutionary models for cool brown dwarfs and extrasolar giant planets. The case of HD 20945
• Accurate masses and radii of normal stars: modern results and applications
• A Non-Gray Theory of Extrasolar Giant Planets and Brown Dwarfs
• Improved Age Estimation for Solar-Type Dwarfs Using Activity-Rotation Diagnostics
• Evolutionary models for very-low-mass stars and brown dwarfs with dusty atmospheres
• Type Ia Supernova Explosion Models
• Mind your Ps and Qs: the Interrelation between Period (P) and Mass-ratio (Q) Distributions of Binary Stars
• batman: BAsic Transit Model cAlculatioN in Python
• A giant planet imaged in the disk of the young star Beta Pictoris
• EXOFAST: A fast exoplanetary fitting suite in IDL
• On the rotational evolution of solar- and late-type stars, its magnetic origins, and the possibility of stellar gyrochronology
• The Evolution of L and T Dwarfs in Color-Magnitude Diagrams
• Discovery and spectroscopy of the young Jovian planet 51 Eri b with the Gemini Planet Imager
• RadVel: The Radial Velocity Modeling Toolkit
• The LCES HIRES/Keck Precision Radial Velocity Exoplanet Survey
• A $Gaia$ Data Release 2 catalogue of white dwarfs and a comparison with SDSS
• The California Planet Survey I. Four New Giant Exoplanets
• Imaging Extrasolar Giant Planets
• The Initial-Final Mass Relation: Direct Constraints at the Low Mass End
• An Assessment of Dynamical Mass Constraints on Pre-Main Sequence Evolutionary Tracks
• Dynamic temperature selection for parallel-tempering in Markov chain Monte Carlo simulations
• The Hipparcos-Gaia Catalog of Accelerations: Gaia EDR3 Edition
• Direct Imaging of a Cold Jovian Exoplanet in Orbit around the Sun-like Star GJ 504
• New cooling sequences for old white dwarfs
• The Hipparcos-Gaia Catalog of Accelerations
• exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series
• Discovery of a probable 4-5 Jupiter-mass exoplanet to HD 95086 by direct-imaging
• Efficient fitting of multiplanet Keplerian models to radial velocity and astrometry data
• Individual Dynamical Masses of Ultracool Dwarfs
• The Extrasolar Planet epsilon Eridani b - Orbit and Mass
• A Mass for the Extrasolar Planet Gl 876b Determined from Hubble Space Telescope Fine Guidance Sensor 3 Astrometry and High-Precision Radial Velocities
• The binarity of the local white dwarf population
• Chromospheric activity as age indicator
• Evidence of an Upper Bound on the Masses of Planets and its Implications for Giant Planet Formation
• orbitize!: A Comprehensive Orbit-fitting Software Package for the High-contrast Imaging Community
• Precise Dynamical Masses of Directly Imaged Companions from Relative Astrometry, Radial Velocities, and Hipparcos-Gaia DR2 Accelerations
• Abundances in the Local Region II: F, G and K Dwarfs and Subgiants
• Peering into the formation history of beta Pictoris b with VLTI/GRAVITY long baseline interferometry
• Grids of white dwarf evolutionary models with masses from M= 0.1 to 1.2 Ms
• New Praesepe white dwarfs and the initial mass-final mass relation
• Keck Laser Guide Star Adaptive Optics Monitoring of 2MASS J1534-2952AB: First Dynamical Mass Determination of a Binary T Dwarf
• The mass of the young planet Pictoris b through the astrometric motion of its host star
• γ^2 Velorum: Orbital Solution and Fundamental Parameter Determination with SUSI
• Evidence for a high mutual inclination between the cold Jupiter and transiting super Earth orbiting $\pi$ Men
• VLT/SPHERE robust astrometry of the HR8799 planets at milliarcsecond-level accuracy Orbital architecture analysis with PyAstrOFit
• Detection of the nearest Jupiter analog in radial velocity and astrometry data
• Characterizing and correcting the proper motion bias of the bright Gaia EDR3 sources
• A Model-Independent Mass and Moderate Eccentricity for $\beta$ Pic b
• Dynamical Mass of GJ 802B: a brown dwarf in a triple system
• A precise architecture characterization of the $\pi$ Men planetary system
• Derivation of the Mass Distribution of Extrasolar Planets with MAXLIMA - a Maximum Likelihood Algorithm
• Dynamical Mass Estimates for Incomplete Orbits: Young Multiple Stars in Taurus and Ophiuchus
• The TRENDS High-Contrast Imaging Survey. VI. Discovery of a Mass, Age, and Metallicity Benchmark Brown Dwarf
• The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet beta Pictoris b from Combined Direct Imaging and Astrometry
• Spectroscopic evolution of massive stars on the main sequence
• A significant mutual inclination between the planets within the $\pi$ Mensae system
• Galactic Tide and Local Stellar Perturbations on the Oort Cloud: Creation of Interstellar Comets
• Orbital Parameter Determination for Wide Stellar Binary Systems in the Age of Gaia
• ShaneAO: wide science spectrum adaptive optics system for the Lick Observatory
• Orbital inclination and mass of the exoplanet candidate Proxima c
• Improving dynamical mass constraints for intermediate-period substellar companions using Gaia DR2
• The Galactic Center: An Improved Astrometric Reference Frame for Stellar Orbits around the Supermassive Black Hole
• The McDonald Accelerating Stars Survey (MASS): White Dwarf Companions Accelerating the Sun-like Stars 12 Psc and HD 159062
• The UV Emission of Stars in LAMOST Survey I. Catalogs
• A dusty benchmark brown dwarf near the ice line of HD 72946
• Mid-Infrared High-Contrast Imaging of HD 114174 B : An Apparent Age Discrepancy in a "Sirius-Like" Binary System
• Kepler's Goat Herd: An Exact Solution to Kepler's Equation for Elliptical Orbits
• The Exoplanet Simple Orbit Fitting Toolbox (ExoSOFT): An Open-Source Tool for Efficient Fitting of Astrometric and Radial Velocity Data
• Dynamical masses for the triple system HD 28363 in the Hyades cluster
• Swimming with ShARCS: Comparison of On-sky Sensitivity With Model Predictions for ShaneAO on the Lick Observatory 3-meter Telescope
• Discovery of a White Dwarf Companion to HD 159062
• Dynamical masses for the Hyades binary 80 Tauri

Cited by in corpus (67)

• A non-interacting Galactic black hole candidate in a binary system with a main-sequence star
• AF Lep b: the lowest mass planet detected coupling astrometric and direct imaging data
• Direct imaging discovery of a super-Jovian around the young Sun-like star AF Leporis
• Astrometric Accelerations as Dynamical Beacons: A Giant Planet Imaged Inside the Debris Disk of the Young Star AF Lep
• Improved Dynamical Masses for Six Brown Dwarf Companions Using Hipparcos and Gaia EDR3
• A Clear View of a Cloudy Brown Dwarf Companion from High-Resolution Spectroscopy
• Direct Imaging Discovery and Dynamical Mass of a Substellar Companion Orbiting an Accelerating Hyades Sun-like Star with SCExAO/CHARIS
• A temperate super-Jupiter imaged with JWST in the mid-infrared
• Assessing the C/O Ratio Formation Diagnostic: A Potential Trend with Companion Mass
• On the Masses, Age, and Architecture of the VHS J1256-1257AB b System
• Astrometric Accelerations as Dynamical Beacons: Discovery and Characterization of HIP 21152 B, the First T-Dwarf Companion in the Hyades
• Dynamical Mass of the Young Substellar Companion HD 984 B
• The TESS-Keck Survey. VIII. Confirmation of a Transiting Giant Planet on an Eccentric 261 day Orbit with the Automated Planet Finder Telescope
• Speckle Interferometry at SOAR in 2021
• The Masses of a Sample of Radial-Velocity Exoplanets with Astrometric Measurements
• Barium and related stars, and their white-dwarf companions. III. The masses of the white dwarfs
• Limits on the Mass and Initial Entropy of 51 Eri b from Gaia EDR3 Astrometry
• Revised orbits of the two nearest Jupiters
• Precise Dynamical Masses of Epsilon Indi Ba and Bb: Evidence of Slowed Cooling at the L/T Transition
• Octofitter: Fast, Flexible, and Accurate Orbit Modelling to Detect Exoplanets
• Two Warm Neptunes transiting HIP 9618 revealed by TESS & Cheops
• The CORALIE survey for southern extrasolar planets XIX. Brown dwarfs and stellar companions unveiled by radial velocity and astrometry
• Polarised radio pulsations from a new T dwarf binary
• Updated characterization of long-period single companion by combining radial velocity, relative astrometry, and absolute astrometry
• Surveying Nearby Brown Dwarfs with HGCA: Direct Imaging Discovery of a Faint, High-Mass Brown Dwarf Orbiting HD 176535 A
• First Observations of the Brown Dwarf HD 19467 B with JWST
• Exploring the Brown Dwarf Desert with Precision Radial Velocities and Gaia DR3 Astrometric Orbits
• Precise dynamical masses of new directly imaged companions from combining relative astrometry, radial velocities, and Hipparcos-Gaia eDR3 accelerations
• The SOPHIE search for northern extrasolar planets-XIX. A system including a cold sub-Neptune potentially transiting a V = 6.5 star HD88986
• Direct-imaging Discovery of a Substellar Companion Orbiting the Accelerating Variable Star, HIP 39017
• Dynamical Mass of the Young Brown Dwarf Companion PZ Tel B
• Atmospheric characterization of the super-Jupiter HIP 99770 b with KPIC
• Radial velocity homogeneous analysis of M dwarfs observed with HARPS I. Exoplanet detection and candidates
• Multi techniques approach to identify and/or constrain radial velocity sub-stellar companions
• Combining Hipparcos and Gaia data for the study of binaries: the BINARYS tool
• The Gliese 86 Binary System: A Warm Jupiter Formed in a Disk Truncated at $\approx$2 AU
• SCExAO and Keck Direct Imaging Discovery of a Low-Mass Companion Around the Accelerating F5 Star HIP 5319
• Significant mutual inclinations between the stellar spin and the orbits of both planets in the HAT-P-11 system
• Use the 4S (Signal-Safe Speckle Subtraction): Explainable Machine Learning reveals the Giant Exoplanet AF Lep b in High-Contrast Imaging Data from 2011
• The McDonald Accelerating Stars Survey (MASS): Architecture of the Ancient Five-Planet Host System Kepler-444
• HD 83443c: A highly eccentric giant planet on a 22-year orbit
• The HARPS search for southern extra-solar planets. XLVII. Five Jupiter-mass planets in long-period orbits, one highly irradiated Neptune, one brown dwarf, and five stellar binaries
• SCExAO/CHARIS Direct Imaging of A Low-Mass Companion At A Saturn-Like Separation from an Accelerating Young A7 Star
• A dynamic view of V Hydrae. Monitoring of a spectroscopic-binary AGB star with an alkaline jet
• Dynamical Architectures of S-type Transiting Planets in Binaries I: Target Selection using Hipparcos and Gaia proper motion anomalies
• A neutron star candidate in the long-period binary 56 UMa
• Global dynamics and architecture of the Kepler-444 system
• The NEID Earth Twin Survey. II. Dynamical Masses in Seven High-acceleration Star Systems
• An Edge-On Orbit for the Eccentric Long-Period Planet HR 5183 b
• Discovery of a cold giant planet and mass measurement of a hot super-Earth in the multi-planetary system WASP-132
• HD 28185 Revisited: An Outer Planet, Instead of a Brown Dwarf, On a Saturn-like Orbit
• Post-processing CHARIS integral field spectrograph data with PyKLIP
• Dynamical Masses and Ages of Sirius-like Systems
• Astrometric and photometric characterization of $\eta$ Tel B combining two decades of observations
• Distinct Rotational Evolution of Giant Planets and Brown Dwarf Companions
• A multi-technique detection of an eccentric giant planet around accelerating star HD 57625
• SCExAO/CHARIS and Gaia Direct Imaging and Astrometric Discovery of a Superjovian Planet 3--4 lambda/D from the Accelerating Star HIP 54515
• Updated Mutual Inclination Measurement for 14 Her b and c
• Orbital Parameters and Binary Properties of 37 FGK stars in the Cores of Open Clusters NGC 2516 and NGC 2422
• PHL 5038AB: Is the brown dwarf causing pollution of its white dwarf host star?
• Spectral analysis of two directly imaged benchmark L dwarf companions at the stellar-substellar boundary
• True Masses using RV data with Hipparcos and Gaia Astrometry
• Precise Masses and Orbits for Nine Radial Velocity Exoplanets
• Direct Imaging Discovery of Giant Exoplanet $\beta$ Pictoris d: A Decade-Long Game of Hide-and-Seek
• 14 Her: a likely case of planet-planet scattering
• A method for statistical research on binary stars using radial velocities
• HIP 67506 C: MagAO-X Confirmation of a New Low-Mass Stellar Companion to HIP 67506 A