Constraining high-energy neutrinos from tidal disruption events with IceCube high-energy starting events

Resumo Técnico: Restrições de Neutrinos de Alta Energia de Eventos de Disrupção de Maré com IceCube HESE

Definição do Problema
Eventos de Disrupção de Maré (TDEs), onde estrelas são despedaçadas por buracos negros supermassivos, são candidatos propostos como fontes de neutrinos de alta energia. Modelos teóricos sugerem que a produção de neutrinos pode ocorrer em várias regiões, incluindo jatos relativísticos (tanto bem-sucedidos quanto sufocados/choked), discos de acreção, coronas e fluxos de saída (outflows). Embora estudos anteriores tenham identificado potenciais coincidências entre TDEs específicos (ex: AT2019dsg) e eventos de neutrinos, melhorias recentes na reconstrução direcional do IceCube mostraram que essas direções reconstruídas de neutrinos não se alinham com as localizações dos TDEs. Consequentemente, uma abordagem estatística é necessária para determinar se os TDEs, como uma população, contribuem para o fluxo difuso de neutrinos de alta energia. O desafio primário reside na baixa estatística tanto do catálogo de TDEs observado quanto do conjunto de dados de neutrinos de alta energia, exigindo uma análise de empilhamento (stacking analysis) robusta para restringir os parâmetros físicos que regem a produção de neutrinos.

Metodologia
Os autores realizam uma análise de máxima verossimilhança não binária (unbinned maximum likelihood analysis) para buscar correlações espaciais e temporais entre o conjunto de dados IceCube 12.5-year High-Energy Starting Events (HESE) e um catálogo de 89 TDEs.

• Catálogos de Dados:

  • TDEs: Uma amostra de 89 TDEs foi selecionada de um catálogo de 134 eventos [73]. Os critérios de seleção exigiram a disponibilidade de fotometria de série temporal para determinar o tempo de pico de brilho ($t_{peak}$) usando um algoritmo de Bayesian Blocks. A amostra abrange vários desvios para o vermelho (redshifts) e inclui eventos detectados em múltiplas bandas ópticas (ZTF, ATLAS, ASASSN, etc.).
  • Neutrinos: A análise utiliza 164 eventos HESE do conjunto de dados de 12,5 anos do IceCube (junho de 2010 – agosto de 2022). O conjunto de dados inclui direções de chegada reconstruídas modeladas por uma distribuição Fisher-Bingham (FB8) de 8 parâmetros para contabilizar incertezas de localização anisotrópicas.

• Análise de Verossimilhança:

  • A análise emprega uma função de verossimilhança não binária (Eq. 3.1) que empilha associações entre cada TDE e todos os eventos de neutrinos na vizinhança temporal e espacial.
  • PDF de Sinal ($S_i$): Modelado como o produto de componentes temporais e espaciais. O componente temporal assume uma distribuição Gaussiana centrada em um atraso de $\lambda_T = 100$ dias em relação ao pico óptico do TDE, com um desvio de $\sigma_T = 10$ dias, motivado por modelos de lançamento de jatos retardados (ex: estados de Disco Magneticamente Arrestado). O componente espacial utiliza a densidade de probabilidade direcional FB8 fornecida pelo IceCube.
  • PDF de Fundo ($B_i$): Assume uma distribuição uniforme no tempo (dentro de uma janela de $\sim$50 dias) e uma distribuição uniforme pelo céu.
  • Estatística de Teste (TS): Definida como $TS = 2 \log(L(ns)/L(ns=0))$, onde $ns$ é o parâmetro de força do sinal. O TS observado ($TS_{obs}$) é comparado contra uma distribuição de $TS_{max}$ gerada a partir de $10^5$ catálogos fictícios (mock) com posições e tempos de TDE gerados aleatoriamente para estabelecer a significância.

