Next Article in Journal
GEMS Embeddings of Hayward Regular Black Holes in Massless and Massive Gravities
Next Article in Special Issue
The EBLM Project—From False Positives to Benchmark Stars and Circumbinary Exoplanets
Previous Article in Journal
Neutron Star Constraints on Neutron Dark Decays
Previous Article in Special Issue
Transiting Circumbinary Planets in the Era of Space-Based Photometric Surveys
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Observational Detection of Higher-Order Secular Perturbations in Tight Hierarchical Triple Stars

by
Tamás Borkovits
1,2,3,4,5,* and
Tibor Mitnyan
2
1
Baja Astronomical Observatory of University of Szeged, Szegedi út, Kt. 766, H-6500 Baja, Hungary
2
HUN-REN-SZTE Stellar Astrophysics Research Group, Szegedi út, Kt. 766, H-6500 Baja, Hungary
3
Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, HUN-REN Hungarian Research Network, Konkoly Thege Miklós út 15–17, H-1121 Budapest, Hungary
4
Gothard Astrophysical Observatory, ELTE Eötvös Loránd University, Szent Imre h. u. 112, H-9700 Szombathely, Hungary
5
HUN-REN-ELTE Exoplanet Research Group, Szent Imre h. u. 112, H-9700 Szombathely, Hungary
*
Author to whom correspondence should be addressed.
Universe 2023, 9(11), 485; https://doi.org/10.3390/universe9110485
Submission received: 16 October 2023 / Revised: 8 November 2023 / Accepted: 13 November 2023 / Published: 20 November 2023
(This article belongs to the Special Issue The Royal Road: Eclipsing Binaries and Transiting Exoplanets)

Abstract

:
In this work, we search for observational evidence of higher-order secular perturbations in three eclipsing binaries. These are slightly eccentric binaries, and they form the inner pairs of tight, compact, hierarchical triple star systems. Simultaneously, we analyze the high-precision satellite (Kepler and TESS) light curves; eclipse timing variations; combined spectral energy distributions (through catalog passband magnitudes); and, where available, radial velocities of KICs 9714358, 5771589, and TIC 219885468. Besides the determination of the robust astrophysical and dynamical properties of the three systems, we find evidence that the observed unusual eclipse timing variations of KIC 9714358 are a direct consequence of the octupole-order secular eccentricity perturbations forced by unusual, resonant behavior between the lines of the apsides of the inner and outer orbital ellipses. We also show that, despite its evident cyclic eclipse depth variations, KIC 5771589 is an almost perfectly coplanar system (to within 0 . 3 ), and we explain the rapid eclipse depth variations with the grazing nature of the eclipses. Finally, we find that the inner pair of TIC 219885468 consists of two twin stars; hence, in this triple there are no octupole-order three-body perturbations. Moreover, we show that this triple is also coplanar on the same level as the former one, but due to its deep eclipses, it does not exhibit eclipse depth variations. We intend to follow this work up with further analyses and a quantitative comparison of the theoretical and the observed perturbations.

1. Introduction

The new window that was opened up by Kepler spacecraft [1] to the exotic world of the tight, compact hierarchical triple stellar systems (TCHTs) has now become more widely opened and transparent with the ongoing operation of the TESS space telescope [2] mission. While Kepler had inevitable primacy in the discovery and principal characterizations of this rare new class of stellar systems, the multiple repeated revisitations of TESS, over five years, have allowed us to extend our knowledge about the properties of such systems. In a wider perspective, our studies of the stellar three-body problem have enhanced our understanding not only quantitatively but also qualitatively, for the reasons which we discuss briefly in the forthcoming text.
Similar to our previous works (e.g., [3]), we consider hierarchical triple or multiple (stellar) systems to be ‘tight’, if the orbital period ratio of the outer orbit to the inner orbit remains below ∼100 (i.e., P out / P in 100 ), while we call a triple or multiple star system ‘compact’ when the outer period does not exceed, let us say, 1000 days1. As is widely accepted (see, e.g., [4,5,6]), the only possible long-term stable configuration amongst three masses, being of the same order of mass, is the hierarchical configuration, i.e., when one of the three bodies is continuously located much farther from the other two than these latter two from each other. In such a system, the motion, in general, can be well approximated with two binaries, i.e., two Keplerian two-body motions. The inner binary is formed by the two closer objects, while the outer binary consists of the third, most distant component, and the center of mass of the inner binary. Hence, tightness chiefly characterizes the strength of the third-body perturbation of such systems. This is the case because the magnitude of the gravitational perturbations relative to the pure Keplerian two-body motions primarily depends on the ratio of the semi-major axes or, more strictly, on the instantaneous separations of the outer to the inner binaries. The former can readily be converted to such a direct observable as the period ratio.
On the other hand, compactness is directly connected to the physical dimensions of the system as a whole. This characteristic size of the given triple star, naturally, has very crucial astrophysical implications, not only in regard to the early and the late evolutionary phases of such systems (see, e.g., [7,8,9] and references therein) but in regard to how tidal and other non-third-body perturbations relate to the gravitational perturbations of the third mass. From our point of view, compactness primarily has practical significance. This significance has two somewhat counteracting aspects. The first aspect is connected to the largest amplitude, and most interesting, orbital perturbations in tight triple star systems, the so-called long-period or apse-node perturbations—see below—which are of primary interest in this paper. The characteristic periods of these perturbations scale with the ratio of P out 2 / P in , i.e., these are essentially driven by the product of the measures of the tightness and the compactness. Hence, the more compact a tight triple is, the shorter the observational window or the length of the necessary dataset for the direct detection and investigation of such effects. Considering the second aspect, however, such triples can be serendipitously discovered and characterized most easily through the medium-period third-body perturbations (see again, below) that drive eclipse timing variations (ETVs) of such eclipsing binaries (EBs) that form the inner binary of a tight hierarchical triple system. The amplitude of such medium-timescale dynamically driven ETVs, however, scales as P in 2 / P out and, in such a way, for the most frequently occurring and shortest period EBs (some hours to few days), the dynamically driven ETVs may remain under the detection limit even in the tightest triples.
This latter fact explains why TCHTs were almost completely unknown before the four-year-long observations of the prime Kepler mission. The vast majority of the EBs known before Kepler had periods of less than 3–4 days, in which case there was no chance of detecting the gravitational third-body perturbations via eclipse timing variations with the accuracies of the highly inhomogeneous ground-based observations2. The quasi-continuous, almost four-year-long observations of the Kepler space telescope have led to the discovery of more than a thousand new EBs with eclipsing periods longer than 5 days. Amongst them, Borkovits et al. [11,12] have identified 46 tight hierarchical triple star candidates through the manifestations of the medium-term ( P out -period) third-body perturbations in their ETVs. (Note, they found 16 further tight triple star candidates amongst the almost 2000 additional, shorter period Kepler EBs, but this was mainly a consequence of the unprecedentedly precise observations of the space telescope and/or other effects, such as, e.g., third-body eclipses.) From this sample of the newly discovered 62 tight triple systems, 39 can be ranked as TCHT, i.e., a triple star system that is not only tight but compact as well. While these former studies depend primarily on the P out -timescale, the so-called medium period perturbations, the four-year-long ETVs of most of the eccentric TCHTs, also exhibited signals of rapid, dynamically driven apsidal motion. Hence, a simplified, quadrupole approximation of the dynamically forced apsidal motion (and of the nodal regression) was included in the analytic ETV model.
The TESS spacecraft, Kepler’s successor, first reobserved the prime Kepler-field almost exactly a decade after Kepler began its survey. Since then, the vast majority of the original Kepler targets and, hence, the Kepler-discovered TCHTs have been revisited typically 4–7 times, each for a four-week-long observing session. Thus, currently the lengths of the available datasets are nearly one and a half decades. Naturally, 14 years is tiny in the context of most interesting astronomical timescales, but in the case of TCHTs, this interval is long enough for the robust detection and quantitative analysis of the ‘long-period’ (or, secular) third-body perturbations. In this paper, we have selected a few such TCHTs from the original Kepler sample, which are exceptional even amongst the TCHTs. In particular, their special properties make possible the clear detection and investigation of even the higher-order secular perturbations with the exclusive use of the high-precision satellite data. In what follows, first we give a brief review in Section 2 of the analytic theory of the doubly averaged or secular stellar three-body problem, focusing on the observable quadrupole and octupole perturbation terms in the low mutual inclination domain. Then, in Section 3 we introduce the three TCHTs that were selected for the current analysis. The details of the observational materials, the data preparation, and the complex photodynamical analysis are described in Section 4 and Section 5, while the results are discussed in Section 6.

2. Dynamics of TCHTs

The hierarchical three-body problem is sometimes called the ‘stellar three-body problem’; however, the fields where the hierarchical treatment can be applied are much wider than the domain of triple and multiple stars. For example, the hierarchical three-body treatment can be used for investigations of the perturbations of the Earth–Moon system due to the Sun [13,14,15], Jupiter’s effects on the motion of the main belt asteroids and/or comets [16,17], and the perturbations on a binary comprising the Earth and an artificial satellite due to the Moon [18].
To our knowledge, Slavenas [19,20] was the first to investigate the hierarchical three-body problem in the context of a stellar triple ( λ Tauri, the very first and, for a long-time, the only TCHT). Due to the fact, however, that there were almost no known tight, or even compact, triple stellar systems, the hierarchical or stellar three-body problem, apart from a very few exceptions (see, e.g., [21,22,23,24,25,26]), has remained beyond the scope of the scientific community. The situation changed completely and quickly when the Kozai mechanism (nowadays called the ‘Von Zeipel–Lidov–Kozai’—ZLK—mechanism) or cycles were essentially ‘rediscovered’ in regard to the formation and evolution of stellar triples during the late 1990s. At that time, the first observational evidence was found in triple stellar and planetary systems for the effectiveness of secular third-body perturbations (see, e.g., [27,28,29]), and, practically, at the same time, in an epochal paper Kiseleva et al. [30] proposed a combined mechanism of third-body perturbations with tidal interactions to explain the formation of the close binary systems. In our view, this last work transformed the ZLK mechanism into its current, widespread, and acknowledged position.
The mechanism, which is currently known as ZLK theory, is practically nothing more than the secular, doubly averaged theory of the hierarchical or stellar three-body problem. This theory describes the long-term dynamical evolution of such a three-body system, where the evolution is driven, in its pure form, only by the mutual gravitational interactions of the three bodies.
If we slightly modify the original classification of Brown [14], in a hierarchical triple system the perturbations are effective over three different timescales, and they can be classified according to these characteristic timescales as follows:
  • short-period perturbations, for which the typical period is in the order of P in , and the amplitude is related to ( P in / P out ) 2 ,
  • medium-period perturbations, with a characteristic period of P out and an amplitude of P in / P out ,
  • long-period perturbations3 having a period of about P out 2 / P in and an amplitude that may reach unity.
