Bayesian inference in single-line spectroscopic binaries with a visual orbit

Abstract

We present a Bayesian inference methodology for the estimation of orbital parameters on single-line spectroscopic binaries with astrometric data, based on the No-U-Turn sampler Markov chain Monte Carlo algorithm. Our approach is designed to provide a precise and efficient estimation of the joint posterior distribution of the orbital parameters in the presence of partial and heterogeneous observations. This scheme allows us to directly incorporate prior information about the system - in the form of a trigonometric parallax, and an estimation of the mass of the primary component from its spectral type - to constrain the range of solutions, and to estimate orbital parameters that cannot be usually determined (e.g. the individual component masses), due to the lack of observations or imprecise measurements. Our methodology is tested by analyzing the posterior distributions of well-studied double-line spectroscopic binaries treated as single-line binaries by omitting the radial velocity data of the secondary object. Our results show that the system's mass ratio can be estimated with an uncertainty smaller than 10% using our approach. As a proof of concept, the proposed methodology is applied to twelve single-line spectroscopic binaries with astrometric data that lacked a joint astrometric-spectroscopic solution, for which we provide full orbital elements. Our sample-based methodology allows us also to study the impact of different posterior distributions in the corresponding observations space. This novel analysis provides a better understanding of the effect of the different sources of information on the shape and uncertainty of the orbit and radial velocity curve.

arXiv commentAccepted for publication to the Astronomical Journal.39 pages, 7 Tables, 18 Figures, two Appendices

References in corpus (21)

• Array Programming with NumPy
• Gaia Data Release 2. Summary of the contents and survey properties
• The Gaia mission
• The SIMBAD astronomical database
• SB9: The Ninth Catalogue of Spectroscopic Binary Orbits
• Quantifying the Uncertainty in the Orbits of Extrasolar Planets
• Measuring the Recoverability of Close Binaries in Gaia DR2 with the Robo-AO Kepler Survey
• RUN DMC: An efficient, parallel code for analyzing Radial Velocity Observations using N-body Integrations and Differential Evolution Markov chain Monte Carlo
• An Affine-Invariant Sampler for Exoplanet Fitting and Discovery in Radial Velocity Data
• Astrometric orbit of a low-mass companion to an ultracool dwarf
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• Discovery of a Very Low Mass Triple with Late-M and T Dwarf Components: LP 704-48/SDSS J0006-0852AB
• Mass estimates for visual binaries with incomplete orbits
• Orbits for eighteen visual binaries and two double-line spectroscopic binaries observed with HRCAM on the CTIO SOAR 4m telescope, using a new Bayesian orbit code based on Markov Chain Monte Carlo
• Are the orbital poles of binary stars in the solar neighbourhood anisotropically distributed?
• Masses of the components of SB2 binaries observed with Gaia. II. Masses derived from PIONIER interferometric observations for Gaia validation
• Bayesian inference for binary neutron star inspirals using a Hamiltonian Monte Carlo Algorithm
• Sensitivity Analyses of Exoplanet Occurrence Rates from Kepler and Gaia
• Orbits and masses of binaries from Speckle Interferometry at SOAR
• Binary orbits from combined astrometric and spectroscopic data
• Visual binary stars with partially missing data: Introducing multiple imputation in astrometric analysis