Eclipse-timing study of new hierarchical triple star candidates in the northern continuous viewing zone of TESS

Abstract

Aims. We compiled a list of more than 3500 eclipsing binaries located in and near the northern continuous viewing zone (NCVZ) of the TESS space telescope that have sufficient TESS photometry to search for additional hidden components in these systems. In addition to discovering their hierarchical nature, we also determined their orbital parameters and analyzed their distributions. Methods. We obtained the TESS light curves of all targets in an automated way by applying convolution-aided differential photometry on the TESS full-frame images from all available sectors up to sector 60. Using a new self-developed Python GUI, we visually confirmed all of these light curves, determined the eclipsing periods of the objects, and calculated their eclipse-timing variations (ETVs). The ETV curves were used in order to search for nonlinear variations that could be attributed to a light travel-time effect (LTTE) or dynamical perturbations caused by additional components in these systems. We preselected 351 such candidates and modeled their ETVs with the analytic formulae of pure LTTE or with a combination of LTTE and dynamical perturbations. Results. We were able to fit a model solution for the ETVs of 135 hierarchical triple candidates, 10 systems of which were known from the literature, and the remaining 125 systems are new discoveries. These systems include some more noteworthy ones, such as five tight triples that are very close to their dynamical stability limit with a period ratio lower than 20, and three newly discovered triply eclipsing triples. We point out that dynamical perturbations occur in GZ Dra, which we found to be a triple, and that the system is one of the most strongly inclined systems known in the literature, with im = 58° ±7°. We also compared the distributions of some orbital parameters from our solutions with those from a previous Kepler sample. Finally, we verified the correlations between the available parameters for systems that have Gaia non-single star orbital solutions with those from our ETV solutions.