arxiv:astro-ph/ v1 25 Jun 2002

Transcription

1 QZ Serpentis: A Dwarf Nova with a 2-Hour Orbital Period and an Anomalously Hot, Bright Secondary Star 1 arxiv:astro-ph/ v1 25 Jun 2002 John R. Thorstensen, William H. Fenton Department of Physics and Astronomy 6127 Wilder Laboratory, Dartmouth College Hanover, NH ; john.thorstensen@dartmouth.edu, w.h.fenton@dartmouth.edu Joseph O. Patterson, Jonathan Kemp 2, Jules Halpern Department of Astronomy, Columbia University 538 West 120th Street, New York, NY 10027; jop@astro.columbia.edu, j.kemp@jach.hawaii.edu,jules@astro.columbia.edu Isabelle Baraffe École Normale Supérieure, Lyon Cedex 07, France; ibaraffe@ens-lyon.fr ABSTRACT We present spectroscopy and time-series photometry of the dwarf nova QZ Ser. The spectrum shows a rich absorption line spectrum of type K4±2. K-type secondary stars are generally seen in dwarf novae with orbital periods P orb 6 h, but in QZ Ser the absorption radial velocities show an obvious modulation (semi-amplitude 207(5) km s 1 ) at P orb = (2) min, much shorter than typical for such a relatively warm and prominent secondary spectrum. The Hα emission-line velocity is modulated at the same period and roughly opposite phase. Time-series photometry shows flickering superposed on a modulation with two humps per orbit, consistent with ellipsoidal variation of the secondary s light. QZ Ser is a second example of a relatively short-period dwarf nova with a surprisingly warm secondary. Model calculations suggest that the secondary is strongly enhanced in helium, and had already undergone significant nuclear evolution when mass transfer began. Several sodium absorption features in the secondary spectrum are unusually strong, which may indicate that the present-day surface was the site of CNO-cycle hydrogen burning in the past. Subject headings: stars individual; stars binary; stars variable. 1. Introduction Cataclysmic variable stars (CVs) are close binary systems in which a low-mass secondary trans- 2 Also at Joint Astronomy Center, Hilo, Hawaii. 1 Based in part on observations obtained at the MDM Observatory, operated by Dartmouth College, Columbia University, Ohio State University, and the University of Michigan. fers mass onto a white dwarf; Warner (1995) wrote an excellent monograph on CVs. The Roche geometry tightly constrains the secondary star s mass at a given orbital period P orb. Short-period systems have low-mass secondaries, so if the chemical composition is normal (X 0.7), the secondary is a faint M dwarf or brown dwarf and contributes negligibly to the visible- 1