It was Harrington [21,22] who, for the first time, described the equations of the stellar three-body problem in Hamiltonian formalism for any arbitrary values of outer eccentricities and mutual inclinations and gave third-order (sometimes called: ‘octupole’) solutions for the long-period or apse-node time-scale variations of the orbital elements, with the application of the double-averaging method of [31,32,33,34].
While the current discoveries of TCHTs (including the currently analyzed systems, as well) are mainly based on the medium-period perturbations of the ETVs of their inner EBs, the scope of this paper is connected to the long-period perturbations, and thus, in what follows, we concentrate only on the latter. The original ZLK theorem is restricted to the lowest-order perturbative terms of the doubly averaged Hamiltonian (which is quadratic in the small parameter α = a in / a out ; hence, it is called the ‘quadrupole approximation’). At that level, the problem has only one degree of freedom, and thus it is integrable; the solution can be given with the use of elliptic functions of Weierstrass P . Note, although the original analytic solution of Kozai [17] is strictly valid only if the outer orbit coincides with the invariable plane of the triple (i.e., all the angular momentum of the system is stored in the outer orbit), as was shown, e.g., by Harrington [21] and Söderhjelm [24], the same solution may remain valid, with small modifications, even in the case where the orbital angular momentum of the inner orbit is non-negligible but small enough. On the other hand, however, as was shown by Ford et al. [35] via numerical integrations and, later analytically by Naoz et al. [36] (amongst others), in the case of an eccentric outer orbit, the octupole-order terms may substantially alter the long-term behavior of a hierarchical triple system.
As was excellently reviewed by Naoz [37], the vast majority of the studies of the doubly averaged stellar three-body problem (being either analytical, numerical or, both) are directed toward the investigation of the very long-timescale evolution of such systems. Specifically, their interests are mainly focused on configurations where large amplitude eccentricity and (mutual) inclination cycles (including prograde-retrograde flip-flops) may occur. Much less effort has been directed toward the short-term (human timescale), directly observable effects of the stellar three-body problem. In this latter context, besides the above mentioned few exceptions, from the point of view of the present work, the paper of Krymolowski & Mazeh [38] is extremely relevant. These authors compared the outcomes of the quadrupole and octupole-order perturbation theories of the doubly averaged problem to each other (and to the results of direct numerical integrations) for different mutual inclinations (from the coplanar case up to i mut = 60 ). They found that the octupole perturbations may be very significant even in the coplanar (or, nearly coplanar) case. This is so because, in an exactly coplanar configuration, the long-period quadrupole perturbations, for example, in the inner eccentricity completely disappear. Despite this fact, numerical integrations show that there long-period perturbations exist in the inner eccentricity in the coplanar configurations, and this effect is chiefly determined by the octupole-order perturbations.
As those TCHTs that we investigate in this paper are actually nearly coplanar systems, in what follows, first we discuss the octupole effects quantitatively. Following the formulae of Harrington [21] but using partly different notations, the doubly averaged Hamiltonian up to the octupole-order takes the form:
H ( Δ Ω dyn = π ) = β 2 2 + 3 e in 2 ( 3 I 2 1 ) + 15 e in 2 1 I 2 cos 2 ω in dyn + β 3 e in e out 1 2 1 + 3 4 e in 2 1 + 11 I 5 I 2 15 I 3 cos ( ω in dyn ω out dyn ) + 1 2 1 + 3 4 e in 2 1 11 I 5 I 2 + 15 I 3 cos ( ω in dyn + ω out dyn ) 35 8 e in 2 1 + I I 2 I 3 cos ( 3 ω in dyn ω out dyn ) 35 8 e in 2 1 I I 2 + I 3 cos ( 3 ω in dyn + ω out dyn ) ,
where e in , out and ω in , out dyn denote the eccentricities and dynamical arguments of periastron of the inner and outer orbits, respectively. These latter quantities (sometimes denoted as the Delaunay variables g in , out ), which should not be confused with the observable arguments of periastrons (the latter of which will hereafter be denoted as ω in , out ), give the angular distances of the pericenter points of the orbits from the ascending nodes of the intersections of the orbital planes. Furthermore, I = cos i mut stands for the cosine of the mutual (or, relative) inclination of the two orbital planes. Finally,
β 2 = G 2 16 m A 7 m AB 3 m Aa 7 ( m Aa m Ab ) 3 L in 4 L out 3 G out 3 = 2 π 30 A G L in ,
β 3 = 15 4 G 2 16 m A 9 m AB 4 m B 9 ( m Aa m Ab ) ( m Aa m Ab ) 5 L in 6 L out 3 G out 5 = 4 π A G oct L in ,
where we have introduced the parameters A G and A G oct , which are more closely related to observable quantities and can be expressed as
A G = 15 8 q out 1 + q out P in P out 2 1 e out 2 3 / 2 ,
A G oct = 1 q in 1 q out 1 1 + q out 1 / 3 P in P out 2 / 3 A G 1 e out 2 .
Additional, newly introduced quantities in the equations above are the gravitational constant G; the individual masses of the inner, close binary stars, as  m A a , m A b ; the total mass of the inner pair, m A = m Aa + m Ab ; the individual mass of the third component, m B ; and the total mass of the triple system, m AB = m A + m B . Moreover, the mass ratios are denoted by q in , out for the inner and outer binaries, respectively.
In regard to the amplitude-like quantity of A G oct , as one can see, it disappears for q in = 1 , i.e., when the inner pair consists of two equal mass stars4. It is common also to introduce the quantity of
ϵ = e out A G oct A G = 1 q in 1 q out 1 1 + q out 1 / 3 P in P out 2 / 3 e out 1 e out 2 ,
which gives the strength of the octupole-order terms relative to the quadrupole ones.
Finally, note that we use the usual Delaunay action variables as
L in = m Aa m Ab m A G m A a in ,
L out = m A m B m AB G m AB a out ,
G in , out = L in , out 1 e in , out 2 ,
H in , out = G in , out cos i in , out dyn ,
while the conjugate angular variables are the mean anomalies ( l in , out ), the dynamical arguments of pericenters ( ω in dyn —or, traditionally, g in , out ), and the dynamical longitudes of the nodes ( Ω in , out dyn —or, h in , out ) of the inner and outer orbits, respectively. One should also keep in mind that, as was noted by Naoz et al. [36], although Equation (1) formally does not contain the nodes ( Ω in , out dyn ), which would suggest the constancy of the conjugate variables H in , out , this arises only from an incorrect application of the elimination of the nodes. For a correct treatment, one should take it into account that, in the original Hamiltonian, both the inner and outer dynamical nodes are present but only through the sines and cosines of their differences, which is a constant Δ Ω dyn = Ω in dyn Ω in out = π and, therefore, sin Δ Ω dyn = 0 , and  cos Δ Ω dyn = 1 . If one substitutes these values into the original Hamiltonian, it leads to the usual form of Equation (1). Hence, strictly speaking, from the simplified (Equation (1)) form of the Hamiltonian, the constancy of H in , out does not follow. This fact, however, does not affect the calculation of the perturbation equations of the other elements through the usual, formal way, and hence, in what follows, we do not take it into account.
The perturbation equations that are of interest to us take the following form:
d e in d τ = 1 e in 2 e in 1 G in d G in d τ
= 2 π A G e in 1 e in 2 1 / 2 1 I 2 sin 2 ω in dyn + 2 π A G oct e out 1 e in 2 1 / 2 1 + 3 4 e in 2 1 + 11 I 5 I 2 15 I 3 sin ( ω in dyn ω out dyn ) + 1 11 I 5 I 2 + 15 I 3 sin ( ω in dyn + ω out dyn ) 105 4 e in 2 1 + I I 2 I 3 sin ( 3 ω in dyn ω out dyn ) + 1 I I 2 + I 3 sin ( 3 ω in dyn + ω out dyn ) ,
d ω in dyn d τ = H G in = 2 π A G 1 e in 2 1 / 2 I 2 1 5 1 e in 2 + 2 + 3 e in 2 5 G in G out I + 1 e in 2 I 2 e in 2 G in G out I cos 2 ω in dyn + 2 π A G oct 1 e in 2 1 / 2 e in e out 1 e in 2 e in 2 1 + 9 4 e in 2 1 + 11 I 5 I 2 15 I 3 + 1 + 3 4 e in 2 I + G in G out 11 10 I 45 I 2 cos ( ω in dyn ω out dyn ) . + 1 e in 2 e in 2 1 + 9 4 e in 2 1 11 I 5 I 2 + 15 I 3 1 + 3 4 e in 2 I + G in G out 11 + 10 I 45 I 2 cos ( ω in dyn + ω out dyn ) 35 4 e in 2 3 1 e in 2 e in 2 1 + I I 2 I 3 + I + G in G out 1 2 I 3 I 2 cos ( 3 ω in dyn ω out dyn ) 35 4 e in 2 3 1 e in 2 e in 2 1 I I 2 + I 3 I + G in G out 1 + 2 I 3 I 2 cos ( 3 ω in dyn + ω out dyn ) ,
d Ω in dyn d τ = H H in = 2 π 5 A G 1 e in 2 1 / 2 C G out I 2 + 3 e in 2 5 e in 2 cos 2 ω in dyn 2 π A G oct 1 e in 2 1 / 2 C G out e in e out 1 + 3 4 e in 2 11 10 I 45 I 2 cos ( ω in dyn ω out dyn ) 11 + 10 I 45 I 2 cos ( ω in dyn + ω out dyn ) 35 4 e in 2 1 2 I 3 I 2 cos ( 3 ω in dyn ω out dyn ) 1 + 2 I 3 I 2 cos ( 3 ω in dyn + ω out dyn ) ,
and, regarding the perturbations in the observable arguments of periastron:
d ω in d τ = d ω in dyn d τ + d Ω in , out dyn d τ cos i in dyn d Ω in d τ cos i in
= 2 π 5 A G 1 e in 2 1 / 2 3 I 2 1 + 5 1 I 2 cos 2 ω in dyn + 2 π A G oct 1 e in 2 1 / 2 e out e in 1 + 9 4 e in 2 1 + 11 I 5 I 2 15 I 3 cos ( ω in dyn ω out dyn ) + 1 11 I 5 I 2 + 15 I 3 cos ( ω in dyn + ω out dyn ) 105 4 e in 2 1 + I I 2 I 3 cos ( 3 ω in dyn ω out dyn ) + 1 I I 2 + I 3 cos ( 3 ω in dyn + ω out dyn ) d Ω in d τ cos i in
d ω out d τ = d ω out dyn d τ + d Ω in , out d τ cos i out dyn d Ω out d τ cos i out = 2 π 10 A G 1 e in 2 1 / 2 G in G out 2 + 3 e in 2 3 I 2 1 + 15 e in 2 1 I 2 cos 2 ω in dyn + 2 π A G oct 1 e in 2 1 / 2 G in G out e in 1 + 4 e out 2 e out 1 + 11 I 5 I 2 15 I 3 cos ( ω in dyn ω out dyn ) + 1 11 I 5 I 2 + 15 I 3 cos ( ω in dyn + ω out dyn ) 35 4 e in 2 1 + I I 2 I 3 cos ( 3 ω in dyn ω out dyn ) + 1 I I 2 + I 3 cos ( 3 ω in dyn + ω out dyn ) d Ω out d τ cos i out .
In what follows, we consider only the prograde, coplanar scenario, i.e., when i mut = 0 . This restriction can be justified by the fact that, as will be shown in Section 6, the mutual inclination in none of the three considered systems exceeds i mut = 0 . 5 . In the coplanar case, when I = 1 , Equations (12)–(17) become much simpler. Before we give these simpler expressions, we note that, for coplanar orbits, the dynamical arguments of periastron, and the dynamical nodes, lose their meanings. The quantities that continue to remain meaningful in this situation are the dynamical longitudes of periastron ( ϖ in , out dyn = ω in , out dyn + Ω in , out dyn ), as well as the angle between the directions of the two periastron points, i.e., ω in dyn ω out dyn = ϖ in dyn ϖ out dyn π , where in the last expression we took into account the fact that Δ Ω dyn = π . Furthermore, we notice that in the coplanar case ω in , out = ϖ in , out dyn . Thus, as can be easily seen, for the prograde, coplanar case, the perturbation equations up to the octupole-order have the following forms:
d e in d τ = 16 π A G oct e out 1 e in 2 1 / 2 1 + 3 4 e in 2 sin ( ω in ω out ) ,
d ω in d τ = 4 π 5 A G 1 e in 2 1 / 2 + 16 π A G oct 1 e in 2 1 / 2 e out e in 1 + 9 4 e in 2 cos ( ω in ω out )
d ω out d τ = 4 π 5 A G 1 e in 2 1 / 2 G in G out 1 + 3 2 e in 2 + 16 π A G oct 1 e in 2 1 / 2 G in G out e in 1 + 4 e out 2 e out cos ( ω in ω out )
and, finally,
d ( ω in ω out ) d τ = 4 π 5 A G 1 e in 2 1 / 2 1 G in G out 1 + 3 2 e in 2 1 e in 2 + 16 π A G oct 1 e in 2 1 / 2 e in e out × 1 + 9 4 e in 2 e in 2 G in G out 1 1 e in 2 1 + 4 e out 2 e out 2 cos ( ω in ω out ) .
As can be easily seen, in this scenario there are no quadrupole-level long-period perturbations in the inner eccentricity, and the observable dynamical apsidal motion rates of both orbits remain constant (again, at the quadrupole level). In contrast to this, the octupole-order perturbations of these elements depend on the trigonometric functions of the angle between the two periastron directions. One should keep in mind that this angle, in general, varies at a lower rate than the individual apsidal lines themselves; hence, one should expect a longer period cyclic, octupole variation than the usual (quadrupole) dynamically forced apsidal motion period.
A detailed analytic and/or numerical analysis of these perturbation equations is beyond the scope of the current paper. Our aim remains only to identify the consequences of these orbital perturbations directly from the high-precision observations of selected TCHTs; to determine the accurate dynamical (and astrophysical) parameters of these systems; and, finally, to compare our findings at least qualitatively with the predictions of the theory, discussed briefly above. A more detailed study including the quantitative comparison of the theoretical and the observed perturbations will be published later.

3. Selected Systems

We selected two TCHTs from the well-observed primary Kepler-sample and, moreover, a third system, which serves as a counterexample in the sense that, although this system is very similar to the previous two in several aspects, the octupole-order perturbations are almost nulled out due to a near unit mass ratio of the inner binary. Unfortunately, we did not find such an illustrative TCHT amongst the Kepler systems; hence, we took one from the TESS northern continuous viewing zone (NCVZ) sample. The main catalog data for the three systems are collected in Table 1, where, in addition to the J2000.0 coordinates and the catalog magnitudes from near ultraviolet to infrared passbands (i.e., form GALEX NUV to Wise W4 magnitudes), we also tabulate the effective temperatures ( T eff ), photogeometric distances, metallicities ([ M / H ]), interstellar reddenings ( E ( B V ) ), and proper motions ( μ α , δ ) from different catalogs. Finally, we also give the specific R U W E (renormailzed unit weight error) parameters, which were introduced in Gaia DR2 as an indicator of the quality of the astrometric solutions. This parameter is relevant in the context of the current study because it was found that a value greater than ≳1.4 could indicate the multiplicity of the source (see, e.g., [39]). In what follows, we briefly discuss the available basic information on these three TCHTs.
KIC 9714359 = KOI-6073 ( P in = 6.48  d; P out = 103.8  d; P out / P in = 16.0 , ϵ 0.016 ) is the compact, tight triple star system that chiefly inspired this study. This target was identified first as an EB in the first release of the Kepler EB catalog [48]. The hierarchical triple star nature of this source was reported first by Rappaport et al. [49], who analyzed the ETVs of the first 13 quarters of data of Kepler EBs, who also pointed out that the ETV is dominated almost exclusively by the dynamical, third-body perturbation effects over the classic light-travel time effect (LTTE). Then, the system was included in the sample of 26 eccentric EBs that were comprehensively reanalyzed by Borkovits et al. [11] with the use of their most complex analytic ETV model. Despite the extreme tightness of the triple, they found a quite satisfactory analytical, coplanar ETV solution with a slightly eccentric inner orbit, a moderately eccentric outer orbit ( e in = 0.015 ± 0.001 and e out = 0.30 ± 0.01 , respectively), and a dynamically forced apsidal motion period of P apse , in = 30.6  yr. The model also indicated inner and outer mass ratios of q in = 0.45 ± 0.01 and q out = 0.21 ± 0.01 , respectively.
Further studies, however, resulted in slightly discrepant values for some of the parameters above. First, Kjurkchieva et al. [50] found a substantially (∼5–6 times) larger inner eccentricity of e in = 0.085 from the Kepler light curves, and, moreover, based on astrophysical reasons, they estimated a lower inner mass ratio ( q in = 0.33 ). Then, Windemuth et al. [51] also reanalyzed the Kepler light curves, adding SED information and theoretical PARSEC [52] evolutionary tracks to them in an automated way. While their eccentricity of e in = 0.016 ± 0.002 was really in good agreement with that of the dynamically inferred value of Borkovits et al. [11], they obtained an even much more extreme inner mass ratio ( q in = 0.234 ), and, correspondingly, they found that the system was in a pre-MS state with an age of log τ 7.44 ± 0.04 .
While these discrepancies, noted above, might be interesting in and of themselves, what makes this system especially important for us is the unusual behavior of the ETV curves during the TESS observations. KIC 9714358 was reobserved by TESS in Sectors 14, 15, 40, 54, and 55 in the lower cadence full frame image (FFI) mode. As one can readily see in the upper left panel of Figure 1, the new ETV points, determined from TESS-observed eclipses, do not follow the former trend, and the long-term behavior of the red (primary) and blue (secondary) ETV curves can no longer be described with two antiphased cosine curves, as one would expect for an EB with low eccentricity and an (almost constant) apsidal motion rate (the latter of which is a clear prediction of the quadrupole-order hierarchical three-body problem). Hence, our immediate suspicion was that what we are seeing is nothing more than a direct, observational manifestation of the octupole-order perturbation effects. We will return to this question in Section 6.
Finally, note that this is the only system in our sample for which radial velocity (RV) data are available from the literature. The data release 16 of the APOGEE-2 project [53] tabulates 16 individual RV observations, which we used for our analysis (see Figure 2).
KIC 5771589 = KOI-6625 ( P in = 10.79  d; P out = 113.1  d; P out / P in = 10.5 , ϵ 0.004 ). This system was identified as an EB and then as a triple star candidate in the same sequence as for the previous target. It was also included in the detailed analytical ETV studies of Borkovits et al. [11,12]. This is the second tightest known triple in the whole studied Kepler sample, and the only one for which the ETV curves, determined purely from the ∼4 yr-long original Kepler observations, cover more than half of an apsidal cycle5. The analytical ETV studies suggested a large outer mass ratio ( q out = 0.98 ± 0.06 ), which was strongly supported by the very shallow eclipses (∼0.2–0.4%).
The most remarkable feature of the Kepler light curve is, however, the varying depths of the eclipses (see Figure 3), implying orbital plane precession forced by the not-exactly coplanar tertiary. According to our knowledge, there are only two EBs in the whole primary Kepler sample that exhibit a clear reversal of the trend of the eclipse depth variations (EDVs) during the four-year-long observations. The other such target is the triply eclipsing triple star system KIC 6964043 [3], which displayed an opposite trend in its EDV compared to KIC 5771589. In the former case, the eclipse depths increased first and then decreased, while in the currently considered KIC 5771589 the eclipses became shallower and shallower during the first half of the Kepler observations, and then they started to increase again. These two situations are not symmetric with respect to each other. This is because an EB can reach the maximum of its eclipse depths in two different ways during a nodal precession cycle. The EB may either pass through the exact edge-on view (i.e., i = 90 , which, as was pointed out by Borkovits et al. [3], really did happen in KIC 6964043) without reaching an extremum in its observable inclination, or it may reach its nearest position to an edge-on view (i.e., | cos i | has a minimum but not zero value or, of course, the inclination itself has an extremum). In contrast to this, when eclipse depths reach a minimum value, one can be certain that the visible inclination has reached its extremum farthest from i = 90 6. In such a way, the current light curve behavior offers more strict constraints on the orbital configuration.
KIC 5771589 was reobserved with the TESS space telescope in Sectors 14, 26, 40, 41, 53, and 54. In contrast to KIC 9714358, for this target 2-min cadence observations are also available for all sectors. Unfortunately, despite the relative brightness of the system, due to the very shallow eclipse depths, the TESS data have much lower quality. Despite the low S/N ratios, however, very shallow eclipses can still be identified at the expected times, but in the Year 2 data they are almost lost in the noise. By contrast, the eclipses can be seen clearly at the end of the Year 4 observations, indicating that the eclipse depths must have increased again between the Year 2 and Year 4 TESS observations (see the lower right panel of Figure 3).
TIC 219885468 has P in = 7.54  d; P out = 110.8  d; P out / P in = 14.7 ; and  ϵ 0.0004 . This system was not observed with the Kepler spacecraft, but instead it is located near to the northern ecliptic pole, in the NCVZ of the TESS space telescope. But, we did not choose it primarily for its relatively long data train but rather for didactic reasons. In particular, we wanted to illustrate that, despite the clear similarity to the above-mentioned KIC 97154358 (in both their inner and outer periods), due to the twin-nature of the inner binary star, and hence, in the absence of octupole-order perturbations, the short-term orbital evolution of the two triples differs significantly.
TIC 219885468, a previously unknown EB, was identified as a likely triple star candidate during our ongoing study of TESS EBs in the NCVZ, which surveys the ETVs of these objects to identify triple star candidates through their eclipse timing variations (Mitnyan et al., to be submitted soon). The ETVs of this EB display typical timing variations that are dominated by third-body perturbations, including medium-term effects as well as dynamically forced rapid apsidal motion (see bottom row of Figure 1). The primary and secondary eclipses look very similar in depth, which indicates that the inner mass ratio should be near unity—and this nulls out the octupole perturbation terms. Moreover, no eclipse depth variations are observed during the ∼1300 days of TESS observations, which suggests a strong coplanarity (see Figure 4).

