A first detailed study of the colliding wind WR+O binary WR 30a

We present a detailed, extensive investigation of the photometric and spectroscopic behaviour of WR 30a. This star is definitely a binary system with a period around 4.6 d. We propose the value         . The identification of the components as WO4+O5((f)) indicates a massive evolved binary system; the O5 component is a main-sequence or, more likely, a giant star. The radial velocities of the O star yield a circular orbit with an amplitude         and a mass function of 0.013     . The spectrum of WR 30a exhibits strong profile variations of the broad emission lines that are phase-locked with the orbital period. We report the detection of the orbital motion of the WO component with     , but this should be confirmed by further observations. If correct, it implies a mass ratio     . The star exhibits sinusoidal light variations of amplitude 0.024 mag peak-to-peak with the minimum of light occurring slightly after the conjunction with the O star in front. On the basis of the phase-locked profile variations of the C  iv λ 4658 blend in the spectrum of the WO, we conclude that a wind–wind collision phenomenon is present in the system. We discuss some possibilities for the geometry of the interaction region.     

The massive Wolf–Rayet binary SMC WR7

We present a study of optical spectra of the Wolf–Rayet star AzV 336a (=SMC WR7) in the Small Magellanic Cloud. Our study is based on data obtained at several Observatories between 1988 and 2001. We find SMC WR7 to be a double-lined WN+O6 spectroscopic binary with an orbital period of 19.56 d. The radial velocities of the He absorption lines of the O6 component and the strong He  ii emission at λ 4686 Å of the WN component describe anti-phased orbital motions. However, they show a small phase shift of ∼1 d. We discuss possible explanations for this phase shift. The amplitude of the radial velocity variations of He  ii emission is twice that of the absorption lines. The binary components have fairly high minimum masses, ∼18 and 34 M_⊙ for the WN and O6 components, respectively.     

XMM–Newton X-ray observations of the Wolf–Rayet binary system WR 147

We present results of an ≈20-ks X-ray observation of the Wolf–Rayet (WR) binary system WR 147 obtained with XMM–Newton . Previous studies have shown that this system consists of a nitrogen-type WN8 star plus an OB companion whose winds are interacting to produce a colliding wind shock. X-ray spectra from the pn and MOS detectors confirm the high extinction reported from infrared studies and reveal hot plasma including the first detection of the Fe Kα line complex at 6.67 keV. Spectral fits with a constant-temperature plane-parallel shock model give a shock temperature   kT _shock= 2.7  keV (   T _shock≈ 31  MK), close to but slightly hotter than the maximum temperature predicted for a colliding wind shock. Optically thin plasma models suggest even higher temperatures, which are not yet ruled out. The X-ray spectra are harder than can be accounted for using 2D numerical colliding wind shock models based on nominal mass-loss parameters. Possible explanations include: (i) underestimates of the terminal wind speeds or wind abundances, (ii) overly simplistic colliding wind models or (iii) the presence of other X-ray emission mechanisms besides colliding wind shocks. Further improvement of the numerical models to include potentially important physics such as non-equilibrium ionization will be needed to rigorously test the colliding wind interpretation.     

WR 146 – observing the OB-type companion

We present new radio and optical observations of the colliding-wind system WR 146 aimed at understanding the nature of the companion to the Wolf–Rayet (WR) star and the collision of their winds. The radio observations reveal emission from three components: the WR stellar wind, the non-thermal wind–wind interaction region and, for the first time, the stellar wind of the OB companion. This provides the unique possibility of determining the mass-loss rate and terminal wind velocity ratios of the two winds, independent of distance. Respectively, these ratios are 0.20±0.06 and 0.56±0.17 for the OB-companion star relative to the WR star. A new optical spectrum indicates that the system is more luminous than had been believed previously. We deduce that the 'companion' cannot be a single, low-luminosity O8 star as suggested previously, but is either a high-luminosity O8 star, or possibly an O8+WC binary system.     

Spectropolarimetric observations of WR 137

We present high signal-to-noise spectropolarimetry of the WC7+OB star WR 137 (HD 192641) covering the spectral range from 3600Å to 1.1µm. The polarization spectrum shows unambiguous evidence for an intrinsic polarization vector.     

WR 7a: a V Sagittae or a qWR star?

