arxiv: v1 [astro-ph.sr] 27 Feb 2012

Transcription

1 Astronomy & Astrophysics manuscript no c ESO 2012 February 28, 2012 The quest for companions to post-common envelope binaries II. NSVS and HS Beuermann, K. 1, Breitenstein, P. 2, Dẹbski, B. 3, Diese, J. 4, Dubovsky, P. A. 5, Dreizler, S. 1, Hessman, F. V. 1, Hornoch, K. 6, Husser, T.-O. 1, Pojmanski, G. 7, Wolf, M. 8, Woźniak, P. R. 9, Zasche, P. 8, Denk, B. 2, Langer, M. 2, Wagner, C. 2, Wahrenberg, D. 2, Bollmann, T. 4, Habermann, F. N. 4, Haustovich, N. 4, Lauser, M. 4, Liebing, F. 4, and Niederstadt, F. 4 arxiv: v1 [astro-ph.sr] 27 Feb Institut für Astrophysik, Georg-August-Universität, Friedrich-Hund-Platz 1, D Göttingen, Germany 2 Westfalen-Kolleg, Rheinische Straße 67, D Dortmund, Germany 3 Astronomical Observatory, Jagiellonian University, ul. Orla 171, Kraków, Poland 4 Max-Planck-Gymnasium, Theaterplatz 10,D Göttingen, Germany 5 Vihorlat Observatory, Mierova 4, Humenne, Slovakia 6 Academy of Sciences, Fričova 298, CZ Ondřejov, Czech Republic 7 Warsaw University Astronomical Observatory, Al. Ujazdowskie 4, Warszawa, Poland 8 Astronomical Institute, Faculty of Math. and Physics, Charles University, CZ Praha 8, V Holešovičkách 2, Czech Republic 9 Los Alamos National Laboratory, MS D466, Los Alamos, NM 87545, USA Received 16 September 2011; accepted 19 December 2011 ABSTRACT We report new mid-eclipse times of the two close binaries NSVS and HS , harboring an sdb primary and a lowmass main-sequence secondary. Both objects display clear variations in the measured orbital period, which can be explained by the action of a third object orbiting the binary. If this interpretation is correct, the third object in NSVS is a giant planet with a mass of roughly 12 M Jup. For HS , we provide evidence that strengthens the case for the suggested periodic nature of the eclipse time variation and reduces the uncertainties in the parameters of the brown dwarf implied by that model. The derived period is 8.4 yr and the mass is 31 M Jup, if the orbit is coplanar with the binary. This research is part of the PlanetFinders project, an ongoing collaboration between professional astronomers and student groups at high schools. Key words. Stars: binaries: close Stars: binaries: eclipsing Stars: subdwarfs Stars: individual: NSVS , HS Planets and satellites: detection 1. Introduction Some eclipsing close binaries exhibit O C variations in the observed mid-eclipse times relative to those calculated from an underlying linear ephemeris. An early example is the cataclysmic variable UX UMa, which displayed what looked like a periodic variation with an amplitude of 150 s and a period of 29 yr. This modulation was interpreted as being possibly caused by a third body (Nather & Robinson 1974) or by apsidal motion (Africano & Wilson 1976) until Rubenstein et al. (1991) proved its aperiodic nature. Since then an explanation of these variations by processes internal to the binary has been favored, which need not be periodic. One such process (Applegate 1992) involves changes in the internal constitution of the secondary, which give rise to changes in its spin and variations in the orbital period by spin-orbit coupling. With more data becoming available, it was realized, however, that internal processes might be too feeble to explain the observations (e.g. Brinkworth et al. 2006; Chen 2009; Schwarz et al. 2009; Potter et al. 2011), but this notion was opposed by Wittenmyer et al. (2011). It was also realized that the detached post-common envelope binaries (PCEBs) represent more suitable laboratories for studying eclipse time variations than the accreting cataclysmic variables, in particular, the variety to which UX UMa belongs. With this in mind, the interpretation of eclipse time variations in terms of the light-travel time effect produced by a third body was revived in a series of papers by Qian et al. (2009, 2010a,b, and references therein). An intriguing case is the semi-detached cataclysmic variable HU Aqr, which displays multi-periodic O C variations. Qian et al. (2011) attributed these variations to the effects of two planets orbiting the binary, but Horner et al. (2011) and Wittenmyer et al. (2011) demonstrated the secular instability of the proposed planetary system and questioned its existence. While this case remains unresolved, the evidence of systematic and possibly periodic eclipse time variations in detached PCEBs increases. The 14.7 mag star HS (in short HS0705) was discovered by Drechsel et al. (2001) to be an eclipsing close binary with an orbital period of 2.3 h, consisting of an sdb primary and an M-dwarf secondary. Qian et al. (2009, 2010a) summarized all the then available eclipse times and, fitting a sinusoidal to the data, suggested the presence of a brown dwarf orbiting the binary with a period of 7.15 yr. NSVS (in short NSVS1425) was discovered as a 13.2 mag variable star in the Northern Sky Variability Survey (Woźniak et al. 2004) and identified as an eclipsing binary with a 2.6 h orbital period and an sdb primary by Wils et al. (2007). The source seems to be in many respects a twin of HS0705 and the prototype sdb/m-dwarf binary HW Vir (Lee et al. 2009). 2. Observations and data analysis We present here the results of an ongoing search for eclipse time variations in these two sources, which is based on new observa- 1

