Search for Long-Transient Gravitational Waves from Supernova SN2023ixf using GFH-v2 Pipeline

Resumo Técnico: Busca por Ondas Gravitacionais de Transientes Longos de SN 2023ixf usando GFH-v2

Problema e Motivação
Magnetares recém-nascidos que rotacionam rapidamente, formados em supernovas de colapso de núcleo (CCSNe), são previstos para emitir ondas gravitacionais contínuas de transientes longos (tCWs). Diferente dos sinais de curta duração associados ao próprio colapso, esses sinais surgem do desaceleramento (spin-down) de uma estrela de nêutrons ou magnetar recém-formada e não axissimétrica. A emissão é impulsionada pela perda de energia rotacional, produzindo um sinal com uma frequência que diminui de acordo com uma lei de potência e uma amplitude de deformação (strain) que decai ao longo do tempo. A detecção desses sinais forneceria insights diretos sobre o nascimento, a estrutura interna e a evolução do spin de estrelas de nêutrons.

A SN 2023ixf, uma supernova de colapso de núcleo do Tipo II descoberta na vizinha galáxia M101 (distância $\approx$ 6,85 Mpc), apresenta uma oportunidade única para tal busca. Sua proximidade e o tempo preciso de sua detecção eletromagnética (dentro de uma janela de $\sim$2 horas) permitem uma busca direcionada. Crucialmente, a época do evento sobrepõe-se ao Engineering Run 15 (ER15) do LIGO, fornecendo um conjunto de dados de segmentos coincidentes dos detectores Hanford (H1) e Livingston (L1) imediatamente precedendo a quarta rodada de observação (O4).

Metodologia
A análise utiliza o pipeline GFH-v2, um método de busca semi-coerente hierárquica baseado na Transformada de Hough de Frequência Generalizada, especificamente projetado para sinais com evolução de frequência por lei de potência.

• Modelo de Sinal: A busca visa uma estrela de nêutrons rígida, não precessante, que sofre desaceleramento (spin-down) primariamente via emissão de ondas gravitacionais. Isso corresponde a um índice de frenagem $n=5$. A frequência de onda gravitacional $f_{gw}$ evolui como $\dot{f}_{gw} = -k f_{gw}^n$, e a amplitude de deformação $h_0(t)$ decai como $t^{-1/2}$. A busca assume parâmetros canônicos de estrela de nêutrons ($M=1,4 M_\odot$, $R=12$ km) e uma posição de céu fixa baseada na localização eletromagnética.
• Dados: A análise utiliza segmentos de ciência coincidentes de H1 e L1 durante o ER15, totalizando aproximadamente 48.619 segundos (0,56 dias) de dados dentro de uma janela on-source conservadora. A qualidade dos dados foi verificada como comparável à sensibilidade da O4a na faixa de 600–2000 Hz.
• Configuração de Busca: O espaço de parâmetros cobre frequências iniciais ($f_0$) de 600 a 2000 Hz (passo de 10 Hz) e ellipticidades ($\epsilon$) de $3 \times 10^{-4}$ a $3 \times 10^{-3}$ (divididas em cinco subintervalos). Isso resulta em 700 configurações de busca independentes.
• Etapas do Pipeline:
  1. Preparação de Dados: Os dados do Short Fourier Transform Database (SFDB) são passados por filtros de banda (band-pass) e transformados em séries temporais analíticas complexas.
  2. Geração de Peakmap: Máximos locais no espectro de potência equalizado são selecionados para formar um mapa de picos no tempo-frequência.
  3. Transformada de Hough: O peakmap é mapeado de $(t, f)$ para coordenadas transformadas $(t, x)$, onde $x = 1/f^{n-1}$. Neste espaço, a evolução da lei de potência torna-se linear. Uma transformada de Hough mapeia isso para o espaço de parâmetros $(x_0, k)$ (frequência inicial e parâmetro de spin-down).
  4. Seleção de Candidatos: Candidatos são identificados como bins com excesso de contagens em relação ao fundo local. Uma razão crítica (CR) é calculada usando estatísticas de ruído estimadas localmente (mediana e dispersão dentro de intervalos de 10 Hz) para contabilizar o ruído não estacionário.
  5. Coincidência: Candidatos de H1 e L1 são comparados no espaço de parâmetros $(x_0, k)$. Um par é considerado coincidente se a métrica de distância $d_{coin} \le 2$.
  6. Follow-up: Candidatos coincidentes com $CR \ge 3,5$ passam por um follow-up de heterodinação. A evolução de fase esperada é removida dos dados, e o sinal é reanalisado com tempos de coerência mais longos ($3\times$ e $6\times$ o original) para verificar se há um aumento na significância, o que indicaria um sinal astrofísico real.

Principais Resultados

• Seleção de Candidatos: Após aplicar o requisito de coincidência e o limiar de CR ($CR \ge 3,5$), 131 candidatos coincidentes restaram.
• Resultado do Follow-up: Nenhum dos 131 candidatos exibiu o aumento esperado na razão crítica durante as etapas de follow-up de heterodinação. Consequentemente, nenhum candidato foi promovido para a segunda etapa de follow-up.
• Detecção: A análise não encontrou evidências de emissão de ondas gravitacionais de transientes longos associadas à SN 2023ixf nos dados do ER15.
• Limites Superiores: Limites superiores para a distância máxima detectável foram estabelecidos como uma função da frequência inicial e da ellipticidade. Para o intervalo de maior ellipticidade ($\epsilon \in [2,5 \times 10^{-3}, 3,0 \times 10^{-3}]$), os limites superiores de eficiência de 90% variam de aproximadamente 1 a 2,5 Mpc em quase toda a faixa de frequência analisada. Para ellipticidades menores, os limites não puderam ser determinados ou foram restritos a distâncias menores devido à duração limitada dos dados e aos menores rapports sinal-ruído.

Significância e Alegações
Os autores afirmam que, embora os limites superiores derivados (1–2,5 Mpc) estejam abaixo da distância real para M101 ($\approx$ 6,85 Mpc), o trabalho representa a primeira aplicação do pipeline GFH-v2 a uma busca direcionada para uma supernova de colapso de núcleo próxima. O estudo caracteriza com sucesso o desempenho do pipeline em dados reais de comissionamento de detectores (ER15). Demonstra que o método GFH-v2 pode ser efetivamente implantado para buscar tCWs de eventos de formação de magnetar, estabelecendo um framework para buscas futuras utilizando conjuntos de dados mais sensíveis de próximas rodadas de observação (O5) e detectores de próxima geração. O artigo conclui que a busca fornece uma estimativa direta de sensibilidade para buscas de transientes longos nos dados do ER15 e serve como um benchmark para esforços futuros.