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.     

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.     

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.     

Circumstellar shocks: colliding stellar winds

In this paper I will review some recent developments in the field of circumstellar shocks, particularly as they relate to colliding stellar winds. I shall review the basic physics of colliding winds and shocks, and discuss recent developments in hydrodynamic modelling of colliding winds. I shall also report on recent X-ray observations of shock emission in Wolf-Rayet binary systems where high resolution X-ray spectra of colliding wind shock emission is being seen. I will discuss the occurrence of colliding winds to such diverse systems as Wolf-Rayet binaries, pre-main sequence binaries, symbiotic stars as well as the Galactic center object IRS 7, where recent results on interacting winds are yielded insight into the structure of winds in general.     

Light curves of SN 1998A and comparison with similar unusual SNe

Following the proposal by Damineli that the central object of Eta Carinae may be an early-type binary, we perform numerical simulations of the X-ray emission from colliding stellar winds. A synthetic light curve has been generated which qualitatively agrees with the recent X-ray variability, and provides further support for the binary model. In particular, the model predicts a rise in the observed X-ray emission towards periastron, followed by a sharp drop and subsequent recovery. This is indeed what is seen in the RXTE light curve, although some problems concerning the X-ray spectrum at periastron still need to be explained. The simulations suggest that the width of the periastron dip will provide strong constraints on the binary and stellar wind properties of the components of Eta Car.     

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 .     

Observations of structure in the stellar wind of HD 152408

We report on the presence of substantial, low-velocity stellar wind structure in the extreme O supergiant HD 152408, based on optical spectroscopic time-series observations. Systematic variations in the form of migrating optical depth enhancements occur in the absorption trough of the He I 5876 P Cygni profile. These variations start deep in the stellar wind, slowly accelerate bluewards to 0.5v over 1–2 days, and recur at intervals of about 1 day. Sympathetic variations are apparent in the Balmer emission lines. The observations provide constraints on the stability of the low-velocity stellar wind regime, and indicate the presence of large-amplitude perturbations at great depths in the outflow.     

Stellar wind accretion in GX 301−2: evidence for a high-density stream

The X-ray binary system GX 301−2 consists of a neutron star in an eccentric orbit accreting from the massive early-type star Wray 977. It has previously been shown that the X-ray orbital light curve is consistent with the existence of a gas stream flowing out from Wray 977 in addition to its strong stellar wind. Here, X-ray monitoring observations by the Rossi X-ray Timing Explorer ( RXTE )/All-Sky Monitor and pointed observations by the RXTE /Proportional Counter Array over the past decade are analysed. We analyse both the flux and column density dependence on orbital phase. The wind and stream dynamics are calculated for various system inclinations, companion rotation rates and wind velocities, as well as parametrized by the stream width and density. These calculations are used as inputs to determine both the expected accretion luminosity and the column density along the line-of-sight to the neutron star. The model luminosity and column density are compared to observed flux and column density versus orbital phase, to constrain the properties of the stellar wind and the gas stream. We find that the change between bright and medium intensity levels is primarily due to decreased mass loss in the stellar wind, but the change between medium and dim intensity levels is primarily due to decreased stream density. The mass-loss rate in the stream exceeds that in the stellar wind by a factor of ∼2.5. The quality of the model fits is better for lower inclinations, favouring a higher mass for Wray 977 in its allowed range of  40–60 M_⊙  .     

The symbiotic star CH Cygni – I. Non-thermal bipolar jets

Very Large Array surface brightness and spectral index maps of the evolving extended emission of the triple symbiotic star CH Cygni are presented. These are derived from observations at 4.8, 8.4 and 14 GHz between 1985 and 1999. The maps are dominated by thermal emission around the central bright peak of the nebula, but we also find unambiguous non-thermal emission associated with the extended regions. Our observations confirm that this is a jet. The central region has been associated with the stellar components through Hubble Space Telescope imaging. If the jets are the result of ejection events at outburst, expansion velocities are consistent with those from other measurement methods. We propose that the non-thermal emission is caused by material ejected in the bipolar jets interacting with the circumstellar wind envelope. The resulting shocks lead to local enhancements in the magnetic field from the compact component of the order of 3 mG.     

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.     

A comprehensive variability study of the enigmatic WN8 stars: final results

