Aims: Our goal was to increase number of known doubly eclipsing systems such that the resulting dataset would allow to study them via statistical means, as well as prove that they constitute gravitationally bound 2+2 quadruple system. Methods: We analysed photometric data for eclipsing binaries provided by the OGLE survey in the LMC fields. We found a large number of new doubly eclipsing systems (our discoveries are 3x more numerous than previous studies). With a typical orbital period of days for the binaries, we sought eclipse time variations (ETVs) on the timescale of years. In the cases where we were able to detect the ETV period, the difference between the inner and outer periods in the quadruple system is large enough. This allows us to interpret ETVs primarily as the light-time effect, thus providing an interesting constraint on masses of the binaries. Results: In addition to significantly enlarging the database of known doubly eclipsing systems, we performed a thorough analysis of 72 cases. ETVs for 28 of them (39% of the studied cases) showed evidence of relative motion. We note OGLE BLG-ECL-145467 as the most interesting case; it is bright (I=12.6 mag), consists of two detached binaries with periods of about 3.3 d and 4.9 d (making it a candidate for a 3:2 resonant system), mutual period about 1538 d. Distribution of the orbital period ratio P_A/P_B of binaries in 2+2 quadruples shows statistically significant excess at 1 and 1.5. The former is likely a natural statistical preference in weakly interacting systems with periods within the same range. The latter is thought to be evidence of a capture in the 3:2 mean motion resonance of the two binaries. This sets important constraints on evolutionary channels in these systems. The total number of doubly eclipsing systems increased to 146, more than 90% of which are at low declinations on the southern sky.

The eclipsing binaries are still important objects for our understanding of the universe. Especially these ones located within the more complex multiple systems can help us solving the problem of their origin and subsequent evolution of these higher order multiples. Photometry and spectroscopy spanning over more than 25 yr were used for the first complete analysis of the two bright triple systems, namely V348 And and V572 Per. The light curves in photometric filters were combined together with the radial velocities and analyzed simultaneously, yielding the precise physical parameters of the eclipsing components of these multiple systems. The system V348 And consists of two eclipsing components with its orbital period of about 27.7 days. The system is a very detached one, and both eclipses are rather narrow, lasting only about 0.016 of its period. The visual orbit of the wide pair has the period of about 87 yr. All three components of the system are probably of B8-9 spectral type, and the parallax of the system was slightly shifted to the value of 2.92 mas. On the other hand, the system V572 Per shows apsidal motion of its inner orbit, the orbital period being of about 1.2 days, while the apsidal motion of about 48 yr. The components are of A and F spectral types, while the motion with the third component around a common barycenter is only negligible. According to our modeling, this system is not a member of open cluster Alpha Persei.

We present an analysis of the apsidal motion and light curve parameters of 54 never-before-studied galactic Algol-type binaries. This is the first analysis of such a large sample of eccentric eclipsing binaries in our Galaxy, and has enabled us to identify several systems that are worthy of further study. Bringing together data from various databases and surveys, supplemented with new observations, we have been able to trace the long-term evolution of the eccentric orbit over durations extending back up to several decades. Our present study explores a rather different sample of stars to those presented in the previously published catalogue of eccentric eclipsing binaries, sampling to fainter magnitudes, covering later spectral types, sensitive to different orbital periods with more than 50% of our systems having periods longer than six days. The typical apsidal motion in the sample is rather slow (mostly of order of centuries long), although in some cases this is less than 50 yr. All of the systems, except one, have eccentricities less than 0.5, with an average value of 0.23. Several of the stars also show evidence for additional period variability. In particular we can identify three systems in the sample, HD 44093, V611 Pup, and HD 313631, which likely represent relativistic apsidal rotators.

