Intrinsic alignments in the FLAMINGO simulations with two-point statistics

Resumo Técnico: Alinhamentos Intrínsecos nas Simulações FLAMINGO com Estatísticas de Dois Pontos

Definição do Problema
Os alinhamentos intrínsecos (IA) de galáxias constituem um grande contaminante sistemático astrofísico para levantamentos de lentes gravitacionais fracas de próxima geração, como o Euclid e o Legacy Survey of Space and Time (LSST). O IA induz correlações nas formas das galáxias que mimetizam o sinal de cisalhamento cósmico (cosmic shear), podendo enviesar a estimativa de parâmetros cosmológicos se não for modelado corretamente. Embora os dados observacionais tenham melhorado, compreender as limitações das estratégias atuais de modelagem de IA e estabelecer priors sensatos requer simulações hidrodinâmicas de grande escala. Essas simulações devem possuir volume suficiente para corresponder às escalas dos levantamentos e resolução suficiente para resolver efeitos de pequena escala que influenciam o campo de maré. Além disso, a dependência do IA em relação à física bariônica (feedback) permanece uma questão em aberto, uma vez que os modelos de alinhamento existentes não contabilizam atualmente as variações de feedback.

Metodologia
Os autores utilizam a suíte de simulações hidrodinâmicas FLAMINGO, especificamente uma caixa de 2,8 Gpc$^3$ e uma suíte de execuções de 1 Gpc$^3$ com variadas implementações de feedback ( feedback de AGN térmico vs. jato, e variações na força do feedback de supernova). O estudo foca em uma amostra de $\approx$4,9 milhões de análogos de Galáxias Elípticas Vermelhas (LRG) em redshift $z=0$, selecionados de subhalos com massas de matéria escura ligadas $>10^{12} M_\odot$ e $>300$ partículas estelares.

Os autores medem estatísticas de dois pontos:

1. Correlação Posição-Posição ($w_{gg}$): Rastreando o agrupamento de galáxias (clustering).
2. Correlação Posição-Forma ($w_{g+}$): Rastreando o sinal de alinhamento intrínseco.

Essas medições são realizadas utilizando os estimadores generalizados de Landy-Szalay e o código IACorr. Os autores modelam conjuntamente os sinais de agrupamento e alinhamento usando o viés (bias) não linear para $w_{gg}$ e dois modelos de IA para $w_{g+}$:

• Alinhamento Não Linear (NLA): Um modelo de alinhamento linear estendido para escalas quase-lineares usando o espectro de potência da matéria não linear, caracterizado pela amplitude $A_1$ e viés $b_1$.
• Alinhamento de Maré e Torque de Maré (TATT): Um modelo mais complexo que incorpora termos lineares e quadráticos, caracterizado pelas amplitudes $A_1$, $A_2$ e $A_{1\delta}$.

A análise emprega uma abordagem Bayesiana usando o amostrador aninhado (nested sampler) nautilus para restringir os parâmetros. Os autores também investigam a dependência desses parâmetros em relação à massa do halo e testam a sensibilidade do sinal de alinhamento a variações no feedback de AGN e de supernova, comparando diferentes execuções da simulação FLAMINGO.

Principais Contribuições e Resultados

• Desempenho do Modelo: Tanto os modelos NLA quanto o TATT fornecem bons ajustes aos dados da FLAMINGO. O modelo TATT produz resíduos menores em pequenas escalas ($r < 10$ Mpc/$h$) em comparação ao NLA, permitindo um corte de escala mínima inferior ($5$ Mpc/$h$ vs. $10$ Mpc/$h$) mantendo um $\chi^2_{red}$ comparável ao NLA. Os resíduos entre os modelos e os dados estão dentro de alguns por cento nas escalas consideradas ($5-100$ Mpc/$h$).
• Comparação com Observações: A amplitude NLA ($A_1$) medida da amostra FLAMINGO mostra boa concordância com vários estudos observacionais de LRGs (ex: KiDS, DES, SDSS) através de uma gama de luminosidades. A simulação estende o intervalo medido para luminosidades mais altas do que as anteriormente sondadas pelas observações.
• Dependência de Massa: Os autores confirmam que a amplitude de alinhamento $A_1$ escala com a massa do halo como uma lei de potência. Crucialmente, eles apresentam a primeira medição da dependência de massa para os parâmetros de ordem superior do TATT ($A_2$ e $A_{1\delta}$). Eles descobrem que $A_2$ é aproximadamente um múltiplo constante de $A_1$, e a razão $A_{1\delta}/A_1$ segue uma dependência de lei de potência em relação à massa do halo.
• Introdução do TATT-M: Com base nas relações empíricas encontradas entre os parâmetros TATT e a massa do halo, os autores introduzem um modelo de TATT dependente de massa, o TATT-M. Este modelo reduz o espaço de parâmetros do TATT ao expressar $A_2$ e $A_{1\delta}$ como funções de $A_1$ e da massa do halo (derivadas do ajuste da amostra de halos).
  • A comparação de modelos Bayesianos mostra que o TATT-M é fortemente preferido em relação ao modelo NLA padrão para a amostra de galáxias centrais.
  • O uso do TATT-M restringe significativamente as restrições sobre $A_1$ (reduzindo o erro por um fator de $\approx 3,7$) em comparação ao modelo TATT padrão, devido aos priors informativos derivados da amostra de halos da simulação.
• Sensibilidade ao Feedback: Os autores testam a sensibilidade dos parâmetros de IA ao feedback bariônico usando as variações de feedback da FLAMINGO.
  • Comparações iniciais de valores brutos de $A_1$ entre as execuções de feedback mostraram variações, particularmente para execuções com funções de massa estelar deslocadas ('M*-$\sigma$') e feedback de jato forte.
  • No entanto, os autores demonstram que essas variações são primariamente impulsionadas por mudanças na massa média do halo da amostra de galáxias selecionada, resultante do corte fixo de partículas estelares, e não por um efeito físico direto do feedback no mecanismo de alinhamento em si.
  • Quando a análise é feita em bins de massa de halo para remover a dependência de massa, todas as variações de feedback (incluindo térmico, jato e variações de força de supernova) produzem amplitudes e inclinações de alinhamento consistentes dentro das incertezas estatísticas. O artigo conclui que, dentro dos limites de resolução da FLAMINGO, o feedback não altera significativamente o sinal de alinhamento além de seu efeito na distribuição de massa estelar das galáxias.

Significância
Este trabalho fornece orientações críticas para a modelagem de alinhamentos intrínsecos em futuros levantamentos cosmológicos. Ao validar os modelos NLA e TATT contra dados hidrodinâmicos de alta resolução, os autores confirmam que esses modelos são suficientes para análises de estatísticas de dois pontos até escalas de $5$ Mpc/$h$. A introdução do modelo TATT-M oferece um caminho para melhorar o poder de restrição de futuras análises de lentes gravitacionais fracas ao utilizar priors derivados de simulações sobre termos de alinhamento de ordem superior. Além disso, os resultados sugerem que, para amostras do tipo LRG, o impacto do feedback bariônico nos alinhamentos intrínsecos é amplamente mediado por mudanças na massa estelar das galáxias, implicando que as estratégias de modelagem atuais podem ser robustas contra incertezas na física do feedback, desde que as dependências de massa sejam corretamente contabilizadas.