Rapid Oscillations in Cataclysmic Variables. XVI. DW Cancri

Joseph Patterson,1 John R. Thorstensen,2 Tonny Vanmunster,3 Robert E. Fried,4 Brian Martin,5 Tut Campbell,6 Jeff Robertson,6 Jonathan Kemp,7,1 David Messier,8 and Eve Armstrong1  

ABSTRACT

We report photometry and spectroscopy of the nova‐like variable DW Cancri. The spectra show the usual broad H and He emission lines, with an excitation and continuum slope characteristic of a moderately high accretion rate. A radial‐velocity search yields strong detections at two periods, 86.1015(3) minutes and 38.58377(6) minutes. We interpret these as respectively the orbital period of the binary, and the spin period of a magnetic white dwarf. The light curve also shows the spin period, plus an additional strong signal at 69.9133(10) minutes, which coincides with the difference frequency . These periods are stable over the 1 year baseline of measurement. This triply‐periodic structure mimics the behavior of several well‐credentialed members of the “DQ Herculis” (intermediate polar) class of cataclysmic variables. DQ Her membership is also suggested by the mysteriously strong sideband signal (at ), attesting to a strong pulsed flux at X‐ray/EUV/UV wavelengths. DW Cnc is a new member of this class, and would be an excellent target for extended observation at these wavelengths.

Received 2004 February 11; accepted 2004 March 7; published 2004 May 17

1. INTRODUCTION

 

DW Cancri is a little‐studied cataclysmic variable, discovered in the Byurakan blue‐galaxy survey and identified as a CV from its Balmer emission lines (Stepanian 1982; Kopylov et al. 1988). Uemura et al. (2002) reported kilosecond quasi‐periodic oscillations, centered on 37 and 73 minutes. Rodriguez‐Gil et al. (2004, hereafter RGAC) reported a radial‐velocity period in the range 77–86 minutes, and suggested photometric periods near 38 and 77 minutes, although their data were severely plagued by aliases.

The curious light curves in the Uemura et al. study launched DW Cnc as a high‐priority target for our global network of telescopes, the Center for Backyard Astrophysics (CBA; Skillman & Patterson 1993). We invested many nights in 2002–2004, studying the star with CBA photometry, along with contemporaneous spectroscopy. In this paper we report the results, which certify DW Cnc as a new member of the DQ Herculis (intermediate polar) class of cataclysmic variables.

2. PHOTOMETRY

 

2.1. Campaign and Light Curves

We obtained time‐series photometry of DW Cnc for a total of 454 hr over 90 nights in 2002–2004. About half the data were acquired through filters (U, V, I) on the MDM 1.3 m telescope, and about half at several CBA observing stations (with the lion’s share from CBA–Belgium, CBA–Arkansas, CBA–Calgary, and Braeside Observatory). Some nights of CBA observation used a V filter, but most were unfiltered (to maximize signal‐to‐noise on this faint star). The primary comparison stars used were 1 9 east ( , , ) and 2 5 west northwest ( , , ; magnitudes from A. A. Henden 2004, private communication) from DW Cnc. We found that aside from small and measurable zero‐point shifts, the V and unfiltered time series were substantially the same, and could be spliced together for our main purpose, the study of periodic signals. The star also cooperated by remaining within ∼0.1 mag of throughout our observations. This appears to be the normal “bright state” of DW Cnc. The photometric observing log is summarized in Table 1.

TABLE 1
TABLE 1 2002–2004 Photometry

Open New Window

A typical night’s light curve is shown in the upper frame of Figure 1. The light curve seems unremarkable for CVs, except that the apparent flickering is of unusually large amplitude on both long (∼1 hr) and short (seconds) timescales. Power spectra of individual nights show no stable features at high frequency (out to 0.2 Hz), but consistent peaks near 37 and 21 c/d.9 In fact, most of the variance at low frequency comes from the presence of these signals, which turn out to be very stable. The mean nightly power spectrum—the incoherent sum over the 10 longest and best nights—is shown in the lower frame of Figure 1. The continuum power in this figure follows , a rather flat distribution for a CV (exponents of −1.5 to −2.2 are the rule).

Fig. 1.— Top: Sample light curve of DW Cnc. Bottom: Nightly average power spectrum, the incoherent sum over the 10 longest nights.

Open New Window

