Free Access
Issue
A&A
Volume 610, February 2018
Article Number A72
Number of page(s) 9
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201630072
Published online 05 March 2018

© ESO 2018

1 Introduction

The discovery for an extraneous component in a previous known stellar system by finding extraneous eclipses on the light curves has a long history. The earliest binary star system β Persei1, known as Algol, was discovered in this way more than 3000 years ago. Most of the binary systems were found (and also studied) by the eclipse signals, hence the name “eclipsing binaries”. This method also works for the multiple systems. The extraneous eclipses on the light curves of previous known binaries suggest a new component accompanying the binary system, just like the eclipses on a single star’s light curve suggest a companion to the single star. Since the studies on eclipsing binaries have provided numerous fundamental properties for the stellar field, similarly, the eclipsing multiple systems will also provide important information, especially in regards to the star formation environment.

The biggest problem in the study of eclipsing multiple systems is small sample sizes, and this situation should not be caused by the scarcity of multiple systems but caused by the limitations of observation. The emergence odds of the extraneous eclipses are slim, and therefore until the release of the Kepler space telescope it has been hard to find them confidently. The Kepler telescope is devoted to discovering earth-size planets orbiting other stars, and so to monitor photometrically many kinds of stars including binaries. The extraordinary features of the Kepler observations, that is, the long-term (four years) uninterrupted brightness monitoring in a 115 deg2 fixed field of view with unprecedented photometric precision (20 per million) and high limiting magnitude (21st magnitude), are ideal and suitable for searching for extraneous eclipses on the binary light curves. In fact, the Kepler telescope was designed and operated for the purpose of finding extraneous eclipses indicating the planets around stars.

In this paper, three binaries with extraneous eclipses, KIC 5255552, KIC 10091110 and KIC 11495766, are studied. Section 2 illustrates the data reduction. Studies of the three objects are outlined in Sect. 3. A first list of eclipsing multiple candidates is presented in Sect. 4. We summarize our work in the last section.

2 Data reduction

Binary light curves were provided by the Kepler mission, and they were downloaded from Mikulski Archive for Space Telescopes (MAST) at a convenient individual webpage2 on October 27, 2015. Besides the photometric data, the Kepler Eclipsing Binary Catalog (Prša et al. 2011; Slawson et al. 2011; Matijevič et al. 2012; Conroy et al. 2014; Kirk et al. 2016) is another important resource. The catalog provides multiple orbital periods for some binaries, which led us to the extraneous eclipses in these binaries.

2.1 Normalization of Kepler light curves

The Pre-Search Data Conditioning (PDC) fluxes in the Kepler data were selected “to identify and correct flux discontinuities that cannot be attributed to known spacecraft or data anomalies” (Kepler Data Processing Handbook3; Fraquelli et al. Fraquelli et al. 2014). However, there are still unphysical jumps and long-term (week or longer) variations in the PDC fluxes. Almost all the jumps occurred after the light curve gaps (where there is no datum); about half of the gaps were the borders of the observation quarters of the Kepler spacecraft, and the other half were caused by the missing data.

These unphysical jumps and long-term variations can be eliminated simultaneously by the flux normalization on each light curve segment (divided by the gaps). However, the long-term variations contain not only the unphysical but also the physical parts, and the physical parts will be also eliminated by the flux normalization together with the unphysical parts. Fortunately, the elimination of the long-term variations does not affect the analysis on the binary light variations, because the binary light variations are always in short term, whether the orbital period is long or short (for the long-period binaries, their light variations, that is, the eclipses, are also very short). In fact, the elimination is often more conducive to the analysis of binary light curves, because some irrelevant light variations (such as magnetic activities, pulsations) are removed.

The flux normalization steps are: (1) the light curves are divided into dozens of segments depending on the width of the gaps on the light curves. If the width of a gap is larger than one day (according to our experiences, somewhat arbitrary), the light curve will be divided into two segments at that gap; (2) each light segment is fitted by a cubic polynomial model, and their fluxes are divided by the fitting values to be normalized. In the fitting, the large scattered points are removed to get a better fit in an iterative process (up to ten times practically); and (3) all the normalized segments are stitched together to form the whole normalized light curve.

2.2 Data binning

