Predictions from $s$-process AGB models of the isotopic variations of zirconium and neodymium for comparison to bulk meteorites

Probleemstelling

Meteoritische bulk-data tonen isotopische variaties van oorsprong uit het slow-neutron-capture proces (s-process) aan voor diverse elementen zwaarder dan ijzer. Een opvallend kenmerk is dat de lichtere s-process elementen (zoals Zirkonium, Zr, en Molybdeen, Mo) grotere anomalieën vertonen dan de zwaardere s-process elementen (zoals Neodymium, Nd, en Wolfraam, W).
De huidige interpretatie van deze data is vaak gebaseerd op modellen van Asymptotic Giant Branch (AGB)-sterren die de gemiddelde zonnelijke s-process verdeling voorspellen over de levensduur van de Melkweg. Dit is echter onnauwkeurig voor de verklaring van isotopische variabiliteit in het vroege zonnestelsel, omdat de specifieke AGB-sterren die de stardustkorrels leverden, niet noodzakelijk dezelfde zijn als die welke de totale zonnelijke abundantie hebben gevormd.
Het centrale vraagstuk is of de relatieve grootte van de isotopische variabiliteit tussen Zr (eerste s-process piek, neutronmagisch getal 50) en Nd (tweede s-process piek, neutronmagisch getal 82) kan worden verklaard door s-process nucleosynthese in AGB-sterren, en zo ja, onder welke omstandigheden (bijv. metaalrijkheid).

Methodologie

De auteurs hebben een uitgebreide reeks van 80 nieuwe AGB-stermodellen ontwikkeld en geanalyseerd om de isotopische verhoudingen van Zr en Nd te voorspellen.

1. Stermodellen:

   • Massa: 2 tot 3,5 zonsmassa's ($M_\odot$).
   • Metaalrijkheid (Z): Variërend van zonne-metaalrijkheid ($Z=0.014$) tot tweemaal zonne-metaalrijkheid ($Z=0.03$), met aanvullende modellen voor lagere en hogere metaalrijkheid.
   • Convectieve Overshoot: Variatie in de parameter $N_{ov}$ (0 tot 4 schaalhoogtes) om de efficiëntie van de 'third dredge-up' (TDU) te testen.
   • Deel-mengzone ($M_{mixed}$): Een geparametriseerde zone aan de basis van de convectieve mantel waar protonen worden gemengd met $^{12}$C om de neutronbron $^{13}$C te vormen. De massa hiervan varieerde tussen $5 \times 10^{-4}$ en $2 \times 10^{-3} M_\odot$.
   • Nucleosynthese: Berekeningen uitgevoerd met de Monash nucleosynthese-code, inclusief een netwerk van 328 isotopen en 2.351 reacties. De neutronbronnen zijn $^{13}$C($\alpha$,n)$^{16}$O en $^{22}$Ne($\alpha$,n)$^{25}$Mg.

2. Data-verwerking en Vergelijking:

   • De oppervlakte-abundanties van de modellen werden verdund met een factor $10^4$ (representatief voor de abundantie van stardust in meteorieten) en gemengd met standaard zonnelijke abundanties.
   • De isotopische verhoudingen werden genormaliseerd om massafhankelijke fractionering uit te sluiten.
   • De afwijkingen werden uitgedrukt in $\varepsilon$-waarden (afwijkingen per 10.000 ten opzichte van de zon): $\varepsilon^{96}\text{Zr}$ en $\varepsilon^{148}\text{Nd}$.
   • De kernvergelijking was de verhouding $\varepsilon^{96}\text{Zr} / \varepsilon^{148}\text{Nd}$. Een hoge verhouding impliceert een sterke s-process signatuur in de eerste piek (Zr) ten opzichte van de tweede (Nd), wat correleert met een lager totaal neutronenflux (neutron exposure).

Belangrijkste Bijdragen

