Magellanic Cloud WC/WO Wolf–Rayet stars – I. Binary frequency and Roche lobe overflow formation

A nearly complete sample of 24 Magellanic Cloud WC/WO subclass Wolf–Rayet stars is studied spectroscopically and photometrically to determine its binary frequency. Theory predicts the Roche lobe overflow produced Wolf–Rayet binary frequency to be 52±14 per cent in the Large Magellanic Cloud and 100 per cent in the Small Magellanic Cloud, not counting non-Roche lobe overflow Wolf–Rayet binaries. Lower ambient metallicity ( Z ) leads to lower opacity, preventing all but the most massive (hence luminous) single stars from reaching the Wolf–Rayet stage. However, theory predicts that Roche lobe overflow even in binaries of modest mass will lead to Wolf–Rayet stars in binaries with periods below approximately 200 d, for initial periods below approximately 1000 d, independent of Z . By examining their absolute continuum magnitudes, radial velocity variations, emission-line equivalent widths and full widths at half-maximum, a WC/WO binary frequency of only 13 per cent, significantly lower than the prediction, is found in the Large Magellanic Cloud. In the unlikely event that all of the cases with a less certain binary status actually turn out to be binary, current theory and observation would agree. (The Small Magellanic Cloud contains only one WC/WO star, which happens to be a binary.) The three WC+O binaries in the Large Magellanic Cloud all have periods well below 1000 d. The large majority of WC/WO stars in such environments apparently can form without the aid of a binary companion. Current evolutionary scenarios appear to have difficulty explaining either the relatively large number of Wolf–Rayet stars in the Magellanic Clouds, or the formation of Wolf–Rayet stars in general.     

The discovery of a low-mass, pre-main-sequence stellar association around γ Velorum

We report the serendipitous discovery of a population of low-mass, pre-main-sequence (PMS) stars in the direction of the Wolf–Rayet/O-star binary system γ ²  Vel and the Vela OB2 association. We argue that γ ²  Vel and the low-mass stars are truly associated and approximately coeval, and that both are at distances between 360 and 490 pc, disagreeing at the 2 σ level with the recent Hipparcos parallax of γ ²  Vel, but consistent with older distance estimates. Our results clearly have implications for the physical parameters of the γ ²  Vel system, but also offer an exciting opportunity to investigate the influence of high-mass stars on the mass function and circumstellar disc lifetimes of their lower mass PMS siblings.     

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.     

The expansion of the giant filamentary shell, DEM 171, in the Magellanic Bridge

The giant filamentary shell, DEM 171, is found to be expanding at approximately 37 km s⁻¹. A supernova and stellar wind origin are both explored as possible causes for the expanding shell. A stellar wind origin would imply a mass-loss rate of the order of 10⁻⁵ M_⊙ yr⁻¹, indicating that it could be caused by a Wolf–Rayet star. A number of blue stars are found to lie within the shell and one is identified as a Wolf–Rayet candidate.     

Multiwavelength study of M33's giant H ii regions NGC 588 and NGC 592

We present a detailed multiwavelength photometric study of giant H  ii regions NGC 592 and NGC 588 in the nearby small spiral galaxy M33. We use data taken with the Wide Field and Planetary Camera 2 (WFPC2) on board the Hubble Space Telescope ( HST ). We detect several massive stars in both ionizing clusters. Six Wolf–Rayet (WR) stars are known to exist within those regions and we are able to constrain their physical properties by comparing their photometry to the latest grid of model atmospheres for WR stars of the nitrogen sequence (WN subclass). We estimate the age and mass of both regions by fitting our photometry to models of integrated stellar populations.     

M4–18: the planetary nebula and its WC10 central star

We present a detailed analysis of the planetary nebula M4–18 (G146.7+07.6) and its WC10-type Wolf–Rayet (WR) central star, based on high‐quality optical spectroscopy (WHT/UES, INT/IDS, WIYN/DensPak) and imaging ( HST /WFPC2). From a non-LTE model atmosphere analysis of the stellar spectrum, we derive T _eff=31 kK,     v _∞=160 km s⁻¹ and abundance number ratios of H/He<0.5, C/He=0.60 and O/He=0.10. These parameters are remarkably similar to those of He 2–113 ([WC10]). Assuming an identical stellar mass to that determined by De Marco et al. for He 2–113, we obtain a distance of 6.8 kpc to M4–18 [ E ( B−V )=0.55 mag from nebular and stellar techniques]. This implies that the planetary nebula of M4–18 has a dynamical age of ∼3100 yr, in contrast to ≥270 yr for He 2–113. This is supported by the much higher electron density of the latter. These observations may be reconciled with evolutionary predictions only if [WC]-type stars exhibit a range in stellar masses.
Photoionization modelling of M4–18 is carried out using our stellar WR flux distribution, together with blackbody and Kurucz energy distributions obtained from Zanstra analyses. We conclude that the ionizing energy distribution from the WR model provides the best consistency with the observed nebular properties, although discrepancies remain.     