Before the light curve analysis, the barycentric Julian dates (BJD) need to be transformed to phase. The orbital periods for calculating the phases were taken from the Kepler Eclipsing Binary Catalog, but the reference minimum times were re-measured, not because the values in the catalog are not accurate enough, but because the mean reference minimum times4 (instead of a single minimum time) are needed. After the phase calculation, the dense phase points are binned to reduce the data amount, and also to get rid of the large scattered points.

The binning phase data points were made to be not evenly distributed along the phases but inclined to place more points at eclipses. In other words, the points at the eclipsing parts will be more dense along the abscissa. It is believed that this kind of data binning is more suitable for the binary analysis, because more points (so more weight) should be given to the important eclipsing parts. The final binning phase data are analyzed in the following section.

3 Three Kepler binaries with extraneous eclipses

3.1 KIC 5255552

3.1.1 The extraneous eclipses on KIC 5255552

The object KIC 5255552 was discovered and analyzed as an eclipsing binary with 32.448635 d period by Prša et al. (2011) and Slawson et al. (2011) based on the Kepler data. Borkovits et al. (2015, 2016) measured and analyzed its O−C curve with a new and robust tool which is capable of dealing with eccentric binaries, and determined the reliable parameters. Borkovits et al. (2016) obtained that the 32.448635 d period inner binary has a 0.7 (± 0.1) M tertiary body revolving around it with 862.1 (± 0.1) d period at only 6.4 (± 0.1) degrees of mutual inclination, and also noted three outer eclipsing events (Borkovits et al. 2015).

In this paper, four groups of extraneous eclipses are reported, illustrated in Figs. 1 and 2. The extraneous eclipses, emerging in groups, are generated by the mutual blocking between the binary components and the outer companion. When the outer companion passes in front of (or behind) the binary, it will block (or be blocked by) the two components of the binary. The blocking may generate several eclipses for one passage, as shown in Fig. 2.

The first, second, and fourth groups of extraneous eclipses (see the red points in Fig. 1) are consistent with the 862.1(±0.1) d period of the third companion concluded by Borkovits et al. (2016). These two mutual confirmations are strong evidence forthe triplicity of this object. However, a much shallower and narrower fading was also found at ~1279 d, which happens to constitute a period of ~589 d with another three groups (see the gray arrows in Fig. 1). Considering that the ~1279 fading is not analogous to any other three groups, so the apparent ~589 d period may simply be a coincidence rather than a real period. If this small fading is physical, it indicates another new companion that makes the KIC 5255552 a quadruple system.

3.1.2 The light analysis of the inner binary

The inner 32.448635 d period binary was changing dramatically on eclipses (see the minimum position changes in lower-right panels of Fig. 1), so the analysis on the whole light curves is unfeasible. In view of this, four light curve parts with lengths of about one period were selected and analyzed independently. The results are listed in Table 1 and shown in Fig. 15.

The mass ratio was taken from Borkovits et al. (2015)6, and was fixed in the fitting process. This is because the photometric determination on the mass ratio for a long-period detached binary tends to have considerably great uncertainty. The periods in Table 1 were calculated separately by subtracting two successive deeper minimum times within the selected time span, and they differ somewhat, which can be seen numerically. Similarly, the other parameters for different time spans have changes of a few percent.

The changing parameters with time indicate that the orbit of the 32.448635 d period binary may differ from Keplerian orbit by some degree. This is likely to be due to the disturbance by the massive (MouterMinner = 0.29), close (PouterPinner = 26.2), and coplanar (im = 6.4) outer companions (Borkovits et al. 2016). The KIC 5255552 is not the only system with a non-Keplerian orbit, the object KIC 2856960 is a known triple system whose extraneous eclipses cannot be fitted while satisfying the Kepler’s laws (Marsh et al. 2014).

In summary, the object KIC 5255552 is a hierarchical eclipsing multiple system consisting of three (or four) component stars. The orbit of the inner binary is changing visibly with time, which makes this object a good sample for studying the interaction in the multiplesystems.

thumbnail Fig. 1

Light curves of KIC 5255552. Upper: the whole light curves. The red points are the extraneous eclipses that can been seen more clearly in Fig. 2, and the gray arrows are used to illustrate the time intervals between the extraneous eclipses. Lower eight panels: the four selected light curve parts in phase form, showing the eclipses at phase zero (left) and phase around 0.43 (right) with their fittings (black line).

thumbnail Fig. 2