The star WR 7a, also known as SPH 2, has a spectrum that resembles that of V Sagittae stars although no O  vi emission has been reported. The Temporal Variance Spectrum – TVS – analysis of our data shows weak but strongly variable emission of O  vi lines which is below the noise level in the intensity spectrum.
Contrary to what is seen in V Sagittae stars, optical photometric monitoring shows very little, if any, flickering. We found evidence of periodic variability. The most likely photometric period is   P _phot= 0.227(±14) d  , while radial velocities suggest a period of   P _spec= 0.204(±13) d  . One-day aliases of these periods can not be ruled out. We call attention to similarities with HD 45166 and DI Cru (= WR 46), where multiple periods are present. They may be associated to the binary motion or to non-radial oscillations.
In contrast to a previous conclusion by Pereira et al., we show that WR 7a contains hydrogen. The spectrum of the primary star seems to be detectable as the N  v 4604 Å  absorption line is visible. If so, it means that the wind is optically thin in the continuum and that it is likely to be a helium main sequence star.
Given the similarity to HD 45166, we suggests that WR 7a may be a qWR – quasi Wolf–Rayet – star. Its classification is WN4h/CE in the Smith, Shara & Moffat three-dimensional classification system.     

Colliding stellar wind models with non-equilibrium ionization: X-rays from WR 147

The effects of non-equilibrium ionization are explicitly taken into account in a numerical model which describes colliding stellar winds (CSW) in massive binary systems. This new model is used to analyse the most recent X-ray spectra of the WR+OB binary system WR 147. The basic result is that it can adequately reproduce the observed X-ray emission (spectral shape, observed flux) but some adjustment in the stellar wind parameters is required. Namely (i) the stellar wind velocities must be higher by a factor of 1.4–1.6 and (ii) the mass loss must be reduced by a factor of ∼2. The reduction factor for the mass loss is well within the uncertainties for this parameter in massive stars, but given the fact that the orbital parameters (e.g. inclination angle and eccentricity) are not well constrained for WR 147, even smaller corrections to the mass loss might be sufficient. Only CSW models with non-equilibrium ionization and equal (or nearly equal) electron and ion post-shock temperature are successful. Therefore, the analysis of the X-ray spectra of WR 147 provides evidence that the CSW shocks in this object must be collisionless .     

Massive stars

We describe the present state of massive star research seen from the viewpoint of stellar evolution, with special emphasis on close binaries. Statistics of massive close binaries are reasonably complete for the Solar neighbourhood. We defend the thesis that within our knowledge, many scientific results where the effects of binaries are not included, have an academic value, but may be far from reality. In chapter I, we summarize general observations of massive stars where we focus on the HR diagram, stellar wind mass loss rates, the stellar surface chemistry, rotation, circumstellar environments, supernovae. Close binaries can not be studied separately from single stars and vice versa. First, the evolution of single stars is discussed (chapter I). We refer to new calculations with updated stellar wind mass loss rate formalisms and conclusions are proposed resulting from a comparison with representative observations. Massive binaries are considered in chapter II. Basic processes are briefly described, i.e. the Roche lobe overflow and mass transfer, the common envelope process, the spiral-in process in binaries with extreme mass ratio, the effects of mass accretion and the merging process, the implications of the (asymmetric) supernova explosion of one of the components on the orbital parameters of the binary. Evolutionary computations of interacting close binaries are discussed and general conclusions are drawn. The enormous amount of observational data of massive binaries is summarized. We separately consider the non-evolved and evolved systems. The latter class includes the semi-detached and contact binaries, the WR binaries, the X-ray binaries, the runaways, the single and binary pulsars. A general comparison between theoretical evolution and observations is combined with a discussion of specially interesting binaries: the evolved binaries HD 163181, HD 12323, HD 14633, HD 193516, HD 25638, HD 209481, Per and Sgr; the WR+OB binary V444 Cyg; the high mass X-ray binaries Vela X-1, Wray 977, Cyg X-1; the low mass X-ray binaries Her X-1 and those with a black hole candidate; the runaway Pup, the WR+compact companion candidates Cyg X-3, HD 50896 and HD 197406. We finally propose an overall evolutionary model of massive close binaries as a function of primary mass, mass ratio and orbital period. Chapter III deals with massive star population synthesis with a realistic population of binaries. We discuss the massive close binary frequency, mass ratio and period distribution, the observations that allow to constrain possible asymmetries during the supernova explosion of a massive star. We focuss on the comparison between observed star numbers (as a function of metallicity) and theoretically predicted numbers of stellar populations in regions of continuous star formation and in starburst regions. Special attention is given to the O-type star/WR star/red supergiant star population, the pulsar and binary pulsar population, the supernova rates. Received 17 July 1998     

