Short gamma-ray burst progenitors have short delay times

Probleemstelling
Korte gamma-ray bursts (SGRBs) worden breed beschouwd als afkomstig van de coalescentie van compacte binaire systemen, specif으로 BNS (binair neutronenster) of NSBH (neutronenster-zwart gat) fusies. Een cruciale parameter voor het begrijpen van de evolutie van de SGRB-populatie is de delay time distribution (DTD), die het tijdsinterval beschrijft tussen de vorming van de progenitor massieve ster binaire en de uiteindelijke fusie. Eerdere studies hebben tegenstrijdige resultaten opgeleverd met betrekking tot de vorm van de DTD en de grootte van de delay tijden. Sommige analyses suggereren lange delay tijden op de orde van gigajaren (Gyr), wat wijst op een significante evolutie van de kosmische reactiedichtheid die afwijkt van de kosmische stervormingsgeschiedenis (CSFH). Daarentegen suggereren andere theoretische overwegingen en lokale populatiestudies het bestaan van een "snelle" fusie-kanaal met vertragingen van slechts enkele honderden megajaren (Myr). De auteurs merken op dat discrepanties in eerdere bevindingen kunnen voortkomen uit de incorrecte behandeling van selectie-effecten en de specifieke luminositeitsfunctie-modellen die gehanteerd worden.

Methodologie
De auteurs maken gebruik van een hiërarchisch Bayesiaans inferentiekader om een grote steekproef van SGRBs te analyseren. De methodologie omvat de volgende kerncomponenten:

• Statistisch Kader: Volgens de aanpak van Mandel et al. (2019) en Salafia et al. (2023, S23), modelleren de auteurs de populatie met hyperparameters ($\lambda'_{pop}$) die de intrinsieke distributies van bronparameters definiëren. De posterieure waarschijnlijkheidsverdeling wordt berekend door de likelihood van de geobserveerde data te combineren met een selectiefunctie ($P_{det}$) die rekening houdt met detectiedrempels en instrumentele biases. Het kader wordt bijgewerkt door het geobserveerde aantal gebeurtenissen toe te voegen via een Poissoniaanse term om de lokale reactiedichtheid ($R_0$) te beperken.
• DTD Modellen: Twee verschillende modellen voor de delay time distribution worden getest:
  1. Een machtswetverdeling met een minimale delay tijd afkapwaarde ($P(\tau_d) \propto \tau_d^{-\alpha_\tau}$ voor $\tau_d \ge \tau_{min}^d$).
  2. Een log-normale verdeling gedefinieerd door een mediaan ($\mu_\tau$) en dispersie ($\sigma_\tau$).
• Luminositeitsfunctie Modellen: Om de robuustheid van de resultaten te testen tegen aannames over intrinsieke helderheid, worden twee modellen gebruikt:
  1. Empirische Luminositeitsfunctie (ELF): Een dubbel gebroken machtswetfunctie die eerdere empirische modellen uitbreidt naar lagere luminositeiten, waarbij rekening wordt gehouden met beperkingen uit GRB 170817A.
  2. Quasi-Universele Gestructureerde Jet (QUSJ): Een fysisch gemotiveerd model waarbij luminositeit en piekfotonenergie gecorreleerd zijn via de kijkhoek ten opzichte van een gestructureerde jet-kern.
• Data Samples: Twee primaire steekproeven worden geconstrueerd uit Fermi/GBM en Swift/BAT data:
  1. Een observer-frame sample (210 gebeurtenissen) geselecteerd op basis van flux volledigheid en spectrale piekfotonenergie-snedes.
  2. Een rest-frame sample (18 gebeurtenissen, 16 met roodverschuivingen) geselecteerd om de roodverschuivings volledigheid te maximaliseren.
• Reproductie van Vorig Werk: Om de bron van de discrepanties met de literatuur te isoleren, reproduceren de auteurs de analyse van Wanderman & Piran (2015, W15) met dezelfde steekproef en selectiesneden, maar passen zij hun eigen hiërarchische Bayesiaanse inferentiemethode toe (verwezen als W15*).

Belangrijkste Resultaten
De analyse levert consistente resultaten op voor zowel de luminositeitsfunctie modellen (ELF en QUSJ) als beide DTD parametrisaties (machtswet en log-normaal):

• Korte Delay Tijden: De studie vindt gemiddelde delay tijden ($\langle \tau_d \rangle$) in de reeks van $10 \lesssim \langle \tau_d \rangle / \text{Myr} \lesssim 900$. Specifiek worden mediane waarden gevonden van ongeveer $77$ Myr (ELF, machtswet) en $105$ Myr (QUSJ, machtswet), met 90% geloofwaardigheidsintervallen die oplopen tot ongeveer 1 Gyr.
• Minimale Delay: De minimale delay tijd wordt beperkt tot $\tau_{d,min} \lesssim 370$ Myr (90% geloofwaardigheidslimiet).
• Machtswet Index: Voor de machtswet DTD is de index gevonden om steil te zijn, met $\alpha_\tau > 1.5$, wat duidt op een voorkeur voor kortere fusietijden.
• Reactiedichtheid Evolutie: De geïnferentieerde kosmische reactiedichtheid evolutie volgt nauwgezet de kosmische stervormingsgeschiedenis (CSFH), in contrast met modellen die een significante vertraging of evolutie voorspellen door lange delay tijden.
• Lokale Reactiedichtheid: De geïnferieerde lokale reactiedichtheden ($R_0$) voor de ELF (machtswet) en QUSJ (log-normaal) modellen zijn grotendeels consistent met de BNS-fusie reactiedichtheidsbeperkingen afgeleid van zwaartekrachtgolf-observaties (GWTC-4).
• Discrepantie met W15: Wanneer de auteurs de W15 analyse (W15*) reproduceren met hun eigen inferentiekader, verkrijgen zij significant langere delay tijden ($\langle \tau_d \rangle \approx 3-4$ Gyr). Dit bevestigt dat het verschil in resultaten niet te wijten is aan de data zelf, maar aan de methodologie.

Betekenis en Claims
Het artikel claimt dat de primaire oorzaak van de langere delay tijden die in eerdere studies werden geïnferieerd, een incorrecte behandeling van selectie-effecten is. Door een rigoureus hiërarchisch Bayesiaans kader te gebruiken dat expliciet de selectiefunctie modelleert en een breder scala aan luminositeitsfuncties incorporeert, demonstreren de auteurs dat de SGRB-populatie wordt gedomineerd door progenitors met korte delay tijden.

De auteurs concluderen dat SGRBs primair geassocieerd zijn met een "snel" fusie-kanaal van binaire compacte objecten, met gemiddelde vertragingen van minder dan een gigajaar. Deze bevinding brengt de evolutie van de SGRB-populatie dichter bij de kosmische stervormingsgeschiedenis en ondersteunt het bestaan van een dominante snelle fusie BNS-populatie, wat consistent is met sommige lokale Galactische BNS-studies en chemische verrijkingsbeperkingen. De studie benadrukt dat robuuste inferentie van populatieparameters zorgvuldige modellering van detectiebiases en de intrinsieke luminositeitsverdeling vereist.