Impact of initial mass function on the chemical evolution of high-redshift galaxies

Probleemstelling

Recente waarnemingen van de James Webb Space Telescope (JWST) hebben een opmerkelijke invariantie blootgelegd in de relatie tussen sterrenmassa, metaliciteit en het sterformatie-tempo (SFR) tot een roodverschuiving van $z \sim 8$. Boven deze roodverschuiving lijkt deze relatie te breken, waarbij sterrenstelsels systematisch metalen-rijker zijn dan verwacht. Het onderliggende fysieke mechanisme dat deze correlaties drijft, is nog niet volledig begrepen. Een cruciale onzekerheid is de Initiële Massafunctie (IMF): de verdeling van sterrenmassa's bij hun geboorte. Traditionele modellen gaan vaak uit van een universele IMF (zoals de Kroupa-IMF) met een bovengrens van $\sim 100 M_\odot$. Echter, de overvloed aan UV-heldere sterrenstelsels en specifieke chemische handtekeningen (zoals verhoogde N- en C-abundanties) in het vroege heelal suggereren dat de IMF in het vroege heelal "top-heavy" kan zijn, met een aanzienlijk aandeel van zeer zware sterren ($>100 M_\odot$) die via paar-instabiliteit supernova's (PISNe) enorme hoeveelheden metalen produceren.

Methodologie

De auteurs gebruiken de semi-analytische code a-sloth om de chemische evolutie van sterrenstelsels op hoge roodverschuiving te modelleren. De studie omvat de volgende methodologische verbeteringen en instellingen:

• Simulatieopzet: De simulaties lopen op merger-bomen (samenvoegingsbomen) afgeleid van een hoog-resolutie kosmologische simulatie (Ishiyama et al. 2016) met een resolutie die de kleinste ster-vormende halo's ($M_h \sim 10^6 M_\odot$) oplost. Er zijn in totaal 165 simulaties uitgevoerd over een parameterrooster.

• Geüpdatete Voorschriften:

  • Sterformatie: De straal van het sterrenvormende gebied ($R_\star$) is aangepast van de halo-schaalstraal naar een schaal gebaseerd op de halve-massastraal van sterrenstelsels ($R_\star = 2.5 R_{50}$), wat beter overeenkomt met observaties en de spreiding in het hoofdreeks van sterformatie verbetert.
  • Sterfeedback en Metalen: Er is gebruikgemaakt van een uitgebreid raster van ster-evolutiemodellen (PARSEC v2.0) die de volledige metaliciteitsrange ($Z \sim 10^{-11} - 0.03$) en massabereik ($2 - 600 M_\odot$) bestrijken. Dit omvat bijdragen van sterwinden, kerninstortings-supernova's (CCSNe), pulsatoire paar-instabiliteit supernova's (PPISNe), paar-instabiliteit supernova's (PISNe) en Type Ia supernova's.
  • IMF Variatie: De IMF voor Populatie II-sterren wordt gemodelleerd als Kroupa-achtig, maar met een variabele bovengrens $m_{max}$ variërend van 100 tot 600 $M_\odot$. De slope van de IMF bij hoge massa's is vastgehouden op $\alpha_{II} = 2.3$.
  • Galactische Uitstromen: De efficiëntie van uitstromen ($\gamma_{out}$) wordt gemodelleerd als afhankelijk van de halo-massa en de supernova-energie, met vrije parameters voor de normalisatie ($M_{out0}$) en de helling ($\alpha_{out}$).

• Observatiebeperkingen: De modellen worden getoetst tegen drie hoofdobservaties:

  1. Cosmische Sterformatiegeschiedenis (CSFH): De dichtheid van de sterformatie (SFRD) en de totale sterrenmassadichtheid (CSMD) van $z \sim 4-16$.
  2. Chemische Evolutie: De invariantie van de Massa-Metaliciteit-SFR relatie (MZSFR) en de Fundamentele Metaliciteitsrelatie (FMR) tot $z \sim 8$, en de afwijking daarvan bij hogere $z$.
  3. Sterren-Halo Relatie: De relatie tussen sterrenmassa en halo-massa (SHMR).

Kernbijdragen

