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.     

Magellanic Cloud WC/WO Wolf–Rayet stars – II. Colliding winds in binaries

A search for evidence of colliding winds is undertaken among the four certain Magellanic Cloud WC/WO spectroscopic binaries found in the companion Paper I, as well as among two Galactic WC/WO binaries of very similar subtype. Two methods of analysis, which allow the determination of orbital inclination and parameters relating to the shock cone from spectroscopic studies of colliding winds, are attempted. In the first method, Lührs' spectroscopic model is fitted to the moderately strong C  iii 5696-Å excess line emission arising in the shock cone for the stars Br22 and WR 9. The four other systems show only very weak C  iii 5696-Å emission. Lührs' model follows well the mean displacement of the line in velocity space, but is unable to reproduce details in the line profile and fails to give a reliable estimate of the orbital inclination. In the second method, an alternative attempt is also made to fit the variation of more global quantities, full width at half-maximum and radial velocity of the excess emission, with phase. This method also gives satisfactory results in a qualitative way, but shows numerical degeneracy with orbital inclination. Colliding wind effects on the very strong C  iv 5808-Å Wolf–Rayet emission line, present in all six binaries, are also found to behave qualitatively as expected. After allowing for line enhancement in colliding wind binaries, it now appears that all Magellanic Cloud WC/WO stars occupy a very narrow range in spectral subclass: WC4/WO3.     

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.     

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 .     

NLTE models of line-driven stellar winds – III. Influence of X-ray radiation on wind structure of O stars

We study the influence of X-rays on the wind structure of selected O stars. For this purpose we use our non-local thermodynamic equilibrium (NLTE) wind code with inclusion of additional artificial source of X-rays, assumed to originate in the wind shocks.
We show that the influence of shock X-ray emission on wind mass-loss rate is relatively small. Wind terminal velocity may be slightly influenced by the presence of strong X-ray sources, especially for stars cooler than   T _eff≲ 35 000 K  .
We discuss the origin of the   L _X/ L ∼ 10⁻⁷  relation. For stars with thick wind this relation can be explained assuming that the cooling time depends on wind density. Stars with optically thin winds exhibiting the 'weak wind problem' display enhanced X-ray emission which may be connected with large shock cooling length. We propose that this effect can explain the 'weak wind problem'.
Inclusion of X-rays leads to a better agreement of the model ionization structure with observations. However, we do not find any significant influence of X-rays on P  v ionization fraction implying that the presence of X-rays cannot explain the P  v problem.
We study the implications of modified ionization equilibrium due to shock emission on the line transfer in the X-ray region. We conclude that the X-ray line profiles of helium-like ions may be affected by the line absorption within the cool wind.     

X-ray Observations of Binary and Single Wolf-Rayet Stars with XMM-Newton and Chandra

We present an overview of recent X-ray observations of Wolf-Rayet (WR) stars with XMM-Newton and Chandra. These observations are aimed at determining the differences in X-ray properties between massive WR + OB binary systems and putatively single WR stars. A new XMM spectrum of the nearby WN8 + OB binary WR 147 shows hard absorbed X-ray emission (including the Fe Kα line complex), characteristic of colliding wind shock sources. In contrast, sensitive observations of four of the closest known single WC (carbon-rich) WR stars have yielded only non-detections. These results tentatively suggest that single WC stars are X-ray quiet. The presence of a companion may thus be an essential factor in elevating the X-ray emission of WC + OB stars to detectable levels.     

Wolf-Rayet stars

Summary Recent literature on Population I Wolf-Rayet star research extending from the Milky Way to blue compact dwarf galaxies is reviewed, broken down into inventory, basic parameters and galactic distribution, atmospheres, binaries, intrinsic variability, mass loss, enrichment and evolution. Also the incidence of Wolf-Rayet stars with variable non-thermal radio emission, excess X-ray fluxes, and episodic/periodic IR excesses is reviewed. These phenomena appear to be associated with wind-wind interaction in wide long-period WR+OB binaries and with wind-compact object interaction in WR+c binaries, with orbit sizes of the order of magnitude of the WR radio photosphere sizes or larger.     

ASCA spectroscopy of the hard X-ray emission from the colliding wind interaction in γ2 Velorum