2 light spectrum (Fig. 4 of Patterson 2001); However, Thorstensen el al. (2002) found a K4 ±2 secondary in the dwarf nova 1RXS J (hereafter RX ), which has P orb = 64 min. They suggested that the secondary was somewhat evolved at the start of mass transfer, with its core substantially enhanced in helium. In this scenario the portion of the secondary remaining today corresponds to the core of the original star, and the enhanced helium greatly affects the masstemperature relation. Here we present new observations of the dwarf nova QZ Ser, which appears to be a close relative of RX QZ Ser was discovered by Katsumi Haseda in 1998 and designated as HadV04 in the discovery notification (vsnet-obs 18349) 3, in which T. Kato suggested an identification with the ROSAT WGA source 1WGAJ (White, Giommi, & Angelini 1994). The outbursting object was discovered on small-scale patrol films, and its quiescent counterpart was initially uncertain. On 2001 June 28 we obtained a spectrum of the star nearest the position, and it proved to be an ordinary late-type star. Returning to the field on 2002 January 20 we obtained a spectrum of a somewhat fainter star 11 arcsec to the northwest. This showed the distinctive broad Balmer emission characteristic of dwarf novae at minimum light, together with absorption features of a late-type star. On 2002 Feb 3.23 UT, P. Schmeer detected an outburst of QZ Ser (vsnet-campaign 1281), and an examination of an outburst image by H. Yamaoka (vsnet-alert 7172) confirmed that the outbursting object was on the position of the emission-line object, cementing the identification. The position of QZ Ser, derived from a fit to 62 USNO A2.0 stars (Monet et al. 1996) on one of our 1.3m images (described below), is α ICRS = 15 h 56 m 54 s.50, δ ICRS = (±0. 3 estimated uncertainty). Its position in the USNO A2.0 is not significantly different, which sets an upper limit on its proper motion µ < 0.03 yr 1. The Living Edition of the Catalog and Atlas of Cataclysmic Variables (Downes et al. 2001) gives an up-to-date finding chart. 3 The vsnet mailing list archives are available at 2. Observations and Analysis Our time-series photometry is from the 1.3 m Mcgraw-Hill telescope at MDM Observatory on Kitt Peak, during 2002 February. After bias subtraction and flat fielding our I-band CCD pictures, we measured instrumental magnitudes using aperture photometry, and differenced the program star against a nearby comparison object. On 2002 January 21 we obtained UBVI exposures in photometric conditions with the 1.3 m telescope. We measured these with DAOPHOT and transformed to standard UBVI magnitudes using observations of Landolt (1992) standard stars. Our BV observations agreed within a few hundredths of a magnitude with photometry listed by Henden 4. QZ Ser had V = 17.69, B V = 0.72, U B = 0.32, and V I KC = 1.22, rather redder than typical short-period dwarf novae at minimum light. Table 1 gives a journal of our spectroscopy, all of which is from the 2.4 m Hiltner telescope at MDM Observatory on Kitt Peak. The procedures, instrumental setup and reduction were essentially as described in Thorstensen et al. (1998). Nearly all the exposures from the extensive 2002 February data set were 480 s, the total exposure in February being s. Because of the season, the February data span only 3.4 h of hour angle. Four 480-s spectra were taken in 2002 June, after this paper was first submitted, mostly to improve the ephemeris. Fig. 1 shows the mean spectrum from February. The mean spectrum from January, before the outburst, appears somewhat similar, but with a flux level 30 per cent lower in the V band, a marginally redder continuum, and much weaker He I emission, and in June the spectrum resembled the one taken in January. Without simultaneous photometry it is difficult to be certain, but in February the star may have been slightly brighter than its minimum, perhaps because of the outburst earlier that month. Table 1 details the emission lines in the February spectrum. The lines appear double-peaked at the top, with separation 600 km s 1. Late-type absorption features are also present. We measured absorption velocities by cross- 4 see ftp://ftp.nofs.navy.mil/pub/outgoing/aah/sequence/ 2

