arxiv: v1 [astro-ph.im] 17 Apr 2012

Transcription

1 Mon. Not. R. Astron. Soc., 1 7 (21) Printed 23 November 218 (MN LATEX style file v2.2) Radio Pulsar Style Timing of Eclipsing Binary Stars from the ASAS Catalogue arxiv: v1 [astro-ph.im] 17 Apr 212 S. K. Koz lowski 1,2, M. Konacki 1,2 and P. Sybilski 2 1 Adam Mickiewicz University, Astronomical Observatory, Poznań, Poland 2 Nicolaus Copernicus Astronomical Center, Toruń, Poland Received... 1 INTRODUCTION The fact that the velocity of light is finite was not obvious till 1676 when Olaus Roemer carried out precise measurementsofthetimes ofeclipses ofjovianmoons. Henotedthat Io eclipses were early before opposition and late after opposition when compared to the Ephemerides Bononiensis Mediceorum Siderum, a work by Cassini published in It includes tables of times of eclipses of Jovian moons which were used to determine the differential longitude by simultaneous observations of the same eclipse from two places. Roemer s conclusion, though not a quantitative one, became a great discovery contradicting the Aristotelean thought. He provided scientists with the basics of the O C (observed minus calculated) procedure (Sterken 25b) and was the first one to analyze the effects caused by finite light speed, hereafter called the light time effect (LTE). Below we analyze the photometric data from the All Sky Automated Survey, (ASAS; Pojmanski 22). In 2 we present our method for analyzing the timing variations. In 3 stan@ncac.torun.pl maciej@ncac.torun.pl sybilski@ncac.torun.pl ABSTRACT The Light-Time Effect (LTE) is observed whenever the distance between the observer and any kind of periodic event changes in time. The usual cause of this distance change is the reflex motion about the system s barycenter due to the gravitational influence of one or more additional bodies. We analyze 532 eclipsing contact (EC) and detached (ED) binaries from the All Sky Automated Survey (ASAS) catalogue to detect variations in the times of eclipses which possible can be due to the LTE effect. To this end we use an approach known from the radio pulsar timing where a template radio pulse of a pulsar is used as a reference to measure the times of arrivals of the collected pulses. In our analysis as a template for a photometric time series from ASAS, we use a best-fitting trigonometric series representing the light curve of a given EC or ED. Subsequently, an O C diagram is built by comparing the template light curve with light curves obtained from subsets of a given time series. Most of the variations we detected in O Cs correspond to a linear period change. Three show evidence of more than one complete LTE-orbit. For these objects we obtained preliminary orbital solutions. Our results demonstrate that the timing analysis employed in radio pulsar timing can be effectively used to study large data sets from photometric surveys. Key words: binaries: eclipsing methods: numerical. we show the outcome of applying our approach to the photometric series of 532 eclipsing contact (EC) and detached (ED) binaries from ASAS. In 4 we discuss several interesting cases of most likely the LTE effect due to companions to the analyzed systems and conclude in 5. 2 AUTOMATED TIMING OF ECLIPSING BINARIES Our basic concept of detecting timing variations in photometry of ED/EC in an automated way is based on the method used in radio pulsar timing. It consists of six steps which are shown as a block diagram in Fig. 1. (i) Get a raw data set. It is assumed that a raw data set consists of n magnitudes (or fluxes), their errors and the times they were recorded. (ii) Prepare a template light-curve model. A raw datasetisphasedwiththeknownperiodandtheparameters of a template light-curve model are calculated using the least squares method. At this stage the period can be improved or corrected during the fitting process. This is often necessary for the ASAS data. (iii) Divide a raw data set into subsets. A raw data

