The Luminosity of the Darkness: Schechter function in dark sirens

Este artigo demonstra que incorporar a evolução da função de Schechter na modelagem de sirenes escuras é crucial para evitar viéses na medição da constante de Hubble ($H_0$) e dos parâmetros de taxa, especialmente para eventos de ondas gravitacionais distantes com catálogos de galáxias incompletos.

O essencial

Imagine que o Universo é uma cidade gigante e escura, e nós, os astrônomos, somos detetives tentando descobrir o tamanho dessa cidade e quão rápido ela está crescendo.

Para medir o tamanho da cidade, usamos dois tipos de ferramentas:

1. As "Sirenes Escuras" (Dark Sirens): São ondas sonoras do espaço (ondas gravitacionais) que nos dizem quão longe um evento ocorreu, mas não nos dizem onde exatamente. É como ouvir um grito no escuro e saber que vem de longe, mas não ver quem gritou.
2. O Catálogo de Galáxias: É como um mapa de endereços da cidade. Quando conseguimos identificar qual galáxia "hospedou" o grito, podemos medir a velocidade de expansão do universo (a Constante de Hubble, ou $H_0$).

O problema é que nosso mapa (o catálogo de galáxias) é incompleto. Ele é muito bom para as galáxias vizinhas, mas fica "borrado" e cheio de buracos quando olhamos para galáxias muito distantes.

O que os autores descobriram?

Neste artigo, os cientistas Cezary Turski e sua equipe perguntaram: "E se o mapa que estamos usando estiver desatualizado?"

Eles focaram em uma regra matemática chamada Função de Schechter. Pense nela como uma "receita" que diz quantas galáxias existem e quão brilhantes elas são em diferentes distâncias.

• A Velha Receita (Modelo Constante): Por muito tempo, os cientistas assumiram que a receita é a mesma para sempre. Eles achavam que, se você olhasse para o passado (galáxias distantes), a quantidade e o brilho das galáxias seriam os mesmos de hoje. Era como se a cidade tivesse crescido, mas o número de prédios e o tamanho deles tivessem permanecido estáticos ao longo da história.
• A Nova Receita (Modelo Evolutivo): Os autores mostraram que a receita muda com o tempo. No passado (redshifts altos), as galáxias eram diferentes: havia menos delas, mas elas eram mais brilhantes (ou vice-versa, dependendo da métrica). É como se, no passado da cidade, houvesse menos prédios, mas todos fossem arranha-céus gigantes, enquanto hoje temos muitos prédios menores espalhados.

A Analogia do "Filtro de Luz"

Imagine que você está tentando contar quantas lâmpadas existem em uma sala escura usando uma câmera.

• O Modelo Antigo diz: "Todas as lâmpadas têm o mesmo brilho, não importa a distância."
• O Modelo Novo diz: "As lâmpadas antigas (longe) eram mais fortes, e as novas (perto) são mais fracas, e a quantidade delas também mudou."

Se você usar o modelo antigo para contar as lâmpadas distantes, você vai errar a contagem. Você pode achar que há menos lâmpadas do que realmente existem, ou que elas estão mais longe do que estão.

O Que Isso Muda para a Ciência?

Os autores testaram o que acontece quando aplicam essa "receita nova" (evolutiva) aos dados das ondas gravitacionais:

1. Correção de Viés: Se ignorarmos que as galáxias mudam com o tempo, nossa medição da velocidade de expansão do universo ($H_0$) pode ficar levemente torta. É como tentar medir a velocidade de um carro usando um velocímetro que não foi recalibrado para diferentes tipos de pneus.
2. O Efeito é Pequeno (por enquanto): Com os dados atuais, a diferença na resposta final é pequena. Mas, à medida que nossos telescópios ficarem mais potentes e ouvirmos "gritos" de galáxias cada vez mais distantes, esse erro vai aumentar.
3. A Solução: Eles descobriram que, se deixarmos outros parâmetros (como a taxa de fusão de buracos negros) variarem livremente durante o cálculo, o erro na medição da velocidade do universo é corrigido. No entanto, isso significa que precisamos entender muito bem como as galáxias evoluem para não confundir uma coisa com a outra.