2.2. Periodic Signals

Since the campaign is long and the data span a large range in terrestrial longitude, we are able to detect periodic signals with great sensitivity and (substantial) freedom from aliases. So we spliced the data in various convenient clusters, and analyzed them for periodic content by calculating the power spectra of these clusters. The interval of densest coverage was in 2003 February (JD 671–693), and the power spectrum is shown in the upper frame of Figure 2. Very powerful signals are seen at 20.595 and 37.321(6) c/d, and another significant signal at 12.993(7) c/d. Each signal is essentially sinusoidal to within limits of measurement, since no harmonics are detected; indeed, there are no other periodic signals at all (to a semiamplitude upper limit ∼0.01 mag) out to 300 c/d. After “pre‐whitening” the data by subtracting these three best‐fitting sinusoids, we calculated the power spectrum of the residual time series. The result, seen in the lower frame of Figure 2, shows only a weak and broad bump near 16 c/d.

Fig. 2.— Top: Power spectrum of a 23 night light curve in 2003, showing three signals with their frequencies (±0.006 c/d) flagged by arrows. Bottom: Power spectrum of the residuals, after these signals are removed. Only a possible very weak signal near 16 c/d remains; the arrow points to the “orbital” frequency identified by spectroscopy.

Open New Window

We repeated this analysis for five other clusters, with basically similar results. An example is the 2002 December cluster (JD 613–622), seen in Figure 3. After subtraction of the two obvious signals, we obtained the power spectrum seen in the lower frame. In the latter case, however, we “cleaned” the power spectrum by removing the alias pattern of the two candidate signals, at 12.99 and 16.80(2) c/d. These signals are too weak to be considered secure, but the one at the lower frequency acquires some extra credibility by its agreement with the secure detection in Figure 2.

Fig. 3.— Top: Power spectrum of a 10 night light curve in 2002 December, showing strong signals consistent with those of Fig. 2. Bottom: Power spectrum of the residuals, after removing the strong signals. Probable detections at 12.99 and 16.80 (±0.02) c/d.

Open New Window

We calculated times of minimum light in each detected periodicity for each cluster of data. These timings are given in Table 2, and reduced to diagrams in Figure 4 (where each point represents an average over 3–10 nights). The strong periodic signals track the following ephemerides: In Figure 5 we present mean waveforms for each of the detected signals during the 2003 February segment. They appear roughly sinusoidal, as suggested by the power spectra. The semiamplitudes are 0.078(4), 0.084(4), 0.034(7), and 0.023(9) mag for the four signals (respectively 70 minutes, 38 minutes, 110 minutes, and orbital—as we shall refer to these signals at 20.60, 37.32, 12.99, and 16.72 c/d). Other segments gave similar results for the two strong signals. Thus, it appears that these signals maintain essentially constant amplitude, although we know this only for long timescales (weeks). Amplitude change on shorter timescales would probably elude detection in this study.

TABLE 2
TABLE 2 Photometric Minima (HJD 2,452,000+)

Open New Window

Fig. 4.— diagrams of the timings of photometric minima, relative to eqs. (1) and (2). The signals maintain constant phase over the 400 day baseline.

Open New Window

Fig. 5.— Mean waveforms in V light of four photometric signals during 2003 February. HJD times of maximum light and semiamplitudes were as follows: 38 minutes, 671.6444, 0.084 mag; 70 minutes, 671.6750, 0.078 mag; 110 minutes, 671.6624, 0.034 mag; orbital, 671.6537, 0.023 mag.

Open New Window

Since the spectroscopy reported below appears to specify the true orbital period with high precision, we studied the data carefully for its possible presence. A marginal detection in the best cluster is possible (shown in Fig. 5); but this detection weakened considerably when the entire year’s data were considered. So it would be better to say that no stable orbital modulation is present, to a semiamplitude upper limit of 0.02 mag.

2.3. Colors

Like nearly all CVs, DW Cnc is a very blue star, with , . We obtained 10 nights of multicolor photometry in order to study the pulsed fraction in several passbands. The main signals, at 20.61 and 37.31 c/d, were found to have similar amplitudes in V and I, but the 20.61 c/d signal was somewhat weaker in U (0.05 compared to 0.08 mag semiamplitude). The data were not strictly simultaneous, and there are night‐to‐night variations of up to 30% in pulse amplitude, so this constraint is not very strong. But to a first approximation, both pulse colors are fairly blue (similar to the star’s mean color), with the 20 c/d signal somewhat less blue.