Modelling forbidden line emission profiles from colliding wind binaries

This paper presents calculations for forbidden emission-line profile shapes arising from colliding wind binaries. The main application is for systems involving a Wolf–Rayet (WR) star and an OB star companion. The WR wind is assumed to dominate the forbidden line emission. The colliding wind interaction is treated as an Archimedean spiral with an inner boundary. Under the assumptions of the model, the major findings are as follows. (i) The redistribution of the WR wind as a result of the wind collision is not flux conservative but typically produces an excess of line emission; however, this excess is modest at around the 10 per cent level. (ii) Deviations from a flat-toped profile shape for a spherical wind are greatest for viewing inclinations that are more nearly face-on to the orbital plane. At intermediate viewing inclinations, profiles display only mild deviations from a flat-toped shape. (iii) The profile shape can be used to constrain the colliding wind bow shock opening angle. (iv) Structure in the line profile tends to be suppressed in binaries of shorter periods. (v) Obtaining data for multiple forbidden lines is important since different lines probe different characteristic radial scales. Our models are discussed in relation to Infrared Space Observatory data for WR 147 and γ Vel (WR 11). The lines for WR 147 are probably not accurate enough to draw firm conclusions. For γ Vel, individual line morphologies are broadly reproducible but not simultaneously so for the claimed wind and orbital parameters. Overall, the effort demonstrates how lines that are sensitive to the large-scale wind can help to deduce binary system properties and provide new tests of numerical simulations.     

Non-thermal emission in Wolf–Rayet stars: are massive companions required?

We examine the radio spectral indices of 23 Wolf–Rayet (WR) stars to identify the nature of their radio emission. We identify nine systems as non-thermal emitters. In seven of these systems the non-thermal emission dominates the radio spectrum, while in the remaining two it is of comparable strength to the thermal, stellar wind emission, giving 'composite' spectra. Among these nine systems, seven have known spectroscopic or visual binary companions. The companions are all massive O or early B-type stars, strongly supporting a connection between the appearance of non-thermal emission in WR stars and the presence of a massive companion. In three of these binaries, the origin of non-thermal emission in a wind-collision region between the stars has been well established in earlier work. The binary systems that exhibit only thermal emission are all short‐period systems where a wind-collision zone is deep within the opaque region of the stellar wind of the WR star. To detect non-thermal emission in these systems requires optically thin lines of sight to the wind-collision region.     

Dust structure around the Wolf–Rayet star WR104 from lunar occultation observations at 2.2 μm

A rare opportunity of observing a lunar occultation of a Wolf–Rayet star (WR104) in the near-infrared K band (2.2 μm) was utilized to probe the thick dust envelope surrounding the star at a high one-dimensional angular resolution (∼2 mas). Analysis of the occultation light curve shows a dust structure departing significantly from the uniform disc profile. Our results are in good agreement with recent aperture-masking interferometry carried out at the Keck I telescope, which shows a pinwheel structure around WR104. We report additional fine structures in the dust envelope.     

A 3D dynamical model of the colliding winds in binary systems

We present a three-dimensional (3D) dynamical model of the orbital-induced curvature of the wind–wind collision region in binary star systems. Momentum balance equations are used to determine the position and shape of the contact discontinuity between the stars, while further downstream the gas is assumed to behave ballistically. An Archimedean spiral structure is formed by the motion of the stars, with clear resemblance to high-resolution images of the so-called 'pinwheel nebulae'. A key advantage of this approach over grid or smoothed particle hydrodynamic models is its significantly reduced computational cost, while it also allows the study of the structure obtained in an eccentric orbit. The model is relevant to symbiotic systems and γ-ray binaries, as well as systems with O-type and Wolf–Rayet stars.
As an example application, we simulate the X-ray emission from hypothetical O+O and WR+O star binaries, and describe a method of ray tracing through the 3D spiral structure to account for absorption by the circumstellar material in the system. Such calculations may be easily adapted to study observations at wavelengths ranging from the radio to γ-ray.     

Coordinated monitoring of the eccentric O-star binary Iota Orionis: optical spectroscopy and photometry