Extraneous eclipses of KIC 5255552. The red points in all the panels are the extraneous eclipses. The mean times of the extraneous eclipses for each group are marked in each panel. In the top panel, a small additional panel is used to magnify a tiny extraneous eclipse around 1537.7 d.

Table 1

Binary light curve solutions of KIC 5255552.

3.2 KIC 10091110

3.2.1 Light curve analysis

The object KIC 10091110 has two sets of very obvious eclipses with periods of 4.2185174 d and 8.5303353 d, which are shown in Fig. 3, panels (a) and (b)7. The 4.2185174 d period eclipses have secondary minimums that can be explicitly seen in panel (d), and the light curves can be well fitted by the binary model Wilson-Devinney program (Wilson & Devinney 1971; Wilson 1979, 1990, 2008, 2012; Van Hamme & Wilson 2007; Wilson et al. 2010). The atmosphere parameters of KIC 10091110 were Teff = 6185(204), log g = 4.472(0.166) and [FeH] = −0.18(0.27) (Rowe et al. 2015), and they are employed in the light analysis. The fits are shown in panels (c) and (d), and the parameters are listed in Table 2.

Like the 4.2185174 d period eclipses, the 8.5303353 d period eclipses were analyzed with the Wilson-Devinney program. Due to the lack of secondary eclipse (no secondary eclipse was found after the data folding and binning), its solution bears much more uncertainty. In view of this point, the 8.5303353 d period eclipses were also analyzed by the Spot and Transit Modeling Tool (STMT; Sun et al. 2017), and the results are listed in the last column of Table 2 for comparison. The comparison shows that, generally, only the primary mean densities coincide with each other, and the differences of the other parameters reach 50%. For parameters concerned with the secondary star, such as luminosity ratio and radius ratio, their uncertainties are very large.

thumbnail Fig. 3

Light curves of KIC 10091110 with the solutions. Panel a: the whole light curves of KIC 10091110. Five large scattered points are not shown in the visible region. Panel b: a light curve part from 1220 d to 1260 d. The green and red points are the eclipses of 4.2185174 d and 8.5303353 d period, respectively. Panels c and d: the folded and binned phase curves by 4.2185174 d period around the primary and secondary minimums with the black fitting lines. Panel e: the mass ratio search curve of 4.2185174 d period binary. Panels f and g: the same as panels c and e but for 8.5303353 d period.

3.2.2 Orbital structure of KIC 10091110

Based on the available information including the light curve solutions above, we now wish to inquire into the orbital structure of KIC 10091110. Is KIC 10091110 an eclipsing triple system, or a system consisting of two irrelevant eclipsing binaries? Or what is the origin ofthe two sets of eclipses? There are three weak pieces of evidence hinting at an eclipsing triple system: (1) the third light proportions in the light curve solutions of two periods are both marginal (0.0065(48) and 0.023(11)), indicating only one luminous primary star with two faint companions; (2) the mean density figures of star 1 for the two periods are coincidentally close to each other (0.3264 is close to 0.3019 or 0.29058), indicating the primary stars of the two periods may be the same one; (3) the ratio of the two periods is very close to integer 2, that is 8.5303353/4.2185174 = 2.022. This suggests that the two orbits may be in mean-motion resonance, thus belonging to the same triple system. However, these observations are not enough to confirm the triple system structure; the third light proportion is not reliable in the light curve solution, and the same primary star density and the integer period ratio may simply be coincidences.

In the meantime, there is strong evidence, in two separate forms, against the triple system scenario: (1) no decent eclipse timing variation (ETV) was found from the two sets of eclipses. Our measurements show that the amplitudes of ETV (if there is one) will be no more than 0.008 d (for 4.2 d periods) and 0.005 d (for 8.5 d periods). Furthermore, no ETV was reported from the careful sweeping work on the entire Kepler binaries performed by Rappaport et al. (2013) and Borkovits et al. (2015, 2016) who are able to discover binaries with ETV amplitudes as small as 0.0002 d. The absence of measurable ETV is extremely unlikely for a triple system with a period ratio of only 2. If considering the 2:1 mean-motion resonance, the ETV amplitude of 4.2185174 dperiod (by adopting the parameters of Table 2) may reach a maximum of 0.5 d (by Eq. (33) in Agol et al. 2005)that significantly exceeds the detection limit; (2) a triple system with a period ratio of 2:1 cannot be stable. Equation (1) in Chambers et al. (1996) was used to calculate the stability based on the parameters of Table 2. The result shows that the triple structure is far from stable. Compared to the three pieces of evidence supporting the triple structure, these observations to the contrary are much stronger. Therefore, the two sets of eclipses cannot be generated from an eclipsing triple system, and the object KIC 10091110 cannot be a triple system.