2 2 S. K. Kozlowski, M. Konacki and P. Sybilski V[mag] phase Figure 2. An example template light-curve model (solid line) calculated for ASAS raw data. Note that the solid curve approximates the primary minimum very well despite oscillations visible between eclipses. fitting a trigonometric series to a raw photometric data set: f (t) = Figure 1. A block diagram describing our eclipse timing procedure. set is divided into M intervals. The intervals can be equal in terms of the number of data points they include or the time that they span. The second variant was chosen in this paper. The intervals can be overlapping or have a non-repeating content. When implementing the first method, appropriate corrections must be applied when calculating the final formal errors due to a multiple usage of the same data points. (iv) Phase data in each interval. Data points are phased separately in each interval using the new period value calculated during the creation of the model. The zeropoint is retained for each interval. This way a local light-curve is created and the mid-time of each interval is associated with it. (v) Compare with the template light-curve model. The most important step in this procedure is the comparison of the local light-curves with the template light-curve model. A this stage a one-parameter least squares fit is performed in order to find the time shift between the two light curves. (vi) Plot an O C diagram, fit an LTE orbit. Finally, the collected O C values can be plotted against time and, if possible, an orbit can be fitted. A template light-curve model We have tested two representations for a template model light curve a polynomial model and a harmonic model. We have decided to use a harmonic model. Such an approach, apart from providing a good model of the input data, enables us to conveniently adjust the initial period of a binary. The model is based on a Fourier series and involves (1) l=1 N has to be chosen so that the resulting light-curve model defined by the coefficients a... an, b1... bn and the period T approximate the raw data as good as possible. Theoretically, the more harmonics are used (big N ), the better the approximation. N has an upper limit though. Obviously, as a least squares algorithm is used to fit f (t), the number of parameters cannot exceed the number of data points used in a fitted, i.e. m = 2N + 2 < n. Moreover, if m is close to n, the fit starts to approximate the data noise. This is not a desired effect In our analysis we used N = 18 which was found to be the best for the types of curves analyzed. An example template model along with the original light-curve are presented in Fig. 2. Though it might seem that a lower value of N would make the model less sensitive to erroneous data points, at the same time the real eclipses (especially deep and short ones) would not have been modeled well enough. Having the primary eclipse well modeled is crucial when searching for time shifts between the model and the local light curves. If one would have high quality photometry it would then be good to optimize the procedure and select N individually for each object. In order to carry out the least-squares fitting we used the Levenberg-Marquardt algorithm and its Minpack1 implementation N 2π 2π a X al cos + lt + bl sin lt. 2 T T Calculating O C Having the light curve model in form of T and a... an, b1... bn coefficients, it can be compared with local light curves in each interval. This is done by fixing the parameters describing the model in Eqn. (1) and slightly modifying the 1

3 Timing of Eclipsing Binary Stars from the ASAS Catalogue 3 formula by introducing a time shift parameter t O: N f O(t) = a 2 + [ ( ) 2π a l cos T l(t to) l=1 ( )] 2π + b l sin T l(t to). Then, using least squares, f O(t) is fitted to local light curves with t O as the only parameter. Effectively, t O is the value of O C at the given point in time. Collecting these for all intervals allows one to obtain a O C diagram. Since Eqn. (2) is fitted using the Levenberg-Marquardt as before, the formal errors of the obtained O C values are derived from the covariance matrix, which, in this case is a one-element matrix due to the fact that the fit has only one parameter. (2) O-C[d] HJD - 24[d] 2.3 Detection Criterion Inspecting visually every single O C diagram is not practical, hence in order to find binaries with significant timing variations we use the following timing activity parameter r Figure 3. ASAS s timing variations. The parabola corresponds to P from Pilecki et al. (27)..6 r = S > rmin, (3) A where S is the standard deviation of the O C values and A is the average error ( (O C)) of these values: A = 1 M M (O C) l. (4) l=1 Objects having r greater than a certain r min are considered interesting. We applied the above criterion on data sets divided into 5-, 6-, 7- and 8-intervals. If an object passes the criterion at least once, it is considered interesting. Such objects are finally inspected visually. O-C[d] HJD - 24[d] 3 TIMING VARIATIONS OF ED AND EC BINARIES FROM THE ASAS CATALOGUE The ASAS Catalogue of Variable Stars (ACVS) is publicly available for download from the ASAS Project homepage 2. We used its version 1.1 in this paper. The catalogue consists of 5124 objects showing variability in brightness. Among them 2761 are uniquely classified as eclipsing contact (EC) binaries and 2271 as eclipsing detached(ed) binaries. Based on the original ACVS index file, we created two subindex files: one containing only EC binaries and a second one containing only ED binaries. This reduced the number of investigated light-curvesto 532. In theanalysis only an A-rated 3 photometry was taken into account. The results for ED are collected in Table 1. It lists 29 binaries. For this table we used r min = 2.. The results for EC are collected in Table 2. It contains 44 objects. In this case we used r min = 2.5. The most frequent variation in O C in both cases has a parabolic shape indicating a linear change in period. The r min threshold has been adjusted individually for each type of curves (ED and EC). In case of contact binaries the eclipses are The ACVS rates the quality of each brightness measurement on a scale form A to D with A being the best and D the worst quality. Figure 4. ASAS s timing variations. The parabola corresponds to P from Pilecki et al. (27). generally shallower and less sharp than in case of detached binaries - this results in different sensitivities of the algorithm in the two cases. Recently Pilecki et al. (27) have found EC binaries with high period change rates (HPCR) in the ASAS data. They have published a list of 31 stars exhibiting large P. 1 out of 44 objects listed in Tab. 2 have been detected by Pilecki et al. (27). Of the remaining 21 HPCR objects, 12 were classified as other than EC hence they were not analyzed by us. Finally, 9 objects did not satisfy our criterion. In order to verify that the results obtained with our proposed algorithm are consistent with Pilecki et al. (27), we compared objects with high P from the HPCR list with our timing measurements. This is demonstrated in Figures 3-4 by plotting the linear P (i.e. parabolic in O C) trend from Pilecki et al. (27) together with our timing measurements. We have also extended the O C diagram from Pilecki et al. (27) for VY Cet, a contact binary system with the well known O C variations. Qian (23) studied this object and found the period of the third body to be 7.3 years with a minimum mass of.62m. Figure 5 shows the original O C diagram with our results overplotted.