Eclipsing binaries remain crucial objects for our understanding of the universe. In particular, those that are components of multiple systems can help us solve the problem of the formation of these systems. Analysis of the radial velocities together with the light curve produced for the first time precise physical parameters of the components of the multiple systems V773 Cas, QS Aql, and BR Ind. Their visual orbits were also analyzed, which resulted in slightly improved orbital elements. What is typical for all these systems is that their most dominant source is the third distant component. The system V773 Cas consists of two similar G1-2V stars revolving in a circular orbit and a more distant component of the A3V type. Additionally, the improved value of parallax was calculated to be 17.6 mas. Analysis of QS Aql resulted in the following: the inner eclipsing pair is composed of B6V and F1V stars, and the third component is of about the B6 spectral type. The outer orbit has high eccentricity of about 0.95, and observations near its upcoming periastron passage between the years 2038 and 2040 are of high importance. Also, the parallax of the system was derived to be about 2.89 mas, moving the star much closer to the Sun than originally assumed. The system BR Ind was found to be a quadruple star consisting of two eclipsing K dwarfs orbiting each other with a period of 1.786 days; the distant component is a single-lined spectroscopic binary with an orbital period of about 6 days. Both pairs are moving around each other on their 148 year orbit.

Aims: We present a new updated study of the unique quintuple system of two eclipsing binaries, V994 Her. Methods: Based on new obtained photometric observations, we analysed the period variations of the two eclipsing pairs of V994 Her. Results: We found that the orbital period of the two eclipsing binaries around a common barycenter is about 2.9 yrs instead of the original period, which was reported as 6.3 yrs. Moreover, the system now seems to be close to coplanar (i ~ 84deg), and we also revised the individual masses of all four eclipsing components. The two eclipsing pairs show slow apsidal motion with periods of about 116 and 111 yrs, respectively. Pairs A and B are distant by about 17 mas on the sky, which is encouraging for observers because an interferometric detection may be possible in the next years. Conclusions: All of these results are based on the photometric observations of V994 Her and its detailed analysis of the period changes of the two eclipsing binaries.

The systems V819 Her, V2388 Oph, and V1031 Ori are triples comprised of an eclipsing binary orbiting with a distant visual component on a longer orbit. A detailed analysis of these interesting systems, combining the two observational techniques: interferometry and apparent period variation, was performed. The interferometric data for these three systems obtained during the last century determine the visual orbits of the distant components in the systems. The combined analysis of the positional measurements together with the analysis of apparent period changes of the eclipsing binary (using the minima timings) can be used to study these systems in a combined approach, resulting in a set of parameters otherwise unobtainable without the radial velocities. The main advantage of the technique presented here is the fact that one needs no spectroscopic monitoring of the visual orbits, which have rather long periods: 5.5 yr for V819 Her, 9.0 yr for V2388 Oph, and 31.3 yr for V1031 Ori, respectively. The eccentricities of the outer orbits are 0.69, 0.33, and 0.92, respectively. Moreover, the use of minima timings of the eclipsing pairs help us to derive the orientation of the orbit in space with no ambiguity around the celestial plane. And finally, the combined analysis yielded also an independent determination of the distance of V819 Her (68.7 +/- 1.8 pc), and V2388 Oph (70.6 +/- 8.9 pc). We also present a list of similar systems, which would be suitable for a combined analysis like this one.

The minimum timings of eclipsing binaries V418 Aql, SU Boo, RV CVn, CR Cas, GV Cyg, V432 Per, and BD+42 2782 were collected and analyzed. Their long-term behavior was studied via period analysis, revealing a periodic term in eclipse times. We derived 576 new times of minimum. Hence, to describe the periodic variation, a thirdbody hypothesis was proposed and the resulting orbital periods are as follows: 70, 7.4, 53, 37, 27, 53, and 18 yr, respectively. For the system V432 Per an additional 9.5 yr variation was also found. The predicted minimum masses of these distant bodies were calculated and their detectability discussed. The light curves of SU Boo and RV CVn were analyzed using the PHOEBE program, resulting in physical parameters of the components. New variable stars in the field of V418 Aql were discovered.

More than three hundred years ago, Ole Romer measured the speed of light only by observing the periodic shifting of the observed eclipse arrival times of Jupiter's moons arising from the varying Earth-Jupiter distance. The same method of measuring the periodic modulation of delays is still used in astrophysics. The ideal laboratories for this effect are eclipsing binaries. The unique system V994 Her consists of two eclipsing binaries orbiting each other. However, until now it was not certain whether these are gravitationally bound and what their orbital period is. We show that the system is in fact quintuple and the two eclipsing binaries are orbiting each other with period about 6.3 years. This analysis was made only from studying the periodic modulation of the two periods, when during the periastron passage one binary has an apparently shorter period, while the other one longer, exactly as required by a theory. Additionally, it was found that both inner eclipsing pairs orbit with slightly eccentric orbits undergoing a slow apsidal motions with a period of the order of centuries.

