Galaxy formation in modified gravity -- II. galaxy halo connection and assembly bias

Resumo Técnico: Formação de galáxias em gravidade modificada – II. Conexão halo-galáxia e viés de montagem (assembly bias)

Problema
Pesquisas cosmológicas modernas (por exemplo, DESI, Euclid) visam restringir parâmetros cosmológicos e testar teorias da gravidade mapeando a estrutura em grande escala (LSS) do Universo. No entanto, a interpretação desses dados é limitada pelas incertezas na "conexão galáxia-halo" — o mapeamento entre as galáxias observadas e os halos de matéria escura subjacentes. As análises padrão frequentemente dependem do modelo de Distribuição de Ocupação de Halo (HOD), que assume que a ocupação de galáxias depende exclusivamente da massa do halo. Essa suposição negligencia o "viés de montagem" (AB - assembly bias), onde a formação de galáxias também depende de propriedades secundárias do halo, como o histórico de formação, o ambiente e a concentração. Embora o AB seja conhecido por afetar o agrupamento (clustering) no modelo padrão $\Lambda$CDM, seu comportamento em teorias de Gravidade Modificada (MG), particularmente na gravidade $f(R)$, permanece pouco compreendido. A complexa interação entre os mecanismos de blindagem (screening) de MG (por exemplo, blindagem camaleão) e a física de formação de galáxias pode introduzir vieses significativos na estimativa de parâmetros cosmológicos se não forem devidamente contabilizados.

Metodologia
Os autores utilizam simulações hidrodinâmicas de última geração realizadas com o código MG-AREPO, que incorpora toda a física bariônica do modelo IllustrisTNG, mas adaptado para a gravidade $f(R)$. A suíte de simulações inclui uma gama de modelos $f(R)$ de Hu-Sawicki com variações na força da quinta força, denotada pela magnitude do campo escalar de fundo hoje, $|f_{R0}| = 10^{-x}$ (onde $x \in \{4.0, 4.5, 5.0, 5.5, 6.0\}$), juntamente com um controle padrão $\Lambda$CDM ($f_{R0}=0$).

A partir dessas simulações, os autores geram catálogos fictícios (mock catalogues) para duas populações de galáxias relevantes para levantamentos de Estágio-IV:

1. Galáxias Vermelhas Luminosas (LRGs): Selecionadas por classificação de massa estelar.
2. Galáxias de Linha de Emissão (ELGs): Selecionadas com base na taxa específica de formação estelar (sSFR) e um corte de massa estelar ($M_\star > 2.5 \times 10^{10} M_\odot$) para mitigar efeitos de resolução.

O estudo procede em três estágios:

1. Ajuste de HOD: Os autores ajustam modelos HOD padrão (Zheng et al. 2005 para LRGs; Alam et al. 2020 para ELGs) aos dados de simulação para determinar como os parâmetros de HOD evoluem com o redshift e a força da gravidade.
2. Análise de Dependência Secundária: Os catálogos de galáxias são divididos em quartis baseados em propriedades secundárias do halo: densidade ambiental (definida pela sobredensidade em um anel esférico) e concentração do halo. Os parâmetros de HOD são ajustados independentemente para cada quartil para quantificar a dependência da ocupação de galáxias em relação a essas propriedades.
3. Quantificação do Viés de Montagem: Os autores geram catálogos de galáxias fictícios usando três prescrições de HOD: (i) o HOD básico baseado apenas em massa, (ii) um HOD dependente do ambiente e (iii) um HOD dependente da concentração. Eles comparam a função de correlação em grande escala ($\xi_{LS}$) desses mocks contra a "verdade fundamental" (ground truth) medida diretamente das simulações hidrodinâmicas para quantificar o viés de montagem residual.

Principais Contribuições e Resultados