4. Observational Data and Its Preparation for Analysis

Two of our three targets were observed nearly continuously in long cadence (LC; time resolution of 29.4 min) mode during the prime Kepler mission. Furthermore, for KIC 5771589 one month of short cadence (SC; time resolution of 58.9 s) data in Quarter 4.1 is also available. For our analysis, we used those preprocessed LC datasets, which can be publicly downloaded from the Villanova website7 of the third revision of the Kepler Eclipsing Binary Catalog [56,57]. Quarters 4, 12, and 13 data of KIC 9714358, however, are not available at this site; hence, these were downloaded directly from the MAST database8.
Regarding the TESS observations, we downloaded the FFIs of all three targets for all sectors, when they were observed. (As is well known, Year 2 FFIs were obtained with 30-min cadence time and Year 4 ones with 10-min cadence time, while Year 5 FFIs, which are available only for TIC 219885468 among our targets, have a 200-s cadence time.) We then processed these FFIs using a convolution-based differential photometric pipeline implemented in the FITSH package [58]. Finally, we detrended these raw light curves, which were affected by differing amounts of scattered light and other effects, with smoothing polynomials of various degrees to each light curve section.
After the detrending processes, as a first step, we determined mid-eclipse times for each individual eclipse, in the very same manner as was described in Borkovits et al. [12]9. Then, preparing the dataset for the light curve analysis section of our complex photodynamical treatment, first, we binned the Year 4–5 data into 1800 s bins for a homogeneous sampling (not only with respect to Year 2 TESS data but also to the Kepler LC data). Then, in order to reduce computational costs and give higher statistical weight to the eclipses (which carry most of the relevant astrophysical and dynamical information amongst the light curve points), we dropped the larger parts of the out-of-eclipse light curves and kept only small light-curve sections around the primary and secondary eclipses.
In the case of KIC 9714358, we downloaded the 16 heliocentric RV data points, published in APOGEE DR 16 [53], from the Vizier site10 and used them.
Finally, note that, similar to our previous works, we included a composite spectral energy distribution (SED) analysis in our joint photodynamical runs (see below), for which we mainly used those catalog magnitudes, which are tabulated in Table 1.

5. Photodynamical Modeling

We have carried out complex photodynamical analyses of our three targets with our own software package Lightcurvefactory, developed by one of us (T. Borkovits, see, e.g., [59,60] and references therein). This code contains (i) a built-in numerical integrator to calculate the gravitationally (three-body effects), tidally, and (optionally) relativistically perturbed Jacobian coordinates and velocities of the three stars in the system; (ii) emulators for multiband light curves, ETV curves, and RV curves; (iii) built-in, tabulated PARSEC grids [52] for interpolating fundamental stellar parameters (e.g., radii and effective temperatures) as well as a composite SED model for the three stars; and (iv) an MCMC-based search routine for fitting the parameters. The latter applies our own implementation of the generic Metropolis–Hastings algorithm (see, e.g., [61]). The use of this software package and the consecutive steps of the entire analysis process have been previously explained in detail in several former papers (see, e.g., [59,60,62,63]). In this regard, we note that, despite the fact that Lightcurvefactory has a built-in subroutine for correcting the emulated light curves for finite cadence times (see [59] for the description of the process), due to the relatively long cadence times of 1800 s, we did not applied such corrections. This decision can be justified a posteriori by the fact that applying cadence time corrections in the case of KIC 9714358 resulted in parameters that remained within the ± 1 σ domains of the uncorrected results (tabulated in Table 2)—even in the case of the most sensitive parameters (the inclinations). However, since the analysis runs with the cadence time corrections result in computation times about four-five times longer (on average), we declined to use this feature routinely.
In the absence of an RV curve for each of the three stars in any of our triple systems (or even for two of the stars), we combined a composite SED analysis and the use of precalculated PARSEC grids as proxies to determine the masses of the constituent stars. Naturally, such a photodynamical model solution is no longer astrophysically model-independent. Nonetheless, in a previous work [64] we have shown that for compact, dynamically interacting triple stars (as is the case for the currently investigated systems), such a solution may result in stellar masses and radii within 5–10% of the results of an astrophysical model-independent solution, the latter of which is based on dynamical masses from RV data. The geometric and dynamical parameters (i.e., the orbital elements and also the mass ratios), as well as dimensionless quantities such as the fractional radii and temperature ratios of the stars, remain practically unchanged between the two kinds of solutions.
In the case of these model-dependent runs, the adjusted parameters were as follows:
(i)
Stars: Three stellar-mass-related parameters: the mass of the most massive component (being either m Aa , the primary of the inner pair, or m B , the tertiary star), and the inner and outer mass ratios ( q in , out ). Moreover, the metallicity of the system ([ M / H ]), the (logarithmic) age of the three coeval stars ( log τ ), and the interstellar reddening E ( B V ) for the given triple were also varied. Additionally, the ‘extra light’ contamination 4 parameter was also freely adjusted for two of the three systems.
(ii)
Orbits: Three of six orbital-element related parameters of the inner orbits, and six parameters of the outer orbits, i.e., the eccentricity vector components of both orbits ( e sin ω ) in , out , ( e cos ω ) in , out ; the inclinations relative to the plane of the sky ( i in , out ); and three other parameters for the outer orbit, including the period ( P out ), the longitude of the node relative to the inner binary’s node ( Ω out ), and the periastron passage times ( τ out ), were adjusted.
A couple of other parameters were constrained instead of being adjusted or held constant during our analyses, as follows:
(i)
Stars: The radii and temperatures of the three stars were calculated with the use of three-linear interpolations from the precomputed 3D (metallicity, logarithmic age, and stellar mass) PARSEC grids. Additionally, the distance of the system (which is necessary for the SED fitting) was calculated a posteriori at the end of each trial step by minimizing the value of χ SED 2 .
(ii)
The atmospheric parameters of the stars: we handled them in a similar manner as in our previous photodynamical studies. We utilized a logarithmic limb-darkening law [65] for which the passband-dependent linear and non-linear coefficients were interpolated at each trial step via the tables from the 0.xx versions of the Phoebe software, developed and maintained by A. Prša [66]. We set the gravity darkening exponents for all late type stars to β = 0.32 in accordance with the classic model of Lucy [67] valid for convective stars and held them constant. Note, however, that the choice of this parameter has only minor consequences since the stars in the present study are close to spheroids.
(iii)
Orbits: The orbital period of the inner binary ( P in ) and its orbital phase (through the time of an arbitrary primary eclipse or, more strictly, the time of the inferior conjunction of the secondary star— T in inf ) were constrained internally through the ETV curves. Finally, in the case of KIC 9714358, the systemic radial velocity ( V γ ) was constrained in a similar way, as was done with the distances, described above, i.e., this parameter was calculated a posteriori with a minimization of the value of χ RV 2 .
The median values of the orbital and physical parameters, as well as some derived quantities, of the three triple systems, computed from the MCMC posteriors and their 1 σ uncertainties, are tabulated in Table 2, Table 3 and Table 4. Furthermore, the observed vs. model lightcurves are plotted in Figure 3, Figure 4 and Figure 5, while the observed vs. model ETV curves are shown in Figure 1. In the case of KIC 9714358, we also plot the APOGEE-2 RV data vs. the best photodynamical model, both in the time and the phase domains in Figure 2.
While the majority of the tabulated parameters are self-explanatory, here we add some notes on the apsidal and nodal-motion-related parameters. We give the theoretically estimated apsidal motion periods both in the observational ( P apse ) and the dynamical ( P apse dyn ) frame of references. For their calculation, we do not exclusively use the theoretical (quadrupole) model of the point-mass third-body perturbations, but the tidal and general relativistic effects are also taken into account. In this regard, we also tabulate the individual contributions of the third-body, tidal, and relativistic effects to the apsidal advance rate ( Δ ω 3 b , tide , GR ). Moreover, besides the apsidal motion parameters, we also give the theoretical, quadrupole, nodal regression period ( P node dyn ), i.e., the time needed for a 360 regression of the parameters Ω in , out dyn . Further details of these derivations are discussed in Section 6.2 of Kostov et al. [68].

6. Discussion

In what follows, we discuss our results for the three investigated systems individually. Though we are primarily interested in the inferred dynamical behavior of our triples, first we briefly consider the astrophysical implications of our results and then discuss the dynamics of the systems in detail.

6.1. KIC 9714358