4 4 S. K. Koz lowski, M. Konacki and P. Sybilski ASAS ID P[d] (ASAS) P[d] (corrected) r 5 r 6 r 7 r 8 Other ID CD Eri CI Pup SV Pyx DU Vel RV Crt NSV NSV NSV DT Lup SW Oph V467 Sco V496 Sco V795 Sgr V1173 Sgr BT Pav U Gru Table 1. LTE search results: ASAS ED objects. The table lists the object s ID, the period value from the ACVS catalogue, a corrected period obtained during reference model computation. Four r values for 5-, 6-, 7- and 8-interval runs are given. The last column contains other IDs, if any. Figure 5. The timing variations of VY Cet. This O C diagram is from Pilecki et al. (27). Our timing variations are denoted with squares. 4 ECLIPSING BINARIES WITH LTE ORBITS ASAS , ASAS AND ASAS During a visual inspection of the detected timing variations we identified three interesting cases most likely corresponding to an LTE effect due to a third companion. We subsequently analyzed them using a modified version of our approach. It differs from the original one described above in the way the intervals are chosen. Rather than setting the number of intervals, this time, the length of the intervals (in days) and the shift between them (in days) are used as parameters. This enables one to create overlapping intervals. E.g. choosing a length of an interval of 3 days and setting the shift between intervals to 1 days produces 3 intervals assuming that the time span is 32 days. An important difference is that data points are used multiple times for an O C computation. This way it is possible to produce an infinite number of points in the O C diagram having the same input data as in the standard algorithm. The statis-

