Some OB stars show variable non-thermal radio emission. The non-thermal emission is due to synchrotron radiation that is emitted by electrons accelerated to high energies. The electron acceleration occurs at strong shocks created by the collision of radiatively-driven stellar winds in binary systems. Here we present results of our modelling of two colliding wind systems: Cyg OB2 No. 8A and Cyg OB2 No. 9.
We study the non-thermal radio emission of the binary Cyg OB2 No. 8A, to see if it is variable and if that variability is locked to the orbital phase. We investigate if the synchrotron emission generated in the colliding-wind region of this binary can explain the observations and we verify that our proposed model is compatible with the X-ray data. We use both new and archive radio data from the Very Large Array (VLA) to construct a light curve as a function of orbital phase. We also present new X-ray data that allow us to improve the X-ray light curve. We develop a numerical model for the colliding-wind region and the synchrotron emission it generates. The model also includes free-free absorption and emission due to the stellar winds of both stars. In this way we construct artificial radio light curves and compare them with the observed one. The observed radio fluxes show phase-locked variability. Our model can explain this variability because the synchrotron emitting region is not completely hidden by the free-free absorption. In order to obtain a better agreement for the phases of minimum and maximum flux we need to use stellar wind parameters for the binary components which are somewhat different from typical values for single stars. We verify that the change in stellar parameters does not influence the interpretation of the X-ray light curve. Our model has trouble explaining the observed radio spectral index. This could indicate the presence of clumping or porosity in the stellar wind, which - through its influence on both the Razin effect and the free-free absorption - can considerably influence the spectral index. Non-thermal radio emitters could therefore open a valuable pathway to investigate the difficult issue of clumping in stellar winds.
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.
In colliding-wind binaries, shocks accelerate a fraction of the electrons up to relativistic speeds. These electrons then emit synchrotron radiation at radio wavelengths. Whether or not we detect this radiation depends on the size of the free-free absorption region in the stellar winds of both components. One expects long-period binaries to be detectable, but not the short-period ones. It was therefore surprising to find that Cyg OB2 No. 8A (P = 21.9 d) does show variability locked with orbital phase. To investigate this, we developed a model for the relativistic electron generation (including cooling and advection) and the radiative transfer of the synchrotron emission through the stellar wind. Using this model, we show that the synchrotron emitting region in Cyg OB2 No. 8A does extend far enough beyond the free-free absorption region to generate orbit-locked variability in the radio flux. This model can also be applied to other non-thermal emitters and will prove useful in interpreting observations from future surveys, such as COBRaS - the Cyg OB2 Radio Survey.
In the colliding-wind region of massive binaries, non-thermal radio emission occurs. This non-thermal radio emission (due to synchrotron radiation) has so far been observed at centimetre wavelengths. At millimetre wavelengths, the stellar winds and the colliding-wind region emit more thermal free-free radiation, and it is expected that any non-thermal contribution will be difficult or impossible to detect. We aim to determine if the material in the colliding-wind region contributes substantially to the observed millimetre fluxes of a colliding-wind binary. We also try to distinguish the synchrotron emission from the free-free emission. We monitored the massive binary Cyg OB2 #8A at 3 mm with the NOrthern Extended Millimeter Array (NOEMA) interferometer of the Institut de Radioastronomie Millimetrique (IRAM). The data were collected in 14 separate observing runs (in 2014 and 2016), and provide good coverage of the orbital period. The observed millimetre fluxes range between 1.1 and 2.3 mJy, and show phase-locked variability, clearly indicating that a large part of the emission is due to the colliding-wind region. A simple synchrotron model gives fluxes with the correct order of magnitude, but with a maximum that is phase-shifted with respect to the observations. Qualitatively this phase shift can be explained by our neglect of orbital motion on the shape of the colliding-wind region. A model using only free-free emission results in only a slightly worse explanation of the observations. Additionally, on the map of our observations we also detect the O6.5 III star Cyg OB2 #8B, for which we determine a 3 mm flux of 0.21 +- 0.033 mJy. The question of whether synchrotron radiation or free-free emission dominates the millimetre fluxes of Cyg OB2 #8A remains open. More detailed modelling of this system, based on solving the hydrodynamical equations, is required to give a definite answer.
The radio emission from the well-studied massive stellar system Cyg OB2 #5 is known to fluctuate with a period of 6.7 years between a low-flux state when the emission is entirely of free-free origin, and a high-flux state when an additional non-thermal component (of hitherto unknown nature) appears. In this paper, we demonstrate that the radio flux of that non-thermal component is steady on timescales of hours, and that its morphology is arc-like. This shows that the non-thermal emission results from the collision between the strong wind driven by the known contact binary in the system, and that of an unseen companion on a somewhat eccentric orbit with a 6.7-yr period and a 5 to 10 mas semi-major axis. Together with the previously reported wind-collision region located about 0.8 arcsec to the north-east of the contact binary, Cyg OB2 #5 appears to be the only multiple system known so far to harbor two radio-imaged wind-collision regions.