We discuss an ASCA observation of the eccentric WC8+O7.5 III binary γ ² Velorum near apastron. The X-ray spectrum is compared with two previous observations obtained when the system was near periastron. All three spectra display a hard-emission component that undergoes strong variability over the orbital cycle. The properties of the hard X-ray emission of γ ² Vel are constrained by taking into account the contribution from contaminating soft X-ray sources in the vicinity of γ ² Vel. We find that the observed variations are in qualitative agreement with the predictions of colliding wind models. We investigate for the first time the effect of uncertainties in the chemical composition of the X-ray emitting plasma on our understanding of the high-energy properties of the wind interaction region. Our results indicate that these uncertainties significantly affect the derived shock temperature and absorption column, but play a smaller role in determining the intrinsic X-ray luminosity of the colliding wind zone. We further find that the intrinsic luminosity from the hard X-ray component in γ ² Vel does not follow the 1/ D distance relation expected from simple models of adiabatic shocks.     

Dissipative models of colliding stellar winds — I. Effects of thermal conduction in wide binary systems

The influence of electron thermal conduction on the 2D gas dynamics of colliding stellar winds is investigated. It is shown that, as a result of the non-linear dependence of the electron thermal flux on the temperature, the pre-heating zones (in which the hot gas in the interaction region heats the cool winds in front of the shocks) have finite sizes. The dependence of the problem of the structure of the flow in the interaction region on the dimensionless parameters is studied, and a simple expression is derived for the size of the pre-heating zones at the axis of symmetry. It is shown that small values of the thermal conductivity do not suppress the Kelvin–Helmholtz instability if the adiabatic flow is subject to it. Further studies, both numerical and analytical, in this direction will be of great interest. The influence of thermal conduction on the X-ray emission from the interaction region is also estimated.     

ASCA observations of colliding stellar winds in gamma velorum

We present new high spectral resolution X-ray observations of the colliding wind binary Vel taken with the ASCA satellite. We find two spectral components, one of which is post-shock emission from the colliding winds. Spectral variability is also seen, consistent with current notions of colliding wind phenomena.     

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.     

Dust Formation around Wolf-Rayet Stars

We discuss the observational evidence in favour of dust formation around massive, evolved carbon-rich Wolf-Rayet stars. Large dust shells containing (amorphous carbon) grains are observed in persistent dust makers, such as the pinwheel system WR104 (WC9+OB) and periodic/episodic dust makers, such asWR140 (WC7+O) close to periastron. Problems with the nucleation and growth of dust in such harsh environments are discussed. It is probable that mixing of carbon-rich material from the WC star plus hydrogen from an OB companion, together with efficient shielding via the shocked region of their colliding winds, provides the necessary ingredients. Many problems with the nucleation and growth remain to be addressed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Colliding Winds in Binaries: Observations

Colliding winds in binaries are discussed mainly from an observational point of view. Collisions are especially energetic in the case of hot, luminous stars, which drive strong, fast winds. Emphasis is therefore devoted to binaries containing Wolf-Rayet stars. The subject is divided up into (1) continuum radiation (X-ray and non-thermal radio from the hot bow shock head, IR from dust formed in some WC + O binaries far downstream in the collision shock cone) and (2) line radiation (optical and UV, both from various regions downstream from the bow shock head). The latter is particularly useful in providing constraints on the velocity field and hence ultimately the geometry of the wind collision and the binary system itself. This revised version was published online in July 2006 with corrections to the Cover Date.     

Observational aspects of wind variability

The observational methods to derive information about the instability and variability of winds of hot stars are discussed. The expected effects of blobs and shells in the winds on the P Cygni profiles of the UV resonance lines and on the wind emission lines such as H are described. The evidence for the variability and instability of the winds is derived from the UV observations of Discrete Absorption Components and their variability, the X-ray flickering in compact binaries and the H variations. The evidence for wind variations in Wolf-Rayet stars and Luminous Blue Variables is discussed briefly, as well as a possible explanation for the V/R variations of Be stars.     

X-ray spectroscopy of stars

Non-degenerate stars of essentially all spectral classes are soft X-ray sources. Their X-ray spectra have been important in constraining physical processes that heat plasma in stellar environments to temperatures exceeding one million degrees. Low-mass stars on the cooler part of the main sequence and their pre-main sequence predecessors define the dominant stellar population in the galaxy by number. Their X-ray spectra are reminiscent, in the broadest sense, of X-ray spectra from the solar corona. The Sun itself as a typical example of a main-sequence cool star has been a pivotal testbed for physical models to be applied to cool stars. X-ray emission from cool stars is indeed ascribed to magnetically trapped hot gas analogous to the solar coronal plasma, although plasma parameters such as temperature, density, and element abundances vary widely. Coronal structure, its thermal stratification and geometric extent can also be interpreted based on various spectral diagnostics. New features have been identified in pre-main sequence stars; some of these may be related to accretion shocks on the stellar surface, fluorescence on circumstellar disks due to X-ray irradiation, or shock heating in stellar outflows. Massive, hot stars clearly dominate the interaction with the galactic interstellar medium: they are the main sources of ionizing radiation, mechanical energy and chemical enrichment in galaxies. High-energy emission permits to probe some of the most important processes at work in these stars, and put constraints on their most peculiar feature: the stellar wind. Medium and high- resolution spectroscopy have shed new light on these objects as well. Here, we review recent advances in our understanding of cool and hot stars through the study of X-ray spectra, in particular high-resolution spectra now available from XMM-Newton and Chandra. We address issues related to coronal structure, flares, the composition of coronal plasma, X-ray production in accretion streams and outflows, X-rays from single OB-type stars, massive binaries, magnetic hot objects and evolved WR stars.     