3 correlating our spectra against a rest-frame sum of late-type IAU velocity standards, using xcsao (Kurtz & Mink 1998). The 6000 to 6500 Å region gave the best results, with mean formal uncertainty 10 km s 1. We discarded a few velocities with large formal errors, leaving four velocities from 2002 January and 69 from 2002 February. A least-squares period search (Thorstensen et al. 1996) showed a strong periodicity near two hours. Because the velocity uncertainties are small compared to the amplitude, the period is determined without ambiguity in the daily cycle count despite the relatively small hour-angle span of the data, and the cycle count between the January and February velocities is also secure. Table 3 gives sinusoidal fit parameters for the velocities, and Fig. 2 shows the velocities folded on the period. Hα emission-line velocities measured using a convolution method sensitive to the line wings (Schneider & Young 1980) gave noisier velocities (also in Fig. 2), but yielded a consistent period. A fit to the emission-line velocities with the period fixed at the more accurate absorption-line value gave parameters listed in Table 3. The phase of the emission line velocities lags the absorption lines by 0.58(2) cycles (formal error), whereas a lag of exactly 0.5 would be expected if the emission accurately traced the white-dwarf motion. We conclude that, as often is the case, the emission lines do not trace the white dwarf motion accurately. The equal periods and roughly opposing phases of the emission and absorption do confirm that they arise in the same system, and that the absorption features are from the secondary and not some interloper. Fig. 3 shows a greyscale representation of the February spectra arranged as a function of phase (Taylor et al. 1998). The sharp absorption features and their modulation are obvious. The HeI λ 5876 line shows evidence of an S-wave, i.e., an emission component apparently arising from the hot spot where the accretion stream strikes the disk. A phase-dependent absorption around λ5910, with 20 Å FWHM is also clearly visible. We measured the feature s equivalent width in the original spectra, and while these measures were noisy they showed a smooth and persistent modulation with binary phase (Fig. 2, lower panel). A sine fit gave a mean EW of 3.8 Å, with a half-amplitude of 2.2 Å. The feature reaches maximum strength around the time of maximum radial velocity of the secondary. It seems likely that this feature is due to sodium, but we have been unable to explain it further. Its diffuseness suggests a velocity spread characteristic of the accretion structures rather than the stellar photosphere. Although phase-dependent absorption calls to mind the superficially similar phenomenon in SW Sex stars, the morphology and phasing of this feature in the single-trailed spectrum is quite different from typical SW Sex absorption events (see e.g. Taylor et al. 1998). To estimate the secondary s spectral type, we used spectra (obtained with the same instrumental setup) of stars classified by Keenan & McNeil (1989). Using two temperature-sensitive line ratios appropriate for our spectral resolution and coverage, we found K4 ± 2 subclasses for QZ Ser s secondary. Most of the absorption features appeared similar to the spectral type standards, but the standards do not show the broad λ5910 feature noted above. A pair of absorption lines at λλ5683, 5688 also appears unusually strong, with equivalent widths of 0.9 and 1.2 Å respectively. These are evidently a Na I triplet at λλ , , and The strength of these lines prompted us to search for other sodium features. The NaD lines are unsuitable because of saturation and confusion with λ5876 (they may also be responsible for the unusual λ5910 feature), but a search of the NIST atomic line database 5 yields other features. One, near at λ , appears enhanced in equivalent width by a factor of about three compared to our K star standards; a neighboring sodium line at λ is unfortunately blended with a strong Ca I line at our resolution. In several other features (a blend of λλ and , and one at λλ and ) the enhancement is not enough to be significant, though the features are clearly detected. All the sodium features considered here arise from the 3p level, further corroborating the line identification. These sodium features move with the secondary s spectrum, so we conclude that there is circumstantial evidence for an enhanced sodium abundance in the secondary s photosphere. Because the spectrum does not appear perfectly normal, we cannot put a strong constraint on the 5 asd 3

4 fraction of the light contributed by the secondary, but it is clearly very substantial. From systematically subtracting the library spectra from the average spectrum, we estimate that the secondary contributes 70 ± 20 percent of the light in the Å range. On face value, our V magnitude then implies V = 18.1 ± 0.4 for the secondary alone. Because the secondary s fractional contribution was probably a little higher in January (when the photometry was taken) than in February (when the secondary contribution was measured), our best estimate for the secondary s magnitude is V = 17.9 ± 0.4. The lower panel of Fig. 2 shows averaged, meansubtracted, differential I-band magnitudes as a function of orbital phase. The time averaging suppresses considerable flickering present in the original time series. The light curve shows two humps per orbit, with mean full amplitude 0.076(9) mag. The minima appear to be asymmetric, differing by A = 0.012(7) mag. This is the waveform characteristic of tidally distorted secondaries within a Roche lobe, as calculated for example by Bochkarev et al. (1979) (hereafter BKS). 3. Inferences Distance. The secondary spectrum and period constrain the distance. Assuming the secondary fills its Roche lobe, we have R 2 /R = 0.234f(q)P 2/3 hr (M 2/M ) 1/3 (Beuermann et al. 1998), where f(q) is very close to unity over the range of interest. Evolutionary models suggest M 2 = M (see below), and because this enters weakly in the present calculation we take this as a guide, assuming an uncertainty of ±0.025 M, which yields R 2 = ± R. K- star surface brightnesses inferred from Table 3 of Beuermann et al. (1999) then imply that the secondary has M V = 9.4 ± 0.4, where nearly all the uncertainty arises from the spectral type. Our photometry and estimated secondary-star contribution then gives m M = 8.5 ± 0.6. QZ Ser lies at b = +47, and Schlegel, Finkbeiner, & Davis (1998) estimate E(B V ) = 0.05 in this location from IRAS 100-micron maps, or A V = Accounting for this yields a distance around 460 (+150, 110) pc. At 460 pc, µ < 0.03 yr 1 corresponds to v T < 65 km s 1, which is not unlikely. Orbital Parameters and Masses. After correcting for the estimated blue-star contamination in the I band, the ellipsoidal light curve has A = 0.098(14) mag and A = 0.015(9) mag. To reproduce these numbers, we interpolated and extrapolated from Table I of BKS, adopting plausible values of limb- and gravity-darkening for a cool star (respectively u = 0.8 and β = 0.6, but the results are not sensitive to these choices). With mass ratios in the range , we obtained the measured A for inclination i = 33 ± 4 degrees (at such low inclinations, A depends mainly on i). In principle we could use the measured A to constrain q(i) further, but our measurement is too crude, and cataclysmic variables can easily produce signals at P orb from effects unrelated to ellipsoidal distortion. The secondary s velocity amplitude gives a mass function of 0.075(5) M. This may be distorted by illumination effects, but we see no evidence of variation of the line features with orbital phase. If we assume for a moment that M 2 = M and a broadly typical white dwarf mass M 1 = 0.7 M, the mass function implies i = 32, essentially identical to the inclination derived from the ellipsoidal variations. We emphasize that we have not measured these masses, but we use them only to show that the data are consistent with our scenario without assuming an unusual mass for the white dwarf. As noted earlier, we believe that the apparent emission line radial velocities do not indicate the white dwarf motion with any interesting accuracy. 4. Discussion QZ Ser is only the second short-period dwarf nova to show a K-star secondary. In all other dwarf novae with P orb < 2 hr, save for a handful of helium systems, the secondaries are late M dwarfs which contribute a small fraction of the visible light. How are we to account for this unusual object? As with RX (Thorstensen el al. 2002), we suggest that the secondary evolved significantly on the main sequence prior to mass transfer, enhancing the core with enough helium to greatly affect the mass-t eff relation. The small, hot secondary we see is the remnant of a once more massive secondary. For illustration, we computed some evolutionary models in the framework of the standard 4

