Population synthesis of Galactic middle-aged pulsar wind nebulae I. Detection prospects for current and future instruments

Enunciato del Problema
Le nebulose a vento di pulsar (PWNe) rappresentano la popolazione più numerosa di sorgenti galattiche di raggi $\gamma$ ad altissima energia (VHE; $E > 100$ GeV) e fungono da laboratori critici per lo studio dell'accelerazione delle particelle e delle interazioni pulsar–resto di supernova (SNR). Sebbene gli strumenti attuali come H.E.S.S., MAGIC e VERITAS abbiano identificato con certezza numerose PWNe, le previsioni realistiche a livello di popolazione hanno storicamente mancato di un trattamento dettagliato della "fase di riverberazione". Questa fase si verifica quando l'onda d'urto inversa dell'SNR viaggia di nuovo verso il centro dell'esplosione, collidendo con la PWN e comprimendola. Questo evento altera significativamente la dinamica della nebulosa, amplifica i campi magnetici, ri-energizza le popolazioni di particelle e aumenta drasticamente le perdite radiative. I precedenti approcci di modellazione hanno affrontato una dicotomia: i semplici modelli a guscio sottile a una zona sono computazionalmente efficienti ma non riescono a catturare la complessa retroazione non lineare della fase di riverberazione, mentre le simulazioni idrodinamiche (HD) e magnetoidrodinamiche (MHD) multidimensionali sono accurate ma proibitivamente costose per studi statistici che coinvolgono migliaia di sorgenti su $10^5$ anni. Di conseguenza, le previsioni per la popolazione galattica di PWN TeV, in particolare per i futuri osservatori, rischiano di essere distorte se questa fase viene trascurata o semplificata eccessivamente.

Metodologia
Per affrontare queste limitazioni, gli autori impiegano il framework ibrido TIDE+L, che combina l'efficienza computazionale di un modello dinamico a guscio sottile (TIDE) per la fase di espansione libera con un trattamento lagrangiano della struttura dell'SNR durante la riverberazione. Questo approccio consente l'evoluzione coerente del sistema PWN–SNR attraverso tutte le fasi, inclusi i complessi episodi di compressione della riverberazione, tracciando lo spettro delle particelle soggetto a perdite radiative e adiabatiche.

Lo studio costruisce una popolazione sintetica di circa 2.400 PWN galattiche, da cui 1.600 sorgenti vengono campionate casualmente in 1.000 realizzazioni per garantire la robustezza statistica. La popolazione è campionata da distribuzioni di:

• Proprietà della pulsar: Periodo di rotazione iniziale (distribuzione gaussiana, troncata a 10 ms), campo magnetico iniziale (log-normale) e indice di frenamento ($n=2.33$).
• Proprietà dell'SNR e dell'ambiente: Energia dell'esplosione (gaussiana troncata), massa degli ejecta (derivata dalla massa della sequenza principale a età zero tramite modelli di evoluzione stellare GENEVA) e densità del mezzo interstellare (ISM) (distribuzione uniforme).
• Parametri spettrali: Spettri di iniezione delle particelle, fattori di contenimento e frazioni magnetiche.

Ogni sistema viene evoluto fino a un'età massima di $10^5$ anni. I flussi di raggi $\gamma$ risultanti vengono quindi confrontati con le soglie di sensibilità e la copertura del cielo degli strumenti attuali e futuri: H.E.S.S., HAWC, l'Osservatorio del Telescopio Cherenkov (CTAO), LHAASO e l'Osservatorio Meridionale Gamma a Grande Campo (SWGO). L'analisi tiene conto degli effetti di estensione della sorgente, che degradano la sensibilità per le sorgenti estese.

Contributi Chiave

1. Modellazione Realistica della Riverberazione: Il documento introduce una sintesi di popolazione statisticamente robusta che incorpora esplicitamente il trattamento lagrangiano della riverberazione, superando le prescrizioni analitiche semplificate dei precedenti modelli a guscio sottile.
2. Applicazione del Framework Ibrido: Dimostra l'applicazione del framework TIDE+L a uno studio di popolazione su larga scala, colmando il divario tra la fisica idrodinamica dettagliata e i requisiti computazionali della sintesi di popolazione.
3. Prospettive di Rilevamento: Lo studio fornisce le prime stime dettagliate della popolazione galattica di PWN rilevabile per gli strumenti di prossima generazione, quantificando specificamente il rendimento di CTAO, SWGO e LHAASO.
4. Fasi Evolutive: Categorizza la popolazione nelle fasi di espansione libera, compressione e post-compressione, analizzando la rilevabilità e le caratteristiche spettrali di ciascuna.

