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.     

Theoretical X-ray line profiles from colliding wind binaries

We present theoretical X-ray line profiles from a range of model colliding wind systems. In particular, we investigate the effects of varying the stellar mass-loss rates, the wind speeds and the viewing orientation. We find that a wide range of theoretical line profile shapes is possible, varying with orbital inclination and phase. At or near conjunction, the lines have approximately Gaussian profiles, with small widths  (HWHM ∼ 0.1 v _∞)  and definite blueshifts or redshifts (depending on whether the star with the weaker wind is in front or behind). When the system is viewed at quadrature, the lines are generally much broader  (HWHM ∼ v _∞)  , flat-topped and unshifted. Local absorption can have a major effect on the observed profiles – in systems with mass-loss rates of a few times  10⁻⁶ M_⊙ yr⁻¹  the lower energy lines  ( E  ≲ 1 keV)  are particularly affected. This generally results in blueward-skewed profiles, especially when the system is viewed through the dense wind of the primary. The orbital variation of the linewidths and shifts is reduced in a low-inclination binary. The extreme case is a binary with   i = 0°  , for which we would expect no line profile variation.     

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.     

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.     

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.     

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.     

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.     

Opening angles, Lorentz factors and confinement of X-ray binary jets

We present a collation of the available data on the opening angles of jets in X-ray binaries, which in most cases are small (≲10°). Under the assumption of no confinement, we calculate the Lorentz factors required to produce such small opening angles via the transverse relativistic Doppler effect. The derived Lorentz factors, which are in most cases lower limits, are found to be large, with a mean >10, comparable to those estimated for active galactic nuclei (AGN) and much higher than the commonly assumed values for X-ray binaries of 2–5. Jet power constraints do not, in most cases, rule out such high Lorentz factors. The upper limits on the opening angles show no evidence for smaller Lorentz factors in the steady jets of Cygnus X-1 and GRS 1915+105. In those sources in which deceleration has been observed (notably  XTE J1550−564  and Cygnus X-3), some confinement of the jets must be occurring, and we briefly discuss possible confinement mechanisms. It is however possible that all the jets could be confined, in which case the requirement for high bulk Lorentz factors can be relaxed.     

Precession and nutation in the η Carinae binary system: evidence from the X-ray light curve

It is believed that η Carinae is actually a massive binary system, with the wind–wind interaction responsible for the strong X-ray emission. Although the overall shape of the X-ray light curve can be explained by the high eccentricity of the binary orbit, other features like the asymmetry near periastron passage and the short quasi-periodic oscillations seen at those epochs have not yet been accounted for. In this paper we explain these features assuming that the rotation axis of η Carinae is not perpendicular to the orbital plane of the binary system. As a consequence, the companion star will face η Carinae on the orbital plane at different latitudes for different orbital phases and, since both the mass-loss rate and the wind velocity are latitude dependent, they would produce the observed asymmetries in the X-ray flux. We were able to reproduce the main features of the X-ray light curve assuming that the rotation axis of η Carinae forms an angle of  29°± 4°  with the axis of the binary orbit. We also explained the short quasi-periodic oscillations by assuming nutation of the rotation axis, with an amplitude of about  5°  and a period of about 22 days. The nutation parameters, as well as the precession of the apsis, with a period of about 274 years, are consistent with what is expected from the torques induced by the companion star.     

A new numerical method for solving radiation driven winds from hot stars

We present a general method for solving the non‐linear differential equation of monotonically increasing steady‐state radiation driven winds. We graphically identify all the singular points before transforming the momentum equation to a system of differential equations with all the gradients explicitly given. This permits a topological classification of all singular points and to calculate the maximum and minimum mass‐loss of the wind. We use our method to analyse for the first time the topology of the non‐rotating frozen‐in ionisation m‐CAK wind, with the inclusion of the finite disk correction factor, and find up to 4 singular points, three of the x‐type and one attractor‐type. The only singular point (and solution passing through) that satisfies the boundary condition at the stellar surface is the standard m‐CAK singular point. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

X-rays from massive OB stars: thermal emission from radiative shocks

Chandra grating spectra of a sample of 15 massive OB stars were analysed under the basic assumption that the X-ray emission is produced in an ensemble of shocks formed in the winds driven by these objects. Shocks develop either as a result of radiation-driven instabilities or due to confinement of the wind by a relatively strong magnetic field, and since they are radiative, a simple model of their X-ray emission was developed that allows a direct comparison with observations. According to our model, the shock structures (clumps, complete or fractional shells) eventually become 'cold' clouds in the X-ray sky of the star. As a result, it is expected that for large covering factors of the hot clumps, there is a high probability for X-ray absorption by the 'cold' clouds, resulting in blueshifted spectral lines. Our analysis has revealed that such a correlation indeed exists for the considered sample of OB stars. As to the temperature characteristics of the X-ray emission plasma, the studied OB stars fall in two groups: (i) one with plasma temperature limited to ∼0.1–0.4 keV and (ii) the other with X-rays produced in plasmas at considerably higher temperatures. We argue that the two groups correspond to different mechanisms for the origin of X-rays: in radiation-driven instability shocks and in magnetically confined wind shocks, respectively.