In accord with the previous results of Windemuth et al. [51], we found this triple to be very young. Our inferred age of τ 46  Myr implies that the two less massive components of the triple must still be in the pre-Main Sequence contraction stage. This conclusion is imposed by the fact that the secondary is too large, while at the same time is too cool relative to the primary. To be more specific, the surface brightness (and, hence, indirectly, the temperature ratio) and the ratio of the stellar radii are strongly constrained through the shape, durations, and depth ratios of the primary and secondary eclipses, and our Lightcurvefactory code was unable to find any reliable coeval PARSEC isochrone MS solutions for the two stars. Moreover, in the present situation the inner mass ratio ( q in = 0.406 ± 0.006 ) is also very robustly constrained dynamically both from the ETV curve (where, however, the inner mass ratio occurs only at the octupole levels of both the medium- and long-term perturbations) and from the RV curve of the primary component of the EB. Strictly speaking, as is well-known, in a single-lined spectroscopic binary (SB1) the amplitude of the RV curve directly contains only the spectroscopic mass function, which can easily be written in the following form:
f ( m Ab ) = m Aa sin 3 i q 3 ( 1 + q ) 2 .
This, coupled with the known (inner) inclination of i in = 87 . 61 ± 0 . 06 and the primary’s mass (at least to within a few percent uncertainty, as in the current situation) of  m Aa = 0.83 ± 0.03 M , uniquely gives the (inner) mass ratio q in = 0.406 ± 0.006 . Thus, this very robust mass ratio clearly excludes any post-MS solutions; hence, insofar as we assume that the triple (or, at least the two inner stars) have formed together at the same time, and without any substantial interactions (in particular, mass exchange) during their prior evolution, we can conclude that this system is at the very beginning of its life.
Comparing our results with former survey results, our SED solution more clearly supports a hotter primary component with T Aa = 5300 ± 100  K than the one listed in either the TIC v8.2 catalog ( 4764 ± 109  K, see in Table 1) or inferred from the APOGEE-2 spectra ( 4807 ± 95  K, Jönsson et al. [53]). In this context, some further caution is needed because our inferred (photometric) distance was found to be d = 666 ± 11  pc, which substantially exceeds the Gaia EDR3 value of d EDR 3 = 577 ± 7  pc [47]. One should keep in mind, however, that Gaia EDR3 parallaxes; hence, the calculated distances have not yet been corrected for multiplicity. Thus, naturally, due to this discrepancy, the astrophysical implications of our results should be considered only with some caution.
On the other hand, these discrepancies do not influence the validity of our dynamical results. Besides the robust inner and outer mass ratios (the latter being q out = 0.265 ± 0.003 ) and the eccentricities of the two orbits ( e in = 0.0286 ± 0.0003 and e out = 0.252 ± 0.002 ), here we mention the complete absence of any eclipse depth variations not only during the four years of the Kepler data but also during the extended 13-year-long interval of the Kepler and TESS observations. This latter point very strongly supports our findings about the nearly exact coplanarity of the inner and outer orbital planes, which was found to be i mut = 0 . 07 ± 0 . 03 . (It will be shown in the next subsection that even a few tenths of a degree departure from exact coplanarity may lead to the robust detection of EDVs for such a tight system).
The most interesting dynamical feature of the current system is, however, the unusual behavior of the lines of the apsides of both orbits, which, in our interpretation, is a clear manifestation of the octupole-order perturbation effects, and it has clearly observable consequences.
The theoretical observable apsidal motion period of the inner EB of KIC 9714358, according to the quadrupole-level perturbation theory, should be P apse , in theo = 26.4 ± 0.2  yr (see the middle of Table 2). As one can readily see in the upper left panel of Figure 1, the ∼5000-day-long observational dataset (which is longer than half the duration of the predicted full apsidal cycle) contradicts the quadrupole theory. In the left panel of Figure 6, we plot the variations of the observable and dynamical arguments of periastron according to our best-fitted solution, extending the numerical integration of Lightcurvefactory to the end of the current century. As one can see, the most characteristic effect is that the lines of the apsides of the inner and outer orbits (red and blue curves) revolve with the same speed, resulting in similar inner and outer apsidal motion periods of P apse meas 78  yr (which is quite close to the theoretical quadrupole outer apsidal motion period, P apse , out theo = 76.8 ± 0.6  yr). However, they do so in such a manner that, while the outer orbital ellipse rotates with a constant rate, the major axis (i.e., the apsidal line) of the inner orbit oscillates or librates around the direction of the outer apsidal line with a period of P apse osc P apse meas / 2.5 and with a half-amplitude of ∼15 . (In other words, the difference of ω in ω out oscillates around 0 ). Naturally, the same is true for the dynamical arguments of periastrons, with the difference here being that ω in dyn ω out dyn oscillates around 180 , as is shown nicely in the right panel of Figure 6.
In the right panel of Figure 6, besides the difference of ω in dyn ω out dyn , we also plot the variations of the inner eccentricity ( e in ) at the same times. As one can see, the inner eccentricity oscillates with the very same period and with a 0.25 phase offset relative to the apsidal oscillations. This is in nice accord with the octupole-order long-timescale perturbation equation for the eccentricity (see above, in Equation (18)). And, one can conclude that the currently observable uneven variations in the ETV curves of KIC 9714358, i.e., the momentarily rapidly diverging primary and secondary ETV curves, are the direct manifestations of the current growth of the eccentricity of the binary orbit due to the octupole-order terms.
At this point, however, one should note an important caveat. As was mentioned above, our analysis has shown that this triple system is quite young, and its components (or, at least, the less massive ones) are most probably in pre-MS stages. Such components may potentially exhibit strong spot activity and may rotate non-synchronously (and perhaps even in a non-aligned way). Both effects can affect the ETV curves either indirectly, through the spot-induced light curve distortions, which may mimic rapid quasi-periodic and anticorrelated variations in the primary and secondary ETV curves [69,70,71], or directly, influencing the apsidal motion rate of the inner binary [72,73]. In the current situation, however, we are convinced that the peculiar behavior of the ETVs cannot be explained with these latter effects. First of all, in all the known cases, the time-scales of the quasi-periodic spot-induced variations range from a few days to a few months (see, again, the plots in [70,71]), and, moreover, their amplitudes do not exceed 5–10 min. Moreover, as these ETVs are virtual and caused only by light curve distortions, which affect the determinations of the mid-eclipse times, these distortions would appear in the residuals of our model light curves and, hence, may become directly detectable11. On the other hand, regarding the possibility of the non-synchronous rotation, this may actually affect the apsidal motion rate and hence the apsidal motion period but will not change the eccentricity of the inner pair. It is the latter that is the main source of that bump in the ETVs, which now looks evident in the TESS measured eclipse times of KIC 9714358.
Here, we note that KIC 9714358 is scheduled to be reobserved by TESS in Sectors 74, 75 (January–February 2024) and 81 (July/August 2024), when the offset between the primary and secondary eclipses is expected to be near the next maximum (see upper right panel of Figure 1). Thus, these findings and the corresponding predictions of our model will be verifiable very soon.

6.2. KIC 5771589

As mentioned above in Section 3, the two noteworthy observables of this triple system are the very rapid apsidal motion and the very quickly varying eclipse depths. Despite our best efforts, we were unable to model the variations of the eclipse depths perfectly. In general, as one can see in Figure 3, Lightcurvefactory returns the main properties of the EDVs in the case of both the Kepler and the TESS data. The problem, however, shows up in the second half of the Kepler data where our models fail to reproduce the correct depth ratios of the secondary vs. primary eclipses. In other words, while the Kepler observations suggest that the ratio of the depths of the secondary and primary eclipses remains near constant during the four years of the observations, our model fits tend to equalize the depths of the two kinds of eclipses toward the end of the Kepler observations.
To understand the origin of this problem, one should keep in mind that for an eccentric ( e 0 ) EB viewed not exactly edge-on ( i 90 ), the ratio of the eclipse depths, in addition to being a function of the ratio of the surface brightnesses, also depends strongly upon the position of the observable argument of periastron. This is so because, for the eclipse that occurs closer to periastron, the projected separation of the two stellar disks will be smaller than for the other kind of eclipse (closer to the apastron). Thus, a larger portion of the surface of the eclipsed star will be occulted. Naturally, the more eccentric the orbit, and the more grazing the eclipse, the more significant this effect may be. And, naturally, one can reach a situation where the eclipse closer to apastron disappears entirely, as was observed recently, e.g., in the case of KIC 5731312 [3]. Note, we graphically illustrate these effects for a hypothetical eccentric EB in Figure 7.
In the case of the current system, both the position of the argument of periastron and the inclination have varied during the four years of Kepler observations. As we mentioned above, which can also be seen directly from the ETV curve (middle left panel of Figure 1), the inner orbital ellipse made slightly more than half a revolution during the Kepler observations. According to our photodynamical model, the observable argument of periastron attained a value of ( ω = 270 ) around the very beginning of the Year 2010 (circa 2,455,200; see also the left panel of Figure 8). This means that during that time, the primary eclipses occurred during the periastron passages, while the secondaries occurred at the apastron points, thereby causing the largest differences in the sizes of the eclipsed areas (bottom row of Figure 7). Then, toward the end of the Kepler era, the situation reversed. During the last days of the perfectly operating three gyroscopes on the Kepler spacecraft, the observable argument of periastron reached ω = 90 , i.e., the last Kepler-observed primary eclipses occurred around the apastron points, while the secondaries were near periastron (second row of Figure 7). This is the reason why the secondary eclipse depths tend to be much closer to the primary ones by the end of the Kepler data in our photodynamical models. We note, however, that the same argument can be made purely from the observations, without the use of our photodynamical model results; hence, one cannot claim this discrepancy as a failure of the model. This is so because the ETV curves show in themselves that the system was seen from the direction of the major axis at the very beginning and the end of the Kepler observations. (This can be deduced from the fact that the primary and secondary curves intersect each other only when the orbital ellipse is seen from the direction of the major axis, i.e., when one of the eclipses occurs at periastron and the other at apastron.) Then, the fact, that the secondary (blue) ETV curve is located below the primary (red) one during the entire Kepler observations reveals that during this interval, the secondary eclipses were “in a hurry” in contrast to the primary ones (or, in other words, the secondary eclipses occurred before photometric phase of 0.5). This means that the periastron passage happened between the primary and the secondary eclipses (i.e., ω had values between 270 [or, 90 ] and 90 , reaching ω = 0 at the mid-time—see the top row of Figure 7). But, even if one assumed retrograde apsidal motion, which would mean that the roles of the apastrons and periastrons were exchanged, this would not explain why the observed ratios of the secondary to primary eclipse depths would have immunity to the apsidal revolution.
In our understanding, the only possibility for explaining the near constant depth ratios during half an apsidal cycle is that the true observable inclination of the system should be closer to i = 90 than our findings of i in = 86 . 05 ± 0 . 09 would suggest (since the more edge-on an eccentric EB is viewed, the smaller the apsidal-phase-dependent difference of the eclipsed disk areas, and, for orbits viewed perfectly edge-on, this effect entirely disappears). This fact, however, would require more contaminating light in the system. The main possibilities to fulfill this requirement are as follows: (i) the system has a fourth, more distant, undetected component; (ii) the tertiary is more luminous than would be inferred from the PARSEC stellar isochrones that we utilized or; (iii) the outer mass ratio ( q out ) could be larger; hence, the tertiary would be more massive (and, consequently brighter) relative to the binary members than our solution suggests. Regarding this last possibility, we made efforts to find an acceptable solution with higher q out , but all of our trials failed. So, in what follows, we discuss our findings in the imperfect case, where the (non-)variation of the eclipse depth ratios is imperfectly modeled. Despite this fact, however, we believe that our findings, to be discussed below, are basically correct.
According to our joint photodynamical solution, in KIC 5771589 the most massive star is the third, distant component, with m B = 1.02 ± 0.01 M , though the primary of the inner binary is only slightly less massive ( m Aa = 0.94 ± 0.01 M ), while its binary companion has a lower mass of m Ab = 0.73 ± 0.01 M . Here, we note that, despite the lack of RV data, our solutions yield masses with surprisingly small uncertainties of only 1–2%. Due to the problematic nature of the light curve solution discussed above, however, we are not convinced that these small uncertainties are realistic. The quantity that is certainly more robust is the outer mass ratio, being q out = 0.612 ± 0.003 . This value is substantially lower than what was inferred from the previous analytic ETV fits of Borkovits et al. [11,12]. One should keep in mind, however, that those former fits were based only on the Kepler data (long before the launch of the TESS space telescope) and, because of the extremely tight nature of the current system, as the authors noted, the analytical ETV models they used were somewhat inadequate.
Our solution suggests a metal-deficient ( [ M / H ] 0.50 ), old ( log τ = 9.83 ± 0.02 ), and slightly evolved system. According to the accepted solution, the massive tertiary is clearly on its way toward the giant branch, having log g B = 3.7 ± 0.1 . One should keep in mind, however, that in the absence of third-body eclipses, the direct effects of the tertiary on the joint, complex solutions (apart from the SED fitting) manifest themselves only in the outer mass ratio (via timings) and the contaminated light (via the eclipse depths). Hence, any further statements about the astrophysical condition of the tertiary depend strongly upon the pre-assumption that the whole system formed and evolved in a coeval manner. Our solution gives a distance of d = 709 55 + 33 pc, which differs substantially from the Gaia EDR3-based value of d EDR 3 = 870 ± 100 pc. This fact might again suggest that the system could have an additional, fourth stellar component (not considered in the photodynamical solution), and the high value of the RUWE parameter (11.34) indeed indicates that the Gaia astrometric solution is probably significantly influenced by the multiplicity of the system.
Turning to the readily observable long-term or secular perturbations in this triple, the rapid, dynamically forced apsidal motion is evident. The doubly averaged secular theory predicts apsidal motion periods of the inner and outer orbits to be ( P apse = 11.22 ± 0.04 and 51.7 ± 0.2 yrs, respectively). From the fact that the ∼4-yr-long Kepler dataset covers slightly more than half of an apsidal cycle, one can infer directly that the quadrupole model, again, gives an incorrect result for this system. This can also be seen in the left panel of Figure 8, which shows that the true apsidal motion period is about P apse meas 6.3 yr. On the other hand, however, note that the outer apsidal motion period according to our numerical integration is P apse meas 58.1 yr, which is, again, in much better agreement with the theoretical quadrupole value.
The significance of the octupole-order perturbations in the variations of the inner eccentricity is, again, nicely demonstrated in the right panel of Figure 8. As one can see, the cyclic variations of e in between ∼0.01 and ∼0.02 have exactly the same period as for the quantity ω in dyn ω out dyn , which, in the current system, varies between 0 and 360 with a period of ∼8 yrs.
The most interesting feature of KIC 5771589, however, is its continuous eclipse depth variations. What may be surprising, therefore, at first sight is the fact that according to our results the inner and outer orbital planes are almost coplanar. The non-zero mutual inclination is the origin of the nodal precession; hence, the eclipse depth variations are very small, being only i mut = 0 . 29 ± 0 . 03 . The question naturally arises as to how is it possible that such a small mutual inclination can cause such a readily observable eclipse depth variation? The answer lies in the extreme grazing nature of the eclipses, as is shown in the left panel of Figure 9.
Finally, in Figure 10 we plot the predicted eclipse depth variations for the present century. As one can see, for the expected near-future observations of TESS (in Sectors 74, 80, 81) the inner inclination will be larger; hence, one can predict somewhat deeper and better detectable eclipses. The numerical integrations give a nodal regression period of P node meas 9.2 yr, which is in perfect accord with the theoretical value tabulated in Table 3.

6.3. TIC 219885468

Our third triple is a good counterexample in regard to both (i) the octupole perturbations driving apsidal motion and eccentricity variations and (ii) the remarkable eclipse depth variations in an almost flat system.
Both the tightness and the compactness of this system is between our previous two example triples. Moreover, similar to the other two triples, the inner orbit has small eccentricity ( e in = 0.0423 ± 0.0006 ), while the outer orbit is moderately eccentric ( e out = 0.3903 ± 0.0007 ). The outer mass ratio, again, makes this triple quite similar to KIC 9714358 ( q out = 0.319 ± 0.002 vs. q out = 0.265 ± 0.003 for the TIC and the KIC systems, respectively), while the mutual inclination is almost identical to that of the EDV system KIC 5771589 ( i mut = 0 . 25 ± 0 . 14 vs. i mut = 0 . 29 ± 0 . 03 , respectively). The similarities, however, end at this point. From a dynamical point of view, the main difference between the present and the two KIC systems is that the inner EB of TIC 219885468 consists of two twin stars ( q in = 0.988 ± 0.017 ); hence, the octupole (and any other even order) perturbation terms can be neglected. This is demonstrated nicely in the two panels of Figure 11, where it can be readily seen that both the inner and outer ellipses rotate with constant rates, as is expected for coplanar systems in the quadrupole theory. Moreover, the theoretically calculated quadrupole apsidal advance rates ( P apse , in theo = 19.49 ± 0.09 yr and P apse , out theo = 51.2 ± 0.1 yr) agree quite well with the apsidal motion periods ‘measured’ via numerical integration ( P apse , in meas 13.6 yr and P apse , out meas 54.5 yr).
The agreement is even better in the nodal regression period, being P node meas 13.9 yr vs. P node theo = 14.09 ± 0.07 yr. Turning to this latter effect, despite the similar amplitude of the precession cone, in the case of TIC 219885468 one cannot detect any EDVs. The two panels of Figure 9 nicely illustrate the substantial differences amongst the eclipse geometries of the two EBs. While in the case of KIC 5771578 there are grazing eclipses, in TIC 219885468 the eclipses are much deeper. Hence, the areas of the eclipsed surfaces are less sensitive to such small inclination variations.
Regarding the astrophysical parameters of the three stars, the twin EB members are found to be late F/early G-type, slightly evolved MS stars ( log τ = 9.49 ± 0.06 ) with masses of m Aa = 1.19 ± 0.03 M and m Ab = 1.18 ± 0.04 M and effective temperatures of T Aa = 6 410 ± 40 K and T Ab = 6 400 ± 40 K. The less massive third component is a K dwarf with m B = 0.75 ± 0.02 M and T B = 4 885 ± 50 K. The inferred distance of the triple star agrees with the Gaia EDR3 distance within two sigma ( d = 1 150 ± 20 pc vs. d EDR 3 = 1 111 ± 13 pc).