• Eerste directe vergelijking: Dit is de eerste studie die de relatieve variabiliteit van Zr en Nd in bulk meteorieten direct vergelijkt met gedetailleerde s-process modellen van AGB-sterren, in plaats van te vertrouwen op gemiddelde galactische chemische evolutiemodellen.
• Rol van Metaalrijkheid: Het onderzoek toont aan dat modellen met super-zonne-metaalrijkheid (hoger dan $Z_\odot$) de beste overeenkomst vertonen met de waargenomen meteoritische data.
• Mechanisme van Overshoot en $M_{mixed}$: De resultaten benadrukken dat een sterkere convectieve overshoot (die leidt tot diepere TDU) en/of een kleinere massa van de $^{13}$C-pockets essentieel zijn om de waargenomen verhoudingen te reproduceren.

Resultaten

1. Verhouding $\varepsilon^{96}\text{Zr} / \varepsilon^{148}\text{Nd}$:

   • Waarnemingen in bulk chondrieten tonen een verhouding van ongeveer 10 (Zr-variabiliteit is ~10x groter dan Nd-variabiliteit).
   • Modellen met zonne-metaalrijkheid ($Z=0.014$) voorspellen over het algemeen lagere verhoudingen (rond 1-3).
   • Modellen met super-zonne-metaalrijkheid ($Z=0.02$ en $0.03$), gecombineerd met hoge $N_{ov}$ en lage $M_{mixed}$, produceren verhoudingen die de waargenomen waarde van ~10 kunnen bereiken (tot wel 13 in extreme gevallen).
   • Fysische verklaring: Bij hogere metaalrijkheid zijn er meer ijzer-kernen (seeds) beschikbaar, maar de hoeveelheid vrije neutronen (voornamelijk uit $^{13}$C, dat onafhankelijk is van metaalrijkheid) blijft gelijk. Hierdoor is de neutronenflux per ijzerkern lager, waardoor het s-process minder ver doorgaat naar de zwaardere elementen (tweede piek, Nd) en sterker blijft bij de lichtere elementen (eerste piek, Zr).

2. Invloed van Parameters:

   • Convectieve Overshoot ($N_{ov}$): Hogere overshoot zorgt voor TDU bij lagere kerntemperaturen. Dit kan leiden tot het inslikken van overgebleven $^{13}$C in de volgende thermische puls, wat de neutronendichtheid tijdelijk verhoogt en de verhouding Zr/Nd beïnvloedt.
   • $M_{mixed}$: Een kleinere $M_{mixed}$ (kleinere $^{13}$C-pocket) resulteert in een hogere verhouding, omdat er minder neutronen beschikbaar zijn om de zwaardere elementen te bereiken.

3. Validatie met andere data:

   • De beste modellen (hoge Z, hoge $N_{ov}$, lage $M_{mixed}$) komen ook overeen met de isotopische samenstelling van grote SiC-stardustkorrels uit AGB-sterren.
   • De voorspellingen zijn consistent met waarnemingen van Ba-sterren (sterren die s-process elementen hebben opgenomen van een AGB-compaan) in het zonnepaneel met hoge metaalrijkheid.

Significantie en Conclusie

De studie concludeert dat nucleosynthese in AGB-sterren met een metaalrijkheid hoger dan die van de zon de waargenomen verschillen in de grootte van de isotopische variaties tussen Zr en Nd in bulk meteorieten kan verklaren.

• Geen nieuwe dragers nodig: Er is geen noodzaak om verschillende stofdragers of complexe chemische fractionering in de protoplanetaire schijf aan te roepen om de data te verklaren; het verschil is inherent aan de nucleosynthese in de bronsterren.
• Oorsprong: De AGB-sterren die deze variaties produceerden, behoorden waarschijnlijk tot een oude populatie ($\sim$5-9 Ga) van sterren met hoge metaalrijkheid die al aanwezig was in het galactische zonnepaneel rond de tijd van de vorming van het zonnestelsel (4,6 Ga geleden).
• Toekomstige werk: Verdere verfijning van de neutronenbron-reacties (via ondergrondse experimenten zoals LUNA en JUNA), nieuwe kruissecties voor Nd-isotopen, en gedetailleerde vergelijkingen met individuele stardustkorrels en andere elementen (zoals Mo, Ru, W) zijn nodig om dit beeld te bevestigen.

Kortom, dit onderzoek biedt een robuust nucleosynthetisch kader voor de isotopische anomalieën in meteorieten, waarbij de metaalrijkheid van de bronsterren de sleutelparameter is.