High-resolution simulations of non-thermal emission from LS 5039

TL;DR

This study uses high-resolution simulations to model non-thermal emission from the LS 5039 system, successfully reproducing observed spectral features and orbital light curves, and highlighting the importance of wind dynamics in particle acceleration.

Contribution

The paper introduces a detailed simulation framework combining wind dynamics and energetic particle transport to accurately model non-thermal emission in LS 5039.

Findings

Model reproduces spectral features of LS 5039

Better matches observed orbital light curves in multiple energy bands

Simulations capture key wind-collision region dynamics

Abstract

In a previous study, we investigated the relativistic wind dynamics in the LS 5039 system. In this work, we analyse energetic-particle transport within this modelling context, where we simulate the high-energy particle distribution and ensuing emission of non-thermal radiation. From these high-resolution simulations, we compute the non-thermal emission from this system and compare it to corresponding observations. We modelled the LS 5039 system assuming a wind-driven scenario. Our numerical model uses a joint simulation of the dynamical wind interaction together with the transport of energetic leptons from the shocked pulsar wind. We computed the non-thermal emission from this system in a post-processing step from the resulting distribution of energetic leptons. In this computation, we took into account the synchrotron and inverse Compton emission, relativistic beaming, and {\gamma}{\gamma}-absorption in the stellar radiation field. We investigated the dynamical variation of the energetic particle spectra on both orbital and on short timescales. Our model successfully reproduces many of the spectral features of LS 5039. We also find a better correspondence between our predicted orbital light curves and the corresponding observations in soft x-rays, low-energy, and high-energy gamma rays than in our previous modelling efforts. We find that our high-resolution and large-scale simulations can successfully capture the relevant parts of the wind-collision region that are related to particle acceleration and emission of non-thermal radiation. The quality of the fit strengthens the wind-driven assumption underlying our model. Desirable extensions for the future include a dynamical magnetic-field model for the synchrotron regime, a revision of our injection parameters, and a consideration of an additional hadronic component that could explain recent observations in the 100~TeV regime.