5 disrupted magnetic braking scenario (see Baraffe & Kolb 2000 and references therein). At the 2 hr orbital period of QZ Ser, a donor which started mass transfer on the zero-age main sequence (ZAMS) would have spectral type M4- M5 and mass M M. Such a standard sequence, based on the models of Baraffe & Kolb (2000), is displayed in Fig. 4 by the solid curve. Also shown in Fig. 4 are some test calculations in which mass transfer begins near the end of central H burning. The tracks shown represent several choices of initial secondary masses M 2 > 1 M, constant mass loss rates Ṁ M yr 1, and different levels of nuclear evolution X c < 0.1. As Fig. 4 shows, these naturally reproduce the observed properties of QZ Ser and RX The choice of our parameters seems reasonable and we do not require particularly extreme assumptions on our evolutionary scenarios to fit these objects 6. As already mentioned in Thorstensen el al. (2002), the models predict altered surface abundances. Fig. 4 displays the surface enrichment of 14 N processed by the CNO cycle in the deeper layers and mixed up to the surface. The surface helium abundance increases also, reaching mass fraction 0.6 at P orb 2h. In RX the strengths of Hα and HeI λ5876 are in the ratio 3.6:1, whereas this ratio is typically 6 and above in SU UMa stars (Thorstensen el al. 2002); in the 2002 February spectra of QZ Ser, the ratio is 5:1, which may indicate some enhancement of He. As noted earlier, though, He features were nearly absent in the 2002 January pre-outburst spectra, so the He:H line ratio is evidently not a consistent measure of abundance in this object. The sodium enhancement noted earlier (if real) may provide an important window on the secondary s nuclear processing history. The temperatures reached in the deepest layers of stars with masses > 1.2 M near the end of central H burning are high enough to allow the production of 23 Na via the 22 Ne(p, γ) 23 Na reaction. During mass transfer, the convective envelope proceeds inward and may reach the Na-enriched layers. At the 2-hour period of QZ Ser, the bottom of the 6 With our assumptions, the system in principle should have initially tranferred mass on a thermal timescale, until M 2 became small enough for the system to reappear as a standard CV (see Baraffe & Kolb 2000 for a discussion). secondary s convective envelope should reach into layers which have attained temperatures of K during prior evolution. Based on the recent NACRE reaction rates (Angulo et al. 1999), this is hot enough to destroy 22 Ne by proton capture, but it is not hot enough for the Na-Ne cycle to contribute significant Na enhancement (Weiss, Denissenkov, & Charbonnel 2000). A combination of deep mixing and Ne-Na cycling has been considered as an explanation for anomalous Na abundance in globular cluster giants (Denisenko & Denisenkova 1990; Kraft et al. 1997; Weiss, Denissenkov, & Charbonnel 2000). The nuclear reaction networks presently implemented in our models unfortunately do not produce a quantitative estimate the sodium enhancement, but work is in progress to remedy this. At present, we can only point to the strong Na lines as a clue that the surface material was processed at relatively high temperatures, mixed upward, and exposed as mass transfer stripped away the overlying layers. The origin of the Na enhancents in the globular cluster stars is uncertain, and as in that case, it is possible that QZ Ser s sodium was already present in the gas from which the system formed. It would therefore be desirable to detect the primary products of CNO processing He and N enrichment, and C and O depletion. Unfortunately, test model atmospheres with T eff = 4300 K and log g = 5.5, kindly computed by F. Allard (private communication) showed essentially no observable effect when He was enhanced to 50% in mass fraction and the CNO abundances were altered as expected. More noticeable effects should appear in models at cooler temperatures (T eff < 4000 K), where molecules involving C or O form. If our interpretation based on nuclearly evolved donors is correct, we may expect that a nonnegligible number of short period CVs with anomalously hot secondaries will be discovered. As already emphasized by Beuermann et al. (1998) and Baraffe & Kolb (2000), such evolved sequences can indeed explain a substantial fraction of the observed CVs with late spectral types and orbital periods P > 6 h. As shown in Fig. 4, the most evolved sequences which provide an explanation for such systems predict as well the existence of early spectral types at shorter periods. The late spectral type systems with P > 6 h represent a non-negligible fraction of systems in the sample of 5

