Oblique Shocks at Supernova Remnants in Massive Star Clusters: A Model for the Cosmic-Ray Knee Observed by LHAASO

Abstract This work establishes oblique shocks in massive star clusters (MSCs) as a primary mechanism for accelerating cosmic rays (CRs) up to the knee of the energy spectrum. We develop a model that incorporates the combined contribution of supernova and collective wind shocks, emphasizing the critical role of the shock obliquity angle in determining the maximum particle energy. We illustrate, within our model, that oblique shocks can significantly enhance acceleration efficiency, allowing particles to reach multi-PeV energies in a rigidity-dependent manner. Our preferred model, which incorporates oblique shocks, reproduces the all-particle spectrum and composition observed by the Large High Altitude Air Shower Observatory (LHAASO), interpreting the knee as arising from a sequence of rigidity-dependent cutoffs. The model also predicts subdominant but detectable gamma-ray and neutrino emissions. This study provides an attempt at building a unified framework connecting MSC particle acceleration to the observed features of the CR knee.