Lepto-Hadronic Origin of gamma-rays from the G54.1+0.3 Pulsar Wind Nebula

G54.1+0.3 is a Crab-like pulsar wind nebula (PWN) with the highest $\gamma$-ray to X-ray luminosity ratio among all the nebulae driven by young rotation-powered pulsars. We model the spectral evolution of the PWN and find it difficult to match the observed multi-band data with leptons alone using reasonable model parameters. In lepton-hadron hybrid model instead, TeV photons come mainly from $\pi^0$ decay in proton-proton interaction and the observed photon spectrum can be well reproduced. The newly discovered infrared loop and molecular cloud in or closely around the PWN can work as the target for the bombardment of the PWN protons.

💡 Research Summary

The paper investigates the unusually high γ‑ray‑to‑X‑ray luminosity ratio of the young, Crab‑like pulsar wind nebula (PWN) G54.1+0.3. Using a time‑dependent spectral evolution model that follows the dynamical expansion of the nebula, the authors first test a pure leptonic scenario in which relativistic electrons (and positrons) are the sole radiating particles. By varying key parameters—pulsar spin‑down power, electron injection efficiency, magnetic field strength, and particle spectral index—they find that the leptonic model can reproduce the observed X‑ray synchrotron emission but dramatically under‑predicts the TeV γ‑ray flux measured by VERITAS and MAGIC. Even extreme choices, such as an electron injection efficiency approaching 100 % or a magnetic field as low as a few microgauss, cannot simultaneously match the X‑ray and γ‑ray data without violating physical constraints (e.g., excessive cooling times or unrealistic energy budgets).

To resolve this discrepancy, the authors introduce a lepto‑hadronic hybrid model. In addition to electrons, they assume that a substantial fraction of the pulsar’s spin‑down energy is transferred to relativistic protons (or ions) that remain confined within the nebula for ≳10³ yr. These protons interact with dense target material surrounding the PWN—specifically, an infrared (IR) loop and a molecular cloud that have been identified in recent Spitzer and CO observations. The target density is estimated to be n ≈ 10–100 cm⁻³, corresponding to a total gas mass of order 10³ M⊙.

Proton–proton (p‑p) collisions produce neutral pions (π⁰), which decay almost instantly into two γ‑rays (π⁰ → 2γ). The authors calculate the p‑p interaction rate τ_pp ≈ n σ_pp c t_int, using a high‑energy cross‑section σ_pp ≈ 3 × 10⁻²⁶ cm² and an interaction time t_int set by the confinement time of the protons. With a proton energy budget of ~10⁴⁹ erg (≈10 % of the total spin‑down energy) and the observed target density, the resulting π⁰‑decay γ‑ray spectrum reproduces both the shape and absolute flux of the measured TeV emission. In the hybrid model, the low‑energy (keV) X‑ray band remains dominated by synchrotron radiation from electrons, the GeV band is a mixture of inverse‑Compton scattering (electrons) and π⁰‑decay photons, while the TeV band is essentially pure π⁰‑decay emission. This naturally explains the high γ‑ray‑to‑X‑ray ratio: the hadronic component provides an efficient channel for converting proton kinetic energy into very‑high‑energy photons without requiring unrealistically low magnetic fields.

The paper also discusses observational implications and future tests. High‑resolution γ‑ray imaging with the Cherenkov Telescope Array (CTA) could reveal a spatial morphology that traces the distribution of the IR loop and molecular cloud, a signature of hadronic interactions. Millimeter interferometry (e.g., ALMA) can refine the density and geometry of the target gas, allowing a more precise calculation of the p‑p interaction rate. Additionally, the same p‑p collisions would generate high‑energy neutrinos (∼TeV), offering a complementary probe with IceCube or KM3NeT. Detection of such neutrinos coincident with G54.1+0.3 would provide decisive evidence for a hadronic component.

In summary, the authors demonstrate that a pure leptonic model cannot account for the observed γ‑ray output of G54.1+0.3 under realistic physical conditions. By incorporating a relativistic proton population that interacts with nearby dense material, the lepto‑hadronic model successfully reproduces the multi‑wavelength spectrum and offers a coherent physical picture linking the nebular dynamics, particle acceleration, and the surrounding environment. This work highlights the potential importance of hadronic processes in young PWNe and suggests that other nebulae with similarly high γ‑ray efficiencies may also harbor significant proton populations.