A Python Code to Determine Orbital Parameters of Spectroscopic Binaries
Abstract
We present the open source Python code BinaryStarSolver that solves for the orbital elements of a spectroscopic binary system. Given a time-series of radial velocity measurements, six orbital parameters are determined: the long-term mean, or systemic, radial velocity, the velocity amplitude, the argument of periastron, the eccentricity, the epoch of periastron, and the orbital period referred to by $\{{\gamma, K, \omega, e, T_0, P}\}$ respectively. Also returned to the user is the projected length of the semi-major axis, $a_{1}\sin(i)$, and the mass function, $f(M)$. The determination of spectroscopic orbits and masses is an example of another important area of astrophysics, once the domain of professional astronomers, to which amateurs can now make significant contributions. This code, available from GitHub, is provided in support of that work, and should be of general use to the amateur and professional astronomical community.
- Publication:
-
arXiv e-prints
- Pub Date:
- November 2020
- DOI:
- 10.48550/arXiv.2011.13914
- arXiv:
- arXiv:2011.13914
- Bibcode:
- 2020arXiv201113914M
- Keywords:
-
- Astrophysics - Solar and Stellar Astrophysics;
- Astrophysics - Instrumentation and Methods for Astrophysics;
- J.2
- E-Print:
- 11 pages, 4 figures, submitted to the Journal of the American Association of Variable Star Observers (JAAVSO)