A Python Code to Determine Orbital Parameters of Spectroscopic Binaries

27 Nov 2020  ·  Nicholas Milson, Caroline Barton, Philip D. Bennett ·

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.

PDF Abstract