A study of late-type low-mass eclipsing binaries provides us with important information about the most common stars in the Universe. We obtain the first light curves and perform period analyses of two neglected eclipsing binaries GK Boo and AE For to reveal their basic physical properties. We performed both a period analysis of the times of the minima and a BVR light curve analysis. Many new times of minima for both the systems were derived and collected from the data obtained by automatic and robotic telescopes. This allowed us to study the long-term period changes in these systems for the first time. From the light curve analysis, we derived the first rough estimates of the physical properties of these systems. We find that the analyzed systems are somewhat similar to each other. Both contain low-mass components of similar types, both are close to the Sun, both have short orbital period, and both contain another low-mass companions on longer orbits of a few years. In the case of GK Boo, both components are probably of K3 spectral type, while the distant companion is probably a late M star. The light curve of GK Boo is asymmetric, which probably causes the shift in the secondary minima in the O-C diagram. System AE For comprises two K7 stars, and the third body is a possible brown dwarf with a minimal mass of only about 47 Jupiter Mass. We succeed in completing period and light curve analyses of both systems, although a more detailed spectroscopic analysis is needed to confirm the physical parameters of the components to a higher accuracy.

Aims: We present the detailed analysis of triple system KR Com with different observational techniques - photometry, interferometry, and period variation. Methods: The use of BVR photometry of the close-contact binary KR Com, which is the primary component of a triple system, helps us to better describe the properties of the components. The interferometric data obtained during the last 30 years sufficiently determine the visual orbit, but the use of minima timings of KR Com for the study of period variation together with the visual orbit is a novel approach in this system. Results: Basic physical parameters resulting from the light curve analysis agree well with the previous results from spectroscopy. The temperatures for the primary and secondary component resulted in 5549 and 6072 K, respectively, and the amount of the third light in all filters is about 1/3 of the total luminosity. The distant third component revolves around the common barycenter on 11 yr orbit with a very high eccentricity (0.934) and this movement is also detectable via the period variation, which is clearly visible in the O-C diagram of times of minima observations. The use of minima times for the combined analysis helps us to independently determine the distance to the system (64.02 +/- 9.42 pc) and also to confirm the orientation of the orbit in space. Conclusions: New minima observations and also spectroscopy would be very profitable, especially during the next periastron passage in the year 2017.

The multiple system V505 Sagittarii is composed of at least three stars: a compact eclipsing pair and a distant component, which orbit is measured directly using speckle interferometry. In order to explain the observed orbit of the third body in V505 Sagittarii and also other observable quantities, namely the minima timings of the eclipsing binary and two different radial velocities in the spectrum, we thoroughly test a fourth-body hypothesis - a perturbation by a dim, yet-unobserved object. We use an N-body numerical integrator to simulate future and past orbital evolution of 3 or 4 components in this system. We construct a suitable chi^2 metric from all available speckle-interferometry, minima-timings and radial-velocity data and we scan a part of a parameter space to get at least some of allowed solutions. In principle, we are able to explain all observable quantities by a presence of a fourth body, but the resulting likelihood of this hypothesis is very low. We also discuss other theoretical explanations of the minima timings variations. Further observations of the minima timings during the next decade or high-resolution spectroscopic data can significantly constrain the model.

A new catalog of visual double systems containing eclipsing binaries as one component is presented. The main purpose of this catalog is to compile a complete list of all known multiples of this variety, both for current analysis and to highlight those in need of additional observations. All available photometric and astrometric data were analyzed, resulting in new orbits for eight systems and new times of minimum light for a number of the eclipsing binaries. Some of the systems in the catalog have acceptable solutions for their visual orbits, although in most cases their orbital periods are too long for simultaneous analysis. Also included, however, are a number of systems which currently lack an orbital solution but which may be suitable for simultaneous analysis in the future.