Photometric variability of WC9 stars

Do some Wolf–Rayet stars owe their strong winds to something else besides radiation pressure? The answer to this question is still not entirely obvious, especially in certain Wolf–Rayet subclasses, mainly WN8 and WC9. Both of these types of Wolf–Rayet stars are thought to be highly variable, as suggested by observations, possibly due to pulsations. However, only the WN8 stars have so far been vigorously and systematically investigated for variability. We present here the results of a systematic survey during three consecutive weeks of 19 Galactic WC9 stars and one WC8 star for photometric variability in two optical bands, V and I . Of particular interest are the correlated variations in brightness and colour index in the context of carbon dust formation, which occurs frequently in WC9 and some WC8 stars. In the most variable case, WR76, we used this information to derive a typical dust grain size of  ∼ 0.1 μm  . However, most photometric variations occur at surprisingly low levels, and in fact almost half of our sample shows no significant variability at all above the instrumental level (  σ∼ 0.005– 0.01  mag).     

The circumstellar environment of Wolf–Rayet stars and gamma-ray burst afterglows

We study the evolution of the circumstellar medium of massive stars. We pay particular attention to Wolf-Rayet stars that are thought to be the progenitors of some long gamma-ray bursts (GRBs). We detail the mass-loss rates we use in our stellar evolution models and how we estimate the stellar wind speeds during different phases. With these details we simulate the interactions between the wind and the interstellar medium to predict the circumstellar environment around the stars at the time of core-collapse. We then investigate how the structure of the environment might affect the GRB afterglow. We find that when the afterglow jet encounters the free-wind/stalled-wind interface, rebrightening occurs and a bump is seen in the afterglow light curve. However, our predicted positions of this interface are too distant from the site of the GRB to reach while the afterglow remains observable. The values of the final wind density,   A _*  , from our stellar models are of the same order (≲1) as some of the values inferred from observed afterglow light curves. We do not reproduce the lowest   A _*  values below 0.5 inferred from afterglow observations. For these cases, we suggest that the progenitors could have been a WO-type Wolf–Rayet (WR) star or a very low-metallicity star. Finally, we turn our attention to the matter of stellar wind material producing absorption lines in the afterglow spectra. We discuss the observational signatures of two WR stellar types, WC and WO, in the afterglow light curve and spectra. We also indicate how it may be possible to constrain the initial mass and metallicity of a GRB progenitor by using the inferred wind density and wind velocity.     

New Wolf–Rayet central stars of planetary nebulae identified on the AAO/UKST Hα Survey

The central stars of two of the new planetary nebulae found during scans of the AAO/UKST H α Survey of the Milky Way have been found to exhibit Wolf–Rayet (WR) emission features. One (PMR 1) is an early-type star of class either [WO4] or [WC4]. The other (PMR 2) is a late [WC] star which, depending on the classification scheme used, is either intermediate in class between [WC9] and [WC10] or the sole member of the [WC10] class. Both stars exhibit unusual spectral features which may be attributed to enhanced nitrogen in their atmospheres and could be indicative of unusual stellar evolution.     

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.     

The dual dust chemistries of planetary nebulae with [WCL] central stars

The rather rare class of central stars of planetary nebulae that show very low-excitation Wolf–Rayet spectra has been a subject of great interest, particularly in the infrared, since its discovery in the late 1960s. Further peculiarities have been found with the advent of infrared spectroscopy from ISO . Notably, these objects simultaneously betray the presence of regions of carbon-rich and oxygen-rich dust chemistry. We compare and contrast complete ISO spectra between 2 and 200 μm of a sample of six [WC8] to [WC11] central stars, finding many similarities. Among this sample, one star provides strong evidence of quasi-periodic light variations, suggestive of a dust cloud orbiting in a plane from which we view the system.     

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.     

Chemical enrichment by Wolf–Rayet stars: non-solar metallicities

