• ### 3D magneto-hydrodynamic models of non-thermal photon emission in the binary system $\gamma^2$ Velorum(1708.05860)

Aug. 19, 2017 astro-ph.HE
Recent reports claiming association of the massive star binary system $\gamma^2$Velorum (WR 11) with a high-energy $\gamma$-ray source observed by \textit{Fermi}-LAT contrast the so-far exclusive role of $\eta$ Carinae as the hitherto only detected $\gamma$-ray emitter in the source class of particle-accelerating colliding-wind binary (CWB) systems. We offer support to this claim of association by providing dedicated model predictions for the nonthermal photon emission spectrum of $\gamma^2$ Velorum. We use three-dimensional magneto-hydrodynamic modeling (MHD) to investigate the structure and conditions of the wind-collision region (WCR) of $\gamma^2$ Velorum including the important effect of radiative braking in the stellar winds. A transport equation is then solved throughout the computational domain to study the propagation of relativistic electrons and protons. The resulting distributions of particles are subsequently used to compute nonthermal photon emission components. In agreement with observation in X-ray spectroscopy, our simulations yield a large shock-cone opening angle. We find the nonthermal $\gamma$-ray emission of $\gamma^2$ Velorum to be of hadronic origin owing to the strong radiation fields in the binary system which inhibit the acceleration of electrons to energies sufficiently high for efficient inverse Compton radiation. We also discuss the strong dependence of a hadronic $\gamma$-ray component on the energy-dependent diffusion used in the simulations. Of two mass-loss rates for the WR star found in literature, only the higher one is able to accommodate the observed $\gamma$-ray spectrum with reasonable values for important simulation parameters such as the injection ratio of high-energy particles within the WCR.
• ### MHD Models of Gamma-ray Emission in WR 11(1701.07284)

Jan. 25, 2017 astro-ph.HE
Recent reports claiming tentative association of the massive star binary system gamma^2 Velorum (WR 11) with a high-energy gamma-ray source observed by Fermi-LAT contrast the so-far exclusive role of Eta Carinae as the hitherto only detected gamma-ray emitter in the source class of particle-accelerating colliding-wind binary systems. We aim to shed light on this claim of association by providing dedicated model predictions for the nonthermal photon emission spectrum of WR 11. We use three-dimensional magneto-hydrodynamic modeling to trace the structure and conditions of the wind-collision region of WR 11 throughout its 78.5 day orbit, including the important effect of radiative braking in the stellar winds. A transport equation is then solved in the wind-collision region to determine the population of relativistic electrons and protons which are subsequently used to compute nonthermal photon emission components. We find that - if WR 11 be indeed confirmed as the responsible object for the observed gamma-ray emission - its radiation will unavoidably be of hadronic origin owing to the strong radiation fields in the binary system which inhibit the acceleration of electrons to energies suffciently high for observable inverse Compton radiation. Different conditions in wind-collision region near the apastron and periastron configuration lead to significant variability on orbital time scales. The bulk of the hadronic gamma-ray emission originates at a 400 solar radii wide region at the apex.
• ### Colliding-wind Binaries with strong magnetic fields(1609.01130)

Sept. 5, 2016 astro-ph.SR, astro-ph.HE
The dynamics of colliding wind binary systems and conditions for efficient particle acceleration therein have attracted multiple numerical studies in the recent years. These numerical models seek an explanation of the thermal and non-thermal emission of these systems as seen by observations. In the non-thermal regime, radio and X-ray emission is observed for several of these colliding-wind binaries, while gamma-ray emission has so far only been found in $\eta$ Carinae and possibly in WR 11. Energetic electrons are deemed responsible for a large fraction of the observed high-energy photons in these systems. Only in the gamma-ray regime there might be, depending on the properties of the stars, a significant contribution of emission from neutral pion decay. Thus, studying the emission from colliding-wind binaries requires detailed models of the acceleration and propagation of energetic electrons. This in turn requires a detailed understanding of the magnetic field, which will not only affect the energy losses of the electrons but in case of synchrotron emission also the directional dependence of the emissivity. In this study we investigate magnetohydrodynamic simulations of different colliding wind binary systems with magnetic fields that are strong enough to have a significant effect on the winds. Such strong fields require a detailed treatment of the near-star wind acceleration zone. We show the implementation of such simulations and discuss results that demonstrate the effect of the magnetic field on the structure of the wind collision region.
• List of contributions from the CTA Consortium presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.
• ### The first full orbit of Eta Carinae seen by Fermi(1503.07637)

March 26, 2015 astro-ph.HE
The binary system eta Carinae has completed its first 5.54y orbit since the beginning of science operation of the Fermi Large Area Telescope (LAT). We are now able to investigate the high-energy gamma-ray source at the position of eta Carinae over its full orbital period. By this, we can address and confirm earlier predictions for temporal and spectral variability. Newer versions of the LAT datasets, instrument response functions and background models allow for a more accurate analysis. Therefore it is important to re-evaluate the previously analyzed time period along with the new data to further constrain location, spectral shape, and flux time history of the gamma-ray source. We confirm earlier predictions of increasing flux values above 10 GeV toward the next periastron passage. For the most recent part of the data sample, flux values as high as those before the first periastron passage in 2008 are recorded. A comparison of spectral energy distributions around periastron and apastron passages reveals strong variation in the high-energy band. This is due to a second spectral component that is present only around periastron. Improved spatial consistency with the gamma-ray source at the position of eta Carinae along with the confirmation of temporal variability above 10 GeV in conjunction with the orbital period strengthens the argument for unambiguous source identification. Spectral variability provides additional constraints for future modeling of the particle acceleration and gamma-ray emission in colliding-wind binary systems.
• ### Simulating three-dimensional nonthermal high-energy photon emission in colliding-wind binaries(1405.6868)

May 27, 2014 astro-ph.HE
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.
• ### High-energy particle transport in 3D hydrodynamic models of colliding-wind binaries(1401.1323)

Jan. 7, 2014 astro-ph.HE
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.
• Compilation of CTA contributions to the proceedings of the 33rd International Cosmic Ray Conference (ICRC2013), which took place in 2-9 July, 2013, in Rio de Janeiro, Brazil
• ### Gamma-ray follow-up studies on Eta Carinae(1203.4939)

Aug. 13, 2012 astro-ph.HE
Observations of high energy gamma rays recently revealed a persistent source in spatial coincidence with the binary system Eta Carinae. Since modulation of the observed gamma-ray flux on orbital time scales has not been reported so far, an unambiguous identification was hitherto not possible. Particularly the observations made by the Fermi Large Area Telescope (LAT) posed additional questions regarding the actual emission scenario owing to the existence of two energetically distinct components in the gamma-ray spectrum of this source, best described by an exponentially cutoff power-law function (CPL) at energies below 10 GeV and a power-law (PL) component dominant at higher energies. The increased exposure in conjunction with the improved instrumental response functions of the LAT now allow us to perform a more detailed investigation of location, spectral shape, and flux time history of the observed gamma-ray emission. For the first time, we are able to report a weak but regular flux decrease over time. This can be understood and interpreted in a colliding-wind binary scenario for orbital modulation of the gamma-ray emission. We find the spectral shape of the gamma-ray signal in agreement with a single emitting particle population in combination with significant absorption by gamma-gamma pair production. Studying the correlation of the flux decrease with the orbital separation of the binary components allows us to predict the behaviour up to the next periastron passage in 2014.