No Arabic abstract
We present a model for the non-thermal emission from a colliding-wind binary. Relativistic protons and electrons are assumed to be accelerated through diffusive shock acceleration (DSA) at the global shocks bounding the wind-wind collision region. The non-linear effects of the back-reaction due to the cosmic ray pressure on the particle acceleration process and the cooling of the non-thermal particles as they flow downstream from the shocks are included. We explore how the non-thermal particle distribution and the keV-GeV emission changes with the stellar separation and the viewing angle of the system, and with the momentum ratio of the winds. We confirm earlier findings that DSA is very efficient when magnetic field amplification is not included, leading to significantly modified shocks. We also find that the non-thermal flux scales with the binary separation in a complicated way and that the anisotropic inverse Compton emission shows only a moderate variation with viewing angle due to the spatial extent of the wind-wind collision.
We present a model for the creation of non-thermal particles via diffusive shock acceleration in a colliding-wind binary. Our model accounts for the oblique nature of the global shocks bounding the wind-wind collision region and the finite velocity of the scattering centres to the gas. It also includes magnetic field amplification by the cosmic ray induced streaming instability and the dynamical back reaction of the amplified field. We assume that the injection of the ions and electrons is independent of the shock obliquity and that the scattering centres move relative to the fluid at the Alfv{e}n velocity (resulting in steeper non-thermal particle distributions). We find that the Mach number, Alfv{e}nic Mach number, and transverse field strength vary strongly along and between the shocks, resulting in significant and non-linear variations in the particle acceleration efficiency and shock nature (turbulent vs. non-turbulent). We find much reduced compression ratios at the oblique shocks in most of our models compared to our earlier work, though total gas compression ratios that exceed 20 can still be obtained in certain situations. We also investigate the dependence of the non-thermal emission on the stellar separation and determine when emission from secondary electrons becomes important. We finish by applying our model to WR 146, one of the brightest colliding wind binaries in the radio band. We are able to match the observed radio emission and find that roughly 30 per cent of the wind power at the shocks is channelled into non-thermal particles.
Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study, we still do not have a complete census of acceleration mechanisms. The collision of strong stellar winds in massive binary systems creates powerful shocks, which have been expected to produce high-energy cosmic-rays through Fermi acceleration at the shock interface. The accelerated particles should collide with stellar photons or ambient material, producing non-thermal emission observable in X-rays and gamma-rays. The supermassive binary star eta Carinae drives the strongest colliding wind shock in the solar neighborhood. Observations with non-focusing high-energy observatories indicate a high energy source near eta Carinae, but have been unable to conclusively identify eta Carinae as the source because of their relatively poor angular resolution. Here we present the first direct focussing observations of the non-thermal source in the extremely hard X-ray band, which is found to be spatially coincident with the star within several arc-seconds. These observations show that the source of non-thermal X-rays varies with the orbital phase of the binary, and that the photon index of the emission is similar to that derived through analysis of the gamma-ray spectrum. This is conclusive evidence that the high-energy emission indeed originates from non-thermal particles accelerated at colliding wind shocks.
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.
The binary stellar system HD 93129A is one of the most massive known binaries in our Galaxy. This system presents non-thermal emission in the radio band, which can be used to infer its physical conditions and predict its emission in the high-energy band. We intend to constrain some of the unknown parameters of HD 93129A through modelling the non-thermal emitter, and also to analyse the detectability of this source in hard X-rays and $gamma$-rays. We develop a broadband radiative model for the wind-collision region taking into account the evolution of the accelerated particles streaming along the shocked region, the emission by different radiative processes, and the attenuation of the emission propagating through the local matter and radiation fields. From the analysis of the radio emission, we find that the binary HD~93129A is more likely to have a low inclination and a high eccentricity. The minimum energy of the non-thermal electrons seems to be between $sim 20 - 100$MeV, depending on the intensity of the magnetic field in the wind-collision region. The latter can be in the range $sim 20 - 1500$ mG. Our model is able to reproduce the observed radio emission, and predicts that the non-thermal radiation from HD~93129A will increase in the near future. With instruments such as textit{NuSTAR}, textit{Fermi}, and CTA, it will be possible to constrain the relativistic particle content of the source, and other parameters such as the magnetic field strength in the wind collision zone, which in turn can be used to obtain upper-limits of the magnetic field on the surface of the very massive stars, thereby inferring whether magnetic field amplification is taking place in the particle acceleration region.
Non-thermal emission has been detected in WR-stars for many years at long wavelengths spectral range, in general attributed to synchrotron emission. Two key ingredients are needed to explain such emissions, namely magnetic fields and relativistic particles. Particles can be accelerated to relativistic speeds by Fermi processes at strong shocks. Therefore, strong synchrotron emission is usually attributed to WR binarity. The magnetic field may also be amplified at shocks, however the actual picture of the magnetic field geometry, intensity, and its role on the acceleration of particles at WR binary systems is still unclear. In this work we discuss the recent developments in MHD modelling of wind-wind collision regions by means of numerical simulations, and the coupled particle acceleration processes related.