Since KIC 10091110 cannot be a triple system, could it be a double binary system? It probably is. The evidence mentioned above is reexamined for the double binary scenario. The absence of ETV (as well as the changes in eclipse depth) and the system stability are consistent with a double binary structure. However, the tiny proportion of the third light on both binary light solutions goes against a double binary structure. Because the third light of one binary is the light of the other binary, the third light cannot be small for the two binaries. A reasonable case is that one third light proportion is small, and the other is large. Considering that the third light is not a reliable parameter in the light solution, especially for the 8.5 d period binary solution (due to the absence of secondary eclipse mentioned above), this is not strong evidence against the double binary structure. As for the close mean density of primary stars and the integer period ratio, although these two points do not contradict the double binary structure, they are amazing coincidences.

In summary, we tend to think that KIC 10091110 is a double binary system but with extraordinary coincidences: the two binaries show huge differences in luminosity, have primary stars with close mean density, the ratio of the two binaries’ periods is close to integer 2, and both have extremely small mass ratio. Whether these coincidences suggest some physical relationship between the two pairs of binaries or whether they belong to one gravitationally bound system is unclear. Further, do they originate in the same environment? Otherwise, why so many coincidences?

Table 2

Binary light curve solutions of KIC 10091110.

3.3 KIC 11495766

3.3.1 The 8.3404432 d period eclipses

The object KIC 11495766 was classified to be a planet candidate with one planet of 4.175 d period by Burke et al. (2014) and Rowe et al. (2015), so it was also named KOI-2630. However, it is suggested that the KIC 11495766 is a triple system containing a stellar binary of ~120.73 d period, and very close to it in celestial sphere there is a blended background stellar binary system with 8.3404432 d period. The background object is close enough to contaminate the light curve of KIC 11495766 in some observation quarters (one quarter is about three months in Kepler mission), because the Kepler spacecraft re-orientates itself each quarter, changing the contamination of the target, so the background object will sometimes be measured together with KIC 11495766 and sometimes not. As a result of this problem, the light curve of the background binary will show up in some quarters and disappear in others. Specifically, the 8.3404432 d period eclipses only appear in one quarter among the 15 quarters (see the blue bar in the upper panel of Fig. 4). Therefore, the 8.3404432 d period eclipses should come from a blended background source, but not from KIC 11495766 itself.

The phase curves of the 8.3404432 d period around eclipses are shown in the lower panels of Fig. 4. It can be seen that the difference between the primary and secondary eclipse is visible but not obvious, and the position of the secondary minimum (see the lower right panel) is clearly deviated from 0.5 phase. These indicate that the 8.3404432 d period eclipses are caused by a stellar binary system in a slightly eccentric orbit.

3.3.2 The extraneous eclipses of ~120.73 d period

Besides the 8.3404432 d period eclipses, seven extraneous eclipses were found and shown in Figs. 4 and 5. The red and magenta points in the upper panel of Fig. 4 represent the extraneous eclipses, which can be seen in detail in Fig. 5. The time intervals between thered eclipses (and equivalently the magenta eclipses) are approximately the same, that is, ~120.73 d. The folded phase curves by 120.73 d period8 were shown in the middle panel of Fig. 4, where all the red eclipses located at phase 1 and the magenta eclipses located around phase 0.4. The existence of both primary and secondary eclipses with sharp shape suggests a stellar binary system rather than a planet system, because in a planet system the secondary eclipse is often flat and almost invisible.

For both the left panels (red eclipses) and the right panels (magenta eclipses) in Fig. 5, the lengths of the abscissa andordinate are the same. Therefore, we can easily estimate the changes of the width and depth of the eclipses by eye. For the red eclipses, they become wider and deeper with time (from bottom to top), and become shallower and narrower at the last eclipse. For the magenta eclipses, they change to wider and shallower by several times. In addition to the changes in the eclipse shape, the changes in the periods are also dramatic. The periods, that is, the time intervals between two successive eclipses, change from 120.95 d down to 120.33 d and then up to 120.93 d for the red eclipses. Meanwhile, the time interval of the two magenta eclipses is 119.3 d, that is, 1.65 d smaller than the period 120.95 d of the nearby red eclipses. Furthermore, the two striking magenta eclipses appear and disappear suddenly after and before (at least) several cycles.

