High-energy particle transport in 3D hydrodynamic models of colliding-wind binaries

Abstract

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 particles 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.Comment: 28 pages, 9 figures / accepted for publication in The Astrophysical Journa