Low-ionization Metal Absorption at $0.7 \lesssim z \lesssim 2$ Confronting Cosmological Simulations with Observations

Probleemstelling

Laag-geïoniseerde metaalabsorptielijnen (zoals Mg I, Mg II en Fe II) in spectra van quasarlicht vormen een cruciale probe voor koud gas in en rondom sterrenstelsels (de circumgalactische medium, CGM) en het intergalactische medium (IGM). Hoewel waarnemingen van deze absorptiesystemen uitgebreid zijn, blijft het lastig om te bepalen hoe goed huidige modellen van galactische vorming deze waarnemingen kunnen verklaren. Bestaande statistieken lijden vaak onder heterogene spectrale resolutie en selectie-effecten, wat de vergelijking met hydrodynamische simulaties bemoeilijkt. Het artikel richt zich op de vraag of de IllustrisTNG simulaties de waargenomen eigenschappen van laag-geïoniseerd metaalgas correct reproduceren, en welke ionisatiemodellen nodig zijn om de waarnemingen te verklaren.

Methodologie

De auteurs gebruiken de hoog-resolutie IllustrisTNG50-1 simulatie en passen een grid-gebaseerd raamwerk toe om ionisatieverhoudingen te berekenen en kolomdichtheden te voorspellen.

1. Ionisatiemodellen: Omdat TNG de ionisatietoestanden van elementen niet expliciet volgt, hebben de auteurs twee ionisatievoorschriften geïmplementeerd:

   • Coll: Een puur collisioneel model (zonder stralingsveld).
   • Coll+UVB: Een model dat collisionele processen combineert met foto-ionisatie door een uniform kosmisch ultraviolet achtergrondveld (UVB).
     De ionisatieverhoudingen worden berekend onder de aanname van ionisatie-evenwicht, gebruikmakend van snelheidscoëfficiënten uit de CHIANTI-atomaire database en het UVB-model van Faucher-Giguère et al. (2009).

2. Berekening van Kolomdichtheden:

   • De ionische massa's in elke gascel worden berekend op basis van de totale gasmassa, de elementaire abundanties en de berekende ionisatiefracties.
   • Deze discrete gascellen worden geprojecteerd op een tweedimensionaal grid met behulp van een Smoothed Particle Hydrodynamics (SPH) kernel om kolomdichtheden ($N$) langs verschillende zichtlijnen te genereren.
   • Er worden twee statistieken gegenereerd: de Kansverdelingsfunctie (PDF) van de kolomdichtheid en de kosmische incidentie ($dN/dz$) van Mg II-systemen.

3. Observatievergelijking:

   • De simulaties worden vergeleken met data van de HIRES/UVES catalogus (Churchill et al. 2020), die hoge resolutie biedt voor kinematische analyse.
   • Er wordt ook een vergelijking gemaakt met de recente DESI (Dark Energy Spectroscopic Instrument) Mg II catalogus om de bruikbaarheid van deze grootschalige, maar lagere-resolutie dataset te beoordelen.

Belangrijkste Resultaten

• Kolomdichtheid PDF's:

  • De berekende PDF's voor Mg II en Fe II reproduceren de waargenomen verdelingen van HIRES/UVES breed over het bereik $10^{11.4} \lesssim N \lesssim 10^{16}$ cm$^{-2}$. Dit suggereert dat de simulatie de dominante fysieke drijfveren van laag-geïoniseerde absorptie correct vastlegt.
  • Voor Mg I toont TNG een onvoldoende steilere helling dan waargenomen. Dit wordt toegeschreven aan de beperkingen van het multiphase ISM-model in TNG, dat gas onder $T \lesssim 10^4$ K niet nauwkeurig oplost.
  • Er is weinig evolutie in de PDF's tussen $z \approx 0.7$ en $z \approx 2$, wat suggereert dat deze absorptiesystemen vroeg zijn gevormd en weinig evolueren.

• Kosmische Incidentie ($dN/dz$):

  • Het Coll+UVB model reproduceert de incidentie van Mg II-systemen nauwkeurig voor zwakke tot matig verzadigde absorbers (equivalente breedte $W_0 < 0.6$ Å). Dit ondersteunt het idee dat deze systemen corresponderen met laag-dicht, foto-geïoniseerd gas in de buitenste CGM.
  • Het puur collisionele model (Coll) overschat het aantal systemen, waarschijnlijk omdat het effect van foto-ionisatie ontbreekt.
  • Voor sterke absorbers ($W_0 > 1$ Å) onderschat TNG de incidentie en faalt het om de waargenomen stijging richting $z \sim 2$ te vangen. Dit kan wijzen op een onderschatting van metaalinhoud buiten sterrenstelsels of onnauwkeurigheden in de conversie van kolomdichtheid naar equivalente breedte voor turbulente systemen.

• Vergelijking met DESI:

  • Hoewel DESI een enorme steekproef biedt, leidt de lagere spectrale resolutie ($R \sim 4500$) tot lijnvervaging (blending). Hierdoor worden meerdere zwakke componenten vaak als één sterke component geteld.
  • Vergelijkingen met de hoge-resolutie HIRES/UVES data tonen aan dat DESI minder afzonderlijke componenten identificeert en mogelijk verkeerde identificaties maakt. Dit beperkt de directe bruikbaarheid van DESI-data voor gedetailleerde vergelijkingen met simulaties zonder aanvullende spectrale synthese.

Bijdragen en Significantie

1. Validatie van Ionisatiemodellen: Het artikel demonstreert dat de inclusie van een uniform UVB-veld essentieel is om de waargenomen populatie van Mg II-absorbers correct te reproduceren, vooral voor de zwakke systemen die het grootste deel van de populatie uitmaken.
2. Superioriteit van $dN/dz$: De auteurs concluderen dat de analyse van de kosmische incidentie ($dN/dz$) een effectievere diagnose is voor ionisatievoorschriften dan de kolomdichtheids-PDF alleen, omdat de incidentie gevoeliger is voor de evolutie en de totale hoeveelheid gas.
3. Efficiënte Modellering: De gepresenteerde grid-gebaseerde methode biedt een rekenkundig efficiënte manier om ionisatie toe te passen op grote cosmologische simulaties, waardoor snelle vergelijkingen met waarnemingen mogelijk zijn zonder volledige spectrale synthese voor elke zichtlijn.
4. Toekomstperspectief: De studie benadrukt dat toekomstige vergelijkingen met grootschalige surveys zoals DESI en Euclid (die ook gevoelig is voor Mg II maar lage resolutie heeft) vereisen dat de simulaties worden aangepast aan de resolutie van de waarnemingen of dat geavanceerde spectrale synthese wordt toegepast om vervagingseffecten te corrigeren.

Samenvattend biedt dit werk een robuust kader om de fysica van het koude gas rondom sterrenstelsels te testen en identificeert het de beperkingen van huidige simulaties bij het reproduceren van de sterkste absorptiesystemen en de evolutie daarvan naar hogere roodverschuivingen.