The reason for such dramatic and rapid changes could not have come from the binary itself because the changes driven by the internal factors within a binary system, such as apsidal motion, stellar evolution, and so on, are very small and take part very slowly. The eclipses and periods could never change visibly within a time span of one to two orbital periods. Therefore, the dramatic changes in the eclipses hint at strong disturbance force from (at least) an extraneous body. That is to say, the ~120.73 d period binary must belong to a triple system; a simple binary system could not change in the manner shown in Fig. 5.

Further information can be gained about the ~120.73 d period binary through the comparison between the red and magenta eclipses. Because the red eclipses are narrower and shallower than the magenta eclipses, the red (magenta) eclipses should be caused by the cooler (hotter) component blocked by the hotter (cooler) one near the periastron (apastron). The changes of the red eclipses may be the result of the movement of the eclipse point toward the periastron, and meanwhile the magenta eclipse point moves toward the apastron causing it to become wider and shallower. Since the red and magenta eclipses have large differences in depth (10 to 40 times difference), the two components should have greatly different temperatures (perhaps an approx. two-fold difference, because the fourth power of temperature ratio is proportional to the luminosity ratio). Furthermore, the difference inmass is probably large (likely more than two-fold difference, considering the stellar mass luminosity relation).

It is suggested that the ~120.73 d period binary consists of a star and a much less massive component, and around them there is another body with strong gravitational perturbation on the binary. Therefore, the object KIC 11495766 is suggested to be a triple system containing an eclipsing stellar binary with rapid changes. It should be noted that very close to KIC 11495766 in celestial sphere, there is a background object of another stellar binary system with an 8.3404432 d period.

thumbnail Fig. 4

Upper: whole light curve of KIC 11495766. The blue segment line indicates the time span when the 8.3404432 d period eclipses show up, and the red and magenta points represent the extraneous eclipses. The gray arrows with values indicate the time intervals between the eclipses. Middle: folded phase curves of the eclipses of 120.73 d period. Thered and magenta points correspond to those in the upper panel. Lower panels: folded phase curves of 8.3404432 d period eclipses at primary (left) and secondary (right) minimums, respectively; the black points are the folded points, and the blue lines are the folded and binned points.

thumbnail Fig. 5

Extraneous eclipses of KIC 11495766. Left panels: the red points are the extraneous secondary eclipses. The time intervals between the red eclipses are, or times of, ~120.73 d. Right panels: the magenta points are the extraneous primary eclipses with a time interval of 119.3 d.

4 A first list of eclipsing multiple candidates

A first list of eclipsing multiple candidates is presented in Table 3, and all the objects were found by the Kepler data. For KIC 7941635, a dramatic huge extraneous fading event was found at 420 d from a version of data from August 1, 2014, which is not present in a later version of the data. The older version of the data has been shown, by the case of KIC 7668648 and KIC 8938628 (Zhang et al. 2017), to be better at reserving authentic extraneous eclipses than the newer version data. Therefore, KIC 7941635 is also listed in Table 3. The disappeared fading is probably caused by the improved PDC which works well for the majority of Kepler targets, but inevitably degrades performance for some individual targets by removing legitimate stellar signals on them.

All the periods of the inner binaries are taken from the catalog of Kepler Eclipsing Binary Stars (Prša et al. 2011; Slawson et al. 2011; Matijevič et al. 2012; Conroy et al. 2014; Kirk et al. 2016). The numbers of extraneous eclipses were counted from our processed data (except for KIC 9632895), and most of the period change amplitudes were estimated from the O−C curves calculated by us. Some O−C curves have obvious variation, and then a general change amplitude is given; some others are very diffuse and chaotic, so only a rough upper value headed witha < symbol is given. For these latter curves, minimum time cannot be measured, and so neither for O−C values, change ranges are null in the table. The object KIC 11495766 is not an authentic eclipsing multiple system, so it is not listed here.

5 Summary and discussion

