# The Orbital Period of HDE226868/Cyg X-1

## 1 Introduction

Cygnus X-1, identified with the bright (V$``$8) star HDE226868 (Bolton 1972; Webster & Murdin 1972), has long been regarded as the best black hole candidate among the high-mass X-ray binaries. As such, it has been an object of extensive observation over the past two and a half decades. Perhaps surprisingly, there remain some important uncertainties and discrepancies in the derived properties of the system. In particular, Ninkov, Walker, & Yang (1987, hereafter NWY) report evidence of possible period variation, discussed below, and additional periodicities on timescales ranging from 39 days to 4.5 years have been suggested by Kemp, Herman, & Barbour (1978), Wilson and Fox (1981), Priedhorsky, Terrell, & Holt (1983), and Walker and Quintanilla (1978).   

HDE226868 is a single-lined spectroscopic binary. The optical spectrum has been classified O9.7Iab by Walborn (1973), with variable emission at HeII $`\lambda `$4686 and, less prominently, at the Balmer lines. The radial velocity period of 5.6 days has been established for over two decades, and the definitive orbital elements and ephemeris were published by Gies and Bolton (1982, hereafter GB). They give a period of 5.59974 $`\pm `$ 0.00008 days, a precision which translates to an uncertainty of only $`\pm `$0.02 cycle in phase at the present time.   

However, NWY, in a subsequent radial-velocity study, reported a somewhat longer period of 5.6017 $`\pm `$ 0.0001 days, and suggest the possibility of a period increase with time. These results have the effect that the orbital phase of the Cyg X-1 system is then highly uncertain at the present epoch. The phase calculated using the ephemeris of NWY differs from the GB phase by 0.5 cycle, and if the period is in fact varying, the phase at the present epoch is completely indeterminate.   

If correct, a time-varying period would have important implications for the mass transfer rate and hence evolution of the binary. Furthermore, with this indeterminacy in phase, it is impossible to ascribe and interpret accurately features in the X-ray light curve such as the dipping behaviour, in particular the distribution of dipping with orbital phase (e.g. Remillard and Canizares 1984). To resolve these ambiguities we therefore undertook a program of observations to re-establish the orbital ephemeris of HDE226868. These data also allowed us to undertake a more detailed investigation of the long-term stability of the orbit, as high quality data now exist over a $``$30 year baseline. Our new ephemeris has thereby allowed us to complete a survey of the distribution of X-ray dipping with orbital phase (Balucinska-Church et al. 1998).   

## 2 Observations

Through the La Palma Service Programme, and with the assistance of regularly scheduled observers, we obtained CCD spectra of HDE226868 using the Intermediate Dispersion Spectrograph on the Isaac Newton Telescope on 17 nights in May and June 1996. These included 11 consecutive nights covering two complete orbital cycles. In all we obtained 37 spectra of the object, each covering a spectral range of 4100 to 4900Å, with a dispersion of approximately 0.8 Å/pixel. Exposures varied from 100 to 200 seconds, and yielded a S/N in excess of 100. In addition, we observed 19 Cep, a radial velocity standard of similar spectral type, on two nights, for reference and calibration purposes.   

The spectra were extracted and reduced using standard IRAF routines. The extracted spectra were wavelength calibrated using a Cu-Ar arc recorded immediately before and after the program spectra, and then normalized. Our wavelength calibration accuracy had rms residuals of $``$0.10Å. Typical HDE226868 and 19 Cep spectra are shown in Figure 1.   

## 3 Period Determination

We determined radial velocities by cross-correlation of the HDE226868 spectra with the 19 Cep spectra, using IRAF routines. Separate determinations were made for the hydrogen and helium lines. We used the HeI absorption lines at 4388, 4472, 4713 and 4921Å. The resulting velocities are given in Table 1, and were reduced to heliocentric velocities using the catalogue value (Hoffleit 1982) for the radial velocity of 19 Cep. We performed period searches on both sets of velocities, as well as on our velocities in combination with previously determined velocities from the literature. We used several different period-search techniques, including Scargle periodograms, Fourier searches, and chi-square minimization to a circular orbit fit. All methods yielded similar results, as follows.   

