While Supernova Remnants (SNRs) are widely considered the primary accelerators of cosmic rays (CRs) up to hundreds of TeV, they struggle to account for the CR flux at PeV energies, suggesting the existence of additional PeVatrons. Observations from LHAASO (Large High Altitude Air Shower Observatory) have identified several PeVatron candidates, including some SNRs, pulsar wind nebulae, TeV halos and young massive star clusters (YMSCs). These objects accelerate particles that interact with the surrounding interstellar medium and radiation fields, producing very-high- energy gamma rays (>100 TeV), a key signature of both leptonic and hadronic PeVatrons.
We simulate and model the emission of TeV halos and YMSCs, adopting radial emission profiles derived from observational data. Given the current angular resolution of gamma-ray instruments, these profiles often appear similar, making it challenging to distinguish between source classes. We explore how next-generation Imaging Atmospheric Cherenkov Telescopes (IACTs), namely the CTAO (Cherenkov Telescope Array Observatory) and the ASTRI Mini-Array (Astrofisica con Specchi a Tecnologia Replicante Italiana), can classify these sources based on their morphology. We test our classification methods, derived from the profile features of known sources, on simulated CTAO and ASTRI Mini-Array observations of unidentified extended sources from the first LHAASO catalog.
We present the results of our analysis to highlight the potential of future IACT observations in identifying the nature of extended gamma-ray sources, refining PeVatron candidate classifications, and improving our understanding of cosmic-ray accelerators.