In this paper, three Kepler binaries with extraneous eclipses on the binary light curves were studied. The object KIC 5255552 is an eclipsing triple system with a rapidly changing inner binary and a hierarchical tertiary body in 862.1(±0.1) d period. The tertiary body is confirmed by both the extraneous eclipses and the O−C analysis by Borkovits et al. (2016). Four segmented inner binary light curves were analyzed independently, so the parameter changes with time can be seen quantitatively. A small suspicious eclipse at 1279 d was found, suggesting an additional possible body, making KIC 5255552 a potential quadruple eclipsing system.

The object KIC 10091110 is suggested to be a quadruple stellar system consisting of two eclipsing binaries with coincidental and interesting features. Based on the light curve analysis using a binary model, both binaries have similar extremely small mass ratios (0.060(13) and 0.0564(18)) and similar mean primary densities (0.3264(42) ρ and 0.3019(28) ρ). Not only that, the orbital period ratio of the two binaries is very close to integer 2 (8.5303353/4.2185174 = 2.022). This could easily cause one to consider orbital resonance, but this is in conflict with the lack of ETV. We wondered if these coincidences suggest some physical relationship between the two pairs of binaries. Follow-up observations, especially of the high-resolution spectra, are suggested to reveal the mystery of KIC 10091110.

Excluding the already known 8.3404432 d period eclipses on KIC 11495766, seven extraneous eclipses with a period of ~120.73 d were discovered. In this paper, the ~120.73 d period eclipses are interpreted as being generated by KIC 11495766, and the 8.3404432 d period eclipses were deemed to be caused by a blended background source near KIC 11495766. It cannot be overlooked that the ~120.73 d eclipses exhibit fast and dramatic changes which cannot be explained by the physics of binary itself. Therefore, there must be additional bodies responsible for the eclipse changes, and the object KIC 11495766 must be (at least) a triple system, despite lacking outer eclipses.

A list of 25 eclipsing multiple candidates is presented for the first time. All the 25 eclipsing multiple candidates come from the Kepler observations, and they are perhaps the most reliable candidates so far. The eclipsing multiple systems are not only the multiple systems, their extraneous eclipses provide a new important resource for research on multiple systems. This list could be greatly expanded by future projects similar to the Kepler survey, and we hope this first small list could be beneficial for studying eclipsing multiples.

Table 3

List of 25 eclipsing multiple candidates.

Acknowledgements

We thank the anonymous referee for constructive criticisms and recommendations which have significantly improved this manuscript. This paper includes data collected by the Kepler mission. Funding for the Kepler mission is provided by the NASA Science Mission directorate. We are indebted to the Kepler Science Working Group on Eclipsing Binaries for all their great pilot works which made these discoveries possible. This work is partly supported by Chinese Natural Science Foundation (Nos. 11703082, 11133007, 11325315, 11403056), the West Light Foundation of The Chinese Academy of Sciences, Yunnan Natural Science Foundation (Y5XB071001), Guangdong Provincial Engineering Technology Research Center for Data Science, the Joint Research Fund in Astronomy (grant number U1631108) under cooperative agreement between the National Natural Science Foundation of China (NSFC) and Chinese Academy of Sciences (CAS), and the research fund of Sichuan University of Science and Engineering (grant number 2015RC42), the Key Research Program of the Chinese Academy of Sciences (grant No. KGZD-EW-603), the Science Foundation of Yunnan Province (No. 2012HC011), and the Strategic Priority Research Program “The Emergence of Cosmological Structures” of the Chinese Academy of Sciences (No. XDB09010202).