The best fit period to the hydrogen line velocities for our data alone is 5.566 $`\pm `$ 0.012 days, while the best fit to our helium line data is 5.629 $`\pm `$ 0.015 days. The apparent discrepancy here is startling, as is the amount of apparent change in period from the canonical value in either case. Both of these issues can be readily resolved.   

It is well established that the hydrogen lines are contaminated by emission whose velocity is in approximate antiphase with the absorption lines (GB; Hutchings, Crampton, & Bolton 1979). This produces variable and irregular line profiles which yield spurious velocities and periods when cross-correlated with standard spectra. So any period determined from the hydrogen line velocities is highly suspect. For this reason we reject all velocity determinations based on the hydrogen lines, and in working with historical data below, we use only published velocities determined for the helium lines.   

The fact that the best-fit period to the helium lines seems to be $`2\sigma `$ larger than previously determined values should also be viewed with some scepticism. In fact, with a dense grid of observations over a relatively short period of time, the quality of fit is good over a much broader period range than suggested by the quoted uncertainty. Fitting our data to the NWY period, or to the GB period, produces $`\chi ^2`$ residuals which are not significantly worse than those for our best fit. In fact, the formal uncertainty calculated by the fitting routines is undoubtedly too small; this is because the scatter of the data is such that even for the best orbital models the resulting $`\chi ^2`$ is too large for the formal uncertainty calculations to be valid. We shall see below that this is not unique to our results but in fact has been a problem for most, and perhaps all, previous published orbital solutions for Cyg X-1; uncertainties have been systematically underestimated, resulting in apparent discrepancies where none really exist.   

To test for the possibility of period variation, we combined our observations with all previously published helium-line velocity data (Brucato and Zappala 1974; GB; NWY) and performed a $`\chi ^2`$ minimization weighted fit to a circular orbit both with a fixed period and with a period which varies linearly with time, using the MINUIT package of function-fitting routines (James 1994). The results are summarized in Table 2 and plotted in figure 2. The best fixed period is 5.5998 $`\pm `$ 0.0001 days, with a $`\chi ^2`$ of 4284; the best variable period is 5.5997 $`\pm `$ 0.0001 days, with $`\dot{P}`$ of 3.8 $`\times `$ 10<sup>-7</sup> and $`\chi ^2`$ of 4231. But with only 216 data points, this means that the $`\chi _\nu ^2`$ is $``$20 in both cases. Our new period agrees well with the photometric period derived recently by Voloshina, Lyutyi & Tarasov (1997) and a solution from HeI $`\lambda `$6678 measurements by Sowers et al (1998). Our values for $`\gamma `$ and K are within 2$`\sigma `$ of previous results (NWY and GB).   

The variable-period model produces a slightly better fit to the data, but, again, with a $`\chi ^2`$ difference of less than 1.5%, the improvement in quality of fit between that and the fixed-period model is not statistically significant. We conclude that there is no compelling evidence for period variation in Cygnus X-1. We have therefore adopted the orbital elements of the best-fit constant period model in all further discussion.   

We computed O–C residuals for all the helium-line velocities and performed period searches on these. No significant periodicity was found in these residuals. In particular, we find no evidence of the suggested periodicities at 39d or 78d (Kemp, Herman, & Barbour 1978), or 4.5 years (Wilson and Fox 1981) or of the 294 day X-ray modulation (Priedhorsky, Terrell, & Holt 1983) or the 91d photometric period (Walker and Quintanilla 1978). The absence of a 294 d signal was also noted by Gies & Bolton (1984).   

## 4 Discussion

Our best-fit period for the combined datasets, 5.5998 $`\pm `$ 0.0001 d is identical to well within the standard errors with the GB period of 5.59974 $`\pm `$ 0.00008 d, and our T<sub>0</sub>, or epoch of inferior conjunction, agrees with the prediction of the GB ephemeris to within 0.014 phase, or 0.078 d. Since the uncertainty in the periods yields an uncertainty in the GB prediction of T<sub>0</sub> of $`\pm `$ 0.12 days at this epoch, this is excellent agreement. Thus our results confirm the ephemeris of GB and refute suggestions of period variation by NWY.   

