Analytical determination of orbital elements using Fourier analysis

Abstract

The ESA global astrometry space mission Gaia has been monitoring the position of a billion stars since 2014. The analysis of such a massive dataset is challenging in terms of the data processing involved. In particular, the blind detection and characterization of single or multiple companions to stars (planets, brown dwarfs, or stars) using Gaia astrometry requires highly efficient algorithms. In this article, we present a set of analytical methods to detect and characterize companions in scanning space astrometric time series as well as via a combination of astrometric and radial velocity time series. We propose a general linear periodogram framework and we derive analytical formulas for the false alarm probability (FAP) of periodogram peaks. Once a significant peak has been identified, we provide analytical estimates of all the orbital elements of the companion based on the Fourier decomposition of the signal. The periodogram, FAP, and orbital elements estimates can be computed for the astrometric and radial velocity time series separately or in tandem. These methods are complementary with more accurate and more computationally intensive numerical algorithms (e.g., least-squares minimization, Markov chain Monte-Carlo, genetic algorithms). In particular, our analytical approximations can be used as an initial condition to accelerate the convergence of numerical algorithms. Our formalism has been partially implemented in the Gaia exoplanet pipeline for the third Gaia data release. Since the Gaia astrometric time series are not yet publicly available, we illustrate our methods on the basis of Hipparcos data, together with on-ground CORALIE radial velocities, for three targets known to host a companion: HD 223636 (HIP 117622), HD 17289 (HIP 12726), and HD 3277 (HIP 2790).