Empirical estimates of how massive galaxies can be in {\Lambda}CDM

Resumo Técnico: Estimativas empíricas de quão massivas podem ser as galáxias em ΛCDM

Declaração do Problema
Observações recentes do Telescópio Espacial James Webb (JWST) revelaram uma população de galáxias inesperadamente massivas em altos desvios para o vermelho ($6 \lesssim z \lesssim 10$), com massas estelares inferidas superiores a $M_* \sim 10^{10} M_\odot$. Essas descobertas reacenderam o debate sobre a compatibilidade dessas galáxias com o modelo cosmológico padrão $\Lambda$CDM. Estudos anteriores sugeriram que tais galáxias poderiam exigir eficiências de formação estelar que se aproximam ou excedem o orçamento de bárions cósmicos disponível, desafiando potencialmente a montagem hierárquica de estruturas prevista pelo $\Lambda$CDM. No entanto, essas comparações frequentemente recorrem às funções de massa estelar observadas (GSMF) sem levar plenamente em conta incertezas sistemáticas, como o viés de Eddington decorrente de erros aleatórios nas estimativas de massa estelar e a dispersão intrínseca na relação massa-estrela/massa-halo (SHMR). Este artigo aborda se a tensão aparente entre as observações do JWST e o $\Lambda$CDM persiste quando esses efeitos estatísticos e observacionais são rigorosamente incorporados.

Metodologia
Os autores empregam um quadro empírico baseado em Estatística de Valores Extremos (EVS) para estimar os limites probabilísticos das galáxias e halos de matéria escura mais massivos em função do desvio para o vermelho e da área do levantamento.

1. Quadro EVS: Em vez de depender de modelos específicos de formação de galáxias, o estudo utiliza a GSMF observada e a Função de Massa de Halos (HMF) para derivar a distribuição de probabilidade dos objetos mais extremos. A função de distribuição acumulada (CDF) para a massa estelar máxima, $M_{*,\text{max}}$, é derivada da CDF da GSMF, $F_*(M_*)$, elevada à potência do número de galáxias independentes, $N$, em um determinado volume de levantamento.
2. Correção para Incertezas: A análise distingue três níveis de GSMF:
   • Observada ($\phi_{*,\text{obs}}$): Os dados brutos.
   • Verdadeira ($\phi_{*,\text{true}}$): Desconvoluída para remover erros aleatórios (viés de Eddington) nas estimativas de massa estelar.
   • Intrínseca ($\phi_{*,\text{int}}$): Further desconvoluída para contabilizar a dispersão intrínseca na SHMR.
     Os autores modelam esses erros usando distribuições Gaussianas, com faixas de incerteza aumentando com o desvio para o vermelho (até $\sim 0,6$ dex em $z \sim 7-8$).
3. Restrições de Halos: O estudo calcula a massa bariônica máxima teórica disponível nos halos mais massivos ($M_{\text{bar,max}} = f_{\text{bar}} M_{\text{vir,max}}$) usando sete HMFs diferentes (por exemplo, Despali et al. 2016, Sheth & Tormen 1999) e parâmetros cosmológicos do Planck.
4. Conexão Galáxia-Halo: Para comparar galáxias e halos, os autores exploram três cenários de correlação:
   • Estocástico: Assumindo que a galáxia mais massiva e o halo mais massivo são variáveis aleatórias independentes.
   • Dispersão Correlacionada Positivamente (PCS): Assumindo uma correspondência de ordem de classificação onde o halo mais massivo hospeda a galáxia mais massiva.
   • Classificação Invertida: O halo mais massivo hospeda a galáxia menos massiva (usado como limite inferior).
5. Reservatórios Bariônicos: O estudo estima a massa total de gás frio ($M_{\text{cold}}$) disponível para formação estelar modelando fluxos frios e fluxos de resfriamento radiativo, distinguindo entre modos de acreção quente e frio.

Principais Contribuições e Resultados

