Population synthesis of Galactic middle-aged pulsar wind nebulae II. Observational signatures of superefficiency

Technical Summary: Population Synthesis of Galactic Middle-Aged Pulsar Wind Nebulae II. Observational Signatures of Superefficiency

Problem Statement
Pulsar wind nebulae (PWNe) interacting with their host supernova remnants (SNRs) can enter a "reverberation phase" driven by the SNR reverse shock. During this phase, the nebula undergoes significant compression, leading to magnetic field amplification and adiabatic heating of particles. A specific outcome of this process is "superefficiency," a state where the integrated radiative luminosity of a PWN in a specific frequency band exceeds the contemporaneous spin-down power ($\dot{E}_{SD}$) of the central pulsar. While previous theoretical work (e.g., Torres & Lin 2018; Torres et al. 2019) established the physical possibility of this phenomenon, its prevalence across the Galactic population and its detailed observational signatures remained largely unconstrained. Furthermore, existing models often rely on simplified "one-zone purely thin-shell" prescriptions that may fail to accurately capture the complex hydrodynamics of reverberation, particle reprocessing, and magnetic field evolution.

Methodology
The authors employ a population synthesis approach to investigate the manifestation of superefficiency across a synthetic Galactic population of middle-aged PWNe (ages $\le 10^5$ yr). The study utilizes the hybrid TIDE+L framework, which combines the computational efficiency of the one-zone TIDE code with Lagrangian modules that self-consistently model the hydrodynamic structure of the surrounding SNR during the reverberation phase. This allows for a realistic treatment of shell thickening, delayed compression, and multiple shock reflections, which are critical for accurate particle and magnetic field evolution.

Key methodological steps include:

• Population Generation: A synthetic population of 1,600 sources (drawn from a total of 2,400 simulated) is generated per realization, consistent with Galactic core-collapse SN rates, progenitor masses (Salpeter IMF), and pulsar birth properties.
• Comparative Modeling: The study compares results from the TIDE+L framework against the standard TIDE model (a purely thin-shell approach without Lagrangian reverberation modules) to isolate the impact of detailed reverberation physics.
• Evolutionary Tracking: The population is tracked through three evolutionary stages: free expansion, reverberation-compressing, and reverberation-post-compression.
• Spectral Analysis: Superefficiency is evaluated across 12 frequency bands (Radio to TeV) by calculating the ratio $L_{band}/\dot{E}_{SD}$. The analysis includes individual source tracking, ensemble spectral energy distributions (SEDs), and electron spectra for subsets defined by frequency band and compression factor (CF).
• Statistical Robustness: Results are derived from 1,000 random realizations of the population to quantify mean expected numbers and stochastic variances.

Key Contributions and Results

1. Prevalence of Superefficiency:

   • Using the TIDE+L framework, the study predicts a substantial number of superefficient sources at the current Galactic epoch. The Far-Infrared (FIR) band shows the highest prevalence, with approximately 315 sources (19.7% of the population) exhibiting superefficiency.
   • Significant numbers are also predicted in the Mid-Infrared (MIR), Near-Infrared (NIR), Optical, UV, and GeV bands.
   • In contrast, the purely thin-shell TIDE model predicts significantly fewer superefficient sources across all bands (e.g., ~116 in FIR vs. 315 in TIDE+L) and essentially none in the X-ray bands. This discrepancy highlights that simplified models underestimate magnetic field amplification and particle survival, leading to an underestimation of superefficiency.

2. Frequency-Dependent Mechanisms:

   • FIR Superefficiency: This is the most common outcome, driven by the accumulation of long-lived, low-energy electrons. These particles radiate efficiently when the magnetic field is amplified during compression. Superefficiency in this band often persists into the post-compression phase.
   • X-ray Superefficiency: This is rarer and more transient, closely tied to the active compression phase. It requires a balance where magnetic field amplification boosts synchrotron luminosity without causing catastrophic synchrotron losses that deplete high-energy particles too quickly.
   • GeV Superefficiency: Driven by Inverse-Compton (IC) scattering, this phenomenon reflects the long-term survival of accumulated electron populations. It is less dependent on instantaneous compression and more on the history of particle injection relative to the declining current $\dot{E}_{SD}$.

3. Evolutionary Trends:

   • Superefficiency is not a static property but evolves with time. The fraction of superefficient sources increases as the population enters the reverberation phase.
   • Low-frequency (FIR) superefficiency persists long after the initial compression, whereas intermediate-frequency (X-ray) superefficiency is more transient, peaking during active compression and declining as high-energy particles cool.
   • The study identifies that superefficiency can arise from two distinct parameter regimes in the $L_0/L_{ch}$ vs. $\tau_0/t_{ch}$ plane: systems with low initial spin-down and long injection timescales, and systems with high initial spin-down that have significantly declined in power, leaving a "relic" luminosity that exceeds the current $\dot{E}_{SD}$.

4. Compression Factor (CF) Dependence:

   • The strength of reverberation (CF) critically regulates the spectral signature. Mild compression (CF < 10) maximizes the number of FIR-superefficient sources by allowing electron accumulation without severe cooling.
   • Extreme compression (CF > 100) leads to catastrophic synchrotron cooling, depleting the electron reservoir and reducing the number of superefficient sources, particularly in the FIR.

5. Equipartition Stability:

   • The analysis of the particle-to-magnetic energy ratio ($E_p/E_B$) reveals that equipartition is not a stable evolutionary attractor. Systems rapidly evolve away from equipartition, with transient crossings occurring during dynamical reprocessing and reverse-shock interactions.

Significance
The paper establishes the first statistical baseline for interpreting candidate superefficient PWNe within a population synthesis framework. It demonstrates that superefficiency is a natural, albeit extreme, consequence of reverberation physics rather than an anomaly requiring anomalous acceleration mechanisms.

The study claims that accurate modeling of the reverberation phase is essential; simplified thin-shell models fail to capture the true prevalence of superefficiency, particularly in the FIR and X-ray bands. By linking population statistics with detailed hydrodynamic evolution, the work provides a unified physical framework for understanding how PWNe can temporarily exceed their pulsar's spin-down power. This has direct implications for the classification of high-energy sources, suggesting that some luminous PWNe may be misidentified or that their acceleration efficiencies have been overestimated if the reverberation-driven energy transfer is ignored. The results guide targeted multi-wavelength follow-up observations, identifying the FIR and GeV bands as the most promising regimes for detecting these phenomena in the Galactic population.