Assembly bias and the redshift evolution of intrinsic alignments for LRGs

Resumo Técnico: Viés de montagem e a evolução com o redshift dos alinhamentos intrínsecos para LRGs

Declaração do Problema
Os alinhamentos intrínsecos (IA) de galáxias constituem um contaminante sistemático astrofísico primário em levantamentos de lentes gravitacionais fracas, capazes de enviesar a inferência de parâmetros cosmológicos se não forem corretamente modelados. Os esforços de modelagem atuais, como o modelo de Alinhamento Não Linear (NLA), frequentemente assumem que a amplitude de alinhamento depende exclusivamente da massa do halo. No entanto, com a chegada de levantamentos de alta estatística como o Euclid e o LSST, dependências secundárias anteriormente ocultas por incertezas observacionais podem se tornar significativas. Especificamente, o impacto do viés de montagem (assembly bias) — a dependência das propriedades do halo em relação ao histórico de formação além da massa — e a evolução com o redshift do sinal de alinhamento requerem investigação utilizando simulações hidrodinâmicas de grande escala para garantir uma modelagem robusta para futuros estudos de cosmic shear.

Metodologia
Os autores utilizam a execução de caixa de tamanho 2,8 Gpc da suíte de simulações hidrodinâmicas cosmológicas FLAMINGO (Schaye et al. 2023), que inclui física bariônica e uma variação de feedback fiducial. O estudo foca em uma amostra de galáxias centrais do tipo Galáxia Vermelha Luminosa (LRG), definidas por halos com uma massa de matéria escura ligada $M_h > 10^{12} M_\odot$ e galáxias com mais de 300 partículas estelares ($\sim 2,22 \times 10^{11} M_\odot$).

Principais etapas metodológicas incluem:

1. Definição do Redshift de Formação: O redshift de formação ($z_f$) é definido como o redshift no qual um halo continha metade de sua massa atual ($0,5 M_{200\text{mean}}$). Esta definição é validada por sua forte correlação com a concentração do halo ($c_{200\text{crit}}$).
2. Estimativa de Correlação: Funções de correlação posição-posição ($gg$) e posição-forma ($g+$) são computadas usando estimadores generalizados de Landy-Szalay e modeladas usando o código IATheory dentro do arcabouço NLA.
3. Análise de Viés de Montagem: A amostra é agrupada por massa de halo e redshift de formação para isolar a dependência da amplitude de alinhamento ($A_1$) em $z_f$ para uma massa fixa.
4. Modelagem da Evolução com o Redshift: Os autores ajustam uma relação de lei de potência entre a amplitude de alinhamento e a massa do halo ($A_1 \propto M^{\beta_M}$) através de diferentes snapshots da simulação ($z=0$ a $z=1,5$).
5. Rastreamento de Galáxias: Galáxias centrais são rastreadas através de snapshots usando o Track_id do localizador de subhalos HBT-HERONS para desvendar a evolução com o redshift de um simples crescimento de massa. A amostra é selecionada em um redshift inicial ($z_i$) e acompanhada até $z=0$, com objetos removidos se fundirem ou caírem abaixo do limiar de partículas estelares.

Principais Resultados

• Dependência do Viés de Montagem: A amplitude de alinhamento depende não apenas da massa do halo, mas também do redshift de formação. Para uma massa de halo fixa, halos e suas galáxias centrais que se formaram mais cedo (maior $z_f$) exibem uma amplitude de alinhamento sistematicamente mais alta. Este efeito se manifesta como um deslocamento vertical na relação $A_1$–massa, enquanto a inclinação da lei de potência permanece amplamente inalterada.
• Evolução de Parâmetros com o Redshift: A amplitude da lei de potência massa-alinhamento ($\alpha_M$) aumenta constantemente com o redshift. A inclinação ($\beta_M$) diminui com o redshift até $z \sim 1$ e depois estabiliza. Os autores fornecem ajustes quadráticos empíricos para $\alpha_M(z)$ e $\beta_M(z)$ (Eqs. 13 e 14).
• Modelo Massa-Redshift: Um modelo combinado massa-redshift é proposto (Eqs. 12–14) para prever amplitudes de alinhamento para amostras de LRG dados seu redshift e massa média de halo.
• Análise de Rastreamento: Ao rastrear galáxias a partir de redshifts iniciais elevados ($z_i$), a amplitude de alinhamento aumenta conforme o redshift diminui, o que é consistente com o acréscimo de massa. No entanto, mesmo após a normalização pela previsão do modelo massa-redshift, um deslocamento sistemático permanece: galáxias rastreadas de $z_i$ mais altos retêm uma amplitude de alinhamento mais alta do que aquelas rastreadas de $z_i$ mais baixos. Isso confirma que o viés de montagem introduz um espalhamento irredutível na amplitude de alinhamento para uma dada massa e redshift.
• Desvios da Lei de Potência: O ajuste da lei de potência desvia-se em baixas massas ($< 10^{13} M_\odot$) e altas massas ($> 3 \times 10^{14} M_\odot$). Em baixas massas, os autores observam uma diminuição relativa à lei de potência (contrastando com alguns resultados de apenas gravidade), enquanto em altas massas, uma diminuição é observada, atribuída à formação contínua do halo.

Significância e Alegações
O artigo afirma que a modelagem de IA em levantamentos de lentes gravitacionais fracas não pode depender de intervalos de prior arbitrariamente pequenos para a amplitude de alinhamento, pois o viés de montagem introduz um espalhamento fundamental na amplitude para uma massa de halo fixa. Embora a forma funcional dos modelos existentes (como o NLA) possa não exigir modificação para acomodar o viés de montagem, as restrições de prior sobre o parâmetro de amplitude devem ser ampliadas para acomodar esta dependência secundária.

Além disso, os resultados complicam a implementação de modelos de alinhamento baseados em Distribuição de Ocupação de Halo (HOD) para simulações de apenas gravidade. Como os modelos HOD tipicamente atribuem propriedades baseadas apenas na massa, eles falham em capturar a dependência do tempo de formação, a menos que isso seja explicitamente incorporado. O estudo fornece diretrizes baseadas em simulações para um modelo de evolução com o redshift de IA, oferecendo ajustes empíricos para a dependência massa-redshift da amplitude de alinhamento para amostras do tipo LRG.

Limitações
Os autores observam que a resolução da simulação FLAMINGO é insuficiente para resolver a evolução dos discos galácticos (a fração de disco permanece $< 1\%$ até $z \sim 1,5$). Consequentemente, os resultados aplicam-se estritamente a amostras puras de LRG (elípticas) e não consideram a evolução de elípticas para discos ou os sinais de alinhamento de galáxias discoides. O modelo massa-redshift é válido para massas de halo entre $10^{13}$ e $10^{14,5} M_\odot$ e redshifts até $z \sim 1,5$.