A&A 396, 877-884 (2002)
L. Koch-Miramond1 - P. Ábrahám2,3 - Y. Fuchs1,4 - J.-M. Bonnet-Bidaud1 - A. Claret1
1 - DAPNIA/Service d'Astrophysique, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France
2 - Konkoly Observatory, PO Box 67, 1525 Budapest, Hungary
3 - Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany
4 - Université Paris VII, France
Received 5 June 2002 / Accepted 28 June 2002
We present mid-infrared spectrophotometric results obtained with the ISO on the peculiar X-ray binary Cygnus X-3 in quiescence, at orbital phases 0.83 to 1.04. The 2.4-12 m continuum radiation observed with ISOPHOT-S can be explained by thermal free-free emission in an expanding wind with, above 6.5 m, a possible additional black-body component with temperature K and radius at 10 kpc, likely due to thermal emission by circumstellar dust. The observed brightness and continuum spectrum closely match that of the Wolf-Rayet star WR 147, a WN8+B0.5 binary system, when rescaled at the same 10 kpc distance as Cygnus X-3. A rough mass loss estimate assuming a WN wind gives yr-1. A line at 4.3 m with a more than 4.3 detection level, and with a dereddened flux of 126 mJy, is interpreted as the expected He I 3p-3s line at 4.295 m, a prominent line in the WR 147 spectrum. These results are consistent with a Wolf-Rayet-like companion to the compact object in Cyg X-3 of WN8 type, a later type than suggested by earlier works.
Key words: stars: binaries: close - stars: individual: Cyg X-3 - stars: Wolf-Rayet - stars: mass-loss - infrared: stars
Cygnus X-3 has been known as a binary system since its discovery by
Giacconi et al. (1967), but there is still debate about the masses
of the two stars and the morphology of the system (for a review see Bonnet-Bidaud & Chardin 1988).
The distance of
the object is 8-12.5 kpc with an absorption on the line of sight
mag (van Kerkwijk et al. 1996). The flux modulation
at a period of 4.8 hours, first discovered in X-rays
(Parsignault 1972), then at near infrared wavelengths
(Becklin et al. 1973), and observed simultaneously at X-ray
and near-IR wavelengths by Mason et al. (1986), is believed to be the
orbital period of the binary system.
Following infrared spectroscopic measurements (van Kerkwijk et al. 1992),
where WR-like features have been detected in I and K band spectra,
the nature of the mass-donating star is suggested to be a Wolf-Rayet-like star,
but an unambiguous classification, similar
to the other WR stars, is still lacking. Mitra (1996) and Vanbeveren et al. (1998)
pointed out that it is not possible to find
a model that meets all the observed properties of Cygnus X-3 where
the companion star is a normal Population I Wolf-Rayet star with a
spherically symmetric stellar wind. In the evolution model originally
proposed by van den Heuvel & de Loore (1973) a final period of the
order of 4.8 h may result from a system with initial masses
M10 = 15 ,
M20 = 1 ,
P0 = 5d,
the final system being a neutron star
accreting at a limited rate of 10-7 yr-1,
from the wind of a core He burning star of about
Vanbeveren et al. (1998) proposed that the progenitor
of Cygnus X-3 is a 50 + 10
system with P0 = 6 d;
of the black hole into the envelope of the companion, the hydrogen reach
layers are removed, and a 2-2.5
Wolf-Rayet like star remains
with P = 0.2 d. A system containing a black hole and an He core burning
star is also favored by Ergma & van den Heuvel (1998).
In addition, Cygnus X-3 undergoes giant radio bursts and
there is evidence of jet-like structures moving away from Cygnus X-3 at
0.3-0.9 c (Mioduszewski et al. 1998; 2001; Martí et al. 2001).
The main objective of the Infrared Space Observatory (ISO) spectrophotometric measurements in the 2.4-12 m range was to constrain further the nature of the companion star to the compact object: the expected strong He lines as well as the metallic lines in different ionization states are important clues, together with the spectral shape of the continuum in a wavelength range as large as possible. An additional motivation for the imaging photometry with ISOCAM was to provide spatial resolution to a possible extended emission feature as a remnant of the expected high mass loss from the system. The paper is laid out as follows. In Sect. 2 observational aspects are reviewed. Section 3 summarizes the results on the continuum and line emissions from Cygnus X-3 and four Wolf-Rayet stars of WN 6, 7 and 8 types, and reviews the constraints set by the present observations on the wind and on the nature of the companion to the compact object in Cygnux X-3. Finally, Sect. 4 summarizes the conclusions of this paper.
We observed Cygnus X-3 with the Infrared Space Observatory (ISO, see Kessler et al. 1996) on April 7, 1996 corresponding to JD 2 450 180.8033 to 2 450 180.8519. The subsequent observing modes were: ISOCAM imaging photometry at 11.5 m (LW10 filter, bandwith 8 to 15 m), ISOPHOT-S spectrophotometry in the range 2.4-12 m, for 4096 s, covering the orbital phases 0.83 to 1.04 (according to the parabolic ephemeris of Kitamoto et al. 1992); ISOPHOT multi-filter photometry at central wavelengths 3.6, 10, 25 and 60 m. Observing modes and observation times are summarized in Table 1. Preliminary results were presented in Koch-Miramond et al. (2001).
The ISOPHOT-S data were reduced in three steps. We first used the Phot Interactive Analysis (PIA, Gabriel et al. 1997) software (version 8.2) to filter out cosmic glitches in the raw data and to determine signals by performing linear fits to the integration ramps. After a second deglitching step, performed on the signals, a dark current value appropriate to the satellite orbital position of the individual signal was subtracted. Finally we averaged all non-discarded (typically 3) signals in order to derive a signal per chopper step. Due to detector transient effects, at the beginning of the observation the derived signals were systematically lower than those in the consolidated part of the measurement. We then discarded the first 800 s (3 OFF-ON transitions), and determined an average [ON-OFF] signal for the whole measurement by applying a 1-dimensional Fast Fourier Transformation algorithm (for the application of FFT methods for ISOPHOT data reduction see Haas et al. 2000). The [ON-OFF] difference signals were finally calibrated by applying a signal-dependent spectral response function dedicated to chopped ISOPHOT-S observations (Acosta-Pulido & Ábrahám 2001), also implemented in PIA.
In order to verify our data reduction scheme (which is not completely standard due to the application of the FFT algorithm) and to estimate the level of calibration uncertainties, we reduced HD 184400, an ISOPHOT standard star observed in a similar way as Cygnus X-3. The results were very consistent with the model prediction of the star, and we estimate that the systematic uncertainty of our calibration is less than 10.
|Figure 1: Observed spectrum of Cygnus X-3 in the 2.4-12 m range, obtained with the FFT method.|
|Open with DEXTER|
The observed spectral energy distribution is shown in Fig. 1. The observed (not dereddened) continuum flux in the range 2.4-7 m is mJy in good agreement with that observed by Ogley et al. (2001) with ISOCAM on the same day (the dereddened fluxes are shown in Fig. 2) ; the observed flux decreases to about mJy around 9 m.
An unresolved line is observed at about 4.3 m peaking at 57 10 mJy. The linewidth is 0.04 m, consistent with the instrumental response and corresponding to 2500 km s-1. Note that the measured line flux might be underestimated because the ISOPHOT-S pixels are separated by small gaps, and a narrow line might falls into a gap.
The high resolution configuration of the ISOCAM camera has been used to constrain the spatial extension of the infrared source. The measured FWHM for the source is arcsec (mean value of the four individual images composing the final raster map). This can be compared to the ISOCAM catalogued point spread functions at these energy and configuration which show a FWHM mean value of arcsec, including the effects of the satellite jitter and of the pixel sampling. The slightly larger value for Cygnus X-3, though only marginally significant, might therefore indicate an extended source. The deconvolved extension would be arcsec which at a distance of 10 kpc corresponds to a linear extension of cm. The extended infrared source may be the result of the heating of the surrounding medium by the radio jets whose existence have been now clearly demonstrated both at arcsec (Martí et al. 2000, 2001) and sub-arcsec (Mioduszewski et al. 2001) scales, but it clearly deserves confirmation.
|Figure 2: Dereddened spectrum of Cygnus X-3 using either Draine (1989) law (left) or Lutz et al. (1996) law (right). The asterisks (*) in the lower panel represent the Ogley et al. (2001) ISOCAM results, after dereddening by Draine (1989) (left) and Lutz (right).|
|Open with DEXTER|
|Figure 3: Best fitting of the two dereddened spectra of Cygnus X-3 shown in Fig. 2.|
|Open with DEXTER|
Using the ISO archive data we have analysed the SWS
spectra of four Wolf-Rayet stars: WR 147 (WN8+B0.5), WR 136
(WN6b), WR 134 (WN6) and WR 78 (WN7) whose main
characteristics are given in Table 3. We compare them to the
Cygnus X-3 spectrum, after smoothing the SWS spectra to the
resolution of the ISOPHOT-S instrument (using an IDL routine
of B. Schulz dowloaded from the Home Page of the ISO Data
Centre at Vilspa). The observed WR spectra are shown in Fig.
4 on top of the observed Cygnus X-3 spectrum; the
identification of the emission lines is from Morris et al. (2000).
|Object||PL slopea||PL slopeb||Flux at 4.7m (Jy)|
|Cygnus X-3||1.6 0.2c||0.6 0.30.4||0.079 0.011|
|WR 78||1.4 0.2d||1.4 0.2||0.114 0.008|
|WR 134||0.2 0.2d||0.2 0.2||0.081 0.006|
|WR 136||1.0 0.2e||1.0 0.2||0.076 0.006|
|WR 147||1.6 0.1c||1.0 0.1||0.085 0.005|
a Dereddening with Draine (1989) law.
b Dereddening with Lutz et al. (1996) law; fit between 2.4-12 m.
c Fit between 2.4-6.5 m.
d Fit between 2.4-12 m.
e Fit between 2.4-8 m.
The dereddened spectra of the Wolf-Rayet stars, using
either the Draine (1989) law or the Lutz et al. (1996) law with the AV shown in Table 3, have been fitted with power law slopes
given in Table 2. Wolf-Rayet stars emit free-free continuum
radiation from their extended ionized stellar wind envelopes and
the different slopes reflect different conditions in the wind
(Williams et al. 1997). It is noticeable that the mean continuum flux
density of Cygnus X-3 is the same (within a factor 1.5 at
4.7 m as seen in Table 2) as that of the four WR stars when
their flux density is rescaled to a Cygnus X-3 distance of
|WR 78||WN7ha WNL||No||2.0 kpc||1.48-1.87||Crowther et al. (1995a)||45800705|
|WR 134||WN6||possible||2.1 kpc||1.22-1.99||Morel et al. (1999)||17601108|
|WR 136||WN6b(h) WNE-sb||possible||1.8 kpc||1.35-2.25||Stevens & Howarth (1999)||38102211|
|WR 147||WN8(h) WNL||B0.5V at 0.554''||630 pc||11.2||Morris et al. (1999, 2000)||33800415|
a From van der Hucht (2001) except for WR 147.
b From Crowther et al. (1995b).
|Figure 4: Observed spectra of Cygnus X-3 and four Wolf-Rayet stars (not dereddened). An arbitrary vertical offset has been added to the WR spectra for clarity. The identification of the emission lines is from Morris et al. (2000).|
|Open with DEXTER|
The comparison between the Cygnus X-3 spectrum and
that of the Wolf-Rayet WR 147 at 10 kpc is shown in Fig. 5
after dereddening with the Draine (1989) law (left) and with
the Lutz et al. (1996) law (right). The WR 147 spectrum appears
as the closest WR one to the Cygnus X-3 spectrum, with almost
the same mean flux density at 10 kpc and the same power law
slope (whithin the statistical errors).
|Figure 5: Comparison between the spectral energy distribution of Cygnus X-3 and the one of WR 147 when rescaled at 10 kpc: both dereddened with Draine (1989) law (left) and the Lutz et al. (1996) law (right).|
|Open with DEXTER|
We note that WR 147 is known as a colliding-wind binary that has been spatially resolved (Williams et al. 1997; Skinner et al. 1999), with a separation on the sky large enough for the wind-wind collision zone between the stars to be resolved at near-infrared and radio (Williams et al. 1997), and X-ray energies (Pittard et al. 2002). The spectral energy distribution of WR 147 in the 0.5 m to 2 mm wavelength range (including all components) shown by Williams et al. (1997) is dominated by the free-free emission from the stellar wind of the WN8 star; in the 2 to 10 m range these authors find = 1.0, in good agreement with our ISOPHOT-S measurement (when dereddened with the Lutz et al. 1996 law); and in the mid-infrared to radio range they find = 0.66.
The measured 4.3 m line flux above the continuum in the Cygnus X-3 spectrum is 58 11 mJy (dereddening with the Draine 1989 law), and 126 25 mJy (dereddening with the Lutz et al. 1996 law), using respectively = 1.6 and = 0.6, the best fitted continuum slopes as given in Table 3, both detections being at more than 4.3 . This line is interpreted as the HeI (3p-3s) line at 4.295 m, a prominent line in the WR 147 (WN8+B0.5) spectrum as seen in Morris et al. (2000) and in Fig. 4. Again WR 147 appears as the closest WR to Cygnus X-3 as being the only WR in our sample with a HeI emission line at 4.3 m, the only line clearly seen in our Cygnus X-3 data. The other expected He lines at 2.62, 3.73, 4.05, 7.46 and 10.5 m are not detected, probably due to the faintness of the object, at the limits of the instrument's sensitivity. We note (Fig. 4) that the second highest peak in the SS-part of the Cyg X-3 spectrum is at 3.73 m, and there are also local maxima at 4.05 and 10.5 m. These expected lines are all blended with H lines and the absence of H observed by van Kerkwijk et al. (1996) in the I and K band spectra of Cygnus X-3 could explain the weakness or absence of these lines in our data.
We note that the Br+HeI-II line at 4.05 m is not detected in Cygnus X-3 in quiescence, but is present in the four Wolf-Rayet stars. Strong HeI and HeII lines have been previously observed in the K-range in Cygnus X-3 during quiescence (van Kerkwijk et al. 1992; 1996, Fender et al. 1999). These lines have been interpreted (van Kerkwijk et al. 1996; Cherepashchuk & Moffat 1994), as emission from the wind of a massive companion star to the compact object, and Fender et al. (1999) suggest that the best candidate is probably an early WN Wolf-Rayet star. We note that the close match we have found between the mid-infrared luminosity and the spectral energy distribution, the HeI emission line in Cygnus X-3 in quiescence, and that of the WR 147, is consistent with a Wolf-Rayet like companion of WN8 type to the compact object in Cygnus X-3, a later type than suggested by earlier works (van Kerkwijk et al. 1996; Fender et al. 1999; Hanson et al. 2000).
|Figure 6: Flux density in the 2.4-12 m range versus Cygnus X-3 orbital phase.|
|Open with DEXTER|
|Figure 7: Quasi-simultaneous observations of Cygnus X-3 in the radio, infrared, soft and hard X-rays on April 7, 1996, averaged over the orbital phase; the ISOPHOT-S spectrum is dereddened with the Lutz et al. (1996) law and rebinned to the resolution of 0.3 m; note that the GBI was dormant after 1996 April 1 till November 1996, and that the given flux densities (in dash lines) are mean values during the quiescent period March 17 to April 1, 1996.|
|Open with DEXTER|
In the X-ray range, at the same epoch, the XTE/All Sky Monitor count rate was 7.5 count s-1 corresponding to a mean flux of 1 mJy from 2 to 12 keV (see XTE archive and Levine et al. 1996), and the BATSE instrument on board the Compton Gamma Ray Observatory observed a mean photon flux of 0.039 count s-1 corresponding to a flux density of 0.04 mJy in the 20-100 keV range. Thus the mid-infrared continuum spectrum whose shape is explained by thermal free-free emission in an expanding wind has a different origin than the non-thermal radio emission and the hard X-ray emission which are closely coupled (Mioduszewski et al. 2001; Choudhury et al. 2002).
We warmly thank the ISO project and the ISOCAM and ISOPHOT Teams in Villafranca, Saclay and Heidelberg. We express our gratitude to R. Ogley and to R. Fender for helpful comments, to G. Pooley for giving us the Ryle telescope data and to J. L. Starck for very useful discussions on data analysis. We thank the referee J. Martí for helpful comments on the manuscript. This research has made use of data from the Green Bank Interferometer, a facility of the National Science Foundation operated by the NRAO in support of NASA High Energy Astrophysics programs, of data which were generated by the CGRO BATSE Instrument Team at the Marshall Space Flight Center (MSFC) using the Earth occultation technique, and of quick-look results provided by the ASM/RXTE team. P.A. acknowledges the support of a Hungarian science grant.