Why, then, did NWY’s result show such a discrepancy and suggest a change in period? They reported that the period determined from their data alone was 5.60172 $`\pm `$ 0.00003 d, differing by 20 $`\sigma `$ from the period determined from all previous observations. We suggest several sources for this error.

First of all, we performed our own period searches on the data published by NWY; the best fit we find is at 5.6002 d, or only 2 $`\sigma `$ greater than the GB value, so we may be looking at much ado about a typographical error. In addition, there is excellent phase agreement between the NWY data and the GB ephemeris, again suggesting no real period discrepancy.   

Finally, it seems that the error estimates for both the NWY period and their period determined from historical data are too small by as much as an order of magnitude. This is primarily because the formal error calculations used to determine these error estimates are based on the assumption of a very good fit of the model to the data, i.e. that $`\chi _\nu ^2`$1. This condition is not met by any of the data sets and fits used in this work, presumably because the estimates of the uncertainties in the velocities were too small, although perhaps because of variability in the source. We recalculated the fits, trying larger estimates for the velocity uncertainties until the reduced chi-square criterion was satisfied. We found that the period uncertainties then were about an order of magnitude larger than reported by the original authors. For example, treating NWY’s data in this way yields a best-fit period of 5.6002 $`\pm `$ 0.0003 d, and the discrepancy with GB disappears.   

An additional source of underestimated error uncertainty is the inclusion by both NWY and Bolton (1975) of two velocities obtained by Seyfert and Popper (1941), ostensibly to improve the precision of the period determinations. For example, Bolton (1975) uses these points to decrease his uncertainty estimate by a factor of 10; NWY do not discuss the effect of including these points on their uncertainty. Popper (1996, personal communication) suggests that these velocities are “very weak reeds on which to hang significant conclusions”, the velocities having considerable uncertainties and being based on averaging velocities from lines of several different species. Because of the problems with hydrogen-line velocities discussed above, inclusion of these points is therefore unlikely to improve the period uncertainty.   

GB also introduced “velocity corrections” to many of their velocities, shifting spectra to fit a mean interstellar K-line velocity due to instabilities in their spectrograph. It is not unlikely that this introduced uncertainties larger than the mean errors they report.   

## 5 Conclusions

We have carried out a programme of radial velocity determinations for the black hole binary Cyg X-1, and have combined these data with data used previously to determine the period, to provide a new orbital ephemeris for the source. A critical consideration of the errors associated with previous work has shown these to be underestimated. Based on this, our main conclusion is that there is no evidence for a change in orbital period as suggested by Ninkov et al. Finally, our new ephemeris allows the orbital phase calculation for Cyg X-1 with an error that is much reduced compared with the error that would be attached to extrapolating the Gies and Bolton ephemeris with its quoted accuracy to the present epoch.   

## Acknowledgements

We are very grateful to the La Palma support astronomers, and particularly Don Pollacco, who operate the Service Programme, and to those astronomers who allowed the Cyg X-1 observations to be taken during their own time. The Isaac Newton Group of telescopes is operated on the island of La Palma by the Royal Greenwich Observatory in the Spanish Observatorio del Roque de Los Muchachos of the Instituto de Astrofísica de Canarias. JLS wishes gratefully to acknowledge the following sources of partial support for his participation in this work: the Margaret Cullinan Wray Charitable Lead Annuity Trust (through the AAS Small Grants Program); the NASA JOVE program; a Theodore Dunham, Jr., Grant of the Fund for Astrophysical Research; the International Astronomical Union Exchange of Astronomers Program; and the Department of Astrophysics of the University of Oxford.

Figure 1: Typical rectified spectra of HDE226868 (top) and the radial velocity standard 19 Cep (bottom) obtained with the La Palma 2.5m INT. The spectral types are very similar, but note the HeII $`\lambda `$4686 emission in HDE226868.

Figure 2: Best fitting fixed period, circular orbit to all helium radial velocity data (ours plus all previously published material, see text). The lower panel shows the residuals to the fit.