2 Fig. 1. Eclipse light curves of NSVS Left: MONET/N on 5 May 2011, center: Northern Sky Variability Survey, April August 1999, right: All Sky Automated Survey 3 Table 1. New mid-eclipse times of NSVS and O C ell residuals of the elliptic orbit fit displayed in Fig. 2, right panel. Table 2. New mid-eclipse times of HS and O C ell residuals of the elliptic orbit fit displayed in Fig. 3. Cycle E Ecl BJD(TT) Error O C ell Tel (days) (days) (days) I I I I I I I I I II I I I I II I II I I I I I I I I I I I I I I I (1) NSVS, (2) ASAS3, (3) Kolonica Observatory, (4) San Pedro Matir Observatory, (5) Ondřejov Observatory, (6) MONET/North tions extending throughout 2011 and supplemented by archival data reaching back to Part of the 2011 observations were performed within the project PlanetFinders conducted by highschool teachers and their students in collaboration with profes- Cycle E Ecl BJD(TT) Error O C ell Tel (days) (days) (days) I I I I I I II I II I I I I II I I I I I I I I (1) NSVS, (2) Ondřejov Observatory, (3) MONET/North, (4) Kolonica Observatory, sional astronomers. This project provides high school students with first-hand experience in all aspects of an authentic research project. NSVS1425 and HS0705 were observed with the MONET/North telescope at the University of Texas McDonald Observatory via the MONET browser-based remote-observing interface between April and December The photometric data were taken with an Apogee ALTA E47+ 1k 1k CCD camera in white light, with exposure times of usually 20 s. Using MONET/North, high school students carried out part of the observations from their classroom. Additional observations with integration times between 20 s and 100 s were performed with the 1.0-m Vihorlat National Telescope at Kolonica Saddle, 2

