The star Cyg OB2 No. 9 is a well-known non-thermal radio emitter. Recent theoretical work suggests that all such O-stars should be in a binary or a multiple system. However, there is no spectroscopic evidence of a binary component. Re-analysis of radio observations from the VLA of this system over 25 years has revealed that the non-thermal emission varies with a period of 2.35+-0.02 yr. This is interpreted as a strong suggestion of a binary system, with the non-thermal emission arising in a wind-collision region. We derived some preliminary orbital parameters for this putative binary and revised the mass-loss rate of the primary star downward from previous estimates.

We present one-dimensional, time-dependent models of the clumps generated by the line-deshadowing instability. In order to follow the clumps out to distances of more than 1000 stellar radii, we use an efficient moving-box technique. We show that, within the approximations, the wind can remain clumped well into the formation region of the radio continuum.

HD 167971 is a triple system consisting of a 3.3-day eclipsing binary (O5-8 V + O5-8 V) and an O8 supergiant. It is also a well known non-thermal radio emitter. We observed the radio emission of HD 167971 with the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA). By combining these data with VLA archive observations we constructed a radio lightcurve covering a 20-yr time-range. We searched for, but failed to find, the 3.3-day spectroscopic period of the binary in the radio data. This could be due to the absence of intrinsic synchrotron radiation at the colliding-wind region between the two components of the eclipsing binary, or due to the large amount of free-free absorption that blocks the synchrotron radiation. We are able to explain many of the observed characteristics of the radio data if the non-thermal emission is produced in a colliding-wind region between the supergiant and the combined winds of the binary. Furthermore, if the system is gravitationally bound, the orbital motion occurs over a period of ~ 20 years, or longer, as suggested by the long-term variability in the radio data. We argue that the variability is due to the free-free absorption that changes with orbital phase or may also in part be due to changes in separation, should the orbit be eccentric.

We present a model for the non-thermal radio emission from presumably single O stars, in terms of synchrotron emission from relativistic electrons accelerated in wind-embedded shocks. These shocks are associated with an unstable, chaotic wind. The main improvement with respect to earlier models is the inclusion of the radial dependence of the shock velocity jump and compression ratio, based on 1D hydrodynamical simulations. The decrease of the velocity jump and the compression ratio as a function of radius produces a rapidly decreasing synchrotron emissivity. This effectively prohibits the models from reproducing the spectral shape of the observed non-thermal radio emission. We investigate a number of "escape routes" by which the hydrodynamical predictions might be reconciled with the radio observations. Although these escape routes reproduce the observed spectral shape, none of these escape routes are physically plausible. In particular, re-acceleration by feeding an electron distribution through a number of shocks, is in contradiction with current hydrodynamical simulations. These hydrodynamical simulations have their limitations, most notably the use of 1D. At present, it is not feasible to perform 2D simulations of the wind out to the distances required for synchrotron-emission models. Based on the current hydrodynamic models, we suspect that the observed non-thermal radio emission from O stars cannot be explained by wind-embedded shocks associated with the instability of the line-driving mechanism. The most likely alternative mechanism is synchrotron emission from colliding winds. That would imply that all O stars with non-thermal radio emission should be members of binary or multiple systems.

We discuss the results of two XMM-Newton observations of the open cluster NGC6604 obtained in April and September 2002. We concentrate mainly on the multiple system HD167971 (O5-8V + O5-8V + (O8I)). The soft part of the EPIC spectrum of this system is thermal with typical temperatures of about 2 10^6 to 9 10^6 K. The nature (thermal vs non-thermal) of the hard part of the spectrum is not unambiguously revealed by our data. If the emission is thermal, the high temperature of the plasma (~ 2.3 10^7 to 4.6 10^7 K) would be typical of what should be expected from a wind-wind interaction zone within a long period binary system. This emission could arise from an interaction between the combined winds of the O5-8V + O5-8V close binary system and that of the more distant O8I companion. Assuming instead that the hard part of the spectrum is non-thermal, the photon index would be rather steep (~3). Moreover, a marginal variability between our two XMM-Newton pointings could be attributed to an eclipse of the O5-8V + O5-8V system. The overall X-ray luminosity points to a significant X-ray luminosity excess of about a factor 4 possibly due to colliding winds. Considering HD167971 along with several recent X-ray and radio observations, we propose that the simultaneous observation of non-thermal radiation in the X-ray (below 10.0 keV) and radio domains appears rather unlikely. Our investigation of our XMM-Newton data of NGC6604 reveals a rather sparse distribution of X-ray emitters. Including the two bright non-thermal radio emitters HD168112 and HD167971, we present a list of 31 X-ray sources... (see paper for full abstract)

We present a radio lightcurve of the O5.5 III(f+) star HD 168112, based on archive data from the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA). The fluxes show considerable variability and a negative spectral index, thereby confirming that HD 168112 is a non-thermal radio emitter. The non-thermal radio emission is believed to be due to synchrotron radiation from relativistic electrons that have been Fermi accelerated in shocks. For HD 168112, it is not known whether these shocks are due to a wind-wind collision in a binary system or to the intrinsic instability of the stellar wind driving mechanism. Assuming HD 168112 to be a single star, our synchrotron model shows that the velocity jump of the shocks should be very high, or there should be a very large number of shocks in the wind. Neither of these is compatible with time-dependent hydrodynamical calculations of O star winds. If, on the other hand, we assume that HD 168112 is a binary, the high velocity jump is easily explained by ascribing it to the wind-wind collision. By further assuming the star to be an eccentric binary, we can explain the observed radio variability by the colliding-wind region moving in and out of the region where free-free absorption is important. The radio data presented here show that the binary has a period between one and two years. By combining the radio data with X-ray data, we find that the most likely period is about 1.4 yr.