3. SPECTROSCOPY

 

Our spectroscopy is from four observing runs in 2002–2003 (Table 3). We used the Hiltner 2.4 m telescope at MDM, the modular spectrograph, and a thinned 20482 Tektronix CCD yielding 2 Å pixel−1 from 4300 to 7500 Å. The signal‐to‐noise decreased at each end because of vignetting in the camera. The 1 slit gave a spectral resolution of 3.5 Å, and frequent exposure of arc lamps served to maintain an accurate wavelength scale. We observed flux standards and hot stars to allow conversion to absolute flux, and reduced to flux versus wavelength using standard procedures in IRAF. Individual exposures of DW Cnc were 4 to 6 minutes.

TABLE 3
TABLE 3 Journal of Spectroscopy

Open New Window

Figure 6 shows the mean spectrum from 2002 February. The continuum is quite blue, with and emission lines of moderately high excitation. He ii emission at λ4686 and λ5411 is present. Fluxes and equivalent widths (EWs) are contained in Table 4. Double‐peaked emission is evident in the weaker He i lines and in Fe ii λ5169, with separations near ±460 km s−1.

Fig. 6.— Mean spectrum of DW Cnc from 2002 February. Line identifications and strengths are given in Table 4.

Open New Window

TABLE 4
TABLE 4 Emission Lines

Open New Window

We measured radial velocities of Hα by convolving the line with the derivative of a Gaussian, optimized for a 16 Å FWHM (Schneider & Young 1980); these velocities are given in Table 5. We searched for frequencies in the 0–20 c/d range using a “residualgram” method (Thorstensen & Freed 1985), and found two incommensurate frequencies, near 16.7 and 37.3 c/d. These signals and the resultant velocity curves are shown in the upper three frames of Figure 7. They appear to be coherent over the duration of the observations, and Monte Carlo tests show that there is no significant uncertainty in cycle count on a daily, monthly, or annual timescale.

TABLE 5
TABLE 5 Hα Radial Velocities

Open New Window

Fig. 7.— Period searches of Hα emission velocities. The upper three frames show the raw “residualgram” and folded velocity curves for the two obvious signals. The lower four frames show the same after subtraction of the contaminating signals.

Open New Window

To refine this further we fitted the velocities with sinusoids at the two frequencies, then prepared two new time series with the fitted sinusoid subtracted. Each was searched again, with the result seen in the lower frames of Figure 7. The alias structure (“picket fence”) of the subtracted signal disappeared, and no other significant signals were unmasked. Folds of the new time series gave a much cleaner result (lower right frames of Fig. 7) and better‐determined fit parameters (Table 6). We then simultaneously fitted the original data with two sinusoids, with frequencies fixed at the values in Table 6. Figure 8 shows the data for our best six‐night time series, with the two‐frequency fit superposed.

TABLE 6
TABLE 6 Sinusoid Fits

Open New Window

Fig. 8.— Hα radial velocities from the six nights of longest coverage. Each frame covers the same time interval.

Open New Window

Figure 9 shows gray‐scale representations of the spectra averaged in phase, for both ephemerides; the technique is described by Taylor et al. (1999). At the lower frequency (left panels), s‐waves appear in the Hα and He i λ5876 emission lines. (They are more subtly present in other lines, too.) Similar s‐wave features appear in disk‐accreting CVs, and are thought to arise from the hot spot where the mass‐transfer stream strikes the outer disk. Their appearance here strongly suggests that 86 minutes is the true orbital period. The gray‐scale plots at the higher frequency (right panels) show velocity modulation in the envelope of the Hα and He i λ5876 emission lines, but no s‐wave. There is also a striking modulation in the strength of the He ii λ4686 line at this frequency.

Fig. 9.— Montage of gray‐scale representations of phase‐averaged spectra. The vertical axes show phase increasing vertically, and the horizontal axes are wavelength in Å. Left and right panels are respectively folded on the long (86.1 minute) and short (38.6 minute) period. The top panels show the region around Hβ, He ii λ4686, He i λλ4922 and 5015, and Fe i λ5169; the middle panels show He i λ5876; the bottom panels show Hα. Gray scales are selected to maximize visibility of the periodic features.