Em Resumo

Este trabalho é um lembrete importante: para medir o universo com precisão, precisamos entender como ele mudou.

Não basta apenas olhar para o mapa atual; precisamos saber como o mapa era no passado. Se quisermos resolver o mistério da "Tensão de Hubble" (a briga entre diferentes medições da velocidade do universo), precisamos refinar nossa "receita" de como as galáxias nascem, brilham e morrem ao longo do tempo cósmico.

É como se, ao tentar adivinhar a idade de uma pessoa olhando apenas uma foto, você precisasse saber se ela usava lentes de contato ou óculos no passado, caso contrário, sua estimativa de idade pode estar errada.

Resumo técnico

Resumo Técnico

1. O Problema

A astronomia de ondas gravitacionais (OG) oferece uma via independente para medir a constante de Hubble ($H_0$) utilizando "sirenes escuras" (eventos de fusão de objetos compactos sem contraparte eletromagnética identificada). Este método combina a distância inferida a partir do sinal de OG com o redshift de galáxias candidatas extraídas de catálogos.

No entanto, existe um limite fundamental: os catálogos de galáxias atuais (como o GLADE+) são incompletos em altos redshifts (geralmente completos apenas até $z \approx 0.1$), enquanto as detecções de OG ocorrem até $z \approx 1.2$. Para preencher essa lacuna, é necessário modelar a distribuição de luminosidade das galáxias fora do catálogo ("a luminosidade da escuridão").

O problema central abordado neste trabalho é que a maioria das análises assume uma Função de Schechter (SF) constante em relação ao redshift. A literatura, contudo, indica que a função de luminosidade das galáxias evolui com o tempo cósmico (o redshift). Ignorar essa evolução pode introduzir viéses sistemáticos na estimativa de $H_0$, especialmente para eventos distantes com pouco suporte de catálogo, e afetar a inferência dos parâmetros da taxa de fusão de binários compactos.

2. Metodologia

Os autores utilizaram dados de ondas gravitacionais da colaboração LIGO-Virgo-KAGRA (LVK) até a terceira campanha de observação (O3), selecionando 47 eventos (42 BBH, 3 NSBH, 2 BNS) com SNR > 11, tratando-os como sirenes escuras (exceto GW170817).

A abordagem metodológica envolveu:

• Modelo de Evolução da Função de Schechter: Em vez de assumir parâmetros fixos, os autores implementaram uma SF dependente do redshift. Eles parametrizaram a evolução dos parâmetros característicos $M^*$ (magnitude característica) e $\Phi^*$ (densidade normalizada) com base em dados observacionais de Arnouts et al. (2007) para a banda $K$:
  • $M^*(z) = M^*_0 - Qz$
  • $\Phi^*(z) = \Phi^*_0 10^{0.4Pz}$
  • O parâmetro de inclinação $\alpha$ foi mantido constante ($\alpha = -1.1$), devido à dificuldade em medir sua evolução.
• Cálculo do Prior de Redshift (LOS): Utilizaram o pipeline gwcosmo para calcular o prior de redshift na linha de visão ($p(z)$). Este prior é dividido em duas partes:
  1. Dentro do catálogo: Galáxias presentes no GLADE+.
  2. Fora do catálogo: Galáxias não detectadas, modeladas pela Função de Schechter evolutiva ou constante.
• Comparação de Cenários: Realizaram análises bayesianas comparando quatro cenários:
  1. Catálogo GLADE+ com SF constante.
  2. Catálogo GLADE+ com SF evolutiva.
  3. Catálogo vazio (apenas SF constante).
  4. Catálogo vazio (apenas SF evolutiva).