As a conclusion of our all-sky variability survey of the 'enigmatic' variable WN8 stars, we have carried out coordinated multisite photometric and spectroscopic observations of WN8 stars in 1989 and 1994–1995. We confirm the leading role of the stellar core in restructuring the whole wind. This emerges as a statistical trend: the higher the level of the ∼continuum (i.e. ∼core) light variations, the higher the variability of the P Cygni edges of the optical emission lines. However, the form of the correlation between the light and profile variations is generally different for each individual star. The high level of activity of WN8 stars may be supported/induced by pulsational instability.     

Effects of the stellar wind on X-ray spectra of Cygnus X-3

We study X-ray spectra of Cyg X-3 from BeppoSAX , taking into account absorption and emission in the strong stellar wind of its companion. We find the intrinsic X-ray spectra are well modelled by disc blackbody emission, its upscattering by hot electrons with a hybrid distribution, and by Compton reflection. These spectra are strongly modified by absorption and reprocessing in the stellar wind, which we model using the photoionization code cloudy . The form of the observed spectra implies the wind is composed of two phases. A hot tenuous plasma containing most of the wind mass is required to account for the observed features of very strongly ionized Fe. Small dense cool clumps filling ≲0.01 of the volume are required to absorb the soft X-ray excess, which is emitted by the hot phase but not present in the data. The total mass-loss rate is found to be  (0.6–1.6) × 10⁻⁵ M_⊙ yr⁻¹  . We also discuss the feasibility of the continuum model dominated by Compton reflection, which we find to best describe our data. The intrinsic luminosities of our models suggest that the compact object is a black hole.     

Constraints on black hole accretion in V Puppis

In light of the recent suggestion that the nearby eclipsing binary star system V Puppis has a dark companion on a long orbit, we present the results of radio and X-ray observations of it. We find an upper limit on its radio flux of about 300 μJy and a detection of it in the X-rays with a luminosity of about  3 × 10³¹  erg s⁻¹, a value much lower than what had been observed in some of the low angular resolution surveys of the past. These data are in good agreement with the idea that the X-ray emission from V Puppis comes from mass transfer between the two B stars in the system, but can still accommodate the idea that the X-ray emission comes from the black hole accreting stellar wind from one or both of the B stars.     

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.     

Winds from clusters with non-uniform stellar distributions

We present the analytic and numerical models of the 'cluster wind' resulting from the multiple interactions of the winds ejected by the stars of a dense cluster of massive stars. We consider the case in which the distribution of stars (i.e. the number of stars per unit volume) within the cluster is spherically symmetric, has a power-law radial dependence, and drops discontinuously to zero at the outer radius of the cluster. We carry out comparisons between an analytic model (in which the stars are considered in terms of a spatially continuous injection of mass and energy) and 3D gasdynamic simulations (in which we include 100 stars with identical winds, located in 3D space by statistically sampling the stellar distribution function). From the analytic model, we find that for stellar distributions with steep enough radial dependencies, the cluster wind flow develops a very high central density and a non-zero central velocity, and for steeper dependencies, it becomes fully supersonic throughout the volume of the cluster (these properties are partially reproduced by the 3D numerical simulations). Therefore, the wind solutions obtained for stratified clusters can differ dramatically from the case of a homogeneous stellar distribution (which produces a cluster wind with zero central velocity, and a fully subsonic flow within the cluster radius). Finally, from our numerical simulations, we compute predictions of X-ray emission maps and luminosities, which can be directly compared with observations of cluster wind flows.     

Non-thermal emission processes in massive binaries

In this paper, I present a general discussion of several astrophysical processes likely to play a role in the production of non-thermal emission in massive stars, with emphasis on massive binaries. Even though the discussion will start in the radio domain where the non-thermal emission was first detected, the census of physical processes involved in the non-thermal emission from massive stars shows that many spectral domains are concerned, from the radio to the very high energies. First, the theoretical aspects of the non-thermal emission from early-type stars will be addressed. The main topics that will be discussed are respectively the physics of individual stellar winds and their interaction in binary systems, the acceleration of relativistic electrons, the magnetic field of massive stars, and finally the non-thermal emission processes relevant to the case of massive stars. Second, this general qualitative discussion will be followed by a more quantitative one, devoted to the most probable scenario where non-thermal radio emitters are massive binaries. I will show how several stellar, wind and orbital parameters can be combined in order to make some semi-quantitative predictions on the high-energy counterpart to the non-thermal emission detected in the radio domain. These theoretical considerations will be followed by a census of results obtained so far, and related to this topic. These results concern the radio, the visible, the X-ray and the γ-ray domains. Prospects for the very high energy γ-ray emission from massive stars will also be addressed. Two particularly interesting examples—one O-type and one Wolf-Rayet binary—will be considered in details. Finally, strategies for future developments in this field will be discussed.     