Open New Window

4. DISCUSSION

 

4.1. Periods: the Magic Credential

These data prove conclusively the existence of a stable signal at a period much shorter than . This is basically the defining characteristic of the DQ Herculis (intermediate polar) class of cataclysmic variables, in which the fast period is interpreted as the spin of an accreting magnetic white dwarf (Patterson 1994; Chapter 9 of Hellier 2000). In this case, that period is 38.6 minutes. Most of the other periods have a natural place in this picture: 86.1 minutes is , and 69.9 minutes is the low‐frequency orbital sideband (corresponding to ). The latter arises from the frequency with which the rotating searchlight beam from the white dwarf sweeps across the secondary and other structures fixed in the orbital frame.

Our description of the variability appears to be roughly consistent with previous photometry of DW Cnc. Uemura et al. (2002) and RGAC had smaller photometric data sets (respectively 75 and 11 hr, with an unfavorable distribution in time). Neither were well suited for period study, and we suspect that this is why we reached different conclusions. When we simulated their data by selecting subsets of our data with their time distribution, we could (sometimes) produce power spectra resembling Figure 6 of Uemura et al. or Figure 3 of RGAC.10 In detail, these three studies are not in particularly good agreement: Uemura et al. suggested that the signals are quasi‐periodic, and RGAC reported signals at 38.51(2) and 77.37(3) minutes.

The spectroscopy of RGAC is more obviously consistent with ours. In particular, they identify the critical feature of two noncommensurate frequencies, although their data are too heavily aliased to specify those frequencies uniquely. The present data establish this description with long‐term ephemerides.

4.2. Phasing of the Periodic Signals

The basic geometry of the binary is usually quite poorly known in intermediate polars. The secondaries are generally invisible; and the “disk” emission lines are likely to be poor tracers of binary motion, since they are heavily contaminated by Dopper motions in the two columns of infalling gas (the mass‐transfer stream and the magnetically channeled accretion column onto the white dwarf). This is true in DW Cnc also; our ephemeris accurately describes the motion of an emission line, but we do not know what that implies about the true dynamical motions of the stars.

The geometry of the spin pulse can also be complex. Hellier (1997) summarizes the available data. Of the five intermediate polars (IPs) with radial‐velocity variations at , two (AO Psc and FO Aqr) show maximum blueshift at maximum light, as expected in the simplest interpretation of this component (arising from infalling gas onto the “back pole”; see Kim & Beuermann 1995). The other three (EX Hya, BG CMi, PQ Gem) fail to obey this simple relation. DW Cnc appears to join the miscreants in this respect, showing maximum blueshift cycles before maximum light.

Thus, without additional assumptions that we are unwilling to make, the observed phase information does not yet constrain the binary and accretion geometries. A secure measurement of true orbital motion would help a great deal, as would a detection of the spin pulse in X‐rays.

4.3. Distance and Luminosity

Little information is available concerning the star’s distance. Cataclysmic variables tend to follow an empirical law relating to (Figs. 5 and 7 of Patterson 1984), but DW Cnc is unusual because it is magnetic. It is probably wiser to compare with the one star it most closely resembles: EX Hya, the only other short‐period magnetic CV with a long . A recent Hubble Space Telescope parallax puts EX Hya at 65 pc, implying (Beuermann et al. 2003). Since DW Cnc is another magnetic CV of short and long , we might expect a similar , suggesting a distance of ∼200 pc. Assuming a “bolometric” correction appropriate to its color (∼1.2 mag), the luminosity in the high state would be , where is its distance in units of 200 pc.

4.4. He ii Emission and the X‐Ray Searchlight

DW Cnc shows fairly strong He ii λ4686 emission ( Å), which is sometimes used as an identifier of magnetic CVs. But caution is needed on this point. He ii emission is rare in dwarf novae, and common in nova‐like variables (say, at the level of Å; see Fig. 5 of Patterson & Raymond 1985). Stronger emission (say, Å) is generally limited to magnetic CVs and stars of very high (the supersofts). It is worth exploring why. The origin of the line is very likely photoionization, and thus it mainly testifies to an abundance of X‐ray/EUV photons (with energy above 54 eV). This implies a high‐temperature region ( K). Now, in the blackbody approximation, T scales as , where L is the accretion luminosity and f is the fraction of white‐dwarf area that is radiating. So high T can be achieved with high luminosity (high ) or with a very small f. Supersofts choose the former route. A small f is most readily achieved with a magnetic field that channels accretion flow to a very small area around the magnetic pole. This is the main reason11 that magnetic CVs consistently sport strong He ii emission even when L and are quite modest.

