From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations

Resumo Técnico: De Simulações Cosmológicas a Fusões de Buracos Negros Binários

Definição do Problema
Observações de fusões de buracos negros binários (BBH) fornecem insights críticos sobre a evolução de estrelas massivas ao longo do tempo cósmico. No entanto, conectar as propriedades observadas das fusões aos seus progenitores estelares requer modelos precisos da história da formação estelar cósmica (SFRH) e do enriquecimento químico. Atualmente, a área frequentemente recorre a ajustes analíticos para a densidade da taxa de formação estelar cósmica dependente da metalicidade, $S(Z, z)$, como proxies para histórias subjacentes complexas derivadas de simulações cosmológicas ou observações. Essas prescrições analíticas são frequentemente não restringidas e seu impacto nas populações de BBH inferidas permanece insuficientemente compreendido. Além disso, os efeitos da resolução e do tamanho da caixa das simulações cosmológicas no $S(Z, z)$ derivado e, consequentemente, nas taxas de fusão de BBH, não foram sistematicamente quantificados.

Metodologia
Os autores utilizam a suíte de simulações cosmológicas IllustrisTNG (especificamente TNG50-1, TNG100-1 e TNG300-1), que variam em resolução e tamanho de caixa, para derivar o $S(Z, z)$ baseado em simulação diretamente dos dados de células de gás. Eles comparam esses resultados diretos da simulação contra um modelo de ajuste analítico amplamente utilizado (baseado em L. A. C. van Son et al. 2023), que combina uma forma funcional de Madau & Dickinson (2014) para a densidade total de SFR, $S(z)$, com uma distribuição de metalicidade log-normal assimétrica.

Para avaliar o impacto desses diferentes inputs de $S(Z, z)$ nas populações de BBH, os autores acoplam tanto os modelos de $S(Z, z)$ derivados de simulação quanto o analítico ao código de síntese de população binária rápida COMPAS (Compact Object Mergers: Population Astrophysics and Statistics). Eles empregam um framework de "integração cósmica" para calcular as taxas de fusão e distribuições de massa através do redshift. O estudo investiga especificamente:

1. A fidelidade do ajuste analítico em reproduzir a forma e a normalização do $S(Z, z)$ de diferentes resoluções de TNG.
2. Os desvios resultantes nas taxas de fusão de BBH e distribuições de massa primária ($M_{BH,1}$) ao usar o ajuste analítico versus os dados completos da simulação.
3. A influência da resolução da simulação (variando de alta resolução TNG50-1 para baixa resolução TNG300-1) sobre esses resultados.

Principais Resultados
O estudo descobre que prescrições analíticas para $S(Z, z)$ introduzem erros sistemáticos significativos ao modelar populações de BBH:

• Discrepâncias em $S(Z, z)$: O ajuste analítico falha em capturar características estruturais fundamentais da distribuição de metalicidade baseada em simulação. Especificamente, ele perde um "bump" secundário no regime de alta metalicidade ($Z \gtrsim 3 Z_\odot$) e uma cauda de baixa metalicidade achatada presentes nas simulações. O ajuste tende a superestimar a taxa de formação estelar em redshifts altos ($z \gtrsim 6$) e altas metalicidades, enquanto subestima contribuições em baixas metalicidades e redshifts muito altos, dependendo da resolução da simulação.
• Impacto nas Taxas de Fusão: Em redshifts altos ($z \gtrsim 6$), o ajuste analítico superestima as taxas de fusão de BBH por fatores que variam de 10 a $10^4$ em comparação com as taxas baseadas em simulação, sendo que a magnitude do erro aumenta conforme a resolução da simulação diminui. Em baixos redshifts ($z \approx 0$), as taxas de fusão locais derivadas do ajuste geralmente concordam com as taxas baseadas em simulação dentro de um fator de $\sim 1,5$, embora o TNG50-1 mostre uma superestimação de $\sim 15$ Gpc$^{-3}$ yr$^{-1}$.
• Impacto nas Distribuições de Massa: O uso de ajustes analíticos introduz características espúrias na distribuição de massa primária de BBH. Notavelmente, em $z \lesssim 2$, o ajuste produz um excesso artificial (bump) próximo a $8 M_\odot$ na distribuição de massa primária, o qual está ausente quando se usa o $S(Z, z)$ completo baseado em simulação. Essa característica é atribuída ao canal de transferência de massa estável e é causada pela incapacidade do ajuste em modelar corretamente o pico de formação estelar de alta metalicidade. Por outro lado, o ajuste suprime características de alta massa (próximo a $33 M_\odot$ e $45 M_\odot$) por fatores de 3–5 em $M_{BH,1} \gtrsim 35 M_\odot$.
• Dependência de Resolução: A resolução da simulação afeta significativamente a normalização e a localização do pico de $S(z)$. Simulações de maior resolução (TNG50-1) produzem taxas de formação estelar mais altas e atingem o pico em redshifts ligeiramente menores ($z \approx 2$) comparados a execuções de menor resolução ($z \approx 2,5$). A precisão do ajuste analítico degrada-se com a menor resolução da simulação, particularmente ao reproduzir a inclinação de alto redshift.

Significância e Alegações
Os autores concluem que, embora os ajustes analíticos ofereçam conveniência computacional, eles são proxies não confiáveis para o complexo $S(Z, z)$ subjacente derivado de simulações cosmológicas, particularmente ao interpretar observações de ondas gravitacionais. O estudo demonstra que modelos analíticos simples falham em capturar características físicas críticas — como bumps de alta metalicidade e caudas específicas de baixa metalicidade — que influenciam diretamente as taxas de fusão e distribuições de massa de BBH previstas.

O artigo enfatiza que as discrepâncias surgem porque as formas funcionais dos atuais ajustes analíticos não conseguem acomodar a natureza multimodal das distribuições de metalicidade encontradas em simulações. Consequentemente, confiar nesses ajustes pode levar a inferências incorretas sobre os canais de formação e a história evolutiva de BBHs. Os autores argumentam pela integração cuidadosa de simulações cosmológicas, ajustes analíticos e síntese de população para evitar a introdução de características espúrias na interpretação de dados de ondas gravitacionais. Eles não alegam que os ajustes analíticos sejam inteiramente inúteis, mas destacam suas limitações e a necessidade de validá-los contra dados de simulação de alta fidelidade, especialmente à medida que detectores de próxima geração sondam redshifts mais altos.