Modelling DSA, FAST and CRAFT surveys in a z-DM analysis and constraining a minimum FRB energy

Sintesi Tecnica: Modellazione dei Sondaggi DSA, FAST e CRAFT in un'Analisi z-DM e Vincolamento dell'Energia Minima degli FRB

Enunciato del Problema
La scienza dei Fast Radio Burst (FRB) mira a chiarire le origini di questi burst e a utilizzarli per studi cosmologici, in particolare tracciando la distribuzione dei barioni ionizzati e vincolando la costante di Hubble ($H_0$). Le analisi precedenti che hanno utilizzato il framework $z$–DM (redshift–misura di dispersione), come quelle di James, Ghosh, et al. (2022) e Baptista et al. (2023), hanno modellato con successo i bias strumentali e adattato i parametri della popolazione utilizzando dati dai sondaggi Parkes e ASKAP (CRAFT). Tuttavia, questi studi erano limitati dalle statistiche di piccoli numeri e dall'esclusione di sondaggi più recenti e altamente sensibili come il Deep Synoptic Array (DSA) e il Five-hundred-meter Aperture Spherical radio Telescope (FAST). Inoltre, i modelli precedenti spesso trattavano i contributi galattici alla misura di dispersione (DM) come valori fissi, potenzialmente sovrastimando la precisione dei parametri cosmologici e della popolazione derivati.

Metodologia
Questo lavoro estende il codice $z$DM per incorporare dati da 119 FRB, inclusi 29 con redshift associati, modellando i bias osservativi dei sondaggi DSA e FAST insieme ai dati aggiornati CRAFT/ICS.

• Modellazione dei Sondaggi:
  • FAST: Modellato come un telescopio a singola antenna sferico con un ricevitore a 19 fasci. Il pattern del fascio è stato approssimato utilizzando 19 fasci gaussiani con Half-Power Beamwidths (HPBW) derivati da Jiang et al. (2020). L'analisi ha incluso 9 FRB non localizzati rilevati nel Commensal Radio Astronomy FAST Survey (CRAFTS) e nel Galactic Plane Pulsar Snapshot (GPPS).
  • DSA: Modellato utilizzando dati di commissioning da DSA-110. Il pattern del fascio è stato approssimato come un singolo gaussiano con un Full-Width Half-Maximum (FWHM) di $2.6^\circ$. Per evitare un bias di rilevamento contro gli FRB ad alto redshift (che appaiono in galassie ospiti più deboli), i redshift sono stati utilizzati solo per gli FRB in cui la DM extragalattica ($DMEG$) era inferiore a una soglia ($183 \text{ pc cm}^{-3}$), garantendo un campione completo di ospiti localizzati.
  • CRAFT: Ha incluso FRB aggiuntivi rilevati fino alla fine del 2023, escludendo FRB 20171216A poiché il suo Rapporto Segnale-Rumore (SNR) è sceso sotto la soglia aggiornata.
• Miglioramenti Statistici:
  • Incertezza nella DM Galattica: A differenza degli studi precedenti che utilizzavano valori fissi per il Mezzo Interstellare ($DM_{ISM}$) e l'alone della Via Lattea ($DM_{halo}$), questo lavoro modella $DM_{MW}$ come una distribuzione di probabilità. $DM_{ISM}$ è modellato come una distribuzione normale troncata basata sul modello NE2001 con una deviazione standard del 50%, e $DM_{halo}$ è modellato come una distribuzione normale con una media di $50 \text{ pc cm}^{-3}$ e $\sigma = 15 \text{ pc cm}^{-3}$.
  • Implementazione MCMC: A causa dell'aumento del carico computazionale derivante da sondaggi e parametri aggiuntivi, l'analisi è passata da un calcolo della verosimiglianza basato su griglia a un campionatore Markov Chain Monte Carlo (MCMC) (EMCEE) eseguito sul supercomputer OzSTAR.
  • Parametro Libero $E_{min}$: È stato introdotto un cutoff rigido per l'energia minima degli FRB ($E_{min}$) come parametro libero, permettendo all'analisi di sondare il turnover a bassa energia della funzione di luminosità.

Risultati Chiave

• Vincolo sull'Energia Minima: L'analisi produce un vincolo sull'energia minima degli FRB di $\log E_{min} (\text{erg}) = 39.47^{+0.54}_{-1.28}$. Questo valore è significativamente più alto delle energie rilevate nei repeater forti (ad es. FRB 20121102A e FRB 20201124A), che sono spesso osservati fino a $\sim 10^{37}$ erg.
• Parametri della Popolazione: L'indice spettrale $\alpha$ preferisce un valore piatto di $0.11^{+0.66}_{-0.60}$, differendo dalla pendenza negativa preferita nelle precedenti analisi basate solo su ASKAP. Tuttavia, quando $H_0$ è consentito di variare liberamente, la preferenza torna verso una pendenza negativa ($\alpha \approx -0.92$).
• Propagazione dell'Incertezza: L'inclusione dell'incertezza in $DM_{MW}$ porta a incertezze più grandi nei parametri del modello derivati rispetto alle stime precedenti, nonostante l'aumento delle dimensioni del dataset.
• Previsioni per Rilevamenti Futuri:
  • FAST: Previsto rilevare il 25–41% dei suoi FRB a $z \gtrsim 2$ e il 6–20% a $z \gtrsim 3$. La distribuzione raggiunge il picco a $z \approx 1$, con il 73–77% dei rilevamenti oltre $z \gtrsim 1$.
  • DSA: Previsto rilevare il 2–12% degli FRB a $z \gtrsim 1$.
  • Limiti dei Sondaggi: La DM di ricerca massima attuale per FAST (3700–5000 $\text{ pc cm}^{-3}$) è prevista per perdere il 7–18% dei rilevamenti possibili, mentre il limite di DSA di $1500 \text{ pc cm}^{-3}$ è previsto per perdere il 3–8%.

Significato e Affermazioni
Gli autori affermano che questo lavoro fornisce il primo vincolo su $E_{min}$ derivato da un'analisi $z$–DM. L'alto valore di $E_{min}$ suggerisce due possibilità:

1. Esiste un turnover a bassa energia (appiattimento) nella funzione di luminosità della popolazione generale degli FRB, il che significa che i sondaggi non sono sufficientemente sensibili da rilevare la coda a bassa energia.
2. I repeater forti possiedono una funzione di luminosità fondamentalmente diversa rispetto ai rilevamenti di burst singoli (o burst singoli apparenti).

Il documento sottolinea che, sebbene i nuovi dati espandano lo spazio dei parametri accessibile verso redshift più alti, l'incapacità di localizzare in modo robusto gli FRB ad alto $z$ con FAST limita attualmente la possibilità di ottenere redshift associati. Di conseguenza, gli autori affermano che misurazioni precise di $DM_{ISM}$ saranno critiche per i futuri vincoli sui parametri cosmologici come $H_0$. Il lavoro funge da raffinamento metodologico, dimostrando che incorporare le incertezze nei contributi galattici è necessario per evitare di distorcere artificialmente il potere vincolante degli FRB a bassa DM.