• Sequência Evolutiva dos Efeitos de MG: O estudo identifica uma sequência evolutiva distinta em como a MG afeta a conexão galáxia-halo. À medida que a quinta força "desblinda" (unscreens) os halos ao longo do tempo, os modelos progridem através de regimes de blindagem inicial, acompanhamento (catch-up) e supressão (quenching). Essa sequência explica o comportamento não monotônico dos parâmetros de HOD (ex: $M_c$, $p_{max}$, $M_1$) através de diferentes modelos de gravidade e redshifts. Por exemplo, em modelos de MG mais fracos, a transição do estado blindado para o desblindado ocorre mais tarde, levando a comportamentos de HOD diferentes em comparação com modelos mais fortes no mesmo redshift.
• Efeito de Saturação: Os autores confirmam um efeito de "saturação" onde o impacto da MG na formação de galáxias e no agrupamento atinge um limite. Modelos de MG mais fortes (ex: F4.0, F4.5) eventualmente "alcançam" uns aos outros em termos de parâmetros de HOD, sugerindo que, para um critério de seleção de galáxias fixo, o universo não pode produzir indefinidamente mais galáxias; a MG principalmente acelera a linha do tempo de formação em vez de criar um excesso infinito de galáxias.
• Dependências Secundárias:
  • Ambiente: Existe uma forte dependência dos parâmetros de HOD na densidade ambiental tanto em modelos $\Lambda$CDM quanto em modelos MG. Em MG, isso é exacerbado pelo mecanismo de blindagem camaleão, que faz com que a desblindagem ocorra primeiro em ambientes de baixa densidade.
  • Concentração: Embora os parâmetros de HOD mostrem uma dependência significativa da concentração do halo, essa dependência não se traduz em uma redução significativa do viés de montagem quando usada para gerar mocks.
• Magnitude do Viés de Montagem:
  • Negligenciar o viés de montagem (usando um HOD básico baseado apenas em massa) subestima a força do agrupamento de galáxias em 10–20% em $\Lambda$CDM para $z \lesssim 1$.
  • Em modelos de MG, o efeito de AB é não monotônico e pode ser significativamente diferente de $\Lambda$CDM. Por exemplo, no modelo F5.5, o efeito de AB é mais forte do que em GR, enquanto em F4.0 e F4.5, pode ser mais fraco.
  • Crucialmente, a direção e a magnitude do efeito de AB em MG podem diferir de $\Lambda$CDM, potencialmente levando a uma identificação incorreta do modelo de gravidade se as linhas de processamento (pipelines) de análise padrão forem utilizadas.
• Estratégia de Mitigação: Estender o modelo HOD básico para incluir a densidade ambiental como uma variável secundária reduz o efeito de viés de montagem para 2–3% para $z \lesssim 0.5$ em todos os modelos (incluindo $\Lambda$CDM e vários cenários de $f(R)$). Em contraste, estender o modelo com a concentração do halo não reduz significativamente o efeito de AB.

Significância e Alegações
O artigo afirma que o modelo HOD básico é insuficiente para testes de precisão cosmológica envolvendo Gravidade Modificada, particularmente para Galáxias de Linha de Emissão (ELGs). Os autores argumentam que a complexa interação entre a blindagem de MG e a formação de galáxias cria assinaturas de viés de montagem que diferem qualitativa e quantitativamente da gravidade padrão.

A principal significância reside em demonstrar que incorporar a densidade ambiental no modelo HOD fornece um ponto de partida bem motivado e eficaz para minimizar os sistemas de viés de montagem. Isso é particularmente relevante para levantamentos futuros como DESI e Euclid, que sondam escalas onde o viés de montagem é forte. Os autores alertam que falhar em contabilizar esses efeitos poderia enviesar a inferência de parâmetros cosmológicos ou levar a conclusões incorretas sobre a natureza da gravidade. Eles concluem que, embora o HOD dependente do ambiente reduza o AB para um nível residual ($\lesssim 3\%$ em baixo $z$), trabalhos adicionais com simulações de maior resolução e escala são necessários para eliminar totalmente esses sistemas e testar extensões de HOD mais sofisticadas.