7. Conclusions

We investigated three such EBs that form the inner pairs of TCHT stellar systems. Two systems (KICs 9714358 and 5771589) were observed with the Kepler space telescope in its 4-yr-long prime mission, while the third target (TIC 219885468) is located in the NCVZ of TESS. Besides the evident medium-period third-body perturbations, the ETVs of all three EBs show rapid dynamically driven apsidal motion, and one of the three targets (KIC 5771589) also exhibits cyclic EDVs. We carried out complex photodynamical analyses for the three systems, jointly analyzing their Kepler and TESS light curves (naturally, for TIC 219885468 only the TESS light curve was used); ETV curves; composite SEDs; and, in the case of KIC 9714358, APOGEE-2 RV data as well. We determined reliable, robust astrophysical and orbital parameters for all three systems and their constituent stars.
In our study, however, we focused mainly on the dynamical properties and higher-order third-body perturbations of the systems. We found clear evidence that in the case of the almost perfectly coplanar KIC 9714358, the inner and outer orbital ellipses are oriented in the same directions, and while, on average, they rotate evenly with the same apsidal revolution rate (with a period of P apse meas 78 yr), the axis of the inner orbit librates around this average direction with an oscillation period of P apse , in osc 30 yr and with an amplitude of ∼15 . The variations of the inner eccentricity of KIC 9714358 show the same period as the apsidal oscillations but with a phase shift of ∼0.25 with respect to the libration. We connect this behavior to the octupole perturbations; however, we leave the quantitative investigation of this issue for a future study. We expect, however, that the near-future TESS observations will provide us with additional observational evidence for this behavior.
Regarding KIC 5771589, despite the evident cyclic EDVs of this system (detected both in the Kepler and TESS eras), we found, somewhat surprisingly, that the triple was almost coplanar, with a mutual inclination of i mut = 0 . 29 ± 0 . 03 . We explain this unusual finding with the grazing nature of the eclipses as the depths of grazing eclipses are much more sensitive to even very small variations in the observable inclination. This triple is the second tightest Kepler-triple, and this EB has the second shortest apsidal motion period. Despite this fact, we did not find observational evidence for the operation of the octupole-order perturbations; however, our numerical integration revealed that they are present. The weakness of the octupole effects in the present system could be in accord with the fact that the relative strength of the octupole-order effects relative to the quadrupole ones is much lower than in KIC 9714358.
Finally, TIC 219885468 was selected as a counterexample to the two KIC sources. This triple, in both compactness and tightness, is very similar to the other two systems, but the nearly equal primary and secondary eclipse depths suggested that the two stars of the inner EB should be similar in mass; hence, one cannot expect octupole-order third-body perturbations in this system. While our analysis justified this pre-assumption, finding q in = 0.988 ± 0.017 revealed another aspect in which TIC 219885468 serves as counterexample. The mutual inclination of the system was found to be almost identical to that of the former system ( i mut = 0 . 25 ± 0 . 14 vs. 0 . 29 ± 0 . 03 ). In the current system, however, the eclipses are deep; hence, they are insensitive to such small variations in inclination.
As our primary aim was to observationally detect the higher-order secular perturbations, we did not investigate the past and future evolution of these three triple systems as this question is beyond the focus of the present paper. Currently, all three systems are stable in the sense of the semi-empirical dynamical stability limit of Mardling & Aarseth [6]. This fact, however, does not imply that the systems also remain stable in the distant future. Such investigations may also be part of a future study.

Author Contributions

Conceptualization, T.B.; methodology, T.B.; software, T.B.; validation, T.B. and T.M.; formal analysis, T.B.; investigation, T.B. and T.M.; resources, T.B. and T.M.; writing—original draft preparation, T.B.; writing—review and editing, T.B. and T.M.; visualization, T.B. All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the HUN-REN Hungarian Research Network.

Data Availability Statement

The Kepler and TESS observations were provided by NASA and obtained either from the Kepler Eclipsing Binary Catalog at http://keplerebs.villanova.edu/ (last accessed on 7 November 2023) or the Mikulski Archive for Space Telescopes (MAST) at https://mast.stsci.edu/portal/Mashup/Clients/Mast/Portal.html (last accessed on 7 November 2023). The RV data used for the analysis of KIC 9714358 were provided by the DR 16 of the APOGEE-2 project, which was downloaded from the Vizier site at http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=III/284 (last accessed on 7 November 2023).

Acknowledgments

The authors are very grateful to S. A. Rappaport for his extensive linguistic corrections. This paper includes data collected by the Kepler and TESS missions that were obtained from the MAST data archive at the Space Telescope Science Institute (STScI). Funding for the Kepler mission is provided by the NASA Science Mission Directorate. Funding for the TESS mission is provided by the NASA Explorer Program. STScI is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS5–26555. We used the Simbad service operated by the Centre des Données Stellaires (Strasbourg, France) and the ESO Science Archive Facility services (data obtained under request number 396301). This research has made use of NASA’s Astrophysics Data System.

Conflicts of Interest

The authors declare no conflict of interest.

Notes

1
In what follows, for simplicity, we will talk about only hierarchical triple star systems, but our discussion can be extended very simply to hierarchical multiple systems as well.
2
Strictly speaking, Earth-based timing observations are less effective in the detection of any short outer period (i.e., compact) triple star systems irrespective of their tightness. This question was discussed in detail in the recent review of Borkovits [10].
3
The latter two classes of medium-period and long-period perturbations in Brown’s original terminology were named ‘long-period’ and ‘apse-node type perturbations’, respectively, but this would have resulted in some confusion with the terminology of the classical planetary perturbation theories, which was avoided with the use of this modified terminology.
4
Note, that, in general, for equal mass inner stars all the higher, even order perturbative terms disappear as well.
5
Note that in the case of the most compact (but less tight) Kepler triple, KOI-126, the full apsidal motion period is shorter than the length of the Kepler dataset [54,55], but in that case no precise eclipse times and, hence, no ETVs are available due to the very shallow and frequently missing eclipses.
6
Note, this does not mean that ‘at the other end’ of the precession cycle (i.e., a half precession period later) it cannot reach another extremum that is even farther from the edge-on view.
7
http://keplerebs.villanova.edu/ (accessed on 7 November 2023).
8
9
For Kepler data, naturally, we did not recalculate the times of minima, as was done already for the Borkovits et al. [11,12] papers.
10
11
At this point, naturally, we are talking about the residuals of the eclipsing sections of the observed vs. model light curves as the majority of the non-eclipse light curve sections were dropped from our analysis.

