A&A 396, 877-884 (2002)
DOI: 10.1051/0004-6361:20021273
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
Abstract
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
3.8
.
Vanbeveren et al. (1998) proposed that the progenitor
of Cygnus X-3 is a 50
+ 10
system with P0 = 6 d;
after spiral-in
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 ![]() |
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.7![]() |
Cygnus X-3 | 1.6 ![]() |
0.6 ![]() |
0.079 ![]() |
WR 78 | 1.4 ![]() |
1.4 ![]() |
0.114 ![]() |
WR 134 | 0.2 ![]() |
0.2 ![]() |
0.081 ![]() |
WR 136 | 1.0 ![]() |
1.0 ![]() |
0.076 ![]() |
WR 147 | 1.6 ![]() |
1.0 ![]() |
0.085 ![]() |
a Dereddening with Draine (1989) law.
b Dereddening with Lutz et al. (1996) law; fit between 2.4-12 ![]() c Fit between 2.4-6.5 ![]() d Fit between 2.4-12 ![]() e Fit between 2.4-8 ![]() |
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
10 kpc.
Star | Type | Binarity | Distance | AVa | Reference | TDTNUM |
WR 78 | WN7ha WNL | No | 2.0 kpc | 1.48-1.87 | Crowther et al. (1995a) | 45800705 |
WR 134 | WN6 | possible | ![]() |
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 ![]() |
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 ![]() |
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 ![]() |
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).
Acknowledgements
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.