3 Fig. 2. O C diagrams for NSVS with previously published and our new data shown as yellow and green dots, respectively. Left: Residuals O C lin,1 for a linear ephemeris fitted to the pre-2010 data, cycle numbers E < 9000 (dotted line). Right: Residuals O C lin,2 for an elliptic orbit fit relative to the underlying linear ephemeris of the binary (Eq. 4). Weighted mean data points for cyles E 0 are displayed (see text). Slovakia, equipped with an FLI PL1001E CCD in white light; with the 84-cm reflector of the Observatorio Astronomico Nacional en San Pedro Martir, Mexico, in Johnson V; and with the 65-cm reflector at Ondřejov Observatory, Czech Republic, with a G CCD Camera and R c filter. All data were subjected to the usual dark current subtraction and flatfielding. Relative photometry was performed using suitable comparison stars. Fig. 1 (left panel) shows a light curve of NSVS1425 taken with the MONET/N telescope on 5 May 2011 along with the model fit described below. Early coverage of HS0705 and NSVS1425 is available from archival data of the Northern Sky Variability Survey (Woźniak et al. 2004) 1, which provides short glimpses of the sky for an extended period of time in 1999, with exposure times between 20 and 80 s. HS0705 was extensively studied from October 2000 on by Drechsel et al. (2001). The source was subsequently monitored by Qian et al. (2009, 2010a, and references therein). For NSVS1425, on the other hand, there is a dearth of data between 1999 and 2007, when Wils et al. (2007) noted its eclipsing nature. Fortunately, the source was covered by the All Sky Automated Survey 3 performed from Chile (Pojmański 2004) 2 with eclipses contained in data from 2003 to Again, the survey provides only glimpses of the source with its short series of usually 180 s exposures in the V-band. By folding the barycentrically and leap-second corrected data over the known period, mean eclipse light curves can be reconstructed from both the NSVS and the ASAS3 data. These are displayed in the center and right panel of Fig. 1. The model light curves employ parameters fixed at the values determined from fits to the MONET/N data (see next paragraph). The NSVS data that contain eclipses of NSVS1425 cover the period from 19 April to 10 August The ASAS3 data can be broken up into four light curves defining representative mid-eclipse times for dates in 2003, 2005, 2007, and We fitted the eclipse light curves obtained with the MONET/N, Kolonica, Ondřejov, and San Pedro Martir telescopes, using a heuristic model that involves no physics, except requiring symmetry about the mid-eclipse time p 1. In the variable τ = t p 1, with t the time, the observed relative flux F(t) is represented by a modified and truncated inverted Gaussian G(τ) asas3/ multiplied by a polynomial P(τ) that accounts for any real or apparent variation in the out-of-eclipse flux, i.e., F(t) = P(τ) min (1, G(τ)), (1) where [ G(τ) = p 2 p 3 exp 1 ( ) p5 ] τ, (2) 2 p 4 P(τ) = p 6 + p 7 τ + p 8 τ 2, (3) and p i, i = , are the fit parameters. The Gaussian is modified by replacing the square in the exponential with a free parameter p 5, which allows the creation of a more peaked or broader shape of the eclipse light curve. With a zero level p 2 exceeding unity, the truncation provided by the prescription of Eq. 1 reproduces the finite base width of the eclipse. At each step of the iterative fitting procedure, F(t) is furthermore averaged over the finite, not necessarily equidistant exposure intervals. Fig. 1 shows that a close to perfect fit to the observed light curve is achieved. For exposure bins of 20 s, the achievable accuracy of the derived mid-eclipse times is 1 s for the primary eclipses and a few seconds for the secondary eclipses. The formal errors in the mid-eclipse times derived from fits to the NSVS and ASAS3 data amount to about 1 minute. 3. Results 3.1. NSVS Table 1 lists the new mid-eclipse times for NSVS1425, along with the formal errors from the light curve fits. These data complement the only previously published ones from Wils et al. (2007). The first five entries are from the NSVS for 1999 and the ASAS3 for time intervals of a few months in the summers of 2003, 2005, 2007, and These eclipse times are obtained from fits of the known eclipse profile to the barycentrically corrected data folded over the orbital period. The quoted Julian day is that of the lowest data point during the eclipse, but the fraction of the day is derived from the fit. Of the remaining 27 entries, 24 represent primary eclipses (I) and three secondary eclipses (II). All times were corrected to the Solar System barycenter with leap seconds added and quoted as BJD(TT). We derived an 3