Following on from our recent Paper I, we present theoretical models of Wolf–Rayet (WR) stars for non-solar metallicities from   Z = 0.03  to 0.0001 by mass fraction with different mass-loss rate assumptions. We find that some single WR stars may still form even at the lowest metallicities, but whether this occurs or not depends critically on the upper cut-off point of the initial mass function used. As at solar metallicity, a population of binaries is required to fully reproduce WR star observations. For most scenarios, these binaries dominate the low-metallicity WR population but probably not the enrichment. We find comparable carbon enrichment from single WR stars to that from asymptotic giant branch stars at all metallicities for which data are available, but which of them is the dominant source of carbon depends strongly on the set of asymptotic giant branch yields adopted and the assumed initial mass function. We find an increase in carbon enrichment with increasing metallicity but a decrease in oxygen enrichment, as confirmed by observation.     

Searching for hidden Wolf–Rayet stars in the Galactic plane – 15 new Wolf–Rayet stars

We report the discovery of 15 previously unknown Wolf–Rayet (WR) stars found as part of an infrared (IR) broad-band study of candidate WR stars in the Galaxy. We have derived an empirically based selection algorithm which has selected ∼5000 WR candidate stars located within the Galactic plane drawn from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (mid-IR) and Two-Micron All-Sky Survey (near-IR) catalogues. Spectroscopic follow-up of 184 of these reveals 11 nitrogen-rich (WN) and four carbon-rich (WC) WR stars. Early WC subtypes are absent from our sample and none shows evidence for circumstellar dust emission. Of the candidates which are not WR stars, ∼120 displayed hydrogen emission-line features in their spectra. Spectral features suggest that the majority of these are in fact B supergiants/hypergiants, ∼40 of these are identified Be/B[e] candidates.
Here, we present the optical spectra for six of the newly detected WR stars, and the near-IR spectra for the remaining nine of our sample. With a WR yield rate of ∼7 per cent and a massive star detection rate of ∼65 per cent, initial results suggest that this method is one of the most successful means for locating evolved, massive stars in the Galaxy.     

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.     

Hubble Space Telescope imaging of the progenitor sites of six nearby core-collapse supernovae

The search for the progenitors of six core-collapse supernovae (CCSNe) in archival Hubble Space Telescope ( HST ) WFPC2 pre-explosion imaging is presented. These SNe are 1999an, 1999br, 1999ev, 2000ds, 2000ew and 2001B. Post-explosion imaging of the SNe, with the HST ACS/WFC, has been utilized with the technique of differential astrometry to identify the progenitor locations on the pre-explosion imaging. SNe 1999br, 1999ev, 2000ew and 2001B are recovered in late-time imaging, and estimates of the progenitor locations on the pre-explosion imaging, with subpixel accuracy, have been made. Only the progenitor of the Type II-P SN 1999ev has been recovered, on pre-explosion F555W imaging, at a 4.8σ significance level. Assuming a red supergiant progenitor, the pre-explosion observation is consistent with   M _ZAMS= 15–18 M_⊙  . The progenitors of the other five SNe were below the 3σ detection threshold of the pre-explosion observations. The detection thresholds were translated to mass limits for the progenitors by comparison with stellar evolution models. Pre-explosion observations of the peculiarly faint SN 1999br limit the mass of a red supergiant progenitor to   M _ZAMS < 12 M_⊙  . Analysis has been extended, from previous studies, to include possible detections of high- T _eff, high-mass stars by conducting synthetic photometry of model Wolf–Rayet star spectra. The mass limits for the Type II-P SNe 1999an and 1999br are consistent with previously determined mass limits for this type of SN. The detection limits for the progenitors of the Type Ibc SNe (2000ds, 2000ew and 2001B) do not permit differentiation between high-mass Wolf–Rayet progenitors or low-mass progenitors in binaries.     

Models of forbidden line emission profiles from axisymmetric stellar winds

A number of strong infrared forbidden lines have been observed in several evolved Wolf–Rayet (WR) star winds, and these are important for deriving metal abundances and testing stellar evolution models. In addition, because these optically thin lines form at large radius in the wind, their resolved profiles carry an imprint of the asymptotic structure of the wind flow. This work presents model forbidden line profile shapes formed in axisymmetric winds. It is well known that an optically thin emission line formed in a spherical wind expanding at constant velocity yields a flat-topped emission profile shape. Simulated forbidden lines are produced for a model stellar wind with an axisymmetric density distribution that treats the latitudinal ionization self-consistently and examines the influence of the ion stage on the profile shape. The resulting line profiles are symmetric about line centre. Within a given atomic species, profile shapes can vary between centrally peaked, doubly peaked, and approximately flat-topped in appearance depending on the ion stage (relative to the dominant ion) and viewing inclination. Although application to WR star winds is emphasized, the concepts are also relevant to other classes of hot stars such as luminous blue variables and Be/B[e] stars.     