References

  1. Borucki, W.J.; Koch, D.; Basri, G.; Batalha, N.; Brown, T.; Caldwell, D.; Caldwell, J.; Christensen-Dalsgaard, J.; Cochran, W.D.; DeVore, E.; et al. Kepler Planet-Detection Mission: Introduction and First Results. Science 2010, 327, 977. [Google Scholar] [CrossRef] [PubMed]
  2. Ricker, G.R.; Winn, J.N.; Vanderspek, R.; Latham, D.W.; Bakos, G.Á.; Bean, J.L.; Berta-Thompson, Z.K.; Brown, T.M.; Buchhave, L.; Butler, N.R.; et al. Transiting Exoplanet Survey Satellite (TESS). J. Astron. Telesc. Instrum. Syst. 2015, 1, 014003. [Google Scholar] [CrossRef]
  3. Borkovits, T.; Rappaport, S.A.; Toonen, S.; Moe, M.; Mitnyan, T.; Csányi, I. Von Zeipel–Lidov–Kozai cycles in action: Kepler triples with eclipse depth variations: KICs 6964043, 5653126, 5731312, 8023317. Mon. Not. R. Astron. Soc. 2022, 515, 3773–3795. [Google Scholar] [CrossRef]
  4. Harrington, R.S. Stability Criteria for Triple Stars. Celest. Mech. 1972, 6, 322–327. [Google Scholar] [CrossRef]
  5. Zare, K. Bifurcation Points in the Planar Problem of Three Bodies. Celest. Mech. 1977, 16, 35–38. [Google Scholar] [CrossRef]
  6. Mardling, R.A.; Aarseth, S.J. Tidal interactions in star cluster simulations. Mon. Not. R. Astron. Soc. 2001, 321, 398–420. [Google Scholar] [CrossRef]
  7. Tokovinin, A. Architecture of Hierarchical Stellar Systems and Their Formation. Universe 2021, 7, 352. [Google Scholar] [CrossRef]
  8. Offner, S.S.R.; Moe, M.; Kratter, K.M.; Sadavoy, S.I.; Jensen, E.L.N.; Tobin, J.J. The Origin and Evolution of Multiple Star Systems. In Protostars Planets VII; Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K., Tamura, M., Eds.; Astronomical Society of the Pacific: San Francisco, CA, USA, 2023; Volume 534, 275p. [Google Scholar]
  9. Toonen, S.; Portegies Zwart, S.; Hamers, A.S.; Bandopadhyay, D. The evolution of stellar triples. The most common evolutionary pathways. Astron. Astrophys. 2020, 640, A16. [Google Scholar] [CrossRef]
  10. Borkovits, T. Eclipsing Binaries in Dynamically Interacting Close, Multiple Systems. Galaxies 2022, 10, 9. [Google Scholar] [CrossRef]
  11. Borkovits, T.; Rappaport, S.; Hajdu, T.; Sztakovics, J. Eclipse timing variation analyses of eccentric binaries with close tertiaries in the Kepler field. Mon. Not. R. Astron. Soc. 2015, 448, 946–993. [Google Scholar] [CrossRef]
  12. Borkovits, T.; Hajdu, T.; Sztakovics, J.; Rappaport, S.; Levine, A.; Bíró, I.B.; Klagyivik, P. A comprehensive study of the Kepler triples via eclipse timing. Mon. Not. R. Astron. Soc. 2016, 455, 4136–4165. [Google Scholar] [CrossRef]
  13. Brown, E.W. The stellar problem of three bodies. I. Application of satellite theory. Mon. Not. R. Astron. Soc. 1936, 97, 56–61. [Google Scholar] [CrossRef]
  14. Brown, E.W. The stellar problem of three bodies. II. The equations of motion with a simplified solution. Mon. Not. R. Astron. Soc. 1936, 97, 62–65. [Google Scholar] [CrossRef]
  15. Brown, E.W. The stellar problem of three bodies. III. The motions of the apse and node with applications to the Moon. Mon. Not. R. Astron. Soc. 1936, 97, 116–126. [Google Scholar] [CrossRef]
  16. Von Zeipel, H. Sur l’application des séries de M. Lindstedt à l’étude du mouvement des comètes périodiques. Astron. Nachr. 1910, 183, 345. [Google Scholar] [CrossRef]
  17. Kozai, Y. Secular perturbations of asteroids with high inclination and eccentricity. Astron. J. 1962, 67, 591–598. [Google Scholar] [CrossRef]
  18. Lidov, M.L. The evolution of orbits of artificial satellites of planets under the action of gravitational perturbations of external bodies. Planet. Space Sci. 1962, 9, 719–759. [Google Scholar] [CrossRef]
  19. Slavenas, P. The Stellar Case of the Problem of Three Bodies. Trans. Astron. Obs. Yale Univ. 1927, 6, 33. [Google Scholar]
  20. Slavenas, P. A Note on the Triple System, λ Tauri. Proc. Nat. Acad. Sci. USA 1927, 13, 833–837. [Google Scholar] [CrossRef]
  21. Harrington, R.S. Dynamical Evolution of Triple Stars. Astron J. 1968, 73, 190–194. [Google Scholar] [CrossRef]
  22. Harrington, R.S. The Stellar Three-Body Problem. Cel. Mech. 1969, 1, 200–209. [Google Scholar] [CrossRef]
  23. Söderhjelm, S. The three-body problem and eclipsing binaries. Application to Algol and lambda Tauri. Astron. Astroph. 1975, 42, 229–236. [Google Scholar]
  24. Söderhjelm, S. Studies of the stellar three-body problem. Astron. Astroph. 1982, 107, 54–60. [Google Scholar]
  25. Söderhjelm, S. Third-order and tidal effects in the stellar three-body problem. Astron. Astroph. 1984, 141, 232–240. [Google Scholar]
  26. Mazeh, T.; Shaham, J. The orbital evolution of close triple systems: The binary eccentricity. Astron. Astroph. 1979, 77, 145–151. [Google Scholar]
  27. Beust, H.; Corporon, P.; Siess, L.; Forestini, M.; Lagrange, A.-M. Dynamics and evolutionary status of the young triple stellar system, TY Coronae Austinae. Astron. Astroph. 1997, 320, 478–496. [Google Scholar]
  28. Innanen, K.A.; Zheng, J.Q.; Mikkola, S.; Valtonen, M.J. The Kozai Mechanism and the Stability of Planetary Orbits in Binary Star Systems. Astron. J. 1997, 113, 1915–1919. [Google Scholar] [CrossRef]
  29. Mikkola, S.; Tanikawa, K. Does Kozai Resonance Drive CH Cygni? Astron. J. 1998, 116, 444–450. [Google Scholar] [CrossRef]
  30. Kiseleva, L.G.; Eggleton, P.P.; Mikkola, S. Tidal friction in triple stars. Mon. Not. R. Astron. Soc. 1998, 300, 292–302. [Google Scholar] [CrossRef]
  31. Von Zeipel, H. Recherches sur le mouvement des petites planétes. Premiére partie. Ark. Mat. Astron. Fys. 1916, 11, 1–58. [Google Scholar]
  32. Von Zeipel, H. Recherches sur le mouvement des petites planétes. Deuxiéme partie. Ark. Mat. Astron. Fys. 1916, 11, 59–62. [Google Scholar]
  33. Von Zeipel, H. Recherches sur le mouvement des petites planétes. Troisiéme partie. Ark. Mat. Astron. Fys. 1917, 12, 1–89. [Google Scholar]
  34. Von Zeipel, H. Recherches sur le mouvement des petites planétes. Quatriéme partie. Ark. Mat. Astron. Fys. 1917, 13, 1–93. [Google Scholar]
  35. Ford, E.B.; Kozinsky, B.; Rasio, F.A. Secular Evolution of Hierarchical Triple Star Systems. Astrophys. J. 2000, 535, 385–401. [Google Scholar] [CrossRef]
  36. Naoz, S.; Farr, W.M.; Lithwick, Y.; Rasio, F.A.; Teyssandier, J. Secular dynamics in hierarchical three-body systems. Mon. Not. R. Astron. Soc. 2013, 431, 2155–2171. [Google Scholar] [CrossRef]
  37. Naoz, S. The Eccentric Kozai-Lidov Effect and Its Applications. Ann. Rev. Astron. Astrophys. 2016, 54, 441–489. [Google Scholar] [CrossRef]
  38. Krymolowski, Y.; Mazeh, T. Studies of multiple stellar systems—II. Second-order averaged Hamiltonian to follow long-term orbital modulations of hierarchical triple systems. Mon. Not. R. Astron. Soc. 1999, 304, 720–732. [Google Scholar] [CrossRef]
  39. Stassun, K.G.; Torres, G. Parallax Systematics and Photocenter Motions of Benchmark Eclipsing Binaries in Gaia EDR3. Astrophys. J. Lett. 2021, 907, L33. [Google Scholar] [CrossRef]
  40. Paegert, M.; Stassun, K.G.; Collins, K.A.; Pepper, J.; Torres, G.; Jenkins, J.; Twicken, J.D.; Latham, D.W. TESS Input Catalog versions 8.1 and 8.2: Phantoms in the 8.0 Catalog and How to Handle Them. arXiv 2021, arXiv:2108.04778. [Google Scholar]
  41. Brown, A.G.A. et al. [Gaia Collaboration] Gaia Early Data Release 3. Summary of the contents and survey properties. Astron. Astroph. 2021, 649, A1. [Google Scholar]
  42. Chambers, K.C.; Magnier, E.A.; Metcalfe, N.; Flewelling, H.A.; Huber, M.E.; Waters, C.Z.; Denneau, L.; Draper, P.W.; Farrow, D.; Finkbeiner, D.P.; et al. The Pan-STARRS1 Surveys. arXiv 2016, arXiv:1612.05560. [Google Scholar]
  43. Greiss, S.; Steeghs, D.; Gänsicke, B.T.; Martín, E.L.; Groot, P.J.; Irwin, M.J.; González-Solares, E.; Greimel, R.; Knigge, C.; Østensen, R.H.; et al. Initial Data Release of the Kepler-INT Survey. Astron. J. 2012, 144, 24. [Google Scholar] [CrossRef]
  44. Henden, A.A.; Levine, S.; Terrell, D.; Welch, D. APASS—The Latest Data Release; American Astronomical Society: Washington, DC, USA, 2015; p. 336.16. [Google Scholar]
  45. Skrutskie, M.F.; Cutri, R.M.; Stiening, R.; Weinberg, M.D.; Schneider, S.; Carpenter, J.M.; Beichman, C.; Capps, R.; Chester, T.; Elias, J.; et al. The Two Micron All Sky Survey (2MASS). Astron. J. 2006, 131, 1163–1183. [Google Scholar] [CrossRef]
  46. Cutri, R.M.; Wright, E.L.; Conrow, T.; Fowler, J.W.; Eisenhardt, P.R.M.; Grillmair, C.; Kirkpatrick, J.D.; Masci, F.; McCallon, H.L.; Wheelock, S.L.; et al. Explanatory Supplement to the AllWISE Data Release Products. 2013, wise.rept, 1. Available online: https://cdsarc.cds.unistra.fr/viz-bin/cat/II/328 (accessed on 7 November 2023).
  47. Bailer-Jones, C.A.L.; Rybizki, J.; Fouesneau, M.; Demleitner, M.; Andrae, R. Estimating Distances from Parallaxes. V. Geometric and Photogeometric Distances to 1.47 Billion Stars in Gaia Early Data Release 3. Astron. J. 2021, 161, 147. [Google Scholar] [CrossRef]
  48. Prša, A.; Batalha, N.; Slawson, R.W.; Doyle, L.R.; Welsh, W.F.; Orosz, J.A.; Seager, S.; Rucker, M.; Mjaseth, K.; Engle, S.G.; et al. Kepler Eclipsing Binary Stars. I. Catalog and Principal Characterization of 1879 Eclipsing Binaries in the First Data Release. Astron. J. 2011, 141, 83. [Google Scholar] [CrossRef]
  49. Rappaport, S.; Deck, K.; Levine, A.; Borkovits, T.; Carter, J.; El Mellah, I.; Sanchis-Ojeda, R.; Kalomeni, B. Triple-star Candidates among the Kepler Binaries. Astrophys. J. 2013, 768, 33. [Google Scholar] [CrossRef]
  50. Kjurkchieva, D.; Vasileva, D.; Atanasova, T. Orbital Parameters of the Eclipsing Detached Kepler Binaries with Eccentric Orbits. Astron. J. 2017, 154, 105. [Google Scholar] [CrossRef]
  51. Windemuth, D.; Agol, E.; Ali, A.; Kiefer, F. Modelling Kepler eclipsing binaries: Homogeneous inference of orbital and stellar properties. Mon. Not. R. Astron. Soc. 2019, 489, 1644–1666. [Google Scholar] [CrossRef]
  52. Bressan, A.; Marigo, P.; Girardi, L.; Salasnich, B.; Cero, C.D.; Rubele, S.; Nanni, A. PARSEC: Stellar tracks and isochrones with the PAdova and TRieste Stellar Evolution Code. Mon. Not. R. Astron. Soc. 2012, 427, 127–145. [Google Scholar] [CrossRef]
  53. Jönsson, H.; Holtzman, J.A.; Prieto, C.A.; Cunha, K.; García-Hernández, D.A.; Hasselquist, S.; Masseron, T.; Osorio, Y.; Shetrone, M.; Smith, V.; et al. APOGEE Data and Spectral Analysis from SDSS Data Release 16: Seven Years of Observations Including First Results from APOGEE-South. Astron. J. 2020, 160, 120. [Google Scholar] [CrossRef]
  54. Carter, J.A.; Fabrycky, D.C.; Ragozzine, D.; Holman, M.J.; Quinn, S.N.; Latham, D.W.; Buchhave, L.A.; Van Cleve, J.; Cochran, W.D.; Cote, M.T.; et al. KOI-126: A Triply Eclipsing Hierarchical Triple with Two Low-Mass Stars. Science 2011, 331, 562. [Google Scholar] [CrossRef]
  55. Yenawine, M.E.; Welsh, W.F.; Orosz, J.A.; Bieryla, A.; Cochran, W.D.; Endl, M.; Latham, D.W.; Quinn, S.N.; Short, D.R.; Windmiller, G. Photodynamical Modeling of the Fascinating Eclipses in the Triple-star System KOI-126. Astrophys. J. 2022, 924, 66. [Google Scholar] [CrossRef]
  56. Kirk, B.; Conroy, K.; Prša, A.; Abdul-Masih, M.; Kochoska, A.; Matijevič, G.; Hambleton, K.; Barclay, T.; Bloemen, S.; Boyajian, T.; et al. Kepler Eclipsing Binary Stars. VII. The Catalog of Eclipsing Binaries Found in the Entire Kepler Data Set. Astron. J. 2016, 151, 68. [Google Scholar] [CrossRef]
  57. Abdul-Masih, M.; Prša, A.; Conroy, K.; Bloemen, S.; Boyajian, T.; Doyle, L.R.; Johnston, C.; Kostov, V.; Latham, D.W.; Matijevič, G.; et al. Kepler Eclipsing Binary Stars. VIII. Identification of False Positive Eclipsing Binaries and Re-extraction of New Light Curves. Astron. J. 2016, 151, 101. [Google Scholar] [CrossRef]
  58. Pál, A. FITSH—A software package for image processing. Mon. Not. R. Astron. Soc. 2012, 421, 1825–1837. [Google Scholar] [CrossRef]
  59. Borkovits, T.; Rappaport, S.; Kaye, T.; Isaacson, H.; Vanderburg, A.; Howard, A.W.; Kristiansen, M.H.; Omohundro, M.R.; Schwengeler, H.M.; Terentev, I.A.; et al. Photodynamical analysis of the triply eclipsing hierarchical triple system EPIC 249432662. Mon. Not. R. Astron. Soc. 2019, 483, 1934–1951. [Google Scholar] [CrossRef]
  60. Borkovits, T.; Rappaport, S.A.; Tan, T.G.; Gagliano, R.; Jacobs, T.; Huang, X.; Mitnyan, T.; Hambsch, F.-J.; Kaye, T.; Maxted, P.F.L.; et al. The compact triply eclipsing triple star TIC 209409435 discovered with TESS. Mon. Not. R. Astron. Soc. 2020, 496, 4624–4636. [Google Scholar] [CrossRef]
  61. Ford, E.B. Quantifying the Uncertainty in the Orbits of Extrasolar Planets. Astron. J. 2005, 129, 1706–1717. [Google Scholar] [CrossRef]
  62. Borkovits, T.; Sperauskas, J.; Tokovinin, A.; Latham, D.W.; Csányi, I.; Hajdu, T.; Molnár, L. The compact multiple system HIP 41431. Mon. Not. R. Astron. Soc. 2019, 487, 4631–4647. [Google Scholar] [CrossRef]
  63. Borkovits, T.; Rappaport, S.A.; Hajdu, T.; Maxted, P.F.L.; Pál, A.; Forgács-Dajka, E.; Klagyivik, P.; Mitnyan, T. TICs 167692429 and 220397947: The first compact hierarchical triple stars discovered with TESS. Mon. Not. R. Astron. Soc. 2020, 493, 5005–5023. [Google Scholar] [CrossRef]
  64. Borkovits, T.; Mitnyan, T.; Rappaport, S.A.; Pribulla, T.; Powell, B.P.; Kostov, V.B.; Bíró, I.B.; Csányi, I.; Garai, Z.; Gary, B.L.; et al. Triply eclipsing triple stars in the northern TESS fields: TICs 193993801, 388459317 and 52041148. Mon. Not. R. Astron. Soc. 2022, 510, 1352–1374. [Google Scholar] [CrossRef]
  65. Klinglesmith, D.A.; Sobieski, S. Nonlinear Limb Darkening for Early-Type Stars. Astron. J. 1970, 75, 175–182. [Google Scholar] [CrossRef]
  66. Prša, A.; Zwitter, T. A Computational Guide to Physics of Eclipsing Binaries. I. Demonstrations and Perspectives. Astroph. J. 2005, 628, 426–438. [Google Scholar] [CrossRef]
  67. Lucy, L.B. Gravity-Darkening for Stars with Convective Envelopes. Z. Astrophys. 1967, 65, 89–92. [Google Scholar]
  68. Kostov, V.B.; Powell, B.P.; Torres, G.; Borkovits, T.; Rappaport, S.A.; Tokovinin, A.; Zasche, P.; Anderson, D.; Barclay, T.; Berlind, P.; et al. TIC 454140642: A Compact, Coplanar, Quadruple-lined Quadruple Star System Consisting of Two Eclipsing Binaries. Astrophys. J. 2021, 917, 93. [Google Scholar] [CrossRef]
  69. Kalimeris, A.; Rovithis-Livaniou, H.; Rovithis, P. Starspots and photometric noise on observed minus calculated (O-C) diagrams. Astron. Astrophys. 2002, 387, 969–976. [Google Scholar] [CrossRef]
  70. Tran, K.; Levine, A.; Rappaport, S.; Borkovits, T.; Csizmadia, S.; Kalomeni, B. The Anticorrelated Nature of the Primary and Secondary Eclipse Timing Variations for the Kepler Contact Binaries. Astrophys. J. 2013, 774, 81. [Google Scholar] [CrossRef]
  71. Balaji, B.; Croll, B.; Levine, A.M.; Rappaport, S. Tracking the stellar longitudes of starspots in short-period Kepler binaries. Mon. Not. R. Astron. Soc. 2015, 448, 429–444. [Google Scholar] [CrossRef]
  72. Shakura, N.I. On the Apsidal Motion in Binary Stars. Sov. Astr. Lett. 1985, 11, 224–226. [Google Scholar]
  73. Hegedüs, T.; Nuspl, J. On the apsidal motion and its photometric observation in certain close binary systems. Acta Astron. 1986, 36, 381–393. [Google Scholar]