But with Å, and with emission ∼10 times stronger, DW Cnc is only of moderately high excitation, not quite in the realm of the well‐established magnetic CVs. The much more salient credential is the stable fast period. And the presence of a large‐amplitude signal at (i.e., the 70 minute period) is perhaps the most compelling credential of all.

Why? The reason is that the 70 minute signal is of large amplitude, nearly as strong as the 38 minute signal that presumably drives it. That can only mean that the 38 minute signal is very powerful indeed. The sideband signal must arise from reprocessing in nonaxisymmetric structures that orbit the white dwarf at exactly . The latter are most simply described as “the secondary,” so let us adopt that assumption. Now, in Roche geometry the secondary subtends a solid angle of ∼3% of sr at the white dwarf, so only a few percent of white‐dwarf luminosity can be reprocessed in the secondary. This inefficiency is basically why “reflection effects” in the secondary are only rarely seen in CV light curves.

Thus, sideband signals are generally expected to be quite weak, contrary to observation in DW Cnc. There are a few rescue strategies that could be considered.

1.

Beaming: maybe the white‐dwarf radiation illuminates the secondary powerfully, but illuminates terrestrial observers quite weakly.

2.

Maybe the absorbing fraction is much higher than 3%. Another discrete structure that orbits with can be the mass‐transfer stream, or more specifically the “hot spot” where the stream impacts the accretion disk. The hot‐spot size is not well constrained, but could possibly subtend as much as 15% of sr.

3.

Maybe the bolometric luminosity of the white dwarf’s 38 minute signal is very high, due to the presence of hitherto unseen EUV and X‐ray flux.

Each of these excuses could play a role. There is not much to say about (1), since none of the relevant angles are known for this binary. Possibly (2) plays some role, but even 15% is still inadequate, as we will see below. The dominant explanation is likely to be (3). And this is very plausible, since accretion onto a magnetic white dwarf is predominantly radial, and high temperatures are generally achieved in radial accretion flow. There must be a powerful UV/EUV/X‐ray searchlight associated with that 38 minute signal.

More specifically, a minimum luminosity in the 70 minute signal can be estimated by assuming a pulsed fraction of 100%, and a bolometric correction of zero. In this case we obtain ergs s−1. More plausibly, let us assume that the 70 minute signal is actually 50% pulsed,12 with a bolometric correction near 0.5 mag. This implies ergs s−1. In order to power this with a prograde searchlight at minutes, with a generous assumption of 15% reprocessing, we probably need a searchlight ∼7 times stronger, or ergs s−1.

Hard X‐rays are a plausible choice for the powering searchlight. The only X‐ray observation is in the ROSAT All‐Sky Survey, in which the star appears to be detected at counts s−1. For an unabsorbed 10 keV spectrum, typical of CVs in general, this corresponds to a 1–10 keV flux of , or ergs s−1. Intermediate polars commonly show a fairly large absorption from cool gas near the emission site; assuming , we obtain ergs s−1. This is somewhat of an upper limit; ignoring uncertainties of source variability, a hard X‐ray source much brighter than this would have been detected by the all‐sky surveys of the 1970s (HEAO‐1, HEAO‐2).

The other strong candidate for the searchlight is the component in the ultraviolet. Again assuming a pulsed fraction of 50%, but with a bolometric correction of 1.5 mag (since the component appears to be slightly bluer than the star’s mean color, whereas the component is slightly redder), we obtain ergs s−1. Thus, we get a likely optical/UV/X‐ray component of ergs s−1—about half of what we need.

This shortfall can easily be hidden in the several uncertainties of bolometric correction and geometry, described above. Even so, we needed to assume that 15% of the searchlight falls on the “secondary” (where we use this term to denote “mass donor plus hot spot”), which could be considered optimistic. It is obvious that the searchlight powers the signal, but we do not yet understand the energetics of this.13

5. SUMMARY

 

1.