Multifrequency study of the ring nebula SG 13

We investigate the morphology and kinematics of the interstellar medium in the environs of the open cluster Mrk 50, which includes the Wolf–Rayet star WR 157 and a number of early B-type stars. The analysis was performed using radio continuum images at 408 and 1420 MHz, and H  i 21-cm line data taken from the Canadian Galactic Plane Survey, molecular observations of the ¹²CO  ( J = 1 → 0)  line at 115 GHz from the Five College Radio Astronomy Observatory and available mid- and far-infrared (FIR) observations obtained with the Midcourse Space Experiment and IRAS satellites, respectively.
This study allowed the identification of the radio continuum and molecular counterpart of the ring nebula SG 13, while no neutral atomic structure was found to be associated. The nebula is also detected in the images in the mid- and FIR, showing the existence of dust well mixed with the ionized gas. We estimate the main physical parameters of the material linked to the nebula.
The interstellar gas distribution in the environs of Mrk 50 is compatible with a stellar wind bubble created by the mass loss from WR 157.
The distribution of young stellar object candidates in the region shows that the stellar formation activity may be present in the molecular shell that encircles the ring nebula.     

Masses and luminosities of O‐ and B‐type stars and red supergiants

Massive stars are of interest as progenitors of supernovae, i.e. neutron stars and black holes, which can be sources of gravitational waves. Recent population synthesis models can predict neutron star and gravitational wave observations but deal with a fixed supernova rate or an assumed initial mass function for the population of massive stars. Here we investigate those massive stars, which are supernova progenitors, i.e. with O‐ and early B‐type stars, and also all supergiants within 3 kpc. We restrict our sample to those massive stars detected both in 2MASS and observed by Hipparcos, i.e. only those stars with parallax and precise photometry. To determine the luminosities we calculated the extinctions from published multi‐colour photometry, spectral types, luminosity class, all corrected for multiplicity and recently revised Hipparcos distances. We use luminosities and temperatures to estimate the masses and ages of these stars using different models from different authors. Having estimated the luminosities of all our stars within 3 kpc, in particular for all O‐ and early B‐type stars, we have determined the median and mean luminosities for all spectral types for luminosity classes I, III, and V. Our luminosity values for supergiants deviate from earlier results: Previous work generally overestimates distances and luminosities compared to our data, this is likely due to Hipparcos parallaxes (generally more accurate and larger than previous ground‐based data) and the fact that many massive stars have recently been resolved into multiples of lower masses and luminosities. From luminosities and effective temperatures we derived masses and ages using mass tracks and isochrones from different authors. From masses and ages we estimated lifetimes and derived a lower limit for the supernova rate of ≈20 events/Myr averaged over the next 10 Myr within 600 pc from the sun. These data are then used to search for areas in the sky with higher likelihood for a supernova or gravitational wave event (like OB associations) (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Massive stars exploding in a He-rich circumstellar medium – I. Type Ibn (SN 2006jc-like) events

We present new spectroscopic and photometric data of the Type Ibn supernovae 2006jc, 2000er and 2002ao. We discuss the general properties of this recently proposed supernova family, which also includes SN 1999cq. The early-time monitoring of SN 2000er traces the evolution of this class of objects during the first few days after the shock breakout. An overall similarity in the photometric and spectroscopic evolution is found among the members of this group, which would be unexpected if the energy in these core-collapse events was dominated by the interaction between supernova ejecta and circumstellar medium. Type Ibn supernovae appear to be rather normal Type Ib/c supernova explosions which occur within a He-rich circumstellar environment. SNe Ibn are therefore likely produced by the explosion of Wolf–Rayet progenitors still embedded in the He-rich material lost by the star in recent mass-loss episodes, which resemble known luminous blue variable eruptions. The evolved Wolf–Rayet star could either result from the evolution of a very massive star or be the more evolved member of a massive binary system. We also suggest that there are a number of arguments in favour of a Type Ibn classification for the historical SN 1885A (S-Andromedae), previously considered as an anomalous Type Ia event with some resemblance to SN 1991bg.