Radio observations of the planetary nebula around the OH/IR star OH 354.88-0.54 (V1018 Sco)

We present radio observations of the unique, recently formed, planetary nebula (PN) associated with a very long-period OH/IR variable star V1018 Sco that is unequivocally still in its asymptotic giant branch phase. Two regions within the optical nebula are clearly detected in non-thermal radio continuum emission, with radio spectral indices comparable to those seen in colliding-wind Wolf–Rayet binaries. We suggest that these represent shocked interactions between the hot, fast stellar wind and the cold nebular shell that represents the PN's slow wind moving away from the central star. This same interface produces both synchrotron radio continuum and the optical PN emission. The fast wind is neither spherical in geometry nor aligned with any obvious optical or radio axis. We also report the detection of transient H₂O maser emission in this nebula.     

Active galactic nuclei jet mass loading and truncation by stellar winds

Active galactic nuclei can produce extremely powerful jets. While tightly collimated, the scale of these jets and the stellar density at galactic centres implies that there will be many jet/star interactions, which can mass load the jet through stellar winds. Previous work employed modest wind mass outflow rates, but this does not apply when mass loading is provided by a small number of high mass-loss stars. We construct a framework for jet mass loading by stellar winds for a broader spectrum of wind mass-loss rates than has previously been considered. Given the observed stellar mass distributions in galactic centres, we find that even highly efficient (0.1 Eddington luminosity) jets from supermassive black holes of masses M _BH≲ 10⁴ M_⊙ are rapidly mass loaded and quenched by stellar winds. For  10⁴ M_⊙ < M _BH < 10⁸ M_⊙  , the quenching length of highly efficient jets is independent of the jet's mechanical luminosity. Stellar wind mass loading is unable to quench efficient jets from more massive engines, but can account for the observed truncation of the inefficient M87 jet, and implies a baryon-dominated composition on scales ≳2 kpc therein even if the jet is initially pair plasma dominated.     

Chemical composition of a sample of bright solar‐metallicity stars

We present a detailed analysis of seven young stars observed with the spectrograph SOPHIE at the Observatoire de Haute‐Provence for which the chemical composition was incomplete or absent in the literature. For five stars, we derived the stellar parameters and chemical compositions using our automatic pipeline optimized for F, G, and K stars, while for the other two stars with high rotational velocity, we derived the stellar parameters by using other information (parallax), and performed a line‐by‐line analysis. Chromospheric emission‐line fluxes from Caii are obtained for all targets. The stellar parameters we derive are generally in good agreement with what is available in the literature. We provide a chemical analysis of two of the stars for the first time. The star HIP 80124 shows a strong Li feature at 670.8 nm implying a high lithium abundance. Its chemical pattern is not consistent with it being a solar sibling, as has been suggested. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The role of non-thermal electrons in the optical continuum of stellar flares

The continuum emission of stellar flares in UV and visible bands can be enhanced by two or even three orders of magnitude relative to the quiescent level and is usually characterized by a blue colour. It is difficult for thermal atmospheric models to reproduce all these spectral features. If the flaring process involves the acceleration of energetic electrons which then precipitate downwards to heat the lower atmosphere, collisional excitation and ionization of ambient hydrogen atoms by these non-thermal electrons could be important in powering the continuum emission. To explore such a possibility, we compute the continuum spectra from an atmospheric model for a dMe star, AD Leo, at its quiescent state, when considering the non-thermal effects by precipitating electron beams. The results show that if the electron beam has an energy flux large enough (for example, ℱ₁∼10¹² erg cm⁻² s⁻¹), the U -band brightening and, in particular, the U − B colour are roughly comparable with observed values for a typical large flare. Moreover, for electron beams with a moderate energy flux ℱ₁≲10¹¹ erg cm⁻² s⁻¹, a decrease of the emission at the Paschen continuum appears. This can explain at least partly the continuum dimming observed in some stellar flares. Adopting an atmospheric model for the flaring state can further raise the continuum flux, but it yields a spectral colour incomparable with observations. This implies that the non-thermal effects may play the chief role in powering the continuum emission in some stellar flares.