We report the discovery and complex analyses of the first two compact hierarchical triple star systems discovered with TESS in or near its southern continuous viewing zone during Year 1. Both TICs 167692429 and 220397947 were previously unknown eclipsing binaries, and the presence of a third companion star was inferred from eclipse timing variations exhibiting signatures of strong 3rd-body perturbations and, in the first system, also from eclipse depth variations. We carried out comprehensive analyses, including the simultaneous photodynamical modelling of TESS and archival ground-based WASP lightcurves, as well as eclipse timing variation curves. Also, for the first time, we included in the simultaneous fits multiple star spectral energy distribution data and theoretical PARSEC stellar isochrones, taking into account Gaia DR2 parallaxes and cataloged metallicities. We find that both systems have twin F-star binaries and a lower mass tertiary star. In the TIC 167692429 system the inner binary is moderately inclined ($i_{mut}=27^o$) with respect to the outer orbit, and the binary vs. outer (triple) orbital periods are 10.3 vs. 331 days, respectively. The mutually inclined orbits cause a driven precession of the binary orbital plane which leads to the disappearance of binary eclipses for long intervals. In the case of TIC 220397947 the two orbital planes are more nearly aligned and the inner vs. outer orbital periods are 3.5 and 77 days, respectively. In the absence of radial velocity observations, we were unable to calculate highly accurate masses and ages for the two systems. According to stellar isochrones TIC 167692429 might be either a pre-main sequence or an older post-MS system. In the case of TIC 220397947 our solution prefers a young, pre-MS scenario.

The nearby (50 pc) K7V dwarf HIP~41431 (EPIC 212096658) is a compact 3-tier hierarchy. Three K7V stars with similar masses, from 0.61 to 0.63 solar, make a triple-lined spectroscopic system where the inner binary with a period of 2.9 days is eclipsing, and the outer companion on a 59-day orbit exerts strong dynamical influence revealed by the eclipse time variation in the {\em Kepler} photometry. Moreover, the centre-of-mass of the triple system moves on a 3.9-year orbit, modulating the proper motion. The mass of the 4-th star is 0.35 solar. The Kepler and ground-based photometry and radial velocities from four different spectrographs are used to adjust the spectro-photodynamical model that accounts for dynamical interaction. The mutual inclination between the two inner orbits is 2\fdg16$\pm$0\fdg11, while the outer orbit is inclined to their common plane by 21\degr$\pm$16\degr. The inner orbit precesses under the influence of both outer orbits, causing observable variation of the eclipse depth. Moreover, the phase of the inner binary is strongly modulated with a 59-day period and its line of apsides precesses. The middle orbit with eccentricity $e=0.28$ also precesses, causing the observed variation of its radial velocity curve. Masses and other parameters of stars in this unique hierarchy are determined. This system is dynamically stable and likely old.

We report a study of the eclipse timing variation (ETV) of short period ($P_1\le 6^d$) eclipsing binaries (EB) monitored during the photometric survey OGLE-IV. From the 425\,193 EBs we selected approximately 80\,000 binaries that we found suitable for further examination. Among them we identified 992 potential hierarchical triple (or multiple ) system candidates exhibiting light-travel-time effect (LTTE). Besides, we obtained the orbital parameters of these systems and carried out statistical analyses on the properties of these candidates. We found that (i) there is a significant lack of triple systems where the outer period is less than 500 days; (ii) the distribution of the outer eccentricities has a maximum around $e_2\approx0.3$; (iii) the outer mass ratio calculated from an estimated minimum mass of the third component is lower than $q_2\sim 0.5$ for the majority of the sample. We also present some systems that deserve special attention. (i) There are four candidates that show double periodic ETV, which we explain by the presence of a fourth companion. (ii) For two systems the perturbations of the third component are also found to be significant therefore we give a combined dynamical and LTTE ETV solution. (iii) For one system the third component is found to be probably in the substellar mass domain.