Risultati

• Rilevabilità: Il modello prevede che CTAO rileverà un ordine di grandezza di PWN in più rispetto a quelle attualmente rilevate con certezza nel range TeV. Nello specifico, si prevede che CTAO rileverà circa 172 sorgenti (10,7% della popolazione), salendo a 226 (14,1%) con una configurazione aggiornata "CTAO+". Questo è significativamente superiore alle ~52 sorgenti previste per H.E.S.S.
• Impatto della Riverberazione: I confronti tra il framework TIDE+L e il modello TIDE semplificato rivelano che trascurare la dinamica realistica della riverberazione porta a una sottostima delle sorgenti rilevabili. Il modello semplificato tende a sovrastimare la compressione e l'amplificazione del campo magnetico, risultando in perdite di sincrotrone eccessive e in una popolazione di particelle ad alta energia impoverita, riducendo di conseguenza il flusso di raggi $\gamma$ previsto.
• Distribuzione Evolutiva: La maggior parte della popolazione sintetica (87,5%) si trova nella fase post-compressione. Tuttavia, le sorgenti in questa fase sono tipicamente più vecchie, più grandi e intrinsecamente più deboli, risultando in frazioni di rilevabilità inferiori rispetto a quelle nelle fasi di espansione libera o compressione.
• Evoluzione Spettrale: Studi di caso individuali illustrano che il fattore di compressione (CF) governa l'intensità e la durata dell'amplificazione della luminosità e dell'indurimento/ammorbidimento spettrale durante la riverberazione. Le sorgenti ad alto CF mostrano cambiamenti rapidi e netti nella distribuzione dell'energia spettrale (SED), mentre le sorgenti a basso CF mostrano un'evoluzione più lieve e ritardata.
• Correlazioni: Lo studio rileva una correlazione tra la luminosità integrata 1–10 TeV e la potenza di frenamento della pulsar ($L_{1-10 \text{ TeV}} \propto \dot{E}^{1.15}$). Gli autori notano che recuperare la pendenza più ripida osservata nei dati HGPS ($L \propto \dot{E}^{0.59}$) richiede l'applicazione di tagli di selezione specifici (ad esempio, escludendo sorgenti con bassa potenza di frenamento o età avanzata), suggerendo che i bias dei sondaggi e le limitazioni strumentali influenzano le correlazioni osservate.
• Estensione della Sorgente: Le estensioni angolari simulate sono generalmente coerenti con i dati HGPS, sebbene il modello preveda una frazione più alta di sorgenti con raggi inferiori al minimo osservato in HGPS. Ciò è attribuito all'assunzione di simmetria sferica e alla negligenza del mixing e del moto proprio della pulsar.

Significato
Il documento afferma che una modellazione realistica della fase di riverberazione è essenziale per prevedere accuratamente la popolazione galattica di PWN TeV. Incorporando la dinamica complessa delle interazioni dell'onda d'urto inversa dell'SNR, lo studio fornisce una base teorica più affidabile per interpretare i cataloghi attuali di sorgenti e anticipare i rendimenti dei futuri sondaggi. I risultati confermano che CTAO giocherà un ruolo dominante nel prossimo censimento della popolazione TeV, potenzialmente risolvendo la natura delle sorgenti non identificate ad alto flusso nell'HGPS. Inoltre, il lavoro suggerisce un legame evolutivo naturale tra le PWN di mezza età e l'emergere di aloni TeV estesi, notando che i sistemi più vecchi rilevati dai futuri strumenti sono probabilmente corrispondenti a morfologie simili ad aloni piuttosto che a nebulose compatte classiche. Questo studio serve come prima parte di una serie, stabilendo la base statistica per la visibilità, con un successivo documento pianificato per focalizzarsi sugli esiti di "super-efficienza" della riverberazione.