• Limites Empíricos sobre a Massa da Galáxia: Os autores fornecem previsões EVS para as galáxias mais massivas em várias áreas de levantamento (de $10 \text{ arcmin}^2$ até o céu inteiro). Para um levantamento de céu inteiro em $z \sim 0$, a galáxia mais massiva esperada a $5\sigma$ tem $M_* \sim 7 \times 10^{12} M_\odot$, diminuindo para $\sim 10^{10} M_\odot$ em $z \sim 16$. Essas estimativas são amplamente consistentes com as observações atuais do JWST quando a área do levantamento é considerada.
• Impacto do Viés de Eddington: O estudo demonstra que erros aleatórios nas estimativas de massa estelar inflacionam significativamente a abundância observada de galáxias massivas. Ao corrigir esses erros (passando de $\phi_{*,\text{obs}}$ para $\phi_{*,\text{true}}$), a massa estelar máxima inferida diminui em até $\sim 0,4-0,5$ dex em $z \sim 5$ e quase $\sim 1$ dex em $z \sim 16$. Incluir a dispersão intrínseca da SHMR reduz ainda mais a massa máxima em um adicional de $\sim 0,1-0,2$ dex em $z \sim 5$ e $\sim 0,5-0,6$ dex em $z \sim 16$.
• Resolução da Tensão da "Galáxia Impossível": Ao comparar a massa estelar máxima intrínseca ($M_{*,\text{max|int}}$) com a massa bariônica máxima disponível ($M_{\text{bar,max}}$), a tensão aparente desaparece em grande parte. Sob a suposição de uma correspondência perfeita de ordem de classificação (PCS), nenhuma galáxia é encontrada excedendo o limite bariônico teórico ($M_* < 0,16 \times M_{\text{vir}}$) em qualquer desvio para o vermelho, uma vez que as incertezas de medição são incluídas.
• Eficiência de Formação Estelar (SFE):
  • Em desvios para o vermelho intermediários ($2 \lesssim z \lesssim 6$), as galáxias mais massivas aproximam-se do reservatório total de gás frio ($M_{*,\text{max|int}} \sim M_{\text{cold}}$), sugerindo eficiências de formação estelar próximas do máximo ($\epsilon \sim 1$) durante essa época, consistente com o pico de galáxias formadoras de estrelas submilimétricas/poeirentas.
  • Em desvios para o vermelho mais altos ($z \gtrsim 8$), a eficiência inferida cai significativamente ($\epsilon \ll 1$) ao contabilizar as incertezas, indicando que os halos não estão convertendo bárions com a eficiência extrema sugerida anteriormente.
• Luminosidade UV e Poeira: O estudo modela as luminosidades UV levando em conta a atenuação por poeira e a fração de gás já convertida em estrelas. Embora o modelo reproduza o brilho UV observado para pequenas áreas de levantamento em altos desvios para o vermelho, ele prevê que as galáxias mais brilhantes em áreas de levantamento maiores devem ser mais fracas do que algumas detecções atuais do JWST, sugerindo uma possível superestimação das luminosidades ou efeitos de seleção nos maiores campos.
• Parâmetros Cosmológicos: Variações nos parâmetros cosmológicos (dentro das incertezas do Planck) têm um impacto negligenciável nas previsões EVS para os halos mais massivos, indicando que a tensão não é impulsionada por incertezas nos parâmetros cosmológicos.

Significado e Alegações
O artigo alega que o conflito aparente entre as observações do JWST de galáxias massivas de alto desvio para o vermelho e o modelo $\Lambda$CDM é em grande parte um artefato de efeitos estatísticos e observacionais, especificamente o viés de Eddington e a dispersão intrínseca na SHMR. Ao aplicar um quadro EVS totalmente empírico que contabiliza rigorosamente essas incertezas e volumes de levantamento, os autores concluem que as observações atuais permanecem amplamente consistentes com o $\Lambda$CDM.

O estudo enfatiza que a área do levantamento é um parâmetro fundamental na definição dos extremos esperados das populações de galáxias; levantamentos de pequena área são intrinsecamente enviesados para extremos de massa mais baixos e não podem sondar a cauda completa da distribuição. Os autores afirmam que o problema da "galáxia impossível" provavelmente decorre de incertezas sistemáticas nas estimativas de massa estelar, e não de uma falha do modelo cosmológico padrão ou da necessidade de nova física (como Energia Escura Primordial). Eles sugerem ainda que os halos mais massivos em $2 \lesssim z \lesssim 6$ podem de fato hospedar galáxias operando com eficiências de formação estelar próximas do máximo, enquanto em desvios para o vermelho mais altos, a eficiência é menor do que alegado anteriormente.