Gamma-ray binaries and related systems

After initial claims and a long hiatus, it is now established that several binary stars emit high- (0.1–100 GeV) and very high-energy (>100 GeV) gamma rays. A new class has emerged called “gamma-ray binaries”, since most of their radiated power is emitted beyond 1 MeV. Accreting X-ray binaries, novae and a colliding wind binary (η Car) have also been detected—“related systems” that confirm the ubiquity of particle acceleration in astrophysical sources. Do these systems have anything in common? What drives their high-energy emission? How do the processes involved compare to those in other sources of gamma rays: pulsars, active galactic nuclei, supernova remnants? I review the wealth of observational and theoretical work that have followed these detections, with an emphasis on gamma-ray binaries. I present the current evidence that gamma-ray binaries are driven by rotation-powered pulsars. Binaries are laboratories giving access to different vantage points or physical conditions on a regular timescale as the components revolve on their orbit. I explain the basic ingredients that models of gamma-ray binaries use, the challenges that they currently face, and how they can bring insights into the physics of pulsars. I discuss how gamma-ray emission from microquasars provides a window into the connection between accretion–ejection and acceleration, while η Car and novae raise new questions on the physics of these objects—or on the theory of diffusive shock acceleration. Indeed, explaining the gamma-ray emission from binaries strains our theories of high-energy astrophysical processes, by testing them on scales and in environments that were generally not foreseen, and this is how these detections are most valuable.     

3D models of radiatively driven colliding winds in massive O+O star binaries – I. Hydrodynamics

The dynamics of the wind–wind collision in massive stellar binaries are investigated using 3D hydrodynamical models which incorporate gravity, the driving of the winds, the orbital motion of the stars and radiative cooling of the shocked plasma. In this first paper, we restrict our study to main-sequence O+O binaries. The nature of the wind–wind collision region is highly dependent on the degree of cooling of the shocked plasma, and the ratio of the flow time-scale of the shocked plasma to the orbital time-scale. The pre-shock wind speeds are lower in close systems as the winds collide prior to their acceleration to terminal speeds. Radiative inhibition may also reduce the pre-shock wind speeds. Together, these effects can lead to rapid cooling of the post-shock gas. Radiative inhibition is less important in wider systems, where the winds are accelerated to higher speeds before they collide, and the resulting collision region can be largely adiabatic. In systems with eccentric orbits, cold gas formed during periastron passage can persist even at apastron, before being ablated and mixed into its surroundings and/or accelerated out of the system.     

Colliding Winds in Binary Star Systems: Theory, Models

The properties of the interaction region of colliding winds in hot star binaries and evolved low mass binaries are discussed in the frame of 2D and 3D hydrodynamical numerical simulations. The emission of such systems and their nebulae, in particular in X-rays, is reviewed and new results regarding orbital line profile variations are presented. The connection between colliding winds, thin shell instabilities and structure formation is outlined. This revised version was published online in July 2006 with corrections to the Cover Date.     

Neutron star retention and millisecond pulsar production in globular clusters

We investigate the conditions by which neutron star retention in globular clusters is favoured. We find that neutron stars formed in massive binaries are far more likely to be retained. Such binaries are likely to then evolve into contact before encountering other stars, possibly producing a single neutron star after a common envelope phase. A large fraction of the single neutron stars in globular clusters are then likely to exchange into binaries containing moderate-mass main-sequence stars, replacing the lower-mass components of the original systems. These binaries will become intermediate-mass X-ray binaries (IMXBs), once the moderate-mass star evolves off the main sequence, as mass is transferred on to the neutron star, possibly spinning it up in the process. Such systems may be responsible for the population of millisecond pulsars (MSPs) that has been observed in globular clusters. Additionally, the period of mass-transfer (and thus X-ray visibility) in the vast majority of such systems will have occurred 5–10 Gyr ago, thus explaining the observed relative paucity of X-ray binaries today, given the MSP population.