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

Resumo Técnico: Modelagem das Levantamentos DSA, FAST e CRAFT em uma Análise z-DM e Restrição de uma Energia Mínima de FRB

Declaração do Problema
A ciência de Explosões de Rádio Rápidas (FRB) visa elucidar as origens desses surtos e utilizá-los para estudos cosmológicos, especificamente traçando a distribuição de bárions ionizados e restringindo a constante de Hubble ($H_0$). Análises anteriores utilizando o framework $z$–DM (desvio para o vermelho–medida de dispersão), como as de James, Ghosh, et al. (2022) e Baptista et al. (2023), modelaram com sucesso vieses instrumentais e ajustaram parâmetros populacionais usando dados dos levantamentos Parkes e ASKAP (CRAFT). No entanto, esses estudos foram limitados por estatísticas de pequeno número e pela exclusão de levantamentos mais novos e altamente sensíveis, como o Deep Synoptic Array (DSA) e o Telescópio de Radioastronomia Esférico de Quinhentos Metros (FAST). Além disso, modelos anteriores frequentemente tratavam as contribuições da medida de dispersão galáctica (DM) como valores fixos, potencialmente superestimando a precisão dos parâmetros cosmológicos e populacionais derivados.

Metodologia
Este trabalho estende o código $z$DM para incorporar dados de 119 FRBs, incluindo 29 com desvios para o vermelho associados, modelando os vieses observacionais dos levantamentos DSA e FAST juntamente com dados atualizados do CRAFT/ICS.

• Modelagem de Levantamentos:
  • FAST: Modelado como um telescópio de disco único esférico com um receptor de 19 feixes. O padrão de feixe foi aproximado usando 19 feixes gaussianos com Larguras de Feixe a Meia Potência (HPBW) derivadas de Jiang et al. (2020). A análise incluiu 9 FRBs não localizados detectados no Levantamento Commensal de Radioastronomia FAST (CRAFTS) e no Instantâneo de Púlsares do Plano Galáctico (GPPS).
  • DSA: Modelado usando dados de comissionamento do DSA-110. O padrão de feixe foi aproximado como um único gaussiano com uma Largura a Meia Altura Total (FWHM) de $2.6^\circ$. Para evitar viés de detecção contra FRBs de alto desvio para o vermelho (que aparecem em galáxias hospedeiras mais fracas), os desvios para o vermelho foram utilizados apenas para FRBs onde a DM extragaláctica ($DMEG$) estava abaixo de um limiar ($183 \text{ pc cm}^{-3}$), garantindo uma amostra completa de hospedeiros localizados.
  • CRAFT: Incluiu FRBs adicionais detectados até o final de 2023, excluindo o FRB 20171216A devido ao seu Relação Sinal-Ruído (SNR) cair abaixo do limiar atualizado.
• Melhorias Estatísticas:
  • Incerteza na DM Galáctica: Ao contrário de estudos anteriores que usaram valores fixos para o Meio Interestelar ($DM_{ISM}$) e o halo da Via Láctea ($DM_{halo}$), este trabalho modela $DM_{MW}$ como uma distribuição de probabilidade. $DM_{ISM}$ é modelado como uma distribuição normal truncada baseada no modelo NE2001 com um desvio padrão de 50%, e $DM_{halo}$ é modelado como uma distribuição normal com uma média de $50 \text{ pc cm}^{-3}$ e $\sigma = 15 \text{ pc cm}^{-3}$.
  • Implementação MCMC: Devido à carga computacional aumentada por levantamentos e parâmetros adicionais, a análise mudou de um cálculo de verossimilhança baseado em grade para um amostrador de Cadeia de Markov Monte Carlo (MCMC) (EMCEE) executado no supercomputador OzSTAR.
  • Parâmetro Livre $E_{min}$: Um corte rígido para a energia mínima de FRB ($E_{min}$) foi introduzido como um parâmetro livre, permitindo que a análise investigue a virada de baixa energia da função de luminosidade.

Principais Resultados

• Restrição de Energia Mínima: A análise produz uma restrição na energia mínima de FRB de $\log E_{min} (\text{erg}) = 39.47^{+0.54}_{-1.28}$. Este valor é significativamente mais alto do que as energias detectadas em repetidores fortes (por exemplo, FRB 20121102A e FRB 20201124A), que são frequentemente observados até $\sim 10^{37}$ erg.
• Parâmetros Populacionais: O índice espectral $\alpha$ prefere um valor plano de $0.11^{+0.66}_{-0.60}$, diferindo da inclinação negativa preferida em análises anteriores apenas com ASKAP. No entanto, quando $H_0$ é permitido variar livremente, a preferência retorna a uma inclinação negativa ($\alpha \approx -0.92$).
• Propagação de Incerteza: A inclusão da incerteza em $DM_{MW}$ leva a maiores incertezas nos parâmetros do modelo derivados em comparação com estimativas anteriores, apesar do aumento no tamanho do conjunto de dados.
• Previsões para Detecções Futuras:
  • FAST: Prevê-se que detecte 25–41% de seus FRBs em $z \gtrsim 2$ e 6–20% em $z \gtrsim 3$. A distribuição atinge o pico em $z \approx 1$, com 73–77% das detecções além de $z \gtrsim 1$.
  • DSA: Prevê-se que detecte 2–12% dos FRBs em $z \gtrsim 1$.
  • Limites dos Levantamentos: O DM máximo de busca atual para o FAST (3700–5000 $\text{ pc cm}^{-3}$) é previsto para perder 7–18% das detecções possíveis, enquanto o limite do DSA de $1500 \text{ pc cm}^{-3}$ é previsto para perder 3–8%.

Significância e Alegações
Os autores afirmam que este trabalho fornece a primeira restrição sobre $E_{min}$ derivada de uma análise $z$–DM. O alto valor de $E_{min}$ sugere duas possibilidades:

1. Há uma virada de baixa energia (achatamento) na função de luminosidade da população geral de FRBs, significando que os levantamentos não são sensíveis o suficiente para detectar a cauda de baixa energia.
2. Repetidores fortes possuem uma função de luminosidade fundamentalmente diferente em comparação com detecções de surtos únicos (ou surtos únicos aparentes).

O artigo enfatiza que, embora os novos dados expandam o espaço de parâmetros acessível para desvios para o vermelho mais altos, a incapacidade de localizar robustamente FRBs de alto $z$ com o FAST limita atualmente a capacidade de obter desvios para o vermelho associados. Consequentemente, os autores afirmam que medições precisas de $DM_{ISM}$ serão críticas para futuras restrições de parâmetros cosmológicos como $H_0$. O trabalho serve como um refinamento metodológico, demonstrando que incorporar incertezas nas contribuições galácticas é necessário para evitar distorcer artificialmente o poder de restrição de FRBs de baixa DM.