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

Probleemstelling
Pulsar-windnevels (PWNe) vertegenwoordigen de grootste populatie van galactische bronnen met zeer hoge energie (VHE; $E > 100$ GeV) $\gamma$-straling en fungeren als cruciale laboratoria voor het bestuderen van deeltjesversnelling en interacties tussen pulsars en supernovaremnants (SNR). Hoewel huidige instrumenten zoals H.E.S.S., MAGIC en VERITAS talrijke PWNe stevig hebben geïdentificeerd, ontbraken realistische voorspellingen op populatieniveau historisch een gedetailleerde behandeling van de "reverberatiefase". Deze fase treedt op wanneer de terugkerende schokgolf van het SNR terugreist naar het explosiecentrum, waarbij deze botst met en de PWN comprimeert. Deze gebeurtenis verandert de dynamiek van de nevel aanzienlijk, versterkt magnetische velden, herlaadt deeltjespopulaties en verhoogt stralingsverliezen drastisch. Eerdere modelleringbenaderingen stonden voor een dichotomie: eenvoudige één-zone dunne-schilmodellen zijn computationeel efficiënt maar falen in het vastleggen van de complexe, niet-lineaire terugkoppeling van de reverberatiefase, terwijl meerdimensionale hydrodynamische (HD) en magnetohydrodynamische (MHD) simulaties accuraat zijn maar prohibitief duur voor statistische studies die duizenden bronnen over $10^5$ jaar omvatten. Bijgevolg lopen voorspellingen voor de galactische TeV PWN-populatie, met name voor toekomstige observatoria, het risico op vertekening als deze fase wordt verwaarloosd of te sterk wordt vereenvoudigd.

Methodologie
Om deze beperkingen aan te pakken, hanteren de auteurs het hybride TIDE+L-kader, dat de computationele efficiëntie combineert van een dynamisch dunne-schilmodel (TIDE) voor de vrije-expansiefase met een Lagrangiaanse behandeling van de SNR-structuur tijdens de reverberatie. Deze aanpak staat toe voor de zelfconsistentie-evolutie van het PWN-SNR-systeem over alle stadia, inclusief de complexe compressie-episoden van de reverberatie, terwijl het het deeltjesspectrum bijhoudt dat onderhevig is aan stralings- en adiabatische verliezen.

De studie construeert een synthetische populatie van ongeveer 2.400 galactische PWNe, waaruit 1.600 bronnen willekeurig worden bemonsterd in 1.000 realisaties om statistische robuustheid te waarborgen. De populatie wordt bemonsterd uit verdelingen van:

• Pulsareigenschappen: Initiële rotatieperiode (Gaussische verdeling, afgesneden bij 10 ms), initieel magnetisch veld (log-normaal) en remmingsindex ($n=2.33$).
• SNR- en omgevings eigenschappen: Explosie-energie (afgesneden Gaussisch), uitwerpmassa (afgeleid van de massa op het nulpunt van de hoofdreeks via GENEVA-ster-evolutiemodellen) en dichtheid van het interstellaire medium (ISM) (uniforme verdeling).
• Spectrale parameters: Deeltjesinjectiespectra, opsluitingsfactoren en magnetische fracties.

Elk systeem wordt geëvolueerd tot een maximale leeftijd van $10^5$ jaar. De resulterende $\gamma$-stralingfluxen worden vervolgens vergeleken met de gevoeligheidsdrempels en hemeldekking van huidige en toekomstige instrumenten: H.E.S.S., HAWC, het Cherenkov Telescope Array Observatory (CTAO), LHAASO en het Southern Wide-field Gamma-ray Observatory (SWGO). De analyse houdt rekening met bronuitbreidingseffecten, die de gevoeligheid voor uitgebreide bronnen verslechteren.

Belangrijkste bijdragen

1. Realistische Reverberatiemodellering: Het artikel introduceert een statistisch robuuste populatiesynthese die expliciet de Lagrangiaanse behandeling van reverberatie incorporateert, verdergaand dan de vereenvoudigde analytische voorschriften van eerdere dunne-schilmodellen.
2. Toepassing van Hybride Kader: Het demonstreert de toepassing van het TIDE+L-kader op een grootschalig populatiestudie, waardoor de kloof wordt overbrugd tussen gedetailleerde hydrodynamische fysica en de computationele eisen van populatiesynthese.
3. Detectievooruitzichten: De studie biedt de eerste gedetailleerde schattingen van de detecteerbare galactische PWN-populatie voor instrumenten van de volgende generatie, specifiek kwantificerend de opbrengst van CTAO, SWGO en LHAASO.
4. Evolutionaire Stadia: Het categoriseert de populatie in vrije-expansie-, comprimerende en post-compressiefases, waarbij de detecteerbaarheid en spectrale kenmerken van elk worden geanalyseerd.

