ABSTRACT
We present optical spectroscopy of nine cataclysmic binary stars, mostly dwarf novae, obtained primarily to determine orbital periods Porb. The stars and their periods are LX And, 0.1509743(5) days; CZ Aql, 0.2005(6) days; LU Cam, 0.1499686(4) days; GZ Cnc, 0.0881(4) days; V632 Cyg, 0.06377(8) days; V1006 Cyg, 0.09903(9) days; BF Eri, 0.2708804(4) days; BI Ori, 0.1915(5) days; and FO Per, for which Porb is either 0.1467(4) or 0.1719(5) days. Several of the stars proved to be especially interesting. In BF Eri, we detect the absorption spectrum of a secondary star of spectral type K3 ± 1 subclass, which leads to a distance estimate of ∼1 kpc. However, BF Eri has a large proper motion (∼100 mas yr−1), and we have a preliminary parallax measurement that confirms the large proper motion and yields only an upper limit for the parallax. BF Eri's space velocity is evidently large, and it appears to belong to the halo population. In CZ Aql, the emission lines have strong wings that move with large velocity amplitude, suggesting a magnetically channeled accretion flow. The orbital period of V1006 Cyg places it squarely within the 2 to 3 hr "gap" in the distribution of cataclysmic binary orbital periods.
Export citation and abstract BibTeX RIS
1. INTRODUCTION
Cataclysmic variables (CVs) are binary star systems in which the secondary, usually a late‐type main‐sequence star, fills its Roche lobe and loses mass to the white dwarf primary (Warner 1995). CVs are long‐lived systems that are stable against mass transfer, so the mass transfer must be driven by gradual changes in the orbit, or in the secondary star, or both. It is commonly believed that the evolution of most CVs is driven by the slow loss of angular momentum from the orbit, most likely through magnetic braking of the corotating secondary star, at least at longer orbital periods Porb, where gravitational radiation is ineffective (Andronov & Pinsonneault 2004 give a recent discussion). The loss of angular momentum constricts the Roche critical lobe around the secondary and causes the system to transfer mass as it evolves toward shorter Porb. In this scenario, Porb serves as a proxy measurement for the system's evolutionary state. Correct and complete orbital period measurements are fundamental to any accurate theory of CV evolution. Given the usefulness of Porb, it is fortunate that it can usually be measured accurately and precisely.
This paper presents optical spectroscopy of the nine CVs listed in Table 1. We took these observations mostly for the purpose of finding orbital periods using radial velocities (none of these systems are known to eclipse). The long cumulative exposures also allowed us to look for any unusual features. The Catalog and Atlas of Cataclysmic Variables, Archival Edition (Downes et al. 2001)2 lists seven of the stars as dwarf novae, one as either a dwarf nova or a DQ Her star, and one simply as a cataclysmic, possibly a dwarf nova similar to U Gem or SS Cygni (type UGSS). Except for CZ Aql, for which we confirm a 4.8 hr candidate period suggested by Cieslinski et al. (1998), all of these objects lacked published orbital periods when we began working on them. Subsequently, Tappert & Bianchini (2003) found Porb = 0.0883 days for GZ Cancri; we had communicated our advance findings to these authors so they could disambiguate their period determination.
![]() |
2. OBSERVATIONS, REDUCTIONS, AND ANALYSIS
2.1. Observations
All our spectra were taken at the MDM Observatory on Kitt Peak, Arizona, using either the 1.3 m McGraw‐Hill Telescope or the 2.4 m Hiltner Telescope. The earliest observations we report here are from 1995, and the latest were obtained 2007 January. Table 2 gives a journal of the observations.
![]() |
![]() |
At the 1.3 m, we used the Mark III spectrograph and a SITe 1024 × 1024 CCD detector. The spectral resolution is 5.0 Å, covering a range of either 4480–6760 Å with 2.2 Å pixel−1 for the 2001 December BF Eri data, or 4646–6970 Å with 2.3 Å pixel−1 for the remaining data.
The 2.4 m spectra, except for those of FO Per, were obtained using the modular spectrograph and a SITe 20482 CCD detector, with 2.0 Å pixel−1 over a range of 4210–7500 Å and with a spectral resolution of 3.5 Å. The relatively small number of 2.4 m spectra of FO Per were taken with a Loral 20482 detector and cover from 4285 to 6870 Å at 1.25 Å pixel−1.
2.2. Reductions
For the most part, we reduced the spectra using standard IRAF3 procedures. The wavelength calibration was based on exposures of Hg, Ne, and Xe lamps. Prior to 2003, we took lamp exposures through the night and whenever the telescope was moved. For the 2.4 m data from 2003 to the present, we used lamp exposures taken in twilight to find the shape of the pixel‐to‐wavelength relation, and set the zero point individually for each nighttime exposure using the O i λ5577 night‐sky feature. The apparent velocity of the telluric OH emission bands at the far red end of the spectrum, found with a cross‐correlation routine, provided a check; although these are far from the feature used to set the zero point, their apparent velocity typically remain within 10 km s−1 of zero. Because of the increased efficiency of this technique, we attempted to use it at the 1.3 m telescope as well during the 2004 June/July observing run. For unknown reasons, the results were unsatisfactory. To salvage the Hα emission velocities from that run, we determined a correction by cross‐correlating the night‐sky emission features in the 6200–6625 Å range with a well‐calibrated night‐sky spectrum obtained with a similar instrument. The correction was calculated for each individual spectrum and then applied to each measured velocity, and it did reduce the scatter somewhat, evidently because the wavelength range used includes the Hα emission line for which we measured velocities.
On all our runs, we observed flux standards during twilight when the sky was clear, and applied the resulting calibration to the data. The reproducibility of these observations suggests that our fluxes are typically accurate to ±20%. We also took short exposures of bright O and B stars in twilight to map the telluric absorption features and divide them out approximately from our program object spectra. Before flux calibration, we divided our program star spectra by a mean hot‐star continuum in order to remove the bulk of the response variation. Table 1 lists V magnitudes synthesized from our mean spectra, using the IRAF sbands task and the passband tabulated by Bessell (1990); clouds, losses at the slit, and calibration errors make these uncertain by a few tenths of a magnitude, but they do give a rough indication of the brightness of each system at the time of our observation.
2.3. Analysis
Except for a few spectra taken in outburst (which show weak emission or absorption on a strong continuum), all of the stars show the prominent emission lines. Figures 1 and 2 show averaged spectra, and Table 3 gives the equivalent width and FWHM of each line measured for each star from its averaged spectrum.
Fig. 1.— Plots of the average flux‐calibrated spectra for eight of the stars studied here. The weak features seen near λ5577 are artifacts caused by imperfect subtraction of the strong [O i] night‐sky emission.
Fig. 2.— Plot of the averaged spectrum of BF Eri (top) and the spectrum after a scaled late‐type (K3 V) star has been subtracted (bottom). The spectra have been shifted into a rest frame before averaging and do not include the 2006 March or 2007 January data.
![]() |
![]() |
Two stars, BF Eri and BI Ori, showed the spectral features of a late‐type star. To quantify the secondary contribution in these objects, we began by preparing averaged flux‐calibrated spectra (in BF Eri's case, the secondary's radial velocity curve was measurable, so we shifted the individual spectra to the secondary's rest frame before averaging). Over time we have used the 2.4 m and modular spectrograph to collect spectra of K and M stars classified by Keenan & McNeil (1989) or Boeshaar (1976). The wavelength coverage and spectral resolution of these data are similar to the 1.3 m data. We applied a range of scaling factors to the library spectra, subtracted them from the averaged spectra, and examined the results by eye to estimate a range of spectral types and scaling factors giving acceptable cancellation of the late‐type features.
We use the spectral type and secondary flux to estimate the distance in the following manner. We begin by finding the surface brightness of the secondary star in V, on the assumption that the surface brightness is similar to that of main‐sequence stars of the same spectral type; the Barnes‐Evans relation for late‐type stars is discussed by Beuermann (2006). Combining the known Porb with the assumption that the secondary fills its Roche critical lobe yields the secondary's radius R2 as a function of its mass, M2. In the relevant range of mass ratio, R2∝M1/32, approximately, and the dependence on M1 is weak enough to ignore. We generally do not know M2, so we guess at a generous allowable range for this parameter using evolutionary simulations by Baraffe & Kolb (2000) as a guideline; the weakness of the dependence of R2 on M2 means that this (rather questionable) step does not dominate the error budget. Combining the surface brightness with R2 yields the absolute magnitude MV. Subtracting this from the apparent magnitude measured for the secondary star gives a distance modulus. The reddening maps of Schlegel et al. (1998) can then be used to estimate the extinction. Note carefully that we do not assume that the secondary is a "normal" main‐sequence star; we assume only that the secondary's surface brightness is similar to field stars of the same spectral type. The normalization of the secondary's contribution also depends on the assumption that the spectral features used to judge the subtraction are similar in strength to those of a normal star.
As noted earlier, the immediate aim of our observations was to find orbital periods from radial velocity time‐series data. The Hα emission line is usually the strongest feature, and it generally gives good results in dwarf novae. All the emission‐line velocities reported here are of Hα.
We measured radial velocities of Hα emission using convolution methods described by Schneider & Young (1980) and Shafter (1983). In this technique, one convolves an antisymmetric function with the line profile and takes the zero of the convolution (where the two sides of the line contribute equally) as the line center. For the antisymmetric function with which the spectrum is convolved, we used either the derivative of a Gaussian with adjustable width, or positive and negative Gaussians of adjustable width offset from each other by an adjustable separation. Uncertainties in the convolution velocities are estimated by propagating forward the counting‐statistics errors in the individual data channels; in practice, these are lower limits to the true uncertainties, since the line profile can vary in ways unrelated to the orbital modulation. The choice of convolution parameters is dictated by the shape and width of the line, and in practice the parameters are adjusted to give the best detection of the orbit. The physical interpretation of CV emission lines is complicated and controversial (see, e.g., Shafter 1983; Marsh 1988; Robinson 1992), but in almost all cases the emission‐line periodicity accurately reflects Porb (although Araujo‐Betancor et al. [2005] describe a noteworthy exception to this rule). A sample of the radial velocities for each object are listed in Table 4, while the full tables can be found online.
![]() |
One of our systems, BF Eri, has a K‐type absorption component in its spectrum. We measured its velocities with the cross‐correlation radial velocity package described by Kurtz & Mink (1998), using the region from 5000 to 6500 Å and excluding the region containing the He i λ5876 emission line and the Na D absorption complex. For a cross‐correlation template spectrum, we used a velocity‐compensated sum of many observations of IAU velocity standards taken with the same instrument, as described in Thorstensen et al. (2004).
We searched for periods in all the velocity time series using the "residualgram" method (Thorstensen et al. 1996); the resulting periodograms are given in Figures 3 and 4. At the best candidate periods, we fitted least‐squares sinusoids of the form v(t) = γ + Ksin [2π(t - T0)/P]. Figure 5 shows the velocities folded on the best‐fitting periods, and Table 5 gives the parameters of these fits. Because of limitations of the sampling (e.g., the need to observe only at night from a single site), a single periodicity generally manifests as a number of alias frequencies. To assess the confidence with which we could assert that the strongest alias is the true period, we used a Monte Carlo test described by Thorstensen & Freed (1985).
Fig. 3.— Periodigrams for most of the stars studied here. The vertical axis in each case is the inverse of χ2 for the least‐squares best‐fitting sinusoid at each trial frequency. When data from more that one observing run are combined, the periodigram can require hundreds of thousands of points to resolve the fine‐scale ringing; in those cases, the curve shown is formed by connecting local maxima of the periodogram with straight lines. In those cases, the right‐hand panel gives a close‐up view of the region around the highest peak, revealing the alias structure resulting from different choices of cycle count between the observing runs. The periodogram of BF Eri is for the absorption‐line velocities.
Fig. 4.— Periodigrams for the remainder of the stars, plotted in the same manner as the previous figure. Because the choice of daily cycle count for FO Per remains ambiguous, we have not chosen to enlarge either peak region.
Fig. 5.— Radial velocities plotted as a function of phase using the adopted orbital periods. For CZ Aql, GZ Cnc, and V1006 Cyg, the number of cycle counts between observing runs is unknown, and the exact period chosen to fold the velocities is one of a number of possibilities. The two plots for FO Per are for different choices of the daily cycle count, and each of these in turn is also an arbitrary choice among many finely spaced periods. In BF Eri, both emission and absorption velocities are plotted; the absorption velocities are shown with error bars.
![]() |
The alias problem can be particularly irksome over longer timescales; in this case, the uncertain number of cycles that have elapsed between observing runs causes fine‐scale "ringing" in the periodogram. The individual periods have tiny error bars because of the large time span covered, but the ambiguity in period means that a realistic error bar—one that covers the range of possibilities—is much larger. In those cases, the period uncertainties given in Table 5 are estimated by analyzing data from the individual observing runs separately. When only two observing runs are available, the allowable fine‐scale frequencies are well described by a fitting formula
![](https://content.cld.iop.org/journals/1538-3873/119/855/494/revision1/df1.jpg)
where t1 and t2 are the epochs of blue‐to‐red velocity crossing observed on the two runs, and n is the integer number of cycles that have passed between t1 and t2. The allowed range of n is determined from the weighted average of the periods derived from separate fits to the two runs' data. When more than two observing runs are available, the situation becomes more complex. In some happy cases, there are enough overlapping constraints that only a single, very precise period remains tenable. We were able to find such precise periods for LX And, LU Cam, and BF Eri.
3. NOTES ON INDIVIDUAL OBJECTS
We discuss the stars in alphabetical order by constellation.
3.1. LX Andromedae
LX And was first identified as a variable star (RR V‐3) in the Lick RR Lyrae search (Kinman et al. 1982). It was classified incorrectly as an RV Tauri star, and its dwarf nova nature was unrecognized until the photometric study by Uemura et al. (2000). Morales‐Rueda & Marsh (2002) obtained spectra of LX And as part of their study of dwarf novae in outburst and determined the equivalent widths and FWHMs of the Balmer and He ii lines. Our mean spectrum appears typical for a dwarf nova at minimum light.
Because of the large hour‐angle span, the radial velocity time series leaves no doubt about the daily cycle count, which is near 6.6 cycles day−1. The several observing runs constrain the fine‐scale period in a more complicated way, but the Monte Carlo test indicates that a precise period of 0.1509743(5) days is preferred, with about 98% confidence. Two other candidate periods separated from this by 1 cycle per 53.2 days in frequency are much less likely.
3.2. CZ Aquilae
Very little has been published on CZ Aql, which is listed in the Archival Edition as a U Gem dwarf nova. Cieslinski et al. (1998) included the star in their spectroscopic study of irregular variables and noted a probable 4.8 hr period and emission lines typical of dwarf novae. Our velocities confirm the suggested 4.8 hr period, but we cannot determine a unique cycle count between our observing runs.
While the spectrum superficially resembles that of a dwarf nova, a closer look reveals interesting behavior. Figure 6, constructed using methods described by Taylor et al. (1999), presents our spectra as a phase‐averaged gray‐scale image. There is a striking broad component in the stronger Balmer and He i lines that shows a large velocity excursion, with the red wing of Hα reaching to +3100 km s−1 at phase 0.3 (where phase 0 corresponds to the blue‐to‐red crossing of the line core). The broad components around Hβ and λ6678 move in phase with those of Hα and range from 900 to 2600 and −2500 to −900 km s−1 and 700 to 2100 and −1000 to −600 km s−1, respectively. The wings of λ5876 also move in phase with the others, but the red edge is difficult to follow at its minimum because of interference from the Na D absorption lines, which are stationary and hence interstellar. The maximum of the red edge is 3200 km s−1, while the blue edge ranges from −1500 to −600 km s−1. The blueward wing of all these lines is noticeably weaker than the redward wing.
Fig. 6.— Phase‐averaged spectra of CZ Aql (top two panels) and BF Eri (bottom two panels), presented as a gray scale. The scale is inverted so that emission is represented by darker shades. The two CZ Aql spectra are scaled differently to show the line cores (top) and the extent of the line wings. Note the Na D lines in CZ Aql remain stationary, indicating an interstellar origin. The feature at 6280 Å is telluric. BF Eri's spectrum is plotted in two overlapping sections; the K‐star's orbital motion is plainly visible.
Other emission lines present include He ii λ4686, He i λλ4713 and 4921, and very weakly, Fe ii λ5169. We also detect unidentified emission lines at 6344 Å, as is also seen in LS Peg (Taylor et al. 1999), and at 5046 Å. The strength of the λ5780 diffuse interstellar band (Jenniskens & Desert 1994) and the Na D lines suggest that a good deal of interstellar material lies along the line of sight, and that the luminosity is relatively high.
High‐velocity wings reminiscent of the ones seen here have been seen in V795 Her (Casares et al. 1996; Dickinson et al. 1997), LS Peg (Taylor et al. 1999), V533 Her (Thorstensen & Taylor 2000), and RX J1643+34 (Patterson et al. 2002), all of which are SW Sex stars. We do not, however, detect another SW Sex characteristic, namely phase‐dependent absorption in the He i lines (Thorstensen et al. 1991). The orbital periods of most SW Sex stars are shorter than 4 hr, so CZ Aql's 4.8 hr period would be unusually long for an SW Sex star.
3.3. LU Camelopardalis
Jiang et al. (2000) obtained the first spectrum of this dwarf nova in a follow‐up study of CV candidates from the ROSAT All‐Sky Survey. We found no other published spectroscopic studies. Our velocities constrain the period to a unique value, 0.1499685(7) days. The averaged spectrum shows a rather strong, blue continuum, which may indicate a state somewhat above true minimum.
3.4. GZ Cancri
Jiang et al. (2000) confirmed the cataclysmic nature of GZ Cnc by obtaining the first spectrum of the object. Kato et al. (2002) suggested that this star, originally labeled as a dwarf nova, could possibly be an intermediate polar (DQ Her star), based on similarities in its long‐term photometric behavior to that of other intermediate polars. Tappert & Bianchini (2003) conducted a photometric and spectroscopic study of the system. Using advance results from the present study to help decide the daily cycle count, they found Porb = 0.08825(28) days, or 2.118(07) hr, placing the system near the lower edge of the so‐called gap in the CV period distribution—a dearth of systems in the period range from roughly 2 to 3 hr. Tappert & Bianchini (2003) also saw characteristics that could indicate an intermediate polar classification but did not claim their evidence was definitive on this point.
Almost all our observations come from two observing runs a year apart. The full set of velocities strongly indicates an orbital frequency near 11.4 cycles day−1, with the Monte Carlo test giving a discriminatory power greater than 0.99 for the choice of daily cycle count. However, the number of cycles between the two observing runs is not determined. Precise periods that fit the combined data set are given by P = [349.785(3) days]/n, where n is the integer number of cycle counts; n = 3972 ± 8 corresponds to roughly 1 standard deviation. While our period agrees well with that of Tappert & Bianchini (2003), our data neither support nor disprove the claim that GZ Cnc may be an intermediate polar.
3.5. V632 Cygni
Liu et al. (1999) offer the only published spectrum of this dwarf nova. They measured the equivalent widths and integrated line fluxes of the Balmer, He i, and He ii emission lines and suggested that the orbital period is likely short, based on the very strong Balmer emission. Our spectrum appears similar to theirs, and our measured flux level is also nearly the same. The periodigram in Figure 3 clearly favors an orbital frequency near 15.7 cycles day−1, with a discriminatory power of 95% and a correctness likelihood near unity. This confirms the suggestion of Liu et al. (1999) that the period is rather short and suggests that it is an SU UMa–type dwarf nova.
3.6. V1006 Cygni
Bruch & Schimpke (1992) present the only published spectrum we know of and characterized it as a "textbook example" of a dwarf nova spectrum. They noted a slightly blue continuum with strong Balmer and He i emission, as well as clear He ii λ4686 and Fe ii emission. Our spectrum (Fig. 1) is similar to theirs in both appearance and normalization, and our line measurements (Table 3) are also comparable.
The periodogram (Fig. 4) indicates a frequency near 10.1 cycles day−1, and the Monte Carlo test confirms that the daily cycle count is securely determined. Most of our data are from 2004 June, but we returned in 2005 June/July to confirm the unusual period indicated in the earlier data. The periods found by analyzing the two runs separately are consistent within their uncertainties. As with GZ Cnc, there are multiple choices for the cycle count between the two observing runs; the best‐fitting periods are given by P = [369.006(4) days]/n, where n = 3726 ± 4 corresponds to 1 standard deviation. Including a few velocities from other observing runs suggests that n is slightly larger, perhaps 3728. In any case, the period amounts to 2.38 hr, which places V1006 Cyg firmly in the period gap (Warner 1995), where there is apparently a true scarcity of dwarf novae (Hellier & Naylor 1998).
3.7. BF Eri
The first evidence that BF Eridani was a cataclysmic variable came when an Einstein X‐ray source, 1ES 0437−046, was matched to the variable (Elvis et al. 1992). Schachter et al. (1996) confirmed this match and presented an optical spectrum. Kato (1999) and the Variable Star Observers League in Japan (VSOLJ) found photometric variability characteristic of a dwarf nova.
The spectrum of BF Eri (Fig. 2) shows a significant contribution from a K star, along with the usual dwarf‐nova emission lines. Normally, this suggests that Porb>6 hr. Nearly all our spectra yielded good cross‐correlation radial velocity measurements as well as emission‐line velocities. The absorption‐ and emission‐line velocities independently give a period near 6.50 hr (Table 5), in accordance with expectation based on the spectrum. There is no ambiguity in cycle count over the 5 yr span of the observations, so the period is precise to a few parts per million. Figure 6 shows a phase‐resolved average of the BF Eri spectra, with the absorption spectrum shifting in antiphase to the emission lines.
If the emission‐line velocities faithfully trace the primary's center‐of‐mass motion, and the absorption‐line velocities also trace the secondary's motion, then the two velocity curves should be exactly one‐half cycle out of phase. In BF Eri, we find a shift of 0.515 ± 0.007 cycles between the two curves, consistent with 0.5 cycles, so we feel emboldened to explore the system dynamics.
Masses can only be derived when the orbital inclination is known, as in eclipsing systems. To see if BF Eri might eclipse, we derived differential magnitudes from images that were taken for astrometry (discussed below) and plotted them as a function of orbital phase. Some images were taken at the phase at which an eclipse would appear, but no evidence for an eclipse was found. Limits on the depth and duration of the eclipse are difficult to quantify, because the data were taken in short bursts in the presence of strong intrinsic variability, so a weak eclipse cannot be ruled out, but the photometry does suggest that the inclination is not close to edge‐on.
Because the system apparently does not eclipse, we cannot derive masses; rather, we find broad constraints on the inclination by assuming astrophysically reasonable masses for the components. Taken at face value, the velocity amplitudes K imply a mass ratio q = M2/M1 = 0.60 ± 0.03. If we arbitrarily choose a white dwarf mass M1 = 0.9 M⊙ (so that M2 = 0.53 M⊙), the observed K velocities imply i = 50°. To find a rough lower limit on the inclination, we consider a massive white dwarf (M1 = 1.2 M⊙) and, ignoring the constraint on q for the moment, take M2 = 0.4 M⊙; this yields i = 40°. For a rough upper limit, we assume M1 = 0.6 M⊙ and M2 = 0.4 M⊙, which gives i = 67°.
The decomposition procedure described earlier yielded a spectral type K3 ± 1 subclass; the result of the subtraction is shown in Figure 2. Using the V passband tabulated by Bessell (1990) and the IRAF task sbands, we find a synthetic V = 16.9 ± 0.3 for the K star's contribution. Taking the range of plausible secondary star masses to be 0.4–0.8 M⊙ yields R2 = 0.7 ± 0.1 R⊙ at this Porb. Combining this with the surface brightness expected at this spectral type yields MV = 6.8 ± 0.4 for the secondary. If there is no significant interstellar extinction, we have m - M = 10.1 ± 0.5, or a distance of approximately 1100 ± 300 pc. The dust maps of Schlegel et al. (1998) give a total E(B - V) = 0.062 in this direction. Assuming that BF Eri is beyond the Galactic dust, and taking AV/E(B - V) = 3.3, gives an extinction‐corrected (m - M)0 = 9.9 and a distance estimate of 950+250-200 pc.
We can also estimate a distance using the relation found by Warner (1987) between Porb, i, and the absolute magnitude at maximum light, MV(max). Using our inclination constraints, the Warner relation predicts MV(max) = 3.9 ± 0.7 at this orbital period. The General Catalog of Variable Stars (GCVS)4 lists mp = 13.2 at maximum light; taking this to be similar to Vmax yields m - M = 9.3, or 9.1 corrected for extinction, which corresponds to 660 pc.
Given these distance estimates, it is surprising that BF Eri has a very substantial proper motion. The Lick proper motion survey (Hanson et al. 2004) gives (μX,μY) = (+34, -97) mas yr−1. We have begun a series of parallax observations with the Hiltner 2.4 m telescope, using the protocols described by Thorstensen (2003); so far, we have five epochs from 2005 November and 2007 January. The proper motion relative to the background stars is (μX,μY) = (32, -111) mas yr−1, and the parallax is not detected, with a nominal value of 1 ± 2 mas. The parallax determination is very preliminary, but given the data so far, we estimate the lower limit on the distance based on the astrometry alone to be ∼200 pc.
At the nominal 950 pc distance derived from the secondary star, a 100 mas yr−1 proper motion corresponds to a transverse velocity vT = 451 km s−1. This is implausibly large, so we are left wondering how we might have overestimated the distance. One effect might be as follows. Our distance is based on the secondary's apparent brightness, and we estimate the secondary's contribution to the total light by searching for the best cancellation of its features. If the secondary's absorption lines are weaker than those in the spectral‐type standards, we would underestimate the secondary's contribution. In our best decomposition, the secondary is about 2.2 mag fainter than the total light in V. Assuming (unrealistically) that all the light is from the secondary would therefore decrease the distance modulus by 2.2 mag, to a distance of 340 pc.
We do not yet have enough information to resolve the conundrum posed by BF Eri's unmistakably large proper motion and its apparently large distance, but a reasonable compromise might be to put it at something like 400–500 pc, with an underluminous, low‐metallicity secondary. The cross‐correlation velocities of the secondary have a zero point determined to ±5 km s−1, more or less, and give a substantial systemic velocity of −72 ± 3 km s−1, or −86 km s−1 in the local standard of rest. If the star is at 450 pc, its space velocity with respect to the local standard of rest is ∼250 km s−1, with Galactic components (U,V,W) = (-180, -180, -3) km s−1; that is, the velocity is mostly parallel to the Galactic plane and lags far behind the rotation of the Galactic disk. This would put BF Eri on a highly eccentric orbit; these are halo‐population kinematics (even though the star remains close to the plane). These kinematics would be qualitatively consistent with the weak‐line conjecture used earlier.
3.8. BI Ori
Szkody (1987) published the first quiescent spectrum of BI Orionis, which showed emission lines typical of dwarf novae. Morales‐Rueda & Marsh (2002) show a spectrum in outburst and note the possible presence of weak He ii λ4686 emission.
Only the 2006 January velocities are extensive enough for period‐finding; they give Porb = 4.6 hr, with no significant ambiguity in the daily cycle count. The average spectrum shows the usual emission lines; M dwarf absorption features are also visible, although the signal‐to‐noise ratio of the individual spectra was not adequate for finding absorption‐line velocities. Using the procedures described earlier, we estimate the secondary's spectral type to be M2.5 ± 1.5, with the secondary alone having V = 20.0 ± 0.4. Assuming that the secondary's mass lies in the broad range from 0.2 to 0.6 M⊙, its radius at this Porb would be 0.35–0.6 R⊙. Combining this with the surface brightness derived from the spectral types gives an absolute magnitude MV = 10.2 ± 1.0. The distance modulus, uncorrected for extinction, is therefore m - M = +9.8 ± 1.1, corresponding to 910+600-360 pc. Schlegel et al. (1998) estimate a total reddening E(B - V) = 0.11 in this direction; assuming that BI Ori lies beyond all the dust, and taking AV/E(B - V) = 3.3, reduces the distance to ∼770 pc. At maximum light, BI Ori has mp = 13.2 (GCVS). Assuming the color is neutral, we find MV = 3.4 ± 1.1 at maximum. At BI Ori's period, the Warner (1987) relation predicts MV>3.6 (with the brightest value corresponding to i = 0). This agrees broadly with our nominal value based on the secondary star's distance, but is a little fainter, suggesting that BI Ori is not too far from face‐on, or a little closer than our nominal distance, or both.
3.9. FO Per
FO Persei was apparently discovered by Morgenroth (1939), but its cataclysmic nature was not immediately recognized. Bruch (1989) obtained spectra and gave equivalent widths for the Balmer lines for two different nights of observations, between which the continuum changed from relatively flat to inclined toward the red.
The emission lines in FO Per are rather narrow (Fig. 1, Table 3). This is often taken to indicate a low orbital inclination. The velocity amplitude K is small, so that K/σ ≈ 1.6 for the best fits (Table 5). Because of this, the daily cycle count remains ambiguous; the orbital frequency is either 5.8 or 6.8 cycles day−1, corresponding to Porb of 3.52 or 4.13 hr. CVs with periods in the 3–4 hr range tend to be nova‐like variables (Shafter 1992), whereas FO Per is a dwarf nova; thus, the 4.13 hr period is more likely a priori.
4. SUMMARY
We have determined the orbital periods of eight CVs without significant daily cycle count ambiguity; for FO Per, the period is narrowed to two choices. For three of the systems, we find high‐precision periods by establishing secure cycle counts over long baselines.
While most of these objects are similar to others already known, three stand out as especially interesting. CZ Aql shows asymmetric, high‐velocity wings around the Balmer and He i λλ5876 and 6678 lines, possibly indicating a magnetic system. BF Eri's proper motion of ∼100 mas yr−1 is surprising in view of the large distance indicated by its secondary spectrum and by the Warner relation; even if it is somewhat nearer than these indicators suggest, its kinematics are not typical of disk stars. Finally, the orbital period of V1006 Cyg places it squarely in the middle of the so‐called period gap between 2 and 3 hr.
We are most grateful for support from the National Science Foundation through grants AST 99‐87334 and AST 03‐07413. Bill Fenton took most of the spectra of GZ Cnc, and J. Cameron Brueckner assisted with the BF Eri spectroscopy. Some of the astrometric images of BF Eri were obtained by Sébastien Lépine and Michael Shara of the American Museum of Natural History. We would like to thank the MDM Observatory staff for their skillful and conscientious support. Finally, we are grateful to the Tohono O'odham for leasing us their mountain for a while, so that we may study the glorious universe in which we all live.
Online Material
- Table 4
Footnotes
- 1
Based on observations obtained at the MDM Observatory, operated by Dartmouth College, Columbia University, Ohio State University, Ohio University, and the University of Michigan.
- 2
Available at http://archive.stsci.edu/prepds/cvcat/index.html; this had been called the "Living Edition" until its author retired and ceased updates.
- 3
IRAF is distributed by the National Optical Astronomy Observatories.
- 4
VizieR Online Data Catalog II/214A (P. N. Kholopov et al., 1999)