4 ephemeris by adding the mid-eclipse times of Wils et al. (2007) converted from HJD to BJD(TT) 3. Using either the data of Wils et al. (2007) alone or all primary mid-eclipse times until September 2009 (cycle number 7322), we obtained the same orbital periods within the errors, P = (88) days and P = (3) days, respectively. There is no evidence of a period change in the pre-2010 data. Fig. 2 (left panel) shows the O C lin,1 values relative to a linear ephemeris with the latter period. The situation changes drastically for the post-2009 data, where a sharp upturn is observed, which corresponds to a rather abrupt period increase by 9 ms. A physical mechanism that can cause such a change in a detached binary is its response to a companion in a highly elliptic orbit hurdling through periastron. Similar variations in QS Vir (Parsons et al. 2010) and DP Leo (Beuermann et al. 2011) have been ascribed to the action of third bodies. While such an interpretation may seem premature for NSVS1425, it fits the available data exceedingly well and provides a definite prediction, which can be tested in the next couple of years. Given the uncertain orbital period of the putative third body, the eccentricity is likewise uncertain, a situation that is similarly met for some comets in the Solar System observed only near perihelion. Fixing the period at P 3 = 20 yr, similar to values found for other PCEBs, we fitted all data with the light-travel time effect produced by a third body in an elliptic orbit. The fit involves seven free parameters, five for the elliptic orbit and two for the underlying linear ephemeris of the binary (Beuermann et al. 2010, 2011). The latter is BJD(TT) = (4) (3) E, (4) where E is the cycle number. The best-fit eccentricity of the 20-year orbit is e 3 = 0.50, but a valley of low χ 2 in the P 3, e 3 plane stretches from 14 yr, 0.30 to beyond 30 yr, While we fitted all individual mid-eclipse times, a clearer graphical presentation is obtained by plotting the weighted mean O C lin,2 values relative to the ephemeris of Eq. 4 for subgroups of data taken close together in time (Fig. 2, right panel). The subgroups combine data taken over intervals of a week to about a month, which is much shorter than the suggested orbital period. The residuals of the individual new timings from the elliptic orbit fit are quoted in Table 1. The amplitude of the O C lin,2 variation, K 3 59 s, suggests a mass of the third body of M 3 12 M Jup, if we adopt a total mass of the binary of 0.62 M, which is similar to those of HS0705 (0.617 M, Drechsel et al. 2001) and HW Vir (0.627 M, Lee et al. 2009). The value of M 3 varies with the uncertain orbital period P 3 as P 2/3 3 and increases with the unknown inclination i 3 of the third body as 1/sin i 3. If the orbit of the companion is roughly co-planar with the binary, it qualifies as a giant planet. For the 20-year orbit, the semi-major axis is 6.3 AU. Our fit yields an argument of periastron for the motion of the center of gravity of the binary measured from the ascending node in the plane of the sky ϖ bin = ϖ Periastron passage occurred near JD (E 5091). The fit suggests that O C lin,2 will reach a maximum of 88 s in 2019 (E 40000), with actual numbers depending on the uncertain value of P HS A large body of mid-eclipse times of HS0705 has been collected by Qian et al. (2009, 2010a) and interpreted in terms of a third body with a mass of about 30 M Jup, orbiting the binary with a 3 Fig. 3. Elliptic orbit fit for HS Shown are the residuals O C lin,2 relative to the underlying linear ephemeris of the binary (Eq. 5). Color coding is as in Fig. 2. period of 7.15 yr. Çamurdan et al. (2012) added a few timings and found a similar period. We obtained 18 more primary and three secondary eclipses of HS0705, significantly extending the coverage. In addition, we recovered an early eclipse timing from the NSVS. The new data are listed in Table 2. We fitted the entire set of 89 primary mid-eclipse times assuming that the variations are caused by a third body in an elliptic orbit. The fit yields χ 2 = 80.6 for 82 degrees of freedom and is shown in Fig. 3. For a clearer presentation, we collected the individual data points again into appropriate subgroups comprising time intervals of up to a week for the more accurate data and one month for the less accurate ones near minimum O C lin,2 (compare Fig. 3 of Qian et al. 2010a, where the individual O C values are displayed). The underlying intrinsic ephemeris of the binary is BJD(TT) = (6) (2) E. (5) The amplitude of the light-time effect produced by the third body is K 3 = 86 ± 1 s, its orbital period is P 3 = 8.41 ± 0.05 yr, and the eccentricity is e 3 = 0.38 ± The argument of periastron is ϖ bin = ϖ =26 ±2 and periastron passage occurred near JD (E 40738) at O C 22 s. For the present model, a maximum of O C 67 s is predicted to occur near E in June With a binary mass of M (Drechsel et al. 2001), the semi-major axis of the third object is 3.52 AU and its mass is M 3 = (31.5 ± 1.0)/sin i 3 M Jup, very similar to the values quoted by Qian et al. (2010a). This mass is significantly lower than the estimate of Çamurdan et al. (2012). The present data provide no evidence of a deviation of the mid-eclipse times from those created by a single object, but the achievable accuracy of better than 1 s opens the prospect to searching for a fourth object, as found, e.g., in the WD/MS binary NN Ser (Beuermann et al. 2010) Secondary eclipses The present observations locate the secondary eclipse in NSVS1425 and HS0705 at orbital phases of φ= ± and φ = ± , respectively. This finding along with the clearly non-sinusoidal form of the eclipse time variations excludes a sizeable contribution of apsidal motion (Todoran 1972) to the observed effect. 4