6 Beuermann et al. (1998). Baraffe & Kolb (2000) were thus concerned by the fact that such evolved sequences may remain active in the 2-3 h period gap, since such extreme secondaries never become fully convective. In order to prevent populating the period gap and predicting too early spectral types at shorter periods, (Baraffe & Kolb 2000, see their 4) suggested an increase of the mean mass transfer rate during the secular evolution of such evolved donors. However, the existence of QZ Ser at P orb = 2.0 h now supports the idea that evolved sequences may remain active in the gap. Our proposed scenario would be supported if QZ Ser and RX were found to have enhanced CNO-process abundances at their surfaces, and a more quantitative study of the abundances of Na, Al, and other light metals may provide this evidence despite the low T eff. In the models we have computed to date, the secondaries in longerperiod systems are significantly less massive than in standard systems. If this tendency proves robust, accurate measures of donor masses in longperiod systems could identify systems destined to evolve into stars like RX and QZ Ser in our scenario. Finally, if our interpretation is correct, one should find similar systems in the period gap. Acknowledgments. We gratefully acknowledge funding by the NSF (AST ), and we thank the MDM staff for their support. We thank especially France Allard for computing the test model atmospheres. 6

7 REFERENCES Angulo, C., et al. 1999, Nucl. Phys. A., 656, 3 Baraffe, I., & Kolb., U. 2000, MNRAS, 318, 354 Beuermann, K., Baraffe, I., Kolb, U., & Weichhold, M. 1998, A&A, 339, 518 Beuermann, K., Baraffe, I., & Hauschildt, P. 1999, A&A, 348, 524 Warner, B. 1995, Cataclysmic Variables (Cambridge University Press) Weiss, A., Denissenkov, P. A., & Charbonnel, C. 2000, A&A, 356, 181 White, N. E., Giommi, P., & Angelini, L. 1994, IAU Circ., No Bochkarev, N. G., Karitskaya, E. A., & Shakura, N. I. 1979, Sov. Astron., 23, 8 (BKS) Denisenko, P. A., & Denisenkova, S. N. 1990, Pis ma v Astronomicheskii Zhurnal, 16, 642 Downes, R. A., Webbink, R. F., Shara, M. M., Ritter, H., Kolb, U., & Duerbeck, H. W. 2001, PASP, 113, 764 Keenan, P. C., & McNeil, R. C. 1989, ApJS, 71, 245 Kraft, R. P., Sneden, C., Smith, G. H., Shetrone, M. D., Langer, G. E., & Pilachowski, C. A. 1997, AJ, 113, 279 Kurtz, M.J., Mink, D.J. 1998, PASP 110, 934 Landolt, A. U. 1992, AJ, 104, 340 Langer, G. E., Hoffman, R., & Sneden, C. 1993, PASP, 105, 301 Monet, D. et al. 1996, USNO-A2.0, (U. S. Naval Observatory, Washington, DC) Patterson, J. 2001, PASP, 113, 736 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 Schneider, D. and Young, P. 1980, ApJ, 238, 946 Taylor, C. J., Thorstensen, J. R., & Patterson, J. O. 1998, PASP, 111, 184 Thorstensen, J. R., Patterson, J., Shambrook, A., and Thomas, G. 1996, PASP 108, 73 Thorstensen, J. R., Taylor, C. J., & Kemp, J. 1998, PASP, 110, 1405 Thorstensen, J. R., Fenton, W. H., Patterson, J. O., Kemp, J., Krajci, T., & Baraffe, I. 2002, ApJL, 567, L49 This 2-column preprint was prepared with the AAS LATEX macros v5.0. 7

