Residual Galactic binary foreground in LISA stochastic gravitational-wave background inference: source power concentration and spectral degeneracy

Resumo Técnico: Resíduo de Primeiro Plano de Binárias Galácticas em Inferência de Fundo Estocástico de Ondas Gravitacionais de LISA

Definição do Problema
O Laser Interferometer Space Antenna (LISA) observará uma população densa de binárias compactas galácticas, particularmente binárias de anãs brancas, na banda de milihertz. Embora uma fração desses processos possa ser resolvida individualmente, o restante formará um primeiro plano (foreground) não resolvido ou um componente de confusão. Mesmo após a subtração de fontes recuperadas com confiança, o poder residual de binárias subtraídas imperfeitamente ou não resolvidas permanece. Este primeiro plano residual pode enviesar a inferência do Fundo de Ondas Gravitacionais Estocásticas (SGWB) ou aumentar sua incerteza. O artigo aborda como a incerteza em um espectro de resíduo de binárias galácticas, construído a partir de dados de catálogo, propaga-se para a amplitude de $\Omega_0$ do SGWB inferido.

Metodologia
Os autores utilizam dados do conjunto de injeção LISA Data Challenge 2A (LDC2A) "Sangria" e da tabela de comparação de ajuste global de binárias galácticas "Erebor".

• Construção do Resíduo: As fontes são classificadas como "confiadamente recuperadas" se atenderem aos critérios de confiança de recuperação ($C_{rec} \ge 0,9$) e sobreposição ($O_{best} \ge 0,9$). As fontes que falham nesses critérios constituem a população "residual".
• Modelagem Espectral: A análise cobre a banda de 0,4–6,0 mHz, dividida em 560 bins de frequência ($\Delta f = 10^{-5}$ Hz). O espectro residual é construído usando potências de fonte de Michelson X de comprimento de onda longo com média orbital, normalizadas para uma convenção de deformação incidente equivalente.
• Concentração de Potência da Fonte: Os autores introduzem uma métrica de concentração de potência da fonte adimensional, $C_2 = \sum w_i^2$, onde $w_i$ são potências de fonte normalizadas. Em um modelo de fase aleatória, essa concentração determina a curtose excessiva do resíduo da soma ($\gamma_{2,res} = -1,5 C_2$).
• Covariância e Inferência: A covariância para a inferência do SGWB é derivada pela soma das variâncias das amplitudes de Fourier gaussianas independentes em cada bin. O modelo de SGWB assume um espectro de lei de potência $\Omega_{GW}(f) = \Omega_0 (f/f_{ref})^\alpha$. O poder residual é escalonado por um fator adimensional $\beta$.
• Estrutura Estatística: Uma matriz de informação de Fisher é usada para quantificar a inferência conjunta da amplitude do SGWB $\Omega_0$ e do fator de escala do resíduo $\beta$. O estudo avalia o aumento na incerteza de $\Omega_0$ ao marginalizar sobre $\beta$ e o viés introduzido se o resíduo for omitido.

Principais Contribuições e Resultados

• Degenerescência e Incerteza: Para um SGWB sem dependência de frequência ($\alpha=0$) com amplitude fiduciária $\Omega_0 = 10^{-11}$, a marginalização sobre o fator de resíduo $\beta$ aumenta a incerteza em $\Omega_0$ em 13,6% quando o poder do resíduo é distribuído uniformemente entre as frequências de Fourier dentro de cada bin.
• Restrições de Prior: Para limitar o aumento da incerteza a 10%, é necessária uma priori Gaussiana em $\beta$ com um desvio padrão de $\sigma_\beta \approx 0,0073$.
• Impacto da Estrutura de Frequência: O aumento da incerteza é sensível a como o poder residual é distribuído entre as frequências de Fourier. Distribuir o poder sobre o drift intrínseco e o intervalo Doppler anual reduz o aumento para 11,2%, enquanto concentrar o poder na frequência de Fourier mais próxima reduz o aumento para 5,1%. Isso indica que a estrutura de frequência não resolvida afeta quantitativamente a priori de $\beta$.
• Viés por Omissão: Se o resíduo fiduciário for omitido do modelo enquanto a covariância for mantida fixa, o melhor ajuste de $\Omega_0$ desloca-se em 119,5 vezes a incerteza de $\beta$ fixo. Os autores esclarecem que este deslocamento representa uma projeção do resíduo sobre o template do SGWB ponderada pelo inverso da variância, não uma significância de detecção posterior.
• Discreticidade da Fonte vs. Degenerescência: O estudo distingue entre a curtose excessiva de fase aleatória (dirigida pela concentração de potência da fonte) e a degenerescência espectral na potência média por bin. Embora o ruído instrumental e o poder do SGWB fiduciário diluam fortemente a curtose excessiva na maioria dos bins (mediana $|\gamma_{2,tot}| \approx 1,4 \times 10^{-2}$), a sobreposição espectral entre o resíduo e o SGWB continua sendo o principal motor da incerteza de inferência.
• Análise de Subconjuntos: Analisar subconjuntos de bins classificados pelo excesso de curtose diluído mostra que a largura da priori necessária e o fator de aumento da incerteza ($I_0$) dependem da seleção específica de frequência e do acúmulo de elementos de informação de Fisher, e não apenas do grau de discreticidade da fonte.

Significância e Alegações
O artigo afirma fornecer uma avaliação quantitativa de como as incertezas de nível de catálogo nos resíduos de binárias galácticas se propagam na inferência do SGWB. Ele demonstra que:

1. Os primeiros planos residuais introduzem uma degenerescência específica entre a amplitude do SGWB e a escala de potência do resíduo, quantificável via informação de Fisher.
2. A "curtose excessiva" do resíduo da soma, embora seja uma medida da discreticidade da fonte, não determina a covariância da potência média por bin; esta última é governada pela variância do modo Gaussiano.
3. Os resultados numéricos são condicionais ao modelo de escala adotado, ao ruído instrumental fixo e à suposição de covariância de potência de modo de Fourier independente.
4. A projeção do espectro do resíduo sobre o espectro do SGWB pode causar viés significativo se o resíduo for ignorado, mesmo que o poder do resíduo seja pequeno em relação ao primeiro plano total.

Os autores enfatizam que seus valores numéricos são condicionais ao modelo de potência escalar adotado e não consideram incertezas posteriores nos parâmetros das fontes recuperadas, erros de subtração de forma de onda ou respostas do detector dependentes do tempo. Eles sugerem que extensões futuras de nível de canal possam propagar os resíduos do ajuste global através da resposta completa dependente do tempo para quantificar melhor esses efeitos.