5 Timing of Eclipsing Binary Stars from the ASAS Catalogue 5 ASAS ID P[d] (ASAS) P[d] (corrected) r 5 r 6 r 7 r 8 Other ID TU Eri NSV ER Ori ASAS BV GZ Pup NSV NSV NSV SV Cen NSV V678 Cen NSV VZ Lib V343 Tel NSV NP Tel Table 2. LTE search results: ASAS EC objects. The table lists the object s ID, the period value from the ACVS catalogue, a corrected period obtained during reference model computation. Four r values for 5-, 6-, 7- and 8-interval runs are given. The last column contains other IDs, if any. tical significance of these points is obviously appropriately lower than in the case of non-overlapping intervals and this must be taken into account when deriving formal errors. The main purpose of such an approach is to investigate how the O C diagram looks between the points calculated using the standard algorithm and also how this influences the bestfitting LTE orbit. ASAS , also known as NSV5654 is an EC binary that has a period of ḍ It has been identified with high values of r reaching 8.81 in a 6-interval run. The O C plot is shown in Fig. 6. A linear trend introduced by an imprecise period was removed by applying a period correction obtained from an orbital fit which included a correction to the period dp as one of the parameters. The O C diagram calculated using the new period shows evidence of an LTE orbit. Three complete cycles seem to be visible. The final orbital parameters are summarized in Tab. 3. Another interesting EC object is ASAS or NSV3836. This binary, having a period of ḍ 39383, reveals periodic variability in the O C diagram. As in the previous case, we used the standard (circles with error bars, solid line) and the overlapping method. Figure 7 shows two sets of O C data points as well as two corresponding orbital solutions. Parameters are shown in Tab. 3. ASAS is the only object among analyzed ED systems that appears to exhibit periodic (O C) variations with a period clearly shorter than the data time span. The eclipsing system has a period of ḍ Figure 8

6 6 S. K. Koz lowski, M. Konacki and P. Sybilski O - C [s] -6-8 O - C [s] HJD [d] HJD [d] Figure 6. ASAS The O C diagram was obtained using the standard (circles with error bars, solid line) and the overlapping methods. Orbital solutions were fitted to both data sets. Figure 7. ASAS The O C diagram was obtained using the standard (circles with error bars, solid line) and the overlapping methods. Orbital solutions were fitted to both data sets. parameter unit standard method overlapping method ASAS P [d] 911 ± a sin I [AU].961 ± e.159 ± ω [deg] 113 ± 8.98 T [HJD] ± f(m 1,2,m 3 ) [M ].143± RMS [s] 67.8 χ k 3 ASAS P [d] 1188 ± a sin I [AU].52 ±.4.57 e.131 ± ω [deg] 84 ± f(m 1,2,m 3 ) [M ].133± T [HJD] ± RMS [s] 4.77 χ k 3 O - C [s] HJD [d] Figure 8. ASAS The O C diagram was obtained using the standard (circles with error bars, solid line) and the overlapping methods. Orbital solutions were fitted to both data sets. ASAS P [d] 1286 ± a sin I [AU].76 ± e..32 ω [deg] T [HJD] ± f(m 1,2,m 3 ) [M ].35± RMS [s] χ k 2 Table 3. LTE orbital parameters: period (P), projection of the semimajor axis a on the line of sight (asini, I is the orbit s inclination), eccentricity (e), argument of periapsis (ω), mass function (f(m 1,2,m 3 )) for ASAS , and computed using data from the standard and the overlapping method. RMS, χ 2 and the degrees of freedom (k) have been calculated for the standard method. shows two sets of the O C points as well as two corresponding orbital solutions. In the case of the orbit fitted to the points obtained with the standard method, a circular orbit was assumed. Table 3 provides the orbital parameters. The above analysis of three interesting cases shows the usefulness of the proposed methods. The standard algorithm is well suited for detection and general O C computation while the second method based on overlapping intervals does the interpolation. It is not an interpolation in a strict mathematical sense though. There is no model (linear, cubic, spline, etc.) O C values are calculated accordingly to the actual trend in the data set. As shown in the examples, such interpolation can influence the shape of the fitted orbit. The eccentricity is the most sensitive parameter. asini and P do not differ much when comparing these two approaches.

