We study the apsidal motion in close eccentric massive binaries. Measuring the rate of apsidal motion in such a system gives insight into the internal structure and evolutionary state of the stars. We focus on CPD-41° 7742, for which independent studies in the past showed large discrepancies in the longitude of periastron of the orbit, hinting at the presence of apsidal motion. We perform a consistent analysis of all observational data to solve this apparent discrepancy and report the first determination of apsidal motion in this system. This study confirms the need for enhanced mixing in the stellar evolution models of the primary star to reproduce the observational properties. This points towards larger convective cores than usually considered.

The young open cluster NGC 6231 hosts a rich population of O-type binary stars. We study several of these eccentric short-period massive eclipsing binaries and assess their fundamental parameters. The properties of these systems make them interesting targets to study tidally induced apsidal motion. The analysis of apsidal motion offers a powerful means to obtain information about the internal structure of the stars. Indeed, since the rate of apsidal motion in a binary system is proportional to the internal structure constants of the stars composing it, its value gives direct insight into the internal structure and evolutionary state of these stars. Stellar evolution models are constructed based on the observationally-determined fundamental parameters and a theoretical rate of apsidal motion is inferred. The results are striking: Adopting standard stellar evolution models yields a theoretical rate of apsidal motion much larger than the observational value. This discrepancy results from the standard models predicting too low an efficiency of internal mixing and thus too homogenous stars in terms of density. By enforcing the theoretical rates of apsidal motion to match the observational values, enhanced mixing is required, through a large overshooting parameter and/or additional turbulent/rotational mixing. Our analysis leads to the conclusion that the chemically mixed cores in those massive stars must be more extended than anticipated from standard models.

The presence of additional bodies orbiting a binary star system can be detected by monitoring the binary's eclipse timing. These so-called circumbinary objects will reveal themselves by i) either introducing a reflex motion of the binary system about the total system's barycenter creating a geometric light-travel time effect (LITE), ii) by gravitational perturbations on the binary orbit, or iii) a combination of the two effects resulting in eclipse timing (ETV) and transit timing (TTV) variations. Motivated by the four recently detected circumbinary planets by the Kepler space telescope (Kepler-16b, Kepler-34b, Kepler-35b and Kepler-38b) we have begun to study their dynamics in the presence of an additional massive perturber. In particular we used Kepler-16b as a test bed case. We are aiming to study the detectability of non-transiting and inclined circumbinary planets using the ETV effect along with the fast indicator MEGNO to quantitatively map the chaotic/quasi-periodic regions of the orbital parameter-space and to determine where the orbit of a circumbinary planet will be stable. We have calculated the amplitudes of TTV and ETV signals for different values of the mass and orbital elements of the planet and binary.

Among the longest known eclipse durations and binary periods is that of the star ϵ Aurigae, which exhibits 2-year long eclipses every 27.1 years. Oddly, the nature of the secondary in the system continues to elude ready identification. In 1965, Huang proposed a massive disk as the eclipsing body, and study of the 1984 eclipse led Lissauer and Backman to suggest an embedded B star binary in the disk to maintain it. A collaboration of observers allows me to present recent optical photometry and spectroscopy, near-IR spectroscopy and Spitzer space telescope IRS and MIPS observations of ϵ Aurigae as it approaches its next eclipse. These data argue for current detectability of the embedded binary, and precession of the disk axis, suggesting a radical change is possible for the next mid-eclipse brightening. An international monitoring campaign for the 2009-2011 is being organized, and participation invited via website http://www.du.edu/~rstencel/epsaur.htm.

One of the most important changes in observational astronomy of the 21st Century is a rapid shift from classical object-by-object observations to extensive automatic surveys. As CCD detectors are getting better and their prices are getting lower, more and more small and medium-size observatories are refocusing their attention to detection of stellar variability through systematic sky-scanning missions. This trend is aditionally powered by the success of pioneering surveys such as ASAS, DENIS, OGLE, TASS, their space counterpart Hipparcos and others. Such surveys produce massive amounts of data and it is not at all clear how these data are to be reduced and analysed. This is especially striking in the eclipsing binary (EB) field, where most frequently used tools are optimized for object-by-object analysis. A clear need for thorough, reliable and fully automated approaches to modeling and analysis of EB data is thus obvious. This task is very difficult because of limited data quality, non-uniform phase coverage and solution degeneracy. This paper reviews recent advancements in putting together semi-automatic and fully automatic pipelines for EB data processing. Automatic procedures have already been used to process Hipparcos data, LMC/SMC observations, OGLE and ASAS catalogs etc. We discuss the advantages and shortcomings of these procedures.