A young SNR illuminating nearby Molecular Clouds with cosmic rays

TL;DR

This paper models a young supernova remnant and nearby molecular clouds to explain observed gamma-ray emissions as resulting from cosmic rays accelerated by the SNR and illuminating the clouds, using detailed diffusion and environmental models.

Contribution

It introduces a comprehensive 3D modeling approach of SNR evolution, cosmic ray diffusion, and molecular cloud structure to explain gamma-ray observations in a young SNR system.

Findings

TeV emission from runaway cosmic rays can explain HESS J1729-345.

SNR shock speed consistent with progenitors of 20-25 solar masses.

TeV spectrum fits with cosmic rays escaping the SNR.

Abstract

The Supernova Remnant (SNR) HESS J1731-347 displays strong non-thermal TeV gamma-ray and X-ray emission, thus the object is at present time accelerating particles to very high energies. A distinctive feature of this young SNR is the nearby (~30 pc in projection) extended source HESS J1729-345, which is currently unidentified but is in spatial projection coinciding with known molecular clouds (MC). We model the SNR evolution to explore if the TeV emission from HESS J1729-345 can be explained as emission from runaway hadronic cosmic rays (CRs) that are illuminating these MCs. The observational data of HESS J1729-345 and HESS J1731-347 can be reproduced using core-collapse SN models for HESS J1731-347. Starting with different progenitor stars and their pre-supernova environment, we model potential SNR evolution histories along with the CR acceleration in the SNR and the diffusion of the CRs. A simplified 3-dimensional structure of the MCs is introduced based on 12CO data, adopting a distance of 3.2 kpc to the source. A Monte Carlo-based diffusion model for the escaping CRs is developed to deal with the inhomogeneous environment. The fast SNR forward shock speed as implied from the X-ray data can easily be explained when employing scenarios with progenitor star masses between 20 and 25 solar masses, where the SNR shock is still expanding inside the main sequence (MS)-bubble at present time. The TeV spectrum of HESS J1729-345 is satisfactorily fitted by the emission from the highest-energy CRs that have escaped the SNR, using a standard galactic CR diffusion coefficient in the inter-clump medium. The TeV image of HESS J1729-345 can be explained with a reasonable 3-dimensional structure of MCs. The TeV emission from the SNR itself is dominated by leptonic emission in this model. We also explore scenarios where the shock is starting to encounter the dense MS progenitor wind bubble shell.