1. Integratie van Geavanceerde Stermodellen: Voor het eerst wordt een semi-analytisch model gekoppeld aan een volledig raster van ster-evolutietracks die metalenopbrengsten van sterwinden en PISNe voor een breed massabereik en metaliciteiten omvatten.
2. Kwantificering van de IMF-IMF Interactie: De studie isoleert systematisch de invloed van de IMF-bovengrens ($m_{max}$) en galactische uitstromen op de chemische samenstelling van het vroege heelal.
3. Beperking van Parameters: Door het combineren van CSFH en MZSFR-observaties, kunnen de auteurs de parameter ruimte van de IMF en uitstromen sterk beperken, wat leidt tot een specifiek "best-fit" scenario.

Resultaten

De analyse van de 165 modellen leidt tot de volgende cruciale bevindingen:

• Noodzaak van een Top-Heavy IMF: Alleen modellen met een IMF-bovengrens van $m_{max} \gtrsim 200 M_\odot$ kunnen gelijktijdelijk de waargenomen MZSFR-relatie en de kosmische sterformatiegeschiedenis bij $z \gtrsim 4$ reproduceren. Modellen met $m_{max} = 100 M_\odot$ falen in het produceren van voldoende metalen om de waarnemingen te verklaren, zelfs bij optimale uitstroomparameters.
• Rol van PISNe: De verhoogde metalenopbrengst komt voornamelijk van Pair-Instability Supernova's (PISNe) van zeer zware sterren ($\sim 140-260 M_\odot$). Deze explosies leveren een significante bijdrage aan de totale metalenbudget, wat nodig is om de metalenrijkdom van de waargenomen sterrenstelsels te verklaren.
• Best-Fit Model: Het model dat het beste overeenkomt met alle observaties heeft de volgende parameters:
  • $m_{max} = 250 M_\odot$
  • Uitstroomhelling $\alpha_{out} = 0$ (uitstroomefficiëntie is onafhankelijk van de halo-massa in de onderzochte schaal).
  • Uitstroomnormalisatie $M_{out0} = 3 \times 10^9 M_\odot$.
• Onzichtbare Bevolking: Hoewel de modellen de heldere, UV-rijke sterrenstelsels goed reproduceren, voorspellen ze een enorme populatie van zwakke, metalen-arme sterrenstelsels ($M_\star \lesssim 10^8 M_\odot$). Deze vormen bij $z \gtrsim 12$ meer dan 90% van de totale sterformatiedichtheid en spelen een cruciale rol bij de re-ionisatie en de productie van zware elementen.
• Transiënten en Zwaartekrachtsgolven: De hoge frequentie van zeer zware sterren impliceert een hoge productiesnelheid van PISNe (piek rond $z \sim 5$ met $\sim 10$ gebeurtenissen per jaar per stergraad) en een aanzienlijke bijdrage van metalen-arme sterren aan de vorming van binair zwarte gaten die samensmelten en zwaartekrachtsgolven uitzenden.

Betekenis en Conclusie

De studie concludeert dat de aannames van een universele IMF met een lage bovengrens onvoldoende zijn om de chemische evolutie van het vroege heelal te verklaren. De waarnemingen van JWST vereisen een IMF die "top-heavy" is met een bovengrens van minstens $200 M_\odot$.

Dit heeft diepgaande implicaties:

1. Chemische Evolutie: PISNe van Populatie II-sterren zijn essentieel voor de vroege verrijking van het interstellaire medium.
2. Re-ionisatie: De dominante bijdrage aan de sterformatie komt van onzichtbare, zwakke sterrenstelsels, wat de huidige interpretatie van re-ionisatie verandert.
3. Multi-messenger Astronomie: De hoge frequentie van PISNe en de overvloed aan metalen-arme binair zwarte gaten bieden nieuwe doelen voor toekomstige elektromagnetische waarnemingen (zoals met de Roman Space Telescope) en zwaartekrachtsgolf-detectoren (3e generatie).

De auteurs benadrukken dat hun resultaten modelafhankelijk zijn en dat verdere observaties van specifieke elementen (buiten zuurstof) nodig zijn om de degeneratie tussen IMF-slope, bovengrens en ster-evolutieparameters volledig op te lossen.