References

  1. Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567 [NASA ADS] [CrossRef] [Google Scholar]
  2. Borkovits, T., Rappaport, S., Hajdu, T., & Sztakovics, J. 2015, MNRAS, 448, 946 [NASA ADS] [CrossRef] [Google Scholar]
  3. Borkovits, T., Hajdu, T., Sztakovics, J., et al. 2016, MNRAS, 455, 4136 [NASA ADS] [CrossRef] [Google Scholar]
  4. Burke, C. J., Bryson, S. T., Mullally, F., et al. 2014, ApJS, 210, 19 [NASA ADS] [CrossRef] [Google Scholar]
  5. Carter, J. A., Fabrycky, D. C., Ragozzine, D., et al. 2011, Science, 331, 562 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  6. Chambers, J. E., Wetherill, G. W., & Boss, A. P. 1996, Icarus, 119, 261 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  7. Conroy, K. E., Prša, A., Stassun, K. G., et al. 2014, AJ, 147, 45 [NASA ADS] [CrossRef] [Google Scholar]
  8. Conroy, K., Prša, A., Stassun, K., & Orosz, J. 2015, IBVS, 6138, 1 [NASA ADS] [Google Scholar]
  9. Derekas, A., Kiss, L. L., Borkovits, T., et al. 2011, Science, 332, 216 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  10. Fraquelli, D., & Thompson, S. E., 2014, Revision 5: Kepler Archive Manual (KDMC-10008-005) [Google Scholar]
  11. Jetsu, L., & Porceddu, S. 2015, PLoS ONE, 10 [Google Scholar]
  12. Jetsu, L., Porceddu, S., Lyytinen, J., et al. 2013, ApJ, 773, 1 [NASA ADS] [CrossRef] [Google Scholar]
  13. Kirk, B., Conroy, K., Prsa, A., et al. 2016, AJ, 151, 68 [NASA ADS] [CrossRef] [Google Scholar]
  14. Kostov, V. B., McCullough, P. R., Hinse, T. C., et al. 2013, ApJ, 770, 52 [NASA ADS] [CrossRef] [Google Scholar]
  15. Kostov, V. B., McCullough, P. R., Carter, J. A., et al. 2014, ApJ, 784, 14 [NASA ADS] [CrossRef] [Google Scholar]
  16. Lehmann, H., Zechmeister, M., Dreizler, S., Schuh, S., & Kanzler, R. 2012, A&A, 541, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Marsh, T. R., Armstrong, D. J., & Carter, P. J. 2014, MNRAS, 445, 309 [NASA ADS] [CrossRef] [Google Scholar]
  18. Masuda, K., Uehara, S., & Kawahara, H. 2015, ApJ, 806, L37 [NASA ADS] [CrossRef] [Google Scholar]
  19. Matijevi c, G., Prša, A., Orosz, J. A., et al. 2012, AJ, 143, 123 [NASA ADS] [CrossRef] [Google Scholar]
  20. Orosz, J. A., Welsh, W. F., Carter, J. A., et al. 2012a, ApJ, 758, 87 [NASA ADS] [CrossRef] [Google Scholar]
  21. Orosz, J. A., Welsh, W. F., Carter, J. A., et al. 2012b, Science, 337, 1511 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  22. Prša, A., Batalha, N., Slawson, R. W., et al. 2011, AJ, 141, 83 [NASA ADS] [CrossRef] [Google Scholar]
  23. Rappaport, S., Deck, K., Levine, A., et al. 2013, ApJ, 768, 33 [NASA ADS] [CrossRef] [Google Scholar]
  24. Rowe, J. F., Coughlin, J. L., Antoci, V., et al. 2015, ApJS, 217, 16 [NASA ADS] [CrossRef] [Google Scholar]
  25. Shibahashi, H., & Kurtz, D. W. 2012, MNRAS, 422, 738 [NASA ADS] [CrossRef] [Google Scholar]
  26. Slawson, R. W., Prša, A., Welsh, W. F., et al. 2011, AJ, 142, 160 [NASA ADS] [CrossRef] [Google Scholar]
  27. Sun, L., Gu, S., Wang, X., et al. 2017, AJ, 153, 28 [NASA ADS] [CrossRef] [Google Scholar]
  28. Van Hamme, W., & Wilson, R. E. 2007, ApJ, 661, 1129 [Google Scholar]
  29. Welsh, W. F., Orosz, J. A., Carter, J. A., et al. 2012, Nature, 481, 475 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  30. Welsh, W. F., Orosz, J. A., Short, D. R., et al. 2015, ApJ, 809, 26 [NASA ADS] [CrossRef] [Google Scholar]
  31. Wilson, R. E. 1979, ApJ, 234, 1054 [NASA ADS] [CrossRef] [Google Scholar]
  32. Wilson, R. E. 1990, ApJ, 356, 613 [NASA ADS] [CrossRef] [Google Scholar]
  33. Wilson, R. E. 2008, ApJ, 672, 575 [NASA ADS] [CrossRef] [Google Scholar]
  34. Wilson, R. E. 2012, AJ, 144, 73 [NASA ADS] [CrossRef] [Google Scholar]
  35. Wilson, R. E., & Devinney, E. J. 1971, ApJ, 166, 605 [NASA ADS] [CrossRef] [Google Scholar]
  36. Wilson, R. E., Van Hamme, W., & Terrell, D. 2010, ApJ, 723, 1469 [NASA ADS] [CrossRef] [Google Scholar]
  37. Zhang, J., Qian, S.-B., & He, J.-D. 2017, Res. Astron. Astrophys., 17, 22 [NASA ADS] [CrossRef] [Google Scholar]

