Bayes-based orbital elements estimation in triple hierarchical stellar systems

Abstract

Under certain rather prevalent conditions (driven by dynamical orbital evolution), a hierarchical triple stellar system can be well approximated, from the standpoint of orbital parameter estimation, as two binary star systems combined. Even under this simplifying approximation, the inference of orbital elements is a challenging technical problem because of the high dimensionality of the parameter space, and the complex relationships between those parameters and the observations (astrometry and radial velocity). In this work we propose a new methodology for the study of triple hierarchical systems using a Bayesian Markov-Chain Monte Carlo-based framework. In particular, graphical models are introduced to describe the probabilistic relationship between parameters and observations in a dynamically self-consistent way. As information sources we consider the cases of isolated astrometry, isolated radial velocity, as well as the joint case with both types of measurements. Graphical models provide a novel way of performing a factorization of the joint distribution (of parameter and observations) in terms of conditional independent components (factors), so that the estimation can be performed in a two-stage process that combines different observations sequentially. Our framework is tested against three well-studied benchmark cases of triple systems, where we determine the inner and outer orbital elements, coupled with the mutual inclination of the orbits, and the individual stellar masses, along with posterior probability (density) distributions for all these parameters. Our results are found to be consistent with previous studies. We also provide a mathematical formalism to reduce the dimensionality in the parameter space for triple hierarchical stellar systems in general.

arXiv comment58 pages, 21 figures, 2 tables, 4 appendices. Accepted for publication in PASP

References in corpus (31)

• Validation of the new Hipparcos reduction
• MCMC using Hamiltonian dynamics
• The SIMBAD astronomical database
• Simulated Tempering: A New Monte Carlo Scheme
• Parallel Tempering: Theory, Applications, and New Perspectives
• The Eccentric Kozai-Lidov Effect and Its Applications
• SB9: The Ninth Catalogue of Spectroscopic Binary Orbits
• Quantifying the Uncertainty in the Orbits of Extrasolar Planets
• Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy
• Efficient fitting of multiplanet Keplerian models to radial velocity and astrometry data
• The Joker: A custom Monte Carlo sampler for binary-star and exoplanet radial velocity data
• Orbits for the Impatient: A Bayesian Rejection Sampling Method for Quickly Fitting the Orbits of Long-Period Exoplanets
• Constraining the masses of microlensing black holes and the mass gap with Gaia DR2
• Spectroscopic binaries in a sample of ROSAT X-ray sources south of the Taurus molecular clouds
• Retired A Stars and Their Companions IV. Seven Jovian Exoplanets from Keck Observatory
• Bayesian Re-analysis of the Gliese 581 Exoplanet System
• Astrometric orbit of a low-mass companion to an ultracool dwarf
• Mass estimates for visual binaries with incomplete orbits
• Relative orbit orientation in several resolved multiple systems
• Dancing twins: stellar hierarchies that formed sequentially?
• 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
• The Architecture of the GW Ori Young Triple Star System and Its Disk: Dynamical Masses, Mutual Inclinations, and Recurrent Eclipses
• Orbits and Masses in the multiple system LHS 1070
• Inner and Outer Orbits in 13 Resolved Hierarchical Stellar Systems
• Dynamical Masses and Stellar Evolutionary Model Predictions of M-Stars
• Kappa Fornaci, a triple radio-star
• Orbits and masses of binaries from Speckle Interferometry at SOAR
• Nearby low-mass triple system GJ795
• Binary orbits from combined astrometric and spectroscopic data
• Visual binary stars with partially missing data: Introducing multiple imputation in astrometric analysis
• Detecting Extra-solar Planets with a Bayesian hybrid MCMC Kepler periodogram

Cited by in corpus (3)

• Orbital elements and individual component masses from joint spectroscopic and astrometric data of double-line spectroscopic binaries
• Chasing extreme planetary architectures: I- HD196885Ab, a super-Jupiter dancing with two stars?
• Optimal observational scheduling framework for binary and multiple stellar systems