We establish strictly periodic signals at 38.6 and 86.1 minutes in the radial velocities of DW Cnc. The presence of stable noncommensurate frequencies is a strong credential—almost a defining credential—for membership in the DQ Her (intermediate polar) class of cataclysmic variables.

2.

Time‐series photometry reveals strictly periodic 38.6 and 69.9 minute periodic signals that maintain a constant phase over the 1 year baseline of measurement. The three frequencies follow a simple relation , to within limits to measurement. This is a well‐known hallmark of the DQ Her class, and suggests that 86.1 minutes should be regarded as , 38.6 minutes as , and 69.9 minutes as the “reprocessing” period (lower orbital sideband). Ephemerides for all three signals have an expiration date around 2006.

3.

The existence of so strong a signal (0.16 mag full amplitude) at 70 minutes implies a strong “searchlight” signal at 38 minutes sweeping over the parts of the binary fixed in the orbital frame (secondary, mass‐transfer stream). This is likely to be powered in part by a pulsed hard X‐ray source (with ergs s−1), but most of the searchlight’s power probably appears at UV and EUV wavelengths. Since all of these bandpasses are unexplored, searches for this component are very desirable!

4.

There is also some evidence for a weak periodic signal at 110.85(9) minutes. This does not seem related to any of the other clocks in the binary, and we leave it as an unsolved problem.

We acknowledge financial support from NSF (AST‐00‐98254) and NASA (NNG04GA83G), and thank Jerry Foote, Thomas Krajci, Donn Starkey, and David Skillman for contributions to the photometry campaign. We dedicate this paper to our dear friend and colleague Bob Fried, who died on 13 November 2003 while piloting his plane to a remote airstrip in Black Canyon, AZ, on one of his many “Angel Flight” trips to transport needy people for medical help. Bob was one of the most charismatic amateur astronomers in the world, and a legend in the world of variable stars. This paper was one of his last major projects, although his data were widely available to the astronomical community, and will surely be appearing in scientific journals for many years to come.

REFERENCES

 
  • 1 Department of Astronomy, Columbia University, 550 West 120th Street, New York, NY 10027; , .

  • 2 Department of Physics and Astronomy, Dartmouth College, 6127 Wilder Laboratory, Hanover, NH 03755; .

  • 3 Center for Backyard Astrophysics (Belgium), Walhostraat 1A, B‐3401 Landen, Belgium; .

  • 4 Center for Backyard Astrophysics (Flagstaff), Braeside Observatory, Post Office Box 906, Flagstaff, AZ 86002; .

  • 5 Department of Physics, King’s University College, 9125 50th Street, Edmonton, AB T5H 2M1, Canada; .

  • 6 Department of Physical Science, Arkansas Tech University, 1701 North Boulder Avenue, Russellville, AR 72801; , .

  • 7 Joint Astronomy Centre, University Park, 660 North Aòhōkū Place, Hilo, HI 96720; .

  • 8 Center for Backyard Astrophysics (Norwich), 35 Sergeants Way, Lisbon, CT 06351; .

  • 9 In this paper we use the old‐style “c/d” notation as a shorthand for cycles day−1. Likewise, we frequently use a truncated HJD date, namely .

  • 10 RGAC report that the signal is “remarkably coherent” over 4 yr; but that was probably a misprint, since their power spectrum spans only 4 days. And indeed, that data cannot substantiate coherence over 4 days, since their Fig. 3 shows dozens of acceptable periods. Effectively, the RGAC photometry tests coherence only for a single 6 hr observation (i.e., just a few cycles).

  • 11 Not the only reason. Hard X‐rays, usually present in equal or greater supply, also contribute; but they are relatively inefficient ionizers of helium (since they waste much more energy per photon).

  • 12 We do not know the true pulse fraction. It must be at least 16%, because that is observed. And it is probably not very close to 100%, because that is hard to achieve with a highly sinusoidal waveform. So ∼50% is a plausible and simple choice.

  • 13 An obvious candidate is a strong EUV component, seen in a few intermediate polars (Haberl & Motch 1995). These are normally very weak in CVs, but stars with magnetic channeling can easily have small f and therefore high T—high enough to ionize helium and thereby avoid the edges that normally quench the emission. This is definitely worth a search.

© 2004. The Astronomical Society of the Pacific. All rights reserved. Printed in U.S.A.