1

Actually a triple star system. It is also the first discovered variable star. Its periodic light changes were recorded before 1163 B.C. by ancient Egyptians (Jetsu et al. 2013; Jetsu & Porceddu 2015). The unrecorded discovery may be even earlier as it is the most easily discovered variable to the naked eye.

4

The mean reference minimum time is calculated in the following way: 1, calculate the phase and fold them to phase 0.5 to 1.5 with an initial reference minimum time; 2, fit the minimum at phase 1 to obtain its precise minimum phase; 3, the difference between this phase value and 1 is used to correct the reference minimum time; 4, iterate these previous three steps until the difference is small enough.

5

The light analysis was carried out by setting the shallower eclipses fixed at phase zero, and so the position of the deeper eclipses vary with time.

6

Unlike the case of circular orbit, in an eccentric orbit like KIC 5255552, the deeper eclipse is not necessary generated by the massive component being blocked by its companion. It is possible that the deeper eclipse is generated by that of the less massive component. So this need to be clarified first: based on the time of periastron passage and the argument of periastron of the less massive component provided by Borkovits et al. (2015), it can be worked out that for KIC 5255552, the deeper eclipse corresponds to the massive component being blocked by its companion.

7

It should be noted that, in panel (a), the deep eclipses at around 530 d, 920 d and 1300 d are not an independent set of eclipses, but are only the superposition of the two sets of eclipses.

8

The phase curves were calculated only from the data larger than 800 d where the red eclipses exists.

All Tables

Table 1

Binary light curve solutions of KIC 5255552.

Table 2

Binary light curve solutions of KIC 10091110.

Table 3

List of 25 eclipsing multiple candidates.

All Figures

thumbnail Fig. 1

Light curves of KIC 5255552. Upper: the whole light curves. The red points are the extraneous eclipses that can been seen more clearly in Fig. 2, and the gray arrows are used to illustrate the time intervals between the extraneous eclipses. Lower eight panels: the four selected light curve parts in phase form, showing the eclipses at phase zero (left) and phase around 0.43 (right) with their fittings (black line).

In the text
thumbnail Fig. 2

Extraneous eclipses of KIC 5255552. The red points in all the panels are the extraneous eclipses. The mean times of the extraneous eclipses for each group are marked in each panel. In the top panel, a small additional panel is used to magnify a tiny extraneous eclipse around 1537.7 d.

In the text
thumbnail Fig. 3

Light curves of KIC 10091110 with the solutions. Panel a: the whole light curves of KIC 10091110. Five large scattered points are not shown in the visible region. Panel b: a light curve part from 1220 d to 1260 d. The green and red points are the eclipses of 4.2185174 d and 8.5303353 d period, respectively. Panels c and d: the folded and binned phase curves by 4.2185174 d period around the primary and secondary minimums with the black fitting lines. Panel e: the mass ratio search curve of 4.2185174 d period binary. Panels f and g: the same as panels c and e but for 8.5303353 d period.

In the text
thumbnail Fig. 4

Upper: whole light curve of KIC 11495766. The blue segment line indicates the time span when the 8.3404432 d period eclipses show up, and the red and magenta points represent the extraneous eclipses. The gray arrows with values indicate the time intervals between the eclipses. Middle: folded phase curves of the eclipses of 120.73 d period. Thered and magenta points correspond to those in the upper panel. Lower panels: folded phase curves of 8.3404432 d period eclipses at primary (left) and secondary (right) minimums, respectively; the black points are the folded points, and the blue lines are the folded and binned points.

In the text
thumbnail Fig. 5

Extraneous eclipses of KIC 11495766. Left panels: the red points are the extraneous secondary eclipses. The time intervals between the red eclipses are, or times of, ~120.73 d. Right panels: the magenta points are the extraneous primary eclipses with a time interval of 119.3 d.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.