• Estimação Conjunta: Além de fixar a taxa de fusão, realizaram uma estimativa conjunta de $H_0$ e dos parâmetros de evolução da taxa de fusão ($\gamma$), permitindo que variassem livremente.

3. Principais Contribuições

• Primeira Aplicação de SF Evolutiva em Sirenes Escuras: Este trabalho é pioneiro em incorporar explicitamente a evolução da função de luminosidade das galáxias com o redshift no contexto de medição de $H_0$ via sirenes escuras.
• Análise de Viés Sistemático: Quantificaram como a suposição de uma SF estática distorce a distribuição de probabilidade de redshift para eventos fora do catálogo, especialmente em altos redshifts.
• Interação entre Parâmetros: Demonstraram a degenerescência entre a evolução da Função de Schechter e os parâmetros de evolução da taxa de fusão de binários compactos.

4. Resultados

• Impacto no Prior de Redshift: A SF evolutiva prevê uma densidade de galáxias diferente em comparação com a SF constante. Em baixos redshifts, a SF evolutiva sugere uma densidade maior (devido ao aumento da luminosidade característica e densidade no passado recente), enquanto em altos redshifts prevê menos galáxias. Isso altera o prior de redshift para eventos distantes.
• Medição de $H_0$:
  • Ao usar apenas a SF constante, há um viés leve para valores mais baixos de $H_0$ em comparação com a SF evolutiva.
  • No entanto, a magnitude desse viés é menor que as incertezas estatísticas atuais e menor que o impacto da incerteza no modelo da população de OG (massas e taxas).
  • Para eventos bem localizados e próximos (como GW190814), a diferença entre os modelos é insignificante ($O(10^{-5})$). Para eventos distantes e mal localizados (como GW150914), a diferença é mais perceptível, mas ainda dentro das margens de erro atuais.
• Parâmetros de Taxa ($\gamma$): Quando se permite que a taxa de fusão varie em conjunto com $H_0$, o viés em $H_0$ é corrigido. Contudo, o parâmetro de taxa $\gamma$ (inclinação da lei de potência em baixos redshifts) torna-se dependente do modelo de SF escolhido. A evolução da SF é degenerada com a evolução da taxa.
• Tabela de Resultados (Tabela 1):
  • $H_0$ com GLADE+ e SF Evolutiva: $66.06^{+12.32}_{-13.33}$ km s$^{-1}$Mpc$^{-1}$.
  • $H_0$ com GLADE+ e SF Constante: $64.85^{+10.91}_{-14.55}$ km s$^{-1}$Mpc$^{-1}$.
  • A diferença é pequena, mas sistematicamente relevante para futuras medições de alta precisão.

5. Significância e Conclusões

O estudo conclui que, embora a evolução da Função de Schechter não altere drasticamente as medições atuais de $H_0$ (devido às grandes incertezas estatísticas e à dominância do modelo de população de OG), sua importância aumentará significativamente no futuro.

• Futuro: Com detectores mais sensíveis (como o Einstein Telescope e o Cosmic Explorer) e catálogos de galáxias mais profundos (como o do Vera Rubin Observatory), as incertezas estatísticas diminuirão. Nesse cenário, os erros sistemáticos, como a modelagem incorreta da evolução da população de galáxias, tornar-se-ão o fator limitante.
• Recomendação: É crucial incorporar modelos de SF evolutiva nas análises de sirenes escuras e realizar uma estimativa conjunta dos parâmetros de $H_0$ e da taxa de fusão para evitar viéses.
• Limitação Atual: O parâmetro de inclinação $\alpha$ da SF ainda é pouco conhecido em altos redshifts. Futuras sondas mais profundas serão necessárias para medir a evolução de $\alpha$ e refinar ainda mais o modelo.

Em suma, o artigo destaca que a "luminosidade da escuridão" (o modelo de galáxias não detectadas) não é estática e deve evoluir com o tempo cósmico para garantir a precisão da cosmologia de ondas gravitacionais na próxima década.