Resultaten

• Detecteerbaarheid: Het model voorspelt dat CTAO een orde van grootte meer PWNe zal detecteren dan die momenteel stevig zijn gedetecteerd in het TeV-bereik. Specifiek wordt verwacht dat CTAO ongeveer 172 bronnen detecteert (10,7% van de populatie), wat stijgt tot 226 (14,1%) met een geüpgradede "CTAO+" configuratie. Dit is significant hoger dan de ~52 bronnen die voor H.E.S.S. worden voorspeld.
• Impact van Reverberatie: Vergelijkingen tussen het TIDE+L-kader en het vereenvoudigde TIDE-model tonen aan dat het verwaarlozen van realistische reverberatiedynamiek leidt tot een onderschatting van detecteerbare bronnen. Het vereenvoudigde model neigt tot het overschatten van compressie en magnetische veldversterking, wat resulteert in excessieve synchrotronverliezen en een uitgeputte populatie van deeltjes met hoge energie, waardoor de voorspelde $\gamma$-stralingflux wordt verminderd.
• Evolutionaire Verdeling: Het merendeel van de synthetische populatie (87,5%) bevindt zich in de post-compressiefase. Bronnen in deze fase zijn echter doorgaans ouder, groter en intrinsiek zwakker, wat resulteert in lagere detecteerbaarheidsfracties in vergelijking met die in de vrije-expansie- of comprimerende fases.
• Spectrale Evolutie: Individuele casestudies illustreren dat de compressiefactor (CF) de sterkte en duur van luminositeitsversterking en spectrale verharding/verzachting tijdens reverberatie bepaalt. Bronnen met een hoge CF vertonen snelle, scherpe veranderingen in de spectrale energieverdeling (SED), terwijl bronnen met een lage CF mildere, vertraagde evolutie tonen.
• Correlaties: De studie vindt een correlatie tussen geïntegreerde 1–10 TeV luminositeit en pulsar-remmingsvermogen ($L_{1-10 \text{ TeV}} \propto \dot{E}^{1.15}$). De auteurs merken op dat het herwinnen van de ondiepere helling die in HGPS-gegevens wordt waargenomen ($L \propto \dot{E}^{0.59}$) vereist dat specifieke selectiecuts worden toegepast (bijvoorbeeld het uitsluiten van bronnen met laag remmingsvermogen of oudere leeftijden), wat suggereert dat survey-bias en instrumentele beperkingen waargenomen correlaties beïnvloeden.
• Bronuitbreiding: De gesimuleerde hoekuitbreidingen zijn over het algemeen consistent met HGPS-gegevens, hoewel het model een hoger aandeel bronnen voorspelt met stralen kleiner dan het minimum dat in HGPS is waargenomen. Dit wordt toegeschreven aan de aanname van sferische symmetrie en het verwaarlozen van menging en eigenbeweging van de pulsar.

Betekenis
Het artikel beweert dat realistische modellering van de reverberatiefase essentieel is voor het nauwkeurig voorspellen van de galactische TeV PWN-populatie. Door de complexe dynamiek van SNR-terugkeerschokinteracties te incorporeren, biedt de studie een betrouwbaardere theoretische basislijn voor het interpreteren van huidige broncatalogi en het anticiperen op de opbrengst van aankomende surveys. De resultaten bevestigen dat CTAO een dominante rol zal spelen in de komende TeV-populatietelling, en potentieel de aard van hoog-flux niet-geïdentificeerde bronnen in de HGPS zal oplossen. Bovendien suggereert het werk een natuurlijke evolutionaire link tussen middeloude PWNe en het ontstaan van uitgebreide TeV-halos, waarbij wordt opgemerkt dat oudere systemen die door toekomstige instrumenten worden gedetecteerd, waarschijnlijk corresponderen met halo-achtige morfologieën in plaats van klassieke compacte nevels. Deze studie dient als het eerste deel van een reeks, waarbij de statistische basislijn voor zichtbaarheid wordt vastgesteld, met een daaropvolgend artikel dat gepland is om zich te richten op "super-efficiëntie"-resultaten van reverberatie.