Massive stars in binary systems have long been regarded as potential sources of high-energy gamma rays.The emission is principally thought to arise in the region where the stellar winds collide and accelerate relativistic particles which ubsequently
emit gamma rays. On the basis of a three-dimensional distribution function of high-energy particles in the wind collision region - as obtained by a numerical hydrodynamics and particle transport model - we present the computation of the three-dimensional nonthermal photon emission for a given line of sight. Anisotropic inverse Compton emission is modelled using the target radiation field of both stars. Photons from relativistic bremsstrahlung and neutral pion decay are computed on the basis of local wind plasma densities. We also consider photon photon opacity effects due to the dense radiation fields of the stars. Results are shown for different stellar separations of a given binary system comprising of a B star and a Wolf-Rayet star. The influence of orbital orientation with respect to the line of sight is also studied by using different orbital viewing angles. For the chosen electron-proton injection ratio of 0.01, we present the ensuing photon emission in terms of two-dimensional projections maps, spectral energy distributions and integrated photon flux values in various energy bands. Here, we find a transition from hadron-dominated to lepton-dominated high-energy emission with increasing stellar separations. In addition, we confirm findings from previous analytic modeling that the spectral energy distribution varies significantly with orbital orientation.
Massive stars in binary systems (as WR140, WR147 or $eta$ Carinae) have long been regarded as potential sources of high-energy $gamma$-rays. The emission is thought to arise in the region where the stellar winds collide and produce relativistic parti
cles which subsequently might be able to emit $gamma$-rays. Detailed numerical hydrodynamic simulations have already offered insight in the complex dynamics of the wind collision region (WCR), while independent analytical studies, albeit with simplified descriptions of the WCR, have shed light on the spectra of charged particles. In this paper, we describe a combination of these two approaches. We present a 3D-hydrodynamical model for colliding stellar winds and compute spectral energy distributions of relativistic particles for the resulting structure of the WCR. The hydrodynamic part of our model incorporates the line-driven acceleration of the winds, gravity, orbital motion and the radiative cooling of the shocked plasma. In our treatment of charged particles we consider diffusive shock acceleration in the WCR and the subsequent cooling via inverse Compton losses (including Klein-Nishina effects), bremsstrahlung, collisions and other energy loss mechanisms.