8 Table 1 Journal of Spectroscopy UT Date N a HA start HA end 2002 Jan : 32 2 : Jan : 58 2 : Feb : 30 2 : Feb : 30 0 : Feb : 40 1 : Feb : 07 1 : Feb : 34 1 : Feb : 37 0 : Jun : : 18 a Number of spectra. Table 2 Emission Features Feature E.W. a Flux b FWHM c (Å) (10 15 erg cm 2 s 1 ) (Å) Hγ HeI λ Hβ HeI λ HeI λ HeI λ Hα HeI λ HeI λ a Emission equivalent widths are counted as positive. b Absolute line fluxes are uncertain by a factor of about 2, but relative fluxes of strong lines are estimated accurate to 10 per cent. c From Gaussian fits. 8

9 Table 3 Fits to Radial Velocities a Data set b T 0 P K γ N σ (d) (km s 1 ) (km s 1 ) (km s 1 ) Absorption (3) (11) 207(5) 9(4) Hα emission (12) [ ] 69(6) 20(4) a Fits are of the form v(t) = γ + K sin[2π(t T 0 )/P]. The number of points used is N and σ is the standard deviation from the best fit. b Blue-to-red crossing, HJD

10 Fig. 1. Mean flux-calibrated spectrum. In the upper trace, the individual exposures have been shifted into the absorption rest frame and the trace has been shifted upward. Fig. 2. Data folded on the binary period. All points are shown twice for continuity. Upper panel: Radial velocities with best-fitting sinusoids. Emission velocities are shown as round dots; absorption velocities are as follows: 2002 January = open squares, 2002 February = filled squares, and 2002 June = open triangles. Middle panel: Phase-averaged, mean-subtracted, I-band differential magnitudes. Lower Panel: Equivalent widths of the diffuse λ5910 absorption feature. Fig. 3. Spectra near λ5900 rectified and shown as greyscale against phase. All data are shown twice for continuity. The emission feature to the left is HeI λ Note the doppler motion of the sharp absorption features, and the diffuse absorption feature around λ5910. (Included as a separate jpg file in astro-ph preprint. The white horizontal line is an artifact of the ps to jpg conversion.) 10

11 K0 K4 M0 M4 M8 Fig. 4. Spectral type, mass, effective temperature and surface 14 N abundance (normalized to the solar abundance) of the secondary versus orbital period. The solid curve corresponds to a standard sequence with an initially unevolved donor, in the framework of the disrupted magnetic braking model (reproducing the 2-3h period gap). The other curves correspond to evolutionary sequences starting mass transfer near the end of H burning, at a central H mass fraction X c. Long-dashed curve: M 2 = 1.2 M, Ṁ = M yr 1, X c = Dash-doted curve: M 2 = 1.3 M, Ṁ = M yr 1, X c = Short-dashed curve: M 2 = 1.5 M, Ṁ = 10 8 M yr 1, X c = The locations of QZ Ser (present work, filled triangle) and RX 2329 (Thorstensen el al. 2002), filled square) are indicated. The filled circles are observations from Beuermann et al. (1998). 11

12 This figure "figure3.jpg" is available in "jpg" format from: