The Orbital Period and Variability of the Dwarf Nova ES Draconis
TThe Orbital Period and Variability of the Dwarf Nova ES Draconis
F. A. Ringwald and Kenia Velasco Department of Physics California State University, Fresno 2345 E. San Ramon Ave., M/S MH37 Fresno, CA 93740-8031, U. S. A. E-mail: [email protected]
ABSTRACT A radial velocity study of the cataclysmic variable ES Dra (PG 1524+622) is presented. ES Dra is found to have an orbital period of 0.17660 ± 0.00006 d (4.2384 ± 0.0014 h). The mass-losing secondary star of ES Dra is detectable in the spectrum, and it has a spectral type of M2 ± 1. From this, we estimate the absolute magnitude of ES Dra during our spectroscopic observations to have been M R = 6.5 ± 0.5, and its distance to be 720 ± 150 pc. The long-term light curve of ES Dra compiled by the American Association of Variable Star Observers (AAVSO) shows that ES Dra is a Z Cam star, which between 1995 and 2009 spent most of its time in standstill. KEYWORDS stars: individual: ES Dra (PG 1524+622) – novae: cataclysmic variables – binaries: spectroscopic – AAVSO 1. INTRODUCTION Cataclysmic variable binary stars (CVs) are close binary star systems, in which typically a late-type, approximately main-sequence star, called the secondary star, orbits a white dwarf. In all CVs, the secondary star fills its Roche lobe, or gravitational equipotential, and spills gas onto the white dwarf. Because the two stars are orbiting each other, the gas stream does not strike the white dwarf directly: it forms a ring around the white dwarf, which viscosity in the gas broadens into an accretion disk, through which the gas accretes onto the white dwarf. Reviews on CVs include those by Robinson (1976), Warner (1995), and Hellier (2001). CV orbital periods ( P orb ) range from about 78 minutes to 12 hours, or longer. Many properties of CVs depend on their orbital periods, including their secular evolution, uminosity, and outbursts (Shafter, Wheeler, & Cannizzo 1986). The catalogs of Patterson (1984) and Ritter & Kolb (2003) include only CVs of known or suspected orbital period, emphasizing this parameter's importance. Radial velocity studies reveal orbital periods in CVs more reliably than do photometric modulations, save for eclipses, but are not without problems (Araujo-Betancor et al. 2005). The strongest features in CV spectra are their emission lines, which originate in the CVs' accretion disks. Orbital amplitudes, or K –velocities, when measured from the emission lines, do not reliably trace the motion of the white dwarf, since they originate in the disk. Neither do γ em , the emission line mean velocity, nor T , the epoch of spectroscopic phase zero. Eclipsing systems sometimes show T lagging the eclipse by over 70º in phase (e.g., Thorstensen et al. 1991). This paper reports a spectroscopic period measurement of the cataclysmic variable ES Draconis. ES Dra was discovered by the Palomar-Green survey (Green, Schmidt, & Liebert 1986) and listed as PG 1524+622. In Section 2, observational procedure, data reduction, and analysis are described. The average optical spectrum of ES Dra is discussed in Section 3. We estimate the absolute magnitude and distance of ES Dra in Section 4. The long-term light curve and outburst behavior of ES Dra are discussed in Section 5, with conclusions in Section 6. 2. OBSERVATIONS Spectra of ES Dra were taken with the Mark IIIa spectrograph on the 1.3-m McGraw-Hill telescope at MDM Observatory. The long-slit spectra permitted accurate sky subtraction, although absolute spectrophotometry may have been accurate only to within 30% because of losses through the 2.2-arcsecond slit that was used for all exposures. Instrument rotation was deemed unnecessary for such a red setup. The detector was a TI-4849 CCD in the BRICC camera (Luppino 1989). A 300 line/mm grism blazed at λ λλ onvolving the H α emission line with two Gaussians that were 2.8 channels (640 km s -1 ) wide, separated by 1300 km s -1 (see Shafter 1983). Table 1 lists the measured velocities. A Lomb-Scargle periodogram (Lomb 1976; Scargle 1982; Press et al. 1992) of the radial velocities is shown in Figure 1. The correct alias is obvious, at a period of 0.17660 ± 0.00006 days. Both statistics from the Monte Carlo analysis of Thorstensen & Freed (1985), the discriminatory power and the correctness likelihood, are over 99.9%. Figure 2 shows a least-squares fit to a sinusoid with the most likely period. Spectrograph flexure caused variations of ± 10 km s -1 , shown by velocities measured from the night-sky line at λ λλ II λ α is in strong emission, with an equivalent width of 40 ± 5 Å and a full-width at zero intensity of 5300 ± 200 km s -1 . There is no sign of the absorption wings that dwarf novae in outburst maximum or nova-like variables often show (see Williams 1983; see also Figure 3.6 of Hellier 2001). Some of the Paschen lines and He I λ I λ
4. ABSOLUTE MAGNITUDE AND DISTANCE From the average spectrum (Figure 3), we estimate a through-the-slit magnitude for ES Dra at the time the spectra were taken to have been R = 15.8 ± 0.2. We also estimate the fraction of the light of the mass-losing secondary star of ES Dra at λ M V = +9.9 and V − R = 1.50 for an M2 dwarf (p. 388 of Cox 2000), and also assuming that E(B − V) = 0 and V − R ≈ M R = +6.5 ± 0.5 at the time of our spectroscopic observations. It also implies that ES Dra has a distance d = 720 ± 150 pc. This does not agree with the estimate of d = 1358 pc of Ak et al. (2008), which was based on their period-luminosity-color relation for infrared magnitudes. It was also based on their estimate of E(B − V) = 0.018 for ES Dra, which was calculated from the maps of interstellar dust of Schlegel, Finkbeiner, & Davis (1998) and the model of interstellar absorption of Bahcall & Soneira (1980) along the line of sight to ES Dra, which is at b = 47º. We calculate that the recalibration of Ramseyer (1994) of the method of Bailey (1981) gives a lower limit for the distance to ES Dra of 530 ± 70 pc, assuming our orbital period and the 2MASS magnitudes for ES Dra of J = 15.458 ± 0.064, H = 14.880 ± 0.077, and K s = 14.889 ± 0.127 (Cutri et al. 2003), with the transformation of Carpenter (2001) to the C. I. T. system (Elias et al. 1982). Our distance estimate is consistent with the best estimate of Godon et al. (2009) of d = 770 pc, which is based on the estimate of 702 pc from the outburst maximum versus orbital period relation of Warner (1987), the estimate of 840 pc from the similar but recalibrated relation of Harrison et al. (2004), and the estimate of 540-1016 pc from infrared magnitudes and models of secondary stars of Knigge (2006, 2007). Godon et al. (2009) also measured E(B − V) = 0.02, from continuum fits to their ultraviolet spectrum of ES Dra. Misselt & Shafter (1995) showed that ES Dra is not eclipsing, but this only constrains the orbital inclination to i ≤ m pg = 13.9 and 16.3, with one variation of 1.3 magnitudes in one day. A long-term light curve of observations by the American Association of Variable Star Observers (AAVSO) (Henden 2009), from 1995 October to 2009 December, is shown in Figure 4. This light curve includes both V magnitudes measured by CCD imaging and visual magnitudes ( m vis ) by eyepiece observers: we expect the difference between V and m vis to be small, relative to the ± 0.1-magnitude errors typical for eyepiece visual magnitude estimates. Figure 4 does ot include apparent magnitudes measured through other filters, nor does it include magnitude limits of non-detections. Dwarf nova outbursts are obvious in Figure 4, but they have relatively small amplitudes of 1-2 magnitudes brighter than the apparent magnitude at which ES Dra spends most of its time, V ≈ V = 17 or fainter, which go considerably fainter than the archival plates examined by Andronov (1991). One interpretation of this long-term light curve may be that ES Dra is a dwarf nova of the U Gem type, which is a dwarf nova that has only normal outbursts. The amplitude of the outbursts is small for a dwarf nova, though, averaging between 1 and 2 magnitudes: more typical is 2-5 magnitudes (Robinson 1976). Honeycutt, Robertson, & Turner (1998) and Honeycutt (2001) pointed out what they called “stunted” outbursts, which are sometimes present in luminous nova-like variables. The outbursts of ES Dra appear to be similar in amplitude. They do not seem to be always present. During the most active period, between JD 2,452,000 to 2,452,600 (2001 March to 2002 November, respectively), the outbursts recurred at a rate of about 7-8 outbursts per year, and lasted between 5 and 20 days, although these statistics are approximate, since coverage was sparse. However, our spectrum of ES Dra (Figure 3) more closely resembles the relatively optically thin spectrum of a dwarf nova between outbursts. It does not show the weak emission lines and line wings in absorption that dwarf novae in outburst and nova-like variables usually show. Also, the decline in the amplitudes of the outbursts and progressive increase of faintest V after 2001 March is similar to that shown by Z Cam stars going into one of their characteristic standstills (see Figure 5; see also Figure 5 of Szkody & Mattei 1984). We also note an outburst, with an amplitude of 3 magnitudes, occurring during a low state, near JD 2,453,800 (see Figure 6). We therefore think that ES Dra is a Z Cam star, and that it spent most of the time between 1995 and 2009 in standstill, not quiescence. If ES Dra is a Z Cam star, and if its distance is 720 ± 150 pc, as shown above, then at outburst maximum, at V ≈ M V (max) = 4.5 and 4.3, respectively. If ES Dra is a Z Cam star, and if its distance is 720 ± 150 pc, then when V ≈ M V = 5.8-6.2. When V ≈
17 in quiescence, M V = 7.7. All these values are consistent with absolute magnitudes of Z Cam stars, in their various outburst states (see Figure 11 of Warner 1987). 6. CONCLUSIONS This paper presents a definitive orbital period of 4.2384 ± 0.0014 hours for ES Dra, measured by a radial velocity study. We find a spectral type for the mass-losing econdary star to be M2 ± 1, and a distance of 720 ± 150 pc. We also argue that ES Dra is a Z Cam star. ACKNOWLEDGMENTS The radial velocity observations were made at Michigan-Dartmouth-MIT Observatory, Kitt Peak, Arizona, which when these observations were made was owned and operated by a consortium of the University of Michigan, Dartmouth College and the Massachusetts Institute of Technology. This publication makes use of data products from the Two Micron All Sky Survey, which was a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. Special thanks to Professor John Thorstensen for the use of his bubble machine, to do the Monte Carlo analysis. We acknowledge with thanks the variable star observations from the AAVSO International Database contributed by observers worldwide and used in this research. REFERENCES Ak T., Bilir S., Ak S., Eker Z., 2008, New Astronomy, 13, 133 Andronov I. L., 1991, I. B. V. S. 3645 Araujo-Betancor S., Gänsicke B. T., Hagen H.-J., Marsh T. R., Harlaftis E. T., Thorstensen J. R., Fried R. E., Schmeer P., Engels D., 2005, A&A, 430, 629 Bahcall J. N., Soneira R. M. 1980, ApJS, 44, 73 Bailey J. 1981, MNRAS, 197, 31 Carpenter J. M., 2001, AJ, 121, 2851 Cox A. N. (ed.), 2000, Allen's Astrophysical Quantities (4th ed.) (New York: Springer) Cutri R. M. et al. 2003, The IRSA 2MASS All-Sky Point Source Catalog, NASA/IPAC Infrared Science Archive. http://irsa.ipac.caltech.edu/applications/Gator/ Downes R., Hoard D. W., Szkody P., Wachter S., 1995, AJ, 110, 1824 Elias J. H., Frogel J. A., Matthews K., Neugebauer G. 1982, AJ, 87, 1029 riend M. T., Martin J. S., Smith R. C., Jones D. H. P., 1988, MNRAS, 233, 451 Godon P., Sion E. M., Barrett P. E., Szkody P., 2009, ApJ, 701, 1091 Green R. F., Schmidt M., Liebert J., 1986, ApJS, 61, 305 Harrison T. E., Johnson J. J., McArthur B. E., Benedict G. F., Szkody P., Howell S. B., Gelino D. M. 2004, AJ, 127, 460 Hellier C., 2001, Cataclysmic Variable Stars: How and Why They Vary (Chichester: Springer-Praxis) Henden A. A., 2009, Observations from the AAVSO International Database, private communication Honeycutt R. K., 2001, PASP, 113, 473 Honeycutt R. K., Robertson J. W., Turner G. W., 1998, AJ, 115, 2527 Knigge C., 2006, MNRAS, 373, 484 Knigge C., 2007, MNRAS, 382, 1982 Lomb N. R., 1976, ApSS, 39, 447 Luppino G. A., 1989, PASP, 101, 931 Misselt K. A., Shafter A. W., 1995, AJ, 109, 1757 Oke J. B., 1974, ApJS, 27, 21 Patterson J., 1984, ApJS, 54, 443 Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B., 1992, Numerical Recipes in C, Second Edition, Cambridge Univ. Press, Cambridge, p. 575 Ramseyer T. F., 1994, ApJ, 425, 243 Ringwald F. A., Thorstensen J. R., Hamwey R. M., 1994, MNRAS, 271, 323 Ritter H., Kolb U., 2003, A&A, 404, 301 Robinson E. L., 1976, Annual Review of Astronomy and Astrophysics, 14, 119 Scargle J. D., 1982, ApJ, 263, 835 chlegel D. J., Finkbeiner D. P., & Davis M., 1998, ApJ, 500, 525 Shafter A. E., 1983, ApJ, 267, 222 Shafter A. W., Wheeler J. C., Cannizzo J. K., 1986, ApJ, 305, 261 Szkody P. & Mattei J. A. 1984, PASP, 96, 988 Thorstensen J. R., Freed, I. W. 1985, AJ, 90, 2082 Thorstensen J. R., Ringwald F. A., Wade R. A., Schmidt G. D., Norsworthy J. E., 1991, AJ, 102, 272 Warner B., 1987, MNRAS, 227, 23 Warner B., 1995, Cataclysmic Variable Stars (Cambridge Univ. Press, Cambridge) Williams G., 1983, ApJS, 53, 523 Table 1. H α Emission Radial Velocities for ES Dra, 1991 May 24-26 a HJD b V HJD b V HJD b V (km s − ) (km s − ) (km s − ) 8401.839 −
59 8401.963 −
173 8402.806 −
30 8401.850 −
31 8402.693 −
205 8402.817 −
129 8401.862 −
76 8402.704 −
193 8402.829 −
191 8401.873 8 8402.716 −
50 8402.840 −
232 8401.884 55 8402.727 −
62 8402.853 −
219 8401.895 −
34 8402.738 23 8402.866 −
176 8401.907 −
46 8402.749 6 8402.878 −
121 8401.918 −
70 8402.761 17 8402.889 −
106 8401.929 −
117 8402.772 −
2 8402.901 −
60 8401.940 −
122 8402.784 43 8402.912 43 8401.952 −
173 8402.795 −
22 8402.924 35 a Double Gaussian algorithm, separation 1300 km s − . b Heliocentric Julian Date of mid-integration, minus 2,440,000. Table 2. Derived Orbital Parameters for ES Dra, from H α velocities a P orb K em γ em T σ (days) (km s − ) (km s − ) (HJD – 2,440,000) (km s − ) 0.17660 ± 0.00006 134 ± 11 −
75 ± 8 8411.024 ± 0.002 48 a Velocities fitted to