Figure 1. Observed and modeled ETVs of KICs 9714358 (upper), 5771578 (middle), and TIC 219885468 (bottom). The left-hand panels cover the time interval approximately from the beginning of the Kepler observations to the end of the second northern ecliptic hemisphere scan of TESS, while the right-hand panels show longer, 100 yr-long intervals, which offer direct predictions for future observations. Larger red and blue symbols represent primary and secondary ETV points determined from the observations, while the smaller points connected with straight lines stand for ETV points calculated according to the best-fitting complex photodynamical models. Finally, the black, almost horizontal lines represent the LTTE contributions to the curves. Note that in the middle right panel we downshifted the originally overlapping secondary ETV curve of KIC 5771589 by 0.2 days, for a better view of the small effect of variations in the inner eccentricity forced by the octupole perturbations. See text for further details.
Figure 1. Observed and modeled ETVs of KICs 9714358 (upper), 5771578 (middle), and TIC 219885468 (bottom). The left-hand panels cover the time interval approximately from the beginning of the Kepler observations to the end of the second northern ecliptic hemisphere scan of TESS, while the right-hand panels show longer, 100 yr-long intervals, which offer direct predictions for future observations. Larger red and blue symbols represent primary and secondary ETV points determined from the observations, while the smaller points connected with straight lines stand for ETV points calculated according to the best-fitting complex photodynamical models. Finally, the black, almost horizontal lines represent the LTTE contributions to the curves. Note that in the middle right panel we downshifted the originally overlapping secondary ETV curve of KIC 5771589 by 0.2 days, for a better view of the small effect of variations in the inner eccentricity forced by the octupole perturbations. See text for further details.
Universe 09 00485 g001
Figure 2. APOGEE-2 RVs of the primary component of KIC 9714358 (red dots) together with the best photodynamical model (blue curves). Upper panel: A characteristic section of the observed RV points and the corresponding model fit in the time domain. Both the inner and the outer orbits are well visible. Lower left panel: Phase-folded RV data according to the inner period, after the removal of the signals of the outer orbit. Lower right panel: Phase-folded RV data according to the outer period, after the removal of the effects of the inner orbital motion. The varying shape of the folded solution, caused primarily by the rapid apsidal motion of the outer orbit is readily visible. Note, the thin horizontal black lines at V = 66.91 km s 1 denote the systemic radial velocity of the triple star.
Figure 2. APOGEE-2 RVs of the primary component of KIC 9714358 (red dots) together with the best photodynamical model (blue curves). Upper panel: A characteristic section of the observed RV points and the corresponding model fit in the time domain. Both the inner and the outer orbits are well visible. Lower left panel: Phase-folded RV data according to the inner period, after the removal of the signals of the outer orbit. Lower right panel: Phase-folded RV data according to the outer period, after the removal of the effects of the inner orbital motion. The varying shape of the folded solution, caused primarily by the rapid apsidal motion of the outer orbit is readily visible. Note, the thin horizontal black lines at V = 66.91 km s 1 denote the systemic radial velocity of the triple star.
Universe 09 00485 g002
Figure 3. Upper panel: The ∼4-year-long Kepler light curve of KIC 5771589 (blue dots) with the best-fitted photodynamical model (red). Alternating gray-white stripes denote the consecutive observing quarters. Note, Q5, Q9 and Q13 data are missing as in these quarters the stellar light had fallen onto a bad CCD camera. Lower left panel: Consecutive primary and secondary eclipses from the beginning, mid-time, and near-end of the prime Kepler mission. Lower right panel: Characteristic primary and secondary eclipses from the 2019, 2021, and 2022 visits of TESS. Each of the sections is 1 day long.
Figure 3. Upper panel: The ∼4-year-long Kepler light curve of KIC 5771589 (blue dots) with the best-fitted photodynamical model (red). Alternating gray-white stripes denote the consecutive observing quarters. Note, Q5, Q9 and Q13 data are missing as in these quarters the stellar light had fallen onto a bad CCD camera. Lower left panel: Consecutive primary and secondary eclipses from the beginning, mid-time, and near-end of the prime Kepler mission. Lower right panel: Characteristic primary and secondary eclipses from the 2019, 2021, and 2022 visits of TESS. Each of the sections is 1 day long.
Universe 09 00485 g003
Figure 4. Primary and secondary eclipses from the beginning of Year 2 and Year 4, as well as from the end of Year 5 TESS-observations of TIC 219885468. Each of the sections is 1 day long.
Figure 4. Primary and secondary eclipses from the beginning of Year 2 and Year 4, as well as from the end of Year 5 TESS-observations of TIC 219885468. Each of the sections is 1 day long.
Universe 09 00485 g004
Figure 5. Some satellite-observed primary and secondary eclipses of KIC 9714358 (blue dots) with the best-fitting photodynamical model (red curves). Left panel: Consecutive primary and secondary eclipses from the beginning, mid-time, and near to the end of the prime Kepler mission. Right panel: Characteristic primary and secondary eclipses from the 2019 and 2022 visits of TESS. Each of the sections is 0.4 day long.
Figure 5. Some satellite-observed primary and secondary eclipses of KIC 9714358 (blue dots) with the best-fitting photodynamical model (red curves). Left panel: Consecutive primary and secondary eclipses from the beginning, mid-time, and near to the end of the prime Kepler mission. Right panel: Characteristic primary and secondary eclipses from the 2019 and 2022 visits of TESS. Each of the sections is 0.4 day long.
Universe 09 00485 g005
Figure 6. The evolution of the orbits in KIC 9714358 from its discovery to the end of the current century. Left panel: The variations of the observable and the dynamical arguments of periastron of the inner and outer orbits ( ω in , out and ω in dyn , out , respectively). Right panel: The cyclic variation of the inner eccentricity, and the very similar variations of the difference of the dynamical arguments of periastrons of the two orbits. The gray shaded area represents the interval of the Kepler observations. See text for further details.
Figure 6. The evolution of the orbits in KIC 9714358 from its discovery to the end of the current century. Left panel: The variations of the observable and the dynamical arguments of periastron of the inner and outer orbits ( ω in , out and ω in dyn , out , respectively). Right panel: The cyclic variation of the inner eccentricity, and the very similar variations of the difference of the dynamical arguments of periastrons of the two orbits. The gray shaded area represents the interval of the Kepler observations. See text for further details.
Universe 09 00485 g006
Figure 7. Illustration of the effect of the orientation of the orbit on the properties of the eclipses. For illustration purposes, we generated light curves for a hypothetical binary with parameters: m Aa = 1.00 M , m Ab = 0.80 M ; R Aa = 1.03 R , m Ab = 0.799 R ; T Aa = 5878 K, T Ab = 5274 K; P = 5.00 d; e = 0.30 ; and i = 85 . The only difference amongst the four rows is the orientation of the apsidal line of the ellipse, which are ω = 0 , 90 , 180 , and 270 , respectively. The left panels represent the synthetic light curves, which in all four cases are phased to the inferior conjunction of the less massive secondary, while the right panels show the corresponding geometry, as seen from the Earth, and as it would be seen from the pole of the orbit (small, inserted figures). (In the case of the insert plots, the direction of the Earthly observer is along the negative Y axis. Note also that the secondary revolves counter-clockwise in all the orbit plots.) As one can see on the light curve plots, the depth ratio of the primary and secondary eclipses depends strongly on the orientation of the ellipse. When the ellipse is seen along the minor axis ( ω = 0 , 180 ), the eclipsed disk areas are similar during the inferior and superior conjunctions; hence, the ratio of the eclipse depths—similar to the circular orbit case—primarily depends only upon the ratio of the surface brightnesses. On the other hand, when the ellipse is seen along its major axis ( ω = 90 and 270 ), the eclipsed surface area is much larger during the periastron-eclipse than during the apastron-eclipse. In the current illustration, due to the relatively large eccentricity of e = 0.3 , the difference in the eclipsed surface areas is so large that the eclipse depths reverse.
Figure 7. Illustration of the effect of the orientation of the orbit on the properties of the eclipses. For illustration purposes, we generated light curves for a hypothetical binary with parameters: m Aa = 1.00 M , m Ab = 0.80 M ; R Aa = 1.03 R , m Ab = 0.799 R ; T Aa = 5878 K, T Ab = 5274 K; P = 5.00 d; e = 0.30 ; and i = 85 . The only difference amongst the four rows is the orientation of the apsidal line of the ellipse, which are ω = 0 , 90 , 180 , and 270 , respectively. The left panels represent the synthetic light curves, which in all four cases are phased to the inferior conjunction of the less massive secondary, while the right panels show the corresponding geometry, as seen from the Earth, and as it would be seen from the pole of the orbit (small, inserted figures). (In the case of the insert plots, the direction of the Earthly observer is along the negative Y axis. Note also that the secondary revolves counter-clockwise in all the orbit plots.) As one can see on the light curve plots, the depth ratio of the primary and secondary eclipses depends strongly on the orientation of the ellipse. When the ellipse is seen along the minor axis ( ω = 0 , 180 ), the eclipsed disk areas are similar during the inferior and superior conjunctions; hence, the ratio of the eclipse depths—similar to the circular orbit case—primarily depends only upon the ratio of the surface brightnesses. On the other hand, when the ellipse is seen along its major axis ( ω = 90 and 270 ), the eclipsed surface area is much larger during the periastron-eclipse than during the apastron-eclipse. In the current illustration, due to the relatively large eccentricity of e = 0.3 , the difference in the eclipsed surface areas is so large that the eclipse depths reverse.
Universe 09 00485 g007aUniverse 09 00485 g007b
Figure 8. The same plots as in Figure 6 but for KIC 5771589. See text for further details.
Figure 8. The same plots as in Figure 6 but for KIC 5771589. See text for further details.
Universe 09 00485 g008
Figure 9. Variations of the observable inclinations of the inner and outer orbital planes ( i in —red; i out —blue, respectively) of KIC 5771589 (left) and TIC 219885468 (right). The green area represents the domain where the inner binary produces at least one eclipse during a revolution, while the cyan area covers that inclination domain where both eclipses are present. Furthermore, the nearly horizontal black lines denote the lower and upper limits of that domain of the outer inclination ( i out ), where third-body eclipses might occur. As one can see, in the case of KIC 5771589, the inner inclination ( i in ) shows small amplitude oscillations near the border of the EB domain, implying that the eclipses are grazing. Moreover, the blue curve ( i out ) remains continuously below the third-body eclipse inclination limit; hence, no third-body eclipses can occur. In contrast to this, in the case of TIC 219885468 the small amplitude oscillations of i in are located quite deeply in the (cyan/green) EB domain, indicating that the eclipses are quite deep and, hence, the eclipse depths are insensitive to such small inclination variations. Note also that this is an almost triply eclipsing system as the outer inclination is just below the triply eclipsing limit.
Figure 9. Variations of the observable inclinations of the inner and outer orbital planes ( i in —red; i out —blue, respectively) of KIC 5771589 (left) and TIC 219885468 (right). The green area represents the domain where the inner binary produces at least one eclipse during a revolution, while the cyan area covers that inclination domain where both eclipses are present. Furthermore, the nearly horizontal black lines denote the lower and upper limits of that domain of the outer inclination ( i out ), where third-body eclipses might occur. As one can see, in the case of KIC 5771589, the inner inclination ( i in ) shows small amplitude oscillations near the border of the EB domain, implying that the eclipses are grazing. Moreover, the blue curve ( i out ) remains continuously below the third-body eclipse inclination limit; hence, no third-body eclipses can occur. In contrast to this, in the case of TIC 219885468 the small amplitude oscillations of i in are located quite deeply in the (cyan/green) EB domain, indicating that the eclipses are quite deep and, hence, the eclipse depths are insensitive to such small inclination variations. Note also that this is an almost triply eclipsing system as the outer inclination is just below the triply eclipsing limit.
Universe 09 00485 g009
Figure 10. Simulated light curve of KIC 5771589 in Kepler’s photometric band for all of the 21st century. The gray, shaded area, as before, represents the interval of the Kepler observations, while the other vertical lines stand for the past and scheduled future observations of TESS.
Figure 10. Simulated light curve of KIC 5771589 in Kepler’s photometric band for all of the 21st century. The gray, shaded area, as before, represents the interval of the Kepler observations, while the other vertical lines stand for the past and scheduled future observations of TESS.
Universe 09 00485 g010
Figure 11. The same plots as in Figure 6, but for TIC 219885468. The filled brown area represents the interval between the beginning of the very first (S14) and the end of the momentarily last (S60) NCVZ observations of TESS. See text for further details.
Figure 11. The same plots as in Figure 6, but for TIC 219885468. The filled brown area represents the interval between the beginning of the very first (S14) and the end of the momentarily last (S60) NCVZ observations of TESS. See text for further details.
Universe 09 00485 g011
Table 1. Main properties of the three systems from different catalogs.
Table 1. Main properties of the three systems from different catalogs.
ParameterKIC 9714358KIC 5771589TIC 219885468
RA (J2000)19:34:09.66718:58:20.52617:28:57.721
Dec (J2000)+46:26:14.02+41:00:34.46+70:42:03.09
G  b 14.9892 ± 0.0007 11.8349 ± 0.0020 13.0798 ± 0.0002
G BP   b 15.4996 ± 0.0031 12.0314 ± 0.0005 13.3429 ± 0.0001
G RP   b 14.3166 ± 0.0024 11.2839 ± 0.0004 12.6629 ± 0.0005
a 16.336 ± 0.067 12.512 ± 0.013 13.729 ± 0.200
a 15.422 ± 0.1114 11.943 ± 0.027 13.200 ± 0.069
g 15.632 ± 0.004   c 12.228 ± 0.001   d 13.442 ± 0.196  e
r 14.989 ± 0.004   c 11.849 ± 0.001   d 13.015 ± 0.173  e
i 14.696 ± 0.002   c 11.691 ± 0.001   d 12.905 ± 0.205  e
f 13.513 ± 0.023 10.787 ± 0.020 12.233 ± 0.021
f 13.001 ± 0.026 10.511 ± 0.021 11.970 ± 0.023
f 12.884 ± 0.025 10.453 ± 0.018 11.942 ± 0.022
W1  g 12.801 ± 0.023 10.379 ± 0.023 11.885 ± 0.023
W2  g 12.794 ± 0.023 10.408 ± 0.020 11.909 ± 0.021
W3  g 12.619 ± 0.368 10.340 ± 0.049 11.973 ± 0.125
W4  g 9.198 8.842 9.577
NUV 22.713 ± 0.324 17.232 ± 0.032
T eff (K)  a 4764 ± 109 6231 6350 ± 133
Distance (pc)  h 577 ± 7 870 ± 100 1111 ± 13
[ M / H ]   a 0.462 ± 0.013 0.083 ± 0.011
E ( B V )   a 0.050 0.030 ± 0.007
μ α (mas yr 1 b 8.77 ± 0.02 3.02 ± 0.15 1.02 ± 0.01
μ δ (mas   yr 1 b 3.14 ± 0.02 1.19 ± 0.17 0.64 ± 0.02
RUWE  b 1.03356 11.34329 0.94989
Notes. a TESS Input Catalog (TIC v8.2) [40]. b Gaia EDR3 [41]. c PanSTARRS, [42]. d The Kepler-INT survey [43]. e AAVSO Photometric All Sky Survey (APASS) DR9, [44], http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=II/336/apass9 (accessed on 7 November 2023). f 2MASS catalog [45]. g WISE point source catalog [46]. h Photogeometric distances from Bailer-Jones et al. [47]. Note also that for the SED analysis in Section 5, the uncertainties of the passband magnitudes were set to σ mag = max ( σ catalog , 0.030 ) to avoid the strong overdominance of the extremely accurate Gaia magnitudes over the other measurements.
Table 2. Orbital and astrophysical parameters of KIC 97143586 from the joint photodynamical lightcurve, ETV, RV, SED, and PARSEC isochrone solution. Meanings of most of the parameters are explained in the text, with the exceptions of i inv and Ω inv , whose quantities give the position of the invariable plane with respect to the tangent plane of the sky. The osculating orbital elements are given for epoch t 0 = 2,454,953.0.
Table 2. Orbital and astrophysical parameters of KIC 97143586 from the joint photodynamical lightcurve, ETV, RV, SED, and PARSEC isochrone solution. Meanings of most of the parameters are explained in the text, with the exceptions of i inv and Ω inv , whose quantities give the position of the invariable plane with respect to the tangent plane of the sky. The osculating orbital elements are given for epoch t 0 = 2,454,953.0.
KIC 9714358
Orbital Elements
Subsystem
Aa–AbA–B
P [days] 6 . 47075 0.00019 + 0.00019 104 . 083 0.010 + 0.010
a [R ] 15 . 32 0.10 + 0.21 105 . 6 0.7 + 1.4
e 0 . 02862 0.00034 + 0.00038 0 . 2524 0.0020 + 0.0019
ω [deg] 120 . 35 0.51 + 0.52 102 . 06 0.61 + 0.65
i [deg] 87 . 594 0.036 + 0.083 87 . 531 0.040 + 0.086
T 0 inf [BJD −2,400,000]54,967.3844 0.0004 + 0.0004
τ [BJD −2,400,000] 54 , 964 . 7250 0.0095 + 0.0097 54,977.087 0.083 + 0.082
Ω [deg] 0.0 0 . 011 0.041 + 0.026
i mut [deg] 0 . 074 0.026 + 0.041
ϖ dyn [deg] 303 . 4 0.5 + 0.5 282 . 1 0.6 + 0.6
i dyn [deg] 0 . 055 0.020 + 0.030 0 . 019 0.007 + 0.010
Ω dyn [deg] 189 24 + 20 9 24 + 20
i inv [deg] 87 . 547 0.037 + 0.083
Ω inv [deg] 0 . 008 0.030 + 0.019
mass ratio [ q = m sec / m pri ] 0 . 408 0.008 + 0.004 0 . 265 0.003 + 0.003
RV amplitude K pri [km s 1 ] 34 . 69 0.15 + 0.16 11 . 12 0.12 + 0.13
RV amplitude K sec [km s 1 ] 85 . 08 0.76 + 1.73 41 . 92 0.31 + 0.59
V γ [km s 1 ] 66 . 91 0.02 + 0.02
Apsidal and Nodal-Motion-Related Parameters
P apse [year] 26 . 43 0.20 + 0.17 76 . 84 0.47 + 0.61
P apse dyn [year] 11 . 27 0.07 + 0.06 15 . 63 0.07 + 0.07
P node dyn [year] 19 . 63 0.10 + 0.10
Δ ω 3 b [arcsec/cycle] 2037 11 + 13 23622 104 + 100
Δ ω GR [arcsec/cycle] 0 . 622 0.008 + 0.017 0 . 122 0.002 + 0.003
Δ ω tide [arcsec/cycle] 0 . 182 0.018 + 0.007 0 . 00072 0.00007 + 0.00003
Stellar Parameters
AaAbB
Relative Quantities
fractional radius [ R / a ] 0 . 0488 0.0004 + 0.0003 0 . 0321 0.0011 + 0.0002 0 . 00447 0.00014 + 0.00004
temperature relative to ( T eff ) Aa 1 0 . 6118 0.0056 + 0.0079 0 . 6040 0.0056 + 0.0075
fractional flux [in Kepler-band] 0 . 9135 0.0293 + 0.0118 0 . 0365 0.0059 + 0.0018 0 . 0314 0.0040 + 0.0016
fractional flux [in TESS-band] 0 . 8934 0.0221 + 0.0201 0 . 0435 0.0070 + 0.0018 0 . 0353 0.0041 + 0.0018
Physical Quantities
m [M ] 0 . 818 0.018 + 0.039 0 . 334 0.005 + 0.009 0 . 306 0.007 + 0.011
R [R ] 0 . 746 0.006 + 0.010 0 . 491 0.009 + 0.005 0 . 471 0.010 + 0.005
T eff [K] 5312 118 + 82 3237 31 + 56 3197 32 + 55
L bol [L ] 0 . 397 0.036 + 0.034 0 . 023 0.001 + 0.001 0 . 021 0.002 + 0.001
M bol 5 . 77 0.09 + 0.10 8 . 84 0.05 + 0.07 8 . 99 0.04 + 0.10
M V 5 . 94 0.11 + 0.14 11 . 02 0.15 + 0.21 11 . 25 0.12 + 0.25
log g [dex] 4 . 605 0.009 + 0.008 4 . 575 0.006 + 0.032 4 . 572 0.006 + 0.031
Global System Parameters
log(age) [dex] 7 . 666 0.123 + 0.048
[ M / H ] [dex] 0 . 129 0.045 + 0.058
E ( B V ) [mag] 0 . 087 0.046 + 0.047
extra light 4 [in Kepler-band] 0 . 017 0.012 + 0.032
extra light 4 [in TESS-band] 0 . 030 0.018 + 0.023
( M V ) tot 5 . 92 0.11 + 0.14
distance [pc] 666 11 + 10
Table 3. Orbital and astrophysical parameters of KIC 5771589 from the joint photodynamical lightcurve, ETV, SED, and PARSEC isochrone solution. The osculating orbital elements are given for epoch t 0 = 2,454,953.0.
Table 3. Orbital and astrophysical parameters of KIC 5771589 from the joint photodynamical lightcurve, ETV, SED, and PARSEC isochrone solution. The osculating orbital elements are given for epoch t 0 = 2,454,953.0.
KIC 5771589
Orbital Elements
subsystem
Aa–AbA–B
P [days] 10 . 6791 0.0033 + 0.0024 113 . 872 0.020 + 0.023
a [R ] 24 . 20 0.10 + 0.10 137 . 5 0.6 + 0.5
e 0 . 00358 0.00025 + 0.00042 0 . 1615 0.0018 + 0.0018
ω [deg] 176 . 3 3.4 + 3.4 91 . 12 0.99 + 1.07
i [deg] 86 . 048 0.085 + 0.073 86 . 154 0.084 + 0.074
T 0 inf [BJD −2,400,000] 54 , 962 . 0758 0.0044 + 0.0029 ...
τ [BJD −2,400,000] 54 , 959 . 308 0.101 + 0.101 54 , 974 . 152 0.174 + 0.186
Ω [deg] 0.0 0 . 266 0.025 + 0.026
i mut [deg] 0 . 286 0.025 + 0.026
ϖ dyn [deg] 356 . 3 3.4 + 3.4 271 . 1 1.0 + 1.1
i dyn [deg] 0 . 235 0.021 + 0.021 0 . 051 0.005 + 0.010
Ω dyn [deg] 68 . 0 2.1 + 2.0 248 . 0 2.1 + 2.0
i inv [deg] 86 . 135 0.084 + 0.073
Ω inv [deg] 0 . 219 0.021 + 0.021
mass ratio [ q = m sec / m pri ] 0 . 776 0.011 + 0.010 0 . 612 0.003 + 0.003
RV amplitude K pri [km s 1 ] 49 . 98 0.46 + 0.50 23 . 47 0.14 + 0.13
RV amplitude K sec [km s 1 ] 64 . 47 0.45 + 0.53 38 . 31 0.17 + 0.19
Apsidal and Nodal-Motion-Related Parameters
P apse [year] 11 . 216 0.036 + 0.036 51 . 73 0.15 + 0.15
P apse dyn [year] 5 . 059 0.014 + 0.015 7 . 823 0.016 + 0.018
P node dyn [year] 9 . 217 0.024 + 0.023
Δ ω 3 b [arcsec/cycle] 7489 23 + 22 51651 113 + 101
Δ ω GR [arcsec/cycle] 0 . 569 0.005 + 0.005 0 . 166 0.001 + 0.001
Δ ω tide [arcsec/cycle] 0 . 095 0.010 + 0.014 0 . 00054 0.00008 + 0.00010
Stellar Parameters
AaAbB
Relative Quantities
fractional radius [ R / a ] 0 . 0475 0.0012 + 0.0016 0 . 0286 0.0003 + 0.0004 0 . 0168 0.0020 + 0.0010
temperature relative to ( T eff ) Aa 1 0 . 8475 0.0081 + 0.0080 0 . 9176 0.0108 + 0.0207
fractional flux [in Kepler-band] 0 . 2230 0.0176 + 0.0200 0 . 0419 0.0037 + 0.0046 0 . 6397 0.0526 + 0.0467
fractional flux [in TESS-band] 0 . 2332 0.0207 + 0.0230 0 . 0472 0.0048 + 0.0068 0 . 6954 0.0479 + 0.0272
Physical Quantities
m [M ] 0 . 938 0.013 + 0.014 0 . 727 0.010 + 0.013 1 . 017 0.013 + 0.013
R [R ] 1 . 150 0.032 + 0.040 0 . 694 0.009 + 0.012 2 . 303 0.259 + 0.146
T eff [K] 6260 77 + 42 5295 70 + 67 5746 42 + 72
L bol [L ] 1 . 811 0.147 + 0.173 0 . 340 0.021 + 0.024 5 . 209 0.946 + 0.531
M bol 4 . 12 0.10 + 0.09 5 . 94 0.07 + 0.07 2 . 98 0.11 + 0.22
M V 4 . 17 0.10 + 0.09 6 . 12 0.09 + 0.09 3 . 08 0.10 + 0.20
log g [dex] 4 . 289 0.026 + 0.020 4 . 617 0.007 + 0.006 3 . 718 0.052 + 0.111
Global System Parameters
log(age) [dex] 9 . 832 0.018 + 0.020
[ M / H ] [dex] 0 . 504 0.014 + 0.138
E ( B V ) [mag] 0 . 004 0.003 + 0.010
extra light 4 [in Kepler-band] 0 . 100 0.051 + 0.038
extra light 4 [in TESS-band] 0 . 023 0.018 + 0.025
( M V ) tot 2 . 70 0.09 + 0.16
distance [pc] 709 55 + 33
Table 4. Orbital and astrophysical parameters of TIC 219885468 from the joint photodynamical lightcurve, ETV, SED, and PARSEC isochrone solution. The osculating orbital elements are given for epoch t 0 = 2,458,683.0.
Table 4. Orbital and astrophysical parameters of TIC 219885468 from the joint photodynamical lightcurve, ETV, SED, and PARSEC isochrone solution. The osculating orbital elements are given for epoch t 0 = 2,458,683.0.
TIC 219885468
Orbital Elements
Subsystem
Aa–AbA–B
P [days] 7 . 51281 0.00043 + 0.00040 111 . 5498 0.0071 + 0.0063
a [R ] 21 . 52 0.14 + 0.21 142 . 6 0.9 + 1.4
e 0 . 04233 0.00058 + 0.00060 0 . 3903 0.0007 + 0.0007
ω [deg] 259 . 33 0.22 + 0.23 256 . 02 0.64 + 0.69
i [deg] 88 . 227 0.051 + 0.050 88 . 125 0.218 + 0.206
T 0 inf [BJD −2,400,000] 58 , 686 . 8073 0.0004 + 0.0004 ...
τ [BJD −2,400,000] 54 , 964 . 7250 0.0095 + 0.0097 58 , 704 . 054 0.050 + 0.052
Ω [deg] 0.0 0 . 073 0.192 + 0.191
i mut [deg] 0 . 228 0.089 + 0.164
ϖ dyn [deg] 79 . 33 0.22 + 0.23 76 . 01 0.64 + 0.69
i dyn [deg] 0 . 166 0.065 + 0.119 0 . 063 0.025 + 0.045
Ω dyn [deg] 145 72 + 59 325 72 + 59
i inv [deg] 88 . 15 0.17 + 0.16
Ω inv [deg] 0 . 053 0.14 + 0.14
mass ratio [ q = m sec / m pri ] 0 . 988 0.017 + 0.016 0 . 319 0.002 + 0.002
RV amplitude K pri [km s 1 ] 72 . 26 1.14 + 0.90 17 . 00 0.14 + 0.19
RV amplitude K sec [km s 1 ] 73 . 15 0.83 + 0.57 53 . 23 0.36 + 0.54
Apsidal and Nodal-Motion-Related Parameters
P apse [year] 19 . 489 0.087 + 0.087 51 . 24 0.10 + 0.10
P apse dyn [year] 8 . 178 0.033 + 0.034 11 . 051 0.038 + 0.041
P node dyn [year] 14 . 090 0.068 + 0.064
Δ ω 3 b [arcsec/cycle] 3258 13 + 13 35 , 816 131 + 123
Δ ω GR [arcsec/cycle] 0 . 911 0.012 + 0.018 0 . 214 0.003 + 0.004
Δ ω tide [arcsec/cycle] 0 . 997 0.027 + 0.027 0 . 0047 0.0001 + 0.0001
Stellar Parameters
AaAbB
Relative Quantities
fractional radius [ R / a ] 0 . 0651 0.0013 + 0.0015 0 . 0633 0.0014 + 0.0013 0 . 00487 0.00004 + 0.00007
temperature relative to ( T eff ) Aa 1 0 . 9989 0.0016 + 0.0014 0 . 7625 0.0051 + 0.0044
fractional flux [in TESS-band] 0 . 4915 0.0195 + 0.0233 0 . 4628 0.0214 + 0.0202 0 . 0449 0.0021 + 0.0025
Physical Quantities
m [M ] 1 . 189 0.021 + 0.033 1 . 178 0.033 + 0.038 0 . 755 0.015 + 0.024
R [R ] 1 . 407 0.032 + 0.024 1 . 367 0.042 + 0.031 0 . 694 0.009 + 0.018
T eff [K] 6405 36 + 46 6399 36 + 41 4885 50 + 53
L bol [L ] 3 . 009 0.174 + 0.115 2 . 814 0.212 + 0.190 0 . 246 0.016 + 0.023
M bol 3 . 57 0.04 + 0.06 3 . 65 0.07 + 0.08 6 . 29 0.10 + 0.07
M V 3 . 57 0.04 + 0.06 3 . 65 0.07 + 0.09 6 . 64 0.12 + 0.10
log g [dex] 4 . 221 0.020 + 0.015 4 . 241 0.016 + 0.013 4 . 632 0.008 + 0.003
Global System Parameters
log(age) [dex] 9 . 498 0.068 + 0.044
[ M / H ] [dex] 0 . 089 0.023 + 0.054
E ( B V ) [mag] 0 . 029 0.012 + 0.020
( M V ) tot 2 . 84 0.05 + 0.04
distance [pc] 1152 13 + 22
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Borkovits, T.; Mitnyan, T. Observational Detection of Higher-Order Secular Perturbations in Tight Hierarchical Triple Stars. Universe 2023, 9, 485. https://doi.org/10.3390/universe9110485

AMA Style

Borkovits T, Mitnyan T. Observational Detection of Higher-Order Secular Perturbations in Tight Hierarchical Triple Stars. Universe. 2023; 9(11):485. https://doi.org/10.3390/universe9110485

Chicago/Turabian Style

Borkovits, Tamás, and Tibor Mitnyan. 2023. "Observational Detection of Higher-Order Secular Perturbations in Tight Hierarchical Triple Stars" Universe 9, no. 11: 485. https://doi.org/10.3390/universe9110485

APA Style

Borkovits, T., & Mitnyan, T. (2023). Observational Detection of Higher-Order Secular Perturbations in Tight Hierarchical Triple Stars. Universe, 9(11), 485. https://doi.org/10.3390/universe9110485

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop