MIGHTEE-HI: The direct detection of neutral hydrogen in galaxies at $z>0.25$

Probleemstelling

Neutraal waterstof (HI) vormt de oorspronkelijke gasvoorraad waaruit moleculair gas en uiteindelijk sterren in sterrenstelsels ontstaan. Het begrijpen van de relatie tussen HI en sterrenstelsels over kosmische tijdschalen is cruciaal voor het model van galactische evolutie. Echter, directe detectie van HI-emissie (de 21-cm lijn) is tot nu toe beperkt gebleven tot het zeer lokale universum ($z \lesssim 0.1$) vanwege de extreme zwakte van de lijn.

• Uitdagingen: De overgangswaarschijnlijkheid is vele ordes van grootte lager dan die van Lyman-$\alpha$, wat leidt tot zeer zwakke signalen.
• Huidige status: Eerdere pogingen om HI bij hogere roodverschuivingen ($z > 0.25$) te detecteren, maakten gebruik van ongeoriënteerde surveys (beperkt tot lage $z$), spectrale stapeling (statistisch gemiddelde, geen individuele detecties) of enkele-schotel telescopen zoals FAST (wat kan leiden tot verwarring door slechte hoekresolutie). Er was een gebrek aan een steekproef van individueel gedetecteerde HI-sterrenstelsels bij $z > 0.25$ met behulp van een interferometer.

Methodologie

Het onderzoek maakt gebruik van data van de MeerKAT International Tiered Extragalactic Exploration (MIGHTEE) survey, specifiek de Data Release 1 (DR1) spectrale kubus over het COSMOS-veld, gecombineerd met optische spectroscopische data.

1. Data Bronnen:

   • MIGHTEE-Hi: Een spectrale kubus met een frequentiebereik van 960–1150 MHz, een kanalenbreedte van 104.5 kHz en een hoekresolutie van $\sim15$ boogseconden (FWHM). Dit dekt een roodverschuiving van $z \approx 0.23$ tot $0.5$.
   • DESI (Dark Energy Spectroscopic Instrument): De public catalogus van DESI (Early Data Release) werd gebruikt om doelwitten te selecteren. De auteurs focusten op emissielijn-sterrenstelsels (actief stervormend) binnen het roodverschuivingsbereik $0.25 < z < 0.5$.

2. Zoekstrategie (Targeted Search):

   • In plaats van een "blind" zoeken, werden de posities en roodverschuivingen van 915 optisch geselecteerde sterrenstelsels gebruikt om specifieke "cubelets" ($80 \times 80$ boogsec) uit de spectrale data te extraheren.
   • Visuele inspectie: De auteurs collapseerden de 3D-data naar 2D-kaarten en analyseerden 1D-spectra om emissielijnen te identificeren. Criteria voor een kandidaat waren:
     • Een prominent signaal binnen 4 boogsec van de optische positie.
     • Signaalruisverhouding (SNR) significant boven de ruis in het spectrale domein.
     • Een lijnbreedte van minimaal 5 spectrale kanalen ($\sim150$ km/s).

3. Validatie en Ruisanalyse:

   • Er werd geen a-priori statistische veronderstelling van Gaussische ruis gemaakt. De ruis werd gemeten door $\sim500$ vergelijkbare regio's rondom de bron te analyseren (zowel ruimtelijk als spectraal).
   • De 16e en 84e percentielen van de flux-dichtheid werden gebruikt om de ruis te bepalen.
   • Controle-experimenten: Om "confirmation bias" (bevestigingsbias) te voorkomen, werd dezelfde zoekmethode toegepast op posities die 50 boogminuten van de optische bronnen verwijderd waren. Geen enkele van deze offset-posities leverde een significante detectie op ($SNR > 4$).
   • Er werd gecontroleerd op radio-continuüm emissie en artefacten van continuüm-subtractie.

4. Bepaling van Eigenschappen:

   • HI-massa: Berekend uit de geïntegreerde flux, aannemende dat de lijn optisch dun is.
   • Sterrenstelsel-eigenschappen: Gebruik van de Deep Extragalactic VIsible Legacy Survey (DEVILS) voor multi-golflengte data (UV tot Far-IR). De Bagpipes-code werd gebruikt voor Spectral Energy Distribution (SED) fitting om sterrenmassa ($M_*$) en sterformatiesnelheid (SFR) te bepalen.
   • Baryonische Tully-Fisher Relatie (bTFr): Gebruik van de lijnbreedte ($W_{50}$), gecorrigeerd voor instrumentele verbreding en de inclinatie van het sterrenstelsel (bepaald uit g-band imaging), als proxy voor rotatiesnelheid.

Belangrijkste Bijdragen

• Eerste Interferometrische Detecties: Dit artikel presenteert de eerste steekproef van individuele HI-detecties bij $z > 0.25$ met een interferometer (MeerKAT).
• Targeted Approach: Het succesvol toepassen van een "targeted search" strategie, waarbij optische spectroscopische roodverschuivingen worden gebruikt om de detectieproblematiek van zwakke lijnen bij hoge $z$ te overwinnen.
• bTFr bij Hoge Roodverschuiving: De eerste directe meting van de baryonische Tully-Fisher relatie voor HI-sterrenstelsels bij $z > 0.25$.

Resultaten

1. Detecties: Er werden 11 sterrenstelsels gedetecteerd met een hoge betrouwbaarheid ($SNR > 4$) in het bereik $0.25 \lesssim z \lesssim 0.39$. De hoogste detectie is bij $z = 0.3841$.

   • Opmerking: Er werden geen detecties gevonden tussen $0.39 < z < 0.5$, wat voornamelijk wordt toegeschreven aan de roodverschuivingsverdeling van de oudersteekproef (DESI) en de diepte van de data, niet aan een fundamenteel gebrek aan HI.

2. HI-Massa's: De gedetecteerde sterrenstelsels hebben aanzienlijk grotere HI-massa's ($\log(M_{HI}/M_\odot) \approx 10.2 - 10.9$) dan hun tegenhangers bij lage roodverschuiving voor een gegeven sterrenmassa.

   • Interpretatie: Dit wordt niet veroorzaakt door evolutie van de HI-massa-functie, maar door het feit dat bij hogere roodverschuivingen een veel groter kosmologisch volume wordt onderzocht. Hierdoor worden de zeldzame, HI-rijke sterrenstelsels aan het "heldere uiteinde" van de verdeling gevonden.

3. Stellare Eigenschappen:

   • De sterrenstelsels liggen over het algemeen op of boven de "main sequence" voor stervorming bij $z \sim 0.35$.
   • De gas-uitputtingstijden ($t_{dep} = M_{HI}/SFR$) zijn gemiddeld lager dan bij lokale sterrenstelsels (gemiddeld 3.1 Gyr vs 9.6 Gyr), wat wijst op een efficiëntere omzetting van gas naar sterren in deze specifieke steekproef.

4. Baryonische Tully-Fisher Relatie (bTFr):

   • De meeste sterrenstelsels vallen binnen de spreiding van de lokale bTFr.
   • Flattening: Er is een voorlopig bewijs voor een "afvlakking" van de relatie bij hoge massa's/hoge rotatiesnelheden. Dit kan wijzen op echte evolutie (in lijn met hydrodynamische simulaties) of erop dat de rotatiesnelheid ($W_{50}$) de maximale rotatiesnelheid ($V_{max}$) meet in plaats van de vlakke snelheid ($V_{flat}$) bij zware systemen.
   • Moleculair Gas: Als men aanneemt dat het moleculaire gas een significant deel uitmaakt van de baryonische massa (vooral bij zware sterrenstelsels), dan vallen de punten weer beter op de lokale relatie.

Betekenis en Conclusie

Dit onderzoek markeert een mijlpaal in de radioastronomie door aan te tonen dat directe detectie van neutraal waterstof in individuele sterrenstelsels mogelijk is over een kosmische tijdschaal van meer dan 5 miljard jaar ($z > 0.25$) met de MeerKAT-telescoop.

• Toekomstperspectief: De studie legt de basis voor toekomstige onderzoeken, zowel via gerichte zoekopdrachten (zoals hier gedaan) als via pure HI-selectie (blinde zoekopdrachten) in de volledige MIGHTEE- en LADUMA-survey-gebieden.
• Evolutie: Hoewel er geen sterke bewijzen zijn voor een fundamentele evolutie van de bTFr tot $z \sim 0.35$, suggereert de data dat er subtiele veranderingen zijn in de verhouding tussen baryonische massa en dynamische massa, of in de verdeling van moleculair versus neutraal gas bij zware sterrenstelsels.
• Voorbereiding op SKA: De resultaten onderstrepen het potentieel van huidige interferometers en bereiden de weg voor de Square Kilometre Array (SKA), die deze metingen tot nog hogere roodverschuivingen zal kunnen uitbreiden.