With the objective of investigating the windwind collision phenomenon and supporting contemporaneous X-ray observations, we have organized a large-scale, coordinated optical monitoring campaign of the massive, highly eccentric O9 III+B1 III binary Iota Orionis. Successfully separating the spectra of the components, we refine the orbital elements and confirm the rapid apsidal motion in the system. We also see strong interaction between the components during periastron passage and detect phase-locked variability in the spectrum of the secondary star. However, we find no unambiguous signs of the bow shock crashing on the surface of the secondary, despite the predictions of hydrodynamic simulations. Combining all available photometric data, we find rapid, phase-locked variations and model them numerically, thus restricting the orbital inclination to 50° i 70°.     

Repetitive structure in the stellar wind of HD 93843: a normal O-type star

We present the results from a 28-day IUE time-series campaign monitoring the stellar wind of the O5-type giant HD 93843. The principal aim was to study variability in the wind of a star with a normal projected rotation velocity. Systematic changes are identified, amidst continuous line-profile variability, in the absorption troughs of the Si  iv and N  v resonance lines. The patterns observed have characteristic time-scales of several days and are mimicked by fluctuations (of several 100 km s⁻¹) in the blue wings of the saturated C  iv P Cygni profile.   Fourier analysis provides support for the repeatability of wind structures in HD 93843 on a 7.1-d 'period'. Power at this frequency is evident only at intermediate and high velocities (i.e., above ∼0.3 of the terminal velocity). The long modulation time-scale suggests that changes in the star itself probably provide the physical source for triggering the onset of wind structure. Unfortunately the rotational, photometric, pulsational and magnetic properties of HD 93843 are too poorly constrained or known to permit a more detailed interpretation of the 7.1-d wind modulation in terms of potential inhomogeneities at the stellar surface. Nevertheless, our study demonstrates that the incidence of cyclic, possibly regular, stellar-wind variability is not restricted to rapid rotators. Comparisons with other OB stars which have exhibited repetitive wind changes on 'periods' of several days suggest that the time-dependent UV properties of HD 93843 are more akin to those of the O4-type supergiant ζ Puppis.     

The Galactic hybrid Wolf–Rayet WN 7o/CE + O7V((f)) binary system WR 145

A spectroscopic study of the binary Wolf–Rayet (WR)+O system WR 145 is performed, in order to determine the radial velocity orbits of the individual stars, the angle of orbital inclination and the stellar masses. The emission and absorption components are separated from the original spectra, allowing us to confirm the spectral classification WN 7o/CE of the hybrid WR component and to derive a spectral classification O7V((f)) for the O star. A study of the wind-collision properties is performed. Fitting the radial velocity and full width at half-maximum of the excess emission with Lührs' model results in an inclination angle of   i = 63°  , leading to estimates of the stellar masses:   M _WR= 18 M_⊙  and   M _O= 31 M_⊙  . Both of these masses are compatible with those of other stars of similar types.     

A first orbital solution for the very massive 30 Dor main-sequence WN6h+O binary R145

We report the results of a spectroscopic and polarimetric study of the massive, hydrogen-rich WN6h stars R144 (HD 38282 = BAT99-118 = Brey 89) and R145 (HDE 269928 = BAT99-119 = Brey 90) in the Large Magellanic Cloud. Both stars have been suspected to be binaries by previous studies (R144: Schnurr et al.; R145: Moffat). We have combined radial-velocity (RV) data from these two studies with previously unpublished polarimetric data. For R145, we were able to establish, for the first time, an orbital period of 158.8 d, along with the full set of orbital parameters, including the inclination angle i , which was found to be   i = 38°± 9°  . By applying a modified version of the shift-and-add method developed by Demers et al., we were able to isolate the spectral signature of the very faint line companion star. With the RV amplitudes of both components in R145, we were thus able to estimate their absolute masses. We find minimum masses   M _WRsin³ i = 116 ± 33 M_⊙  and   M _Osin³ i = 48 ± 20 M_⊙  for the WR and the O component, respectively. Thus, if the low-inclination angle were correct, resulting absolute masses of the components would be at least 300 and  125 M_⊙  , respectively. However, such high masses are not supported by brightness considerations when R145 is compared to systems with known very high masses such as NGC 3603-A1 or WR20a. An inclination angle close to  90°  would remedy the situation, but is excluded by the currently available data. More and better data are thus required to firmly establish the nature of this puzzling, yet potentially very massive and important system. As to R144, however, the combined data sets are not sufficient to find any periodicity.     

The nearest star of spectral type O3: a component of the multiple system HD 150136

