No Arabic abstract
We present the first results of a Kepler survey of 41 eclipsing binaries that we undertook to search for third star companions. Such tertiaries will periodically alter the eclipse timings through light travel time and dynamical effects. We discuss the prevalence of starspots and pulsation among these binaries and how these phenomena influence the eclipse times. There is no evidence of short period companions (P < 700 d) among this sample, but we do find evidence for long term timing variations in 14 targets (34%). We argue that this finding is consistent with the presence of tertiary companions among a significant fraction of the targets, especially if many have orbits measured in decades. This result supports the idea that the formation of close binaries involves the deposition of angular momentum into the orbital motion of a third star.
We present a hierarchical triple star system (KIC 9140402) where a low mass eclipsing binary orbits a more massive third star. The orbital period of the binary (4.98829 Days) is determined by the eclipse times seen in photometry from NASAs Kepler spacecraft. The periodically changing tidal field, due to the eccentric orbit of the binary about the tertiary, causes a change in the orbital period of the binary. The resulting eclipse timing variations provide insight into the dynamics and architecture of this system and allow the inference of the total mass of the binary ($0.424 pm 0.017 text{M}_odot$) and the orbital parameters of the binary about the central star.
We report eclipse timing variation analyses of 26 compact hierarchical triple stars comprised of an eccentric eclipsing (inner) binary and a relatively close tertiary component found in the {em Kepler} field. We simultaneously fit the primary and secondary $O-C$ curves of each system for the light-travel time effect (LTTE), as well as dynamical perturbations caused by the tertiary on different timescales. For the first time, we include those contributions of three-body interactions which originate from the eccentric nature of the inner binary. These effects manifest themselves both on the period of the triple system, $P_2$, and on the longer apse-node timescale. We demonstrate that consideration of the dynamically forced rapid apsidal motion yields an efficient and independent tool for the determination of the binary orbits eccentricity and orientation, as well as the 3D configuration of the triple. Modeling the forced apsidal motion also helps to resolve the degeneracy between the shapes of the LTTE and the dynamical delay terms on the $P_2$ timescale, due to the strong dependence of the apsidal motion period on the triples mass ratio. This can lead to the independent determination of the binary and tertiary masses without the need for independent radial velocity measurements. Through the use of our analytic method for fitting $O-C$ curves we have obtained robust solutions for system parameters for the ten most ideal triples of our sample, and only somewhat less robust, but yet acceptable, fits for the remaining systems. Finally we study the results of our 26 system parameter fits via a set of distributions of various physically important parameters, including mutual inclination angle, and mass and period ratios.
We utilize multi-dimensional simulations of varying equatorial jet strength to predict wavelength dependent variations in the eclipse times of gas-giant planets. A displaced hot-spot introduces an asymmetry in the secondary eclipse light curve that manifests itself as a measured offset in the timing of the center of eclipse. A multi-wavelength observation of secondary eclipse, one probing the timing of barycentric eclipse at short wavelengths and another probing at longer wavelengths, will reveal the longitudinal displacement of the hot-spot and break the degeneracy between this effect and that associated with the asymmetry due to an eccentric orbit. The effect of time offsets was first explored in the IRAC wavebands by Williams et. al (2006). Here we improve upon their methodology, extend to a broad ranges of wavelengths, and demonstrate our technique on a series of multi-dimensional radiative-hydrodynamical simulations of HD 209458b with varying equatorial jet strength and hot-spot displacement. Simulations with the largest hot-spot displacement result in timing offsets of up to 100 seconds in the infrared. Though we utilize a particular radiative hydrodynamical model to demonstrate this effect, the technique is model independent. This technique should allow a much larger survey of hot-spot displacements with JWST then currently accessible with time-intensive phase curves, hopefully shedding light on the physical mechanisms associated with thermal energy advection in irradiated gas-giants.
The Kwee - van Woerden (KvW) method used for the determination of eclipse minimum times has been a staple in eclipsing binary research for decades, due its simplicity and the independence of external input parameters, which also makes it well-suited to obtaining timings of exoplanet transits. However, its estimates of the timing error have been known to have a low reliability. During the analysis of very precise photometry of CM Draconis eclipses from TESS space mission data, KvWs original equation for the timing error estimate produced numerical errors, which evidenced a fundamental problem in this equation. This contribution introduces an improved approach for calculating the timing error with the KvW method. A code that implements this improved method, together with several further updates of the original method, are presented. An example of the application to CM Draconis light curves from TESS is given. The eclipse minimum times are derived with the KvW methods three original light curve folds, but also with five and seven folds. The use of five or more folds produces minimum timings with a substantially better precision. The improved method of error calculation delivers consistent timing errors which are in excellent agreement with error estimates obtained by other means. In the case of TESS data from CM Draconis, minimum times with an average precision of 1.1 seconds are obtained. Reliable timing errors are also a valuable indicator for evaluating if a given scatter in an O-C diagram is caused by measurement errors or by a physical period variation.
The eclipsing white dwarf plus main-sequence binary NN Serpentis provides one of the most convincing cases for the existence of circumbinary planets around evolved binaries. The exquisite timing precision provided by the deep eclipse of the white dwarf has revealed complex variations in the eclipse arrival times over the last few decades. These variations have been interpreted as the influence of two planets in orbit around the binary. Recent studies have proved that such a system is dynamically stable over the current lifetime of the binary. However, the existence of such planets is by no means proven and several alternative mechanisms have been proposed that could drive similar variations. One of these is apsidal precession, which causes the eclipse times of eccentric binaries to vary sinusoidally on many year timescales. In this paper we present timing data for the secondary eclipse of NN Ser and show that they follow the same trend seen in the primary eclipse times, ruling out apsidal precession as a possible cause for the variations. This result leaves no alternatives to the planetary interpretation for the observed period variations, although we still do not consider their existence as proven. Our data limits the eccentricity of NN Ser to e<0.001. We also detect a 3.3+/-1.0 second delay in the arrival times of the secondary eclipses relative to the best planetary model. This delay is consistent with the expected 2.84+/-0.04 second Romer delay of the binary, and is the first time this effect has been detected in a white dwarf plus M dwarf system.