• Limites Superiores e Injeção de Sinal:

  • Para restringir parâmetros físicos, os autores geram catálogos sintéticos de TDE e realizam injeção de sinal artificial. Eles variam dois parâmetros chave: a fração de TDEs que possuem jatos intrínsecos ($f_{jet}$) e a energia de raios cósmicos isotrópica equivalente ($E_{CR}$).
  • Catálogos ficticos são construídos selecionando uma subamostra de TDEs com jato baseada em $f_{jet}$, calculando os rendimentos de neutrinos esperados com base em $E_{CR}$, e injetando pares de sinal no conjunto de dados com correlação espacial perfeita e atrasos temporalmente correlacionados.

Resultos Principais

• Busca de Correlação: A análise dos dados observados produz um valor máximo de estatística de teste de $TS_{max} = 0$ para $ns = 0$. Os dados observados são consistentes com a hipótese de apenas fundo, não mostrando nenhuma correlação espacial ou temporal significativa entre os 89 TDEs e os 164 neutrinos HESE.
• Restrições em Parâmetros Físicos: Usando o resultado nulo, os autores estabelecem limites superiores no espaço de parâmetros de $f_{jet}$ e $E_{CR}$.
  • Mesmo com estatísticas limitadas, a análise restringe a região onde $f_{jet} \gtrsim 0,6$ e $E_{CR} \gtrsim 3 \times 10^{53}$ erg com confiança superior a 90% (C.L.).
  • Para o caso específico onde todos os TDEs são emissores de neutrinos ($f_{jet} = 1$), a energia de raios cósmicos isotrópica equivalente é restringida a $E_{CR} < 2 \times 10^{53}$ erg a 90% C.L.
• Implicações Teóricas:
  • Os limites derivados restringem a eficiência ($\epsilon_{CR}$) de conversão da energia de ligação gravitacional em raios cósmicos. Para uma estrela padrão de $1 M_\odot$ e um buraco negro de $10^8 M_\odot$, os limites implicam $\epsilon_{CR} \lesssim 0,5$ (para $f_{jet}=1$) ou $\epsilon_{CR} \lesssim 0,75 f_b$ (para $f_{jet} < 0,6$).
  • Os resultados sugerem que, embora bem-sucedidos, os jatos rádio-fortes (que constituem uma fração minúscula dos TDEs) dificilmente dominarão o fluxo difuso; cenários envolvendo ventos "escondidos", jatos sufocados ou emissão isotrópica de TDEs não jetted permanecem viáveis, mas são restringidos pelo orçamento de energia.
  • O limite superior calculado para o fluxo difuso de neutrinos de todo o céu proveniente de TDEs é atualmente fraco ($\sim 10^{-7}$ GeV cm$^{-2}$ s$^{-1}$ sr$^{-1}$) comparado ao fluxo observado pelo IceCube, mas os autores observam que isso melhorará com dados futuros.

Significância e Alegações
O artigo afirma ser a primeira busca de empilhamento correlacionando o conjunto de dados IceCube HESE com um catálogo específico de TDEs observados. A principal significância reside em estabelecer que os dados atuais não apoiam uma forte correlação, portanto, descartando modelos onde uma grande fração dos TDEs ($f_{jet} \gtrsim 0,6$) produz raios cósmicos de energia extremamente alta ($E_{CR} \gtrsim 3 \times 10^{53}$ erg).

Os autores enfatizam que sua metodologia fornece um arcabouço para estabelecer restrições rigorosas sobre os parâmetros físicos dos TDEs à medida que as estatísticas melhoram. Eles declaram explicitamente que com o advento do Observatório Vera C. Rubin (LSST), que deve detectar $\sim 1000$ TDEs por ano, e futuros observatórios de neutrinos (IceCube-Gen2, KM3NeT, GRAND, etc.), os limites atuais podem ser significativamente apertados. O artigo conclui modestamente que, embora as estatísticas atuais sejam limitadas, a análise demonstra o potencial de restringir o orçamento de energia e as frações de jato dos TDEs, oferecendo um caminho para compreender o seu papel no céu de neutrinos de alta energia.