The interaction between the strong winds in stellar colliding-wind binary (CWB) systems produces two shock fronts, delimiting the wind collision region (WCR). There, particles are expected to be accelerated mainly via diffusive shock acceleration (DSA). We investigate the injection and the acceleration of protons in typical CWB systems by means of Monte Carlo simulations, with both a test-particle approach and a non-linear method modelling a shock locally modified by the backreaction of the accelerated protons. We use magnetohydrodynamic simulations to determine the background plasma in the WCR and its vicinity. This allows us to consider particle acceleration at both shocks, on either side of the WCR, with a realistic large-scale magnetic field. We highlight the possible effects of particle acceleration on the local shock profiles at the WCR. We include the effect of magnetic field amplification due to resonant streaming instability (RSI), and compare results without and with the backreaction of the accelerated protons. In the latter case we find a lower flux of the non-thermal proton population, and a considerable magnetic field amplification. This would significantly increase the synchrotron losses of relativistic electrons accelerated in CWB systems, lowering the maximal energy they can reach and strongly reducing the inverse Compton fluxes. As a result, $\gamma$-rays from CWBs would be predominantly due to the decay of neutral pions produced in nucleon-nucleon collisions. This might provide a way to explain why, in the vast majority of cases, CWB systems have not been identified as $\gamma$-ray sources, while they emit synchrotron radiation.

The interaction between the strong winds of the stars in colliding-wind binary (CWB) systems produces two shock fronts, delimiting the wind collision region (WCR). There, particles are expected to be accelerated mainly via diffusive shock acceleration, and to produce $\gamma$-rays, in processes involving relativistic electrons and/or protons. We investigate the injection and the acceleration of protons in typical CWB systems by means of Monte Carlo simulations, with a test-particle approach. We use magnetohydrodynamic simulations to determine the background conditions in the wind collision region. This allows us to consider particle acceleration at both shocks, on either side of the WCR, with a self-consistently determined large-scale magnetic field, which has an impact on the shape of the WCR, and the topology of which plays an important role in particle acceleration at collisionless shocks. Such studies may contribute to improve $\gamma$-ray flux predictions for CWB systems.

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.