Deciphering the gamma-ray emission in the Cygnus region

Abstract: The Cygnus region is a vast star-forming complex harbouring a population of powerful objects, including massive star clusters and associations, Wolf-Rayet stars, pulsars, and supernova remnants. The multi-wavelength picture is far from understood, in particular the recent LHAASO detection of multi-degree scale diffuse gamma-ray emission up to PeV energies. We aim to model the broadband gamma-ray data, discriminating plausible scenarios amongst all candidate accelerators. We consider in particular relic hadronic emission from a supernova remnant expanding in a low-density environment and inverse Compton emission from stellar-wind termination shocks in the Cygnus OB2 stellar association. We first estimate the maximum particle energy from a 3D hydrodynamical simulation of the supernova remnant scenario. The transport equation is then solved numerically to determine the radial distribution of non-thermal protons and electrons. In order to compute synthetic gamma-ray spectra and emission maps, we develop a 3D model of the gas distribution. This includes, firstly, a HI component with a low-density superbubble around Cygnus OB2 and, secondly, molecular clouds lying at the edge of the superbubble and in the foreground. We find that a powerful, ~50 kyr-old supernova remnant can account for both the morphology and spectrum from 10 TeV-PeV. At PeV energies, the microquasar Cygnus X-3 and diffuse Galactic cosmic rays might also contribute to the flux. Below about 10 TeV, hadronic models are incompatible with the expected existence of a superbubble centred on Cygnus OB2. Instead, the spectrum is well fitted with inverse Compton emission from electrons accelerated at stellar-wind termination shocks in Cygnus OB2 in line with existing multi-wavelength limits.