From radial velocities determined in high signal-to-noise digital spectra, we report the discovery that the brightest component of the binary system HD 150136 is of spectral type O3. We also present the first double-lined orbital solution for this binary. Our radial velocities confirm the previously published spectroscopic orbital period of 2.6 d. He  ii absorptions appear double at quadratures, but single lines of N  v and N  iv visible in our spectra define a radial velocity orbit of higher semi-amplitude for the primary component than do the He  ii lines. From our orbital analysis, we obtain minimum masses for the binary components of 27 and  18 M_⊙  . The neutral He absorptions apparently do not follow the orbital motion of any of the binary components, thus they most probably arise in a third star in the system.     

The massive double-lined O-type binary HD 165052

We present a new optical spectroscopic study of the O-type binary HD 165052 based on high- and intermediate-resolution CCD observations. We re-investigated the spectral classification of the binary components, obtaining spectral types of O6.5 V and O7.5 V for the primary and secondary, respectively, finding that both stars display weak C  iii λ 5696 emission in their spectra. We also determined a radial-velocity orbit for HD 165052 with a period of  2.95510±0.00001 d  , and semi-amplitudes of 94.8 and  104.7±0.5 km s^-1  , resulting in a mass ratio   Q =0.9  . From a comparison with previous radial-velocity determinations, we found evidence of apsidal motion in the system. Several signatures of wind–wind collision, such as phase-locked variability of the X-ray flux and the Struve–Sahade effect, are also considered. It was also found that the reddening in the region should be normal, in contrast with previous determinations.     

Shock breakout in Type Ibc supernovae and application to GRB 060218/SN 2006aj

Recently, a soft blackbody component was observed in the early X-ray afterglow of GRB 060218, which was interpreted as shock breakout from the thick wind of the progenitor Wolf–Rayet (WR) star of the underlying Type Ic supernova 2006aj. In this paper, we present a simple model for computing the characteristic quantities (including energy, temperature and time duration) for the transient event from the shock breakout in Type Ibc supernovae produced by the core-collapse of WR stars surrounded by dense winds. In contrast to the case of a star without a strong wind, the shock breakout occurs in the wind region rather than inside the star, caused by the large optical depth in the wind. We find that, for the case of a WR star with a dense wind, the total energy of the radiation generated by the supernova shock breakout is larger than that in the case of the same star without a wind by a factor of >10. The temperature can be either hotter or colder, depending on the wind parameters. The time duration is larger caused by the increase in the effective radius of the star due to the presence of a thick wind. Then, we apply the model to GRB 060218/SN 2006aj. We show that, to explain both the temperature and the total energy of the blackbody component observed in GRB 060218 by the shock breakout, the progenitor WR star has to have an unrealistically large core radius (the radius at optical depth of 20), larger than 100 R_⊙. In spite of this disappointing result, our model is expected to have important applications to the observations on Type Ibc supernovae in which the detection of shock breakout will provide important clues to the progenitors of Type Ibc supernovae.     

KPD 0422 + 5421: a new short-period subdwarf B/white dwarf binary

The subdwarf B (sdB) star KPD 0422 + 5421 was discovered to be a single-lined spectroscopic binary with a period of P  = 0.090 1795 ± (3 × 10⁻⁷) d (2 h 10 min). The U B light curves display an ellipsoidal modulation with amplitudes of ≈ 0.02 mag. The sdB star contributes nearly all of the observed flux. This and the absence of any reflection effect suggest that the unseen companion star is small (i.e. R _comp ≈ 0.01 R_⊙) and therefore degenerate. We modelled the U B light curves and derived i  = 78.05° ± 0.50° and a mass ratio of q  =  M _comp/ M _sdB = 0.87 ± 0.15. The sdB star fills 69 per cent of its Roche lobe. These quantities may be combined with the mass function of the companion [ f ( M ) = 0.126 ± 0.028 M_⊙] to derive M _sdB = 0.72 ± 0.26 M_⊙ and M _comp = 0.62 ±  0.18 M_⊙. We used model spectra to derive the effective temperature, surface gravity and helium abundance of the sdB star. We found T _eff = 25 000 ± 1500 K, log g  = 5.4 ± 0.1 and [He/H] = −1.0. With a period of 2 h 10 min, KPD 0422 + 5421 has one of the shortest known orbital periods of a detached binary. This system is also one of only a few known binaries that contain a subdwarf B star and a white dwarf. Thus KPD 0422 + 5421 represents a relatively unobserved, and short-lived, stage of binary star evolution.