5 4. Discussion We have presented evidence of eclipse time variations in NSVS1425 and HS0705, which are closely in line with variations expected from the light-travel time effect of a third body. The existence of companions orbiting PCEBs can no longer be easily refuted, although the cases of UX UMa and perhaps HU Aqr illustrate that there should be no indiscriminate application of the planetary model to all PCEBs showing eclipse time variations. In HS0705, the relatively short period of the third body of 8.4 yr and the large amplitude of 84 s will allow us to test the periodicity of the signal within the next decade, and confirmation will provide support for the third-body model. The often considered alternative of spin orbit coupling, where star cycles change the spin of the secondary star and thereby the angular momentum of the binary orbit (Applegate 1992), remains a contender, but may be too feeble to account for the observed effects (e.g. Brinkworth et al. 2006; Chen 2009; Schwarz et al. 2009; Potter et al. 2011). Apsidal motion, can be safely excluded, at least for the objects studied here. Given the large number of planets detected around single stars, it is unlikely that substellar objects do not exist in orbits around close binaries. The observational problem is to differentiate between the O C variations caused by these companions and any effects intrinsic to the binary, which would require data covering many years. The fundamental questions then relate to the origin of such companions, which need not be the same for planet-sized and brown-dwarf objects, and to the incidence in different subtypes of PCEBs. With HW Vir (Lee et al. 2009), HS0705, and probably NSVS1425, at least three PCEBs with sdb primaries possess companions. For systems with a WD primary, the companions to NN Ser (Beuermann et al. 2010) and DP Leo (Qian et al. 2010b; Beuermann et al. 2011) appear reasonably well-established. Depending on the type of primary, these systems have different evolutionary histories (Zorotovic et al. 2011). For the former, the mass of the sdb star indicates that they are in the He-burning stage, having formed on the first giant branch. On the other hand, NN Ser and DP Leo contain CO-WDs and have formed on the asymptotic giant branch. The mass distribution of white dwarfs from the SDSS (Rebassa-Mansergas et al. 2012) 4 suggests that some PCEBs contain He-WDs and it will be interesting to find out whether these systems, too, contain companions, or more generally planetary systems, a finding that may shed some light on the origin of companions to PCEBs. Çamurdan, C. M., Zengin Çamurdan, D., & Ibanoǧlu, C. 2012, New Astr., 17, 325 Chen, W.-C. 2009, A&A, 499, L1 Drechsel, H., Heber, U., Napiwotzki, R., et al. 2001, A&A, 379, 893 Horner, J., Marshall, J. P., Wittenmyer, R. A., & Tinney, C. G. 2011, MNRAS, 416, L11 Lee, J. W., Kim, S.-L., Kim, C.-H. et al. 2009, AJ 137, 3181 Nather, R. E., & Robinson, E. L. 1974, ApJ, 190, 637 Parsons, S. G., Marsh, T. R., Copperwheat, C. M. et al. 2010b, MNRAS, 407, 2362 Pojmański, G. 2004, Astronomische Nachrichten, 325, 553 Potter, S. B., Romero-Colmenero, E., Ramsay, G., et al. 2011, MNRAS, 1263 Qian, S- B., Zhu, L. Y., Zola, S. et al., 2009, ApJ 695, L1631 Qian, S.-B., Zhu, L.-Y., Liu, L., et al. 2010b, Ap&SS, 329, 113 Qian, S.-B., Dai, Z.-B.,Liao, W.-P. et al. 2010a, ApJ, 708, L66 Qian, S.-B., Liu, L., Liao, W.-P., et al. 2011, MNRAS, 414, L16 Rebassa-Mansergas, A., Nebot Gómez-Morán, A., Schreiber, M. R., et al. 2012, MNRAS, 419, 806 Rubenstein, E. P., Patterson, J., & Africano, J. L. 1991, PASP, 103, 1258 Schwarz, R., Schwope, A. D., Vogel, J., et al. 2009, A&A, 496, 833 Todoran, I. 1972, Ap&SS, 15, 229 Wils, P., di Scala, G., & Otero, S. A. 2007, IBVS, 5800, 1 Wittenmyer, R. A., Horner, J. A., Marshall, J. P., Butters, O. W., & Tinney, C. G. 2011, MNRAS, 1920 Woźniak, P. R., Vestrand, W. T., Akerlof, C. W., et al. 2004, AJ, 127, 2436 Zorotovic, M., Schreiber, M. R., Gänsicke, B. T. 2011, A&A, 536, A42 Acknowledgements. This work is based in part on data obtained with the MOnitoring NEtwork of Telescopes (MONET), funded by the Alfried Krupp von Bohlen und Halbach Foundation, Essen, and operated by the Georg-August- Universität Göttingen, the McDonald Observatory of the University of Texas at Austin, and the South African Astronomical Observatory. The Astronomie & Internet program of the Foundation and the MONET consortium provides a major fraction of the observation time to astronomical projects in high schools worldwide. The work of JD and his students at the Max-Planck-Gymnasium was supported by the Robert Bosch Foundation. The research of MW and PZ was supported by the Research Program MSM of the Ministry of Education of the Czech Republic. References Africano, J., & Wilson, J. 1976, PASP, 88, 8 Applegate, J. H. 1992, ApJ, 385, 621 Beuermann, K., Hessman, F. V., Dreizler, S. et al., 2010, A&A, 521, L60 Beuermann, K., Buhlmann, J., Diese, J., et al. 2011, A&A, 526, A53 Brinkworth, C. S., Marsh, T. R., Dhillon, V. S., & Knigge, C. 2006, MNRAS, 365,