Probabilistic orbits and dynamical masses of emission-line binaries

TL;DR

This paper introduces a new Gaussian process-based method to analyze emission-line binary star systems, improving orbital parameter estimation by accounting for complex systematics, and provides benchmark datasets for validation.

Contribution

The paper develops a marginalised Gaussian process algorithm for emission-line binaries and offers benchmark datasets to standardize model validation in this field.

Findings

The marginalised GP outperforms standard algorithms on synthetic data with systematics.

Application to real binaries shows systematics affect orbital parameters and mass estimates.

Benchmark datasets facilitate future validation of analysis methods.

Abstract

The observed orbits of emission-line stars may be affected by systematics owing to their broad emission lines being formed in complex and extended environments. This is problematic when orbital parameter probability distributions are estimated assuming radial-velocity data are solely comprised of Keplerian motion plus Gaussian white noise, leading to overconfident and inaccurate orbital solutions, with implications for the inferred dynamical masses and hence evolutionary models. We present a framework in which these systems can be meaningfully analysed. We synthesise benchmark datasets, each with a different and challenging noise formulation, for testing the performance of different algorithms. We make these datasets freely available with the aim of making model validation an easy and standardised practice in this field. Next, we develop an application of Gaussian processes to model the radial-velocity systematics of emission-line binaries, named marginalised $\mathcal{GP}$. We benchmark this algorithm, along with current standardised algorithms, on the synthetic datasets and find our marginalised $\mathcal{GP}$ algorithm performs significantly better than the standard algorithms for data contaminated by systematics. Finally, we apply the marginalised $\mathcal{GP}$ algorithm to four prototypical emission-line binaries: Eta Carinae, GG Carinae, WR 140, and WR 133. We find systematics to be present in several of these case studies; and consequently, the predicted orbital parameter distributions, and dynamical masses, are modified from those previously determined.