7 Timing of Eclipsing Binary Stars from the ASAS Catalogue 7 5 CONCLUSIONS Eclipsing detached and eclipsing contact binaries were investigated. Altogether 532 objects have been analyzed in terms of the LTE. Results from 5-, 6-, 7- and 8-interval runs have undergone a test estimating the likelihood of interesting O C variations. 29 detached and 44 contact binaries have passed the final visual tests. Most of the resulting O C plots have a parabolic shape meaning a linear period increase or decrease. A few objects reveal 3rd degree variations suggesting a long period LTE orbit, of which only a short part is visible in the data set. Finally, three diagrams show evidence of LTE orbits that have periods shorter than the data s span. These objects have been precisely analyzed using a modified version of the LTE-search algorithm. It uses overlapping intervals and generates more points on the O C plot than the standard approach, thus revealing the possible shape of the orbit (without increasing the accuracy). Fitted orbits have semi-major axes smaller than 1AU and 3 year periods. Obtained O C plots were compared with known literature data showing compatibility. The proposed method is well suited for automated data pipelines due to its versatility. It can handle practically any long time-base photometry data and point out most irregularities in eclipse timing. One thing worth mentioning is that, in general, dealing with O C requires very precise periods. In our case, the periods came bundled with photometric data from the ASAS Catalogue. If aliases, like ± 1, 2, 3% of the correct period, shall occur, the algorithm used to detect O C variations will give a false signal. It must be therefore used with caution and understanding of the process. Further simulation work regarding the influence of various light-curve parameters on the detection possibility is in progress. This work is supported by the Foundation for Polish Science through a FOCUS grant and fellowship, by the European Research Council through the Starting Independent Researcher Grant and Polish MNiSW grant no. N N Mukadam A. S., Mullally F., Nather R. E., Winget D. E., 24, ApJ, 67, 982 Mullally F., Winget D. E., De Gennaro S., Jeffery E., Thompson S. E., Chandler D., Kepler S. O., 28, ApJ, 676, 573 Nather R. E., Mukadam A. S., 24, ApJ, 65, 846 Paczynski B., 1997, in Ferlet R., Maillard J.-P., Raban B., eds, Variables Stars and the Astrophysical Returns of the Microlensing Surveys The Future of Massive Variability Searches. p. 357 Paczynski B., Szczygiel D., Pilecki B., Pojmanski G., 26, MNRAS., 368, 1311 Pilecki B., Fabrycky D., Poleski R., 27, MNRAS, 378, 757 Pojmanski G., 22, Acta Astronomica, 52, 397 Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P., 1992, Numerical recipes in C. The art of scientific computing. Cambridge: University Press, 1992, 2nd ed. Qian S., 23, MNRAS, 342, 126 Ribas I., 25, in C. Sterken ed., The Light-Time Effect in Astrophysics: Causes and cures of the O-C diagram Vol. 335 of ASP Conf. Ser., Combining Astrometry and Light-time Effect: Low-mass Companions around Eclipsing Systems. p. 55 Richer H. B., Ibata R., Fahlman G. G., Huber M., 23, ApJ, 597, L45 Sigurdsson S., Richer H. B., Hansen B. M., Stairs I. H., Thorsett S. E., 23, Sci, 31, 193 Sterken C., 25a, in C. Sterken ed., The Light-Time Effect in Astrophysics: Causes and cures of the O-C diagram Vol. 335 of ASP Conf. Ser., Binary Pulsars, General Relativity and Light-Time Effects. p. 215 Sterken C., 25b, in C. Sterken ed., The Light-Time Effect in Astrophysics: Causes and cures of the O-C diagram Vol. 335 of ASP Conf. Ser., Ole Roemer and the Light- Time Effect. p. 181 Wolszczan A., Frail D. A., 1992, Nat, 355, 145 Zasche P., 25, Ap&SS, 296, 127 This paperhasbeentypesetfrom atex/l A TEXfileprepared by the author. REFERENCES Cordes J. M., Kramer M., Lazio T. J. W., Stappers B. W., Backer D. C., Johnston S., 24, New Astronomy Review, 48, 1413 Dolez N., Vauclair G., Chevreton M., 1983, A&A, 121, L23 Friedrich S., Zinnecker H., Correia S., Brandner W., Burleigh M., McCaughrean M., 27, in Napiwotzki R., Burleigh M. R., eds, 15th European Workshop on White Dwarfs Vol. 372 of ASP Conf. Ser., Search for Giant Planets around White Dwarfs with HST, Spitzer, and VLT. p. 343 Hulse R. A., Taylor J. H., 1975, ApJ, 195, L51 Konacki M., Wolszczan A., 23, ApJ, 591, L147 Lee J. W., Kim S., Kim C., Koch R. H., Lee C., Kim H., Park J., 29, AJ, 137, 3181