We present a model for the non-thermal radio emission from bright O stars, in terms of synchrotron emission from wind-embedded shocks. The model is an extension of an earlier one, with an improved treatment of the cooling of relativistic electrons. This improvement limits the synchrotron-emitting volume to a series of fairly narrow layers behind the shocks. We show that the width of these layers increases with increasing wavelength, which has important consequences for the shape of the spectrum. We also show that the strongest shocks produce the bulk of the emission, so that the emergent radio flux can be adequately described as coming from a small number of shocks, or even from a single shock. A single shock model is completely determined by four parameters: the position of the shock, the compression ratio and velocity jump of the shock, and the surface magnetic field. Applying a single shock model to the O5 If star Cyg OB2 No. 9 allows a good determination of the compression ratio and shock position and, to a lesser extent, the magnetic field and velocity jump. Our main conclusion is that strong shocks need to survive out to distances of a few hundred stellar radii. Even with multiple shocks, the shocks needed to explain the observed emission are stronger than predictions from time-dependent hydrodynamical simulations.

We report the results of a multiwavelength study of the non-thermal radio emitter HD168112 (O5.5III(f^+)). The detailed analysis of two quasi-simultaneous XMM-Newton and VLA observations reveals strong variability of this star both in the X-ray and radio ranges. The X-ray observations separated by five months reveal a decrease of the X-ray flux of ~30%. The radio emission on the other hand increases by a factor 5-7 between the two observations obtained roughly simultaneously with the XMM-Newton pointings. The X-ray data reveal a hard emission that is most likely produced by a thermal plasma at kT ~2-3 keV while the VLA data confirm the non-thermal status of this star in the radio waveband. Comparison with archive X-ray and radio data confirms the variability of this source in both wavelength ranges over a yet ill defined time scale. The properties of HD168112 in the X-ray and radio domain point towards a binary system with a significant eccentricity and an orbital period of a few years. However, our optical spectra reveal no significant changes of the star's radial velocity suggesting that if HD168112 is indeed a binary, it must be seen under a fairly low inclination.

We present an efficient technique to study the 1D evolution of instability-generated structure in winds of hot stars out to very large distances (more than 1000 stellar radii). This technique makes use of our previous finding that external forces play little role in the outer evolution of structure. Rather than evolving the entire wind, as is traditionally done, the technique focuses on a representative portion of the structure and follows it as it moves out with the flow. This requires the problem to be formulated in a moving reference frame. The lack of Galilean invariance of the spherical equations of hydrodynamics is circumvented by recasting them in a pseudo-planar form. By applying the technique to a number of problems we show that it is fast and accurate, and has considerable conceptual advantages. It is particularly useful to test the dependence of solutions on the Galilean frame in which they were obtained. As an illustration, we show that, in a one-dimensional approximation, the wind can remain structured out to distances of more than 1300 stellar radii from the central star.

We present a theoretical model for the non-thermal radio emission from single hot stars, in terms of synchrotron radiation from electrons accelerated in wind-embedded shocks. The model is described by five independent parameters each with a straightforward physical interpretation. Applying the model to a high-quality observation of Cyg OB2 No. 9 (O5 If), we obtain meaningful constraints on most parameters. The most important result is that the outer boundary of the synchrotron emission region must lie between 500 and 2200 stellar radii. This means that shocks must persist up to that distance. We also find that relatively weak shocks (with a compression ratio < 3) are needed to produce the observed radio spectrum. These results are compatible with current hydrodynamical predictions. Most of our models also show a relativistic electron fraction that increases outwards. This points to an increasing efficiency of the acceleration mechanism, perhaps due to multiple acceleration, or an increase in the strength of the shocks. Implications of our results for non-thermal X-ray emission are discussed.

We present radio and submillimetre observations of the O4I(n)f star zeta Pup, and discuss structure in the outer region of its wind (~ 10-100 R*). The properties of bremsstrahlung, the dominant emission process at these wavelengths, make it sensitive to structure and allow us to study how the amount of structure changes in the wind by comparing the fluxes at different wavelengths. To look for variability, we acquired 3.6 and 6 cm observations with the Australia Telescope Compact Array (ATCA). We supplemented these with archive observations from the NRAO Very Large Array (VLA). We did not find variability at more than the +- 20 % level. The long integration time does allow an accurate determination of the fluxes at 3.6 and 6 cm. Converting these fluxes into a mass loss rate, we find Mdot = 3.5 x 10^{-6} Msun/yr. This value confirms the significant discrepancy with the mass loss rate derived from the Halpha profile, making zeta Pup an exception to the usually good agreement between Halpha and radio mass loss rates. We also observed zeta Pup at 850 mum with the James Clerk Maxwell Telescope (JCMT) and at 20 cm with the VLA. A smooth wind model shows that the millimetre fluxes are too high compared to the radio fluxes. While recombination of helium in the outer wind cannot be discounted as an explanation, the wealth of evidence for structure strongly suggests this as the explanation for the discrepancy. Model calculations show that the structure needs to be present in the inner ~ 70 R* of the wind, but that it decays significantly, or maybe even disappears, beyond that radius.