We report a comprehensive search for hierarchical triple stellar system candidates amongst eclipsing binaries (EB) observed by the CoRoT spacecraft. We calculate and check eclipse timing variation (ETV) diagrams for almost 1500 EBs in an automated manner. We identify five relatively short-period Algol systems for which our combined light curve and complex ETV analyses (including both the light-travel time effect and short-term dynamical third-body perturbations) resulted in consistent third-body solutions. The computed periods of the outer bodies are between 82 and 272 days, (with an alternative solution of 831 days for one of the targets). We find that the inner and outer orbits are near coplanar in all but one case. The dynamical masses of the outer subsystems determined from the ETV analyses are consistent with both the results of our light curve analyses and the spectroscopic information available in the literature. One of our candidate systems exhibits outer eclipsing events as well, the locations of which are in good agreement with the ETV solution. We also report another certain triply eclipsing triple system which, however, is lacking a reliable ETV solution due to the very short time range of the data, and four new blended systems (composite light curves of 2 eclipsing binaries each), where we cannot decide whether the components are gravitationally binded or not. Amongst these blended systems we identify the longest period and highest eccentricity eclipsing binary in the entire CoRoT sample.

We produce and analyze eclipse time variation (ETV) curves for some 2600 Kepler binaries. We find good to excellent evidence for a third body in 222 systems via either the light-travel-time (LTTE) or dynamical effect delays. Approximately half of these systems have been discussed in previous work, while the rest are newly reported here. Via detailed analysis of the ETV curves using high-level analytic approximations, we are able to extract system masses and information about the three-dimensional characteristics of the triple for 62 systems which exhibit both LTTE and dynamical delays; for the remaining 160 systems we give improved LTTE solutions. New techniques of preprocessing the flux time series are applied to eliminate false positive triples and to enhance the ETV curves. The set of triples with outer orbital periods shorter than ~2000 days is now sufficiently numerous for meaningful statistical analysis. We find that (i) as predicted, there is a peak near i_m~40 deg in the distribution of the triple vs. inner binary mutual inclination angles that provides strong confirmation of the operation of Kozai-Lidov cycles with tidal friction; (ii) the median eccentricity of the third-body orbits is e_2=0.35; (iii) there is a deficit of triple systems with binary periods <1 day and outer periods between ~50 and 200 days which might help guide the refinement of theories of the formation and evolution of close binaries; and (iv) the substantial fraction of Kepler binaries which have third-body companions is consistent with a very large fraction of all binaries being part of triples.

We report eclipse timing variation analyses of 26 compact hierarchical triple stars comprised of an eccentric eclipsing ('inner') binary and a relatively close tertiary component found in the {\em Kepler} field. We simultaneously fit the primary and secondary $O-C$ curves of each system for the light-travel time effect (LTTE), as well as dynamical perturbations caused by the tertiary on different timescales. For the first time, we include those contributions of three-body interactions which originate from the eccentric nature of the inner binary. These effects manifest themselves both on the period of the triple system, $P_2$, and on the longer "apse-node" timescale. We demonstrate that consideration of the dynamically forced rapid apsidal motion yields an efficient and independent tool for the determination of the binary orbit's eccentricity and orientation, as well as the 3D configuration of the triple. Modeling the forced apsidal motion also helps to resolve the degeneracy between the shapes of the LTTE and the dynamical delay terms on the $P_2$ timescale, due to the strong dependence of the apsidal motion period on the triple's mass ratio. This can lead to the independent determination of the binary and tertiary masses without the need for independent radial velocity measurements. Through the use of our analytic method for fitting $O-C$ curves we have obtained robust solutions for system parameters for the ten most ideal triples of our sample, and only somewhat less robust, but yet acceptable, fits for the remaining systems. Finally we study the results of our 26 system parameter fits via a set of distributions of various physically important parameters, including mutual inclination angle, and mass and period ratios.