Intrinsic galaxy alignments in CAMELS simulations and the significant impact of baryon model

Resumo Técnico: Alinhamentos Intrínsecos de Galáxias em Simulações CAMELS e o Impacto Significativo do Modelo de Bárions

Definição do Problema
Os alinhamentos intrínsecos (IA) de galáxias servem tanto como uma sonda da formação de galáxias quanto como uma importante fonte de erro sistemático em análises de lentes gravitacionais fracas. A inferência cosmológica atual de levantamentos fotométricos é limitada por incertezas na modelagem de IA, particularmente no que diz respeito a como os processos de feedback bariônico (supernovas e núcleos galácticos ativos - AGN) e os parâmetros cosmológicos influenciam os sinais de alinhamento. Embora existam estruturas teóricas (por exemplo, alinhamento de maré para galáxias do tipo early-type e torque de maré para late-type), a dependência específica da IA em relação a uma ampla gama de parâmetros astrofísicos e cosmológicos permanece mal restringida, especialmente ao distinguir entre diferentes populações de galáxias (formação estelar vs. quiescentes).

Metodologia
Os autores utilizam a suíte CAMELS (Cosmology and Astrophysics with MachinE Learning Simulations), especificamente o conjunto Latin-Hypercube (LH) de 1000 simulações hidrodinâmicas utilizando o código IllustrisTNG. Estas simulações variam seis parâmetros fundamentais: parâmetros cosmológicos ($\Omega_m$, $\sigma_8$) e parâmetros de feedback astrofísico ($A_{SN1}$, $A_{SN2}$ para supernovas; $A_{AGN1}$, $A_{AGN2}$ para AGN).

• Construção de Dados: As formas das galáxias são medidas em $z=0$ a partir de subhalos identificados pelo algoritmo SUBFIND, selecionados com base na massa estelar ($M_{2r} > 10^8 M_\odot$) e contagem de partículas ($N_s > 10$). As elipticidades são calculadas a partir do tensor de segundo momento da densidade estelar em três projeções cartesianas, gerando 3.000 catálogos.
• Medidas Estatísticas: O estudo mede a função de correlação cruzada projetada entre as elipticidades das galáxias e o campo de densidade de matéria escura subjacente, denotada por $w_{m+}(r_\perp)$. Um teste nulo é realizado utilizando a correlação de componente cruzada $w_{m\times}$.
• Subamostras: As galáxias são divididas em populações "em formação estelar" (star-forming) e "quiescentes" com base em um limiar de taxa de formação estelar específica (sSFR) de $10^{-10.5} \text{ yr}^{-1}$.
• Análise de Orientação Apenas: Para isolar o efeito da orientação da magnitude da elipticidade, os autores constroem catálogos de "orientação apenas", onde todas as magnitudes de elipticidade são normalizadas para o valor médio $\langle|\varepsilon|\rangle$, calculando uma função de correlação correspondente $\tilde{w}_{m+}$.
• Dependência de Parâmetros: As amplitudes de alinhamento são definidas como médias ponderadas de $w_{m+}$ através de bins de separação transversal. Coeficientes de correlação de Pearson e p-valores são computados para testar dependências lineares entre estas amplitudes e os parâmetros da simulação. A razão $w_{m+}/w_{mm}$ também é analisada para testar a consistência com modelos de alinhamento linear (NLA/TATT).

Principais Resultados

1. Detecção e Diferenças de População: Um sinal significativo de alinhamento intrínseco é detectado nas simulações CAMELS. As galáxias quiescentes exibem amplitudes de alinhamento aproximadamente uma ordem de magnitude maiores do que as galáxias em formação estelar.
2. Dependências de Parâmetros:
   • Cosmologia: A amplitude de alinhamento $w_{m+}$ mostra correlações estatisticamente significativas com $\Omega_m$ e $\sigma_8$. No entanto, quando o sinal é normalizado pela autocorrelação da densidade de matéria ($w_{mm}$), a dependência de $\sigma_8$ desaparece, o que é consistente com as previsões de modelos de alinhamento linear.
   • Feedback: Correlações significativas são encontradas entre os parâmetros de feedback de supernova ($A_{SN1}$ e $A_{SN2}$) e as amplitudes de alinhamento. Nenhuma dependência significativa em relação aos parâmetros de feedback de AGN é observada, o que os autores atribuem ao volume limitado da simulação ($(25 h^{-1} \text{ Mpc})^3$), que restringe o alcance dos modos de grande escala.
3. Mecanismos de Formação Estelar vs. Quiescentes:
   • A correlação entre a amplitude de alinhamento e $\sigma_8$ é insignificante para galáxias em formação estelar, mas forte para amostras quiescentes e totais. Os autores sugerem que a dependência na amostra total é impulsionada pelo aumento da fração de galáxias quiescentes altamente alinhadas em valores mais altos de $\sigma_8$.
   • O feedback de supernova ($A_{SN1}$) mostra uma correlação mais forte com o alinhamento em galáxias em formação estelar do que em galáxias quiescentes. Inversamente, em galáxias quiescentes, $A_{SN1}$ regula primariamente as taxas de formação estelar específica em vez do alinhamento.
4. Ellipticidade vs. Orientação: A normalização das magnitudes de elipticidade (análise de orientação apenas) mantém a sensibilidade aos parâmetros de feedback de supernova. Isso indica que o feedback de supernova influencia não apenas a magnitude da elipticidade da galáxia, mas também a orientação das formas das galáxias em relação à distribuição de matéria escura.
5. Análise de Razão: A razão $w_{m+}/w_{mm}$ remove a dependência de $\sigma_8$, mas retém correlações significativas com $\Omega_m$ e parâmetros de feedback de supernova. Os autores observam uma surpreendente correlação negativa com $\Omega_m$ nesta razão, o que contradiz as previsões simples de NLA/TATT e requer investigação adicional.

Significância e Alegações
O artigo afirma fornecer a primeira medição direta da dependência dos alinhamentos intrínsecos de galáxias em relação aos parâmetros cosmológicos e astrofísicos sobre uma ampla e realista gama utilizando um grande conjunto de simulações hidrodinâmicas. Os resultados oferecem insights fundamentais para:

• Modelagem: Demonstrar que o feedback bariônico (especificamente supernovas) impacta significativamente os sinais de alinhamento, sugerindo que os modelos devem considerar estes processos para evitar vieses na estimativa de parâmetros cosmológicos.
• Evolução de Galáxias: Destacar mecanismos de alinhamento distintos operando através de tipos de galáxias, onde o feedback de supernova impulsiona mudanças de orientação em galáxias em formação estelar enquanto regula as taxas de formação estelar em galáxias quiescentes.
• Levantamentos Futuros: Fornecer restrições necessárias para a validação de pipelines e inferência baseada em simulações para levantamentos de lentes gravitacionais de próxima geração, enfatizando que desajustes entre os modelos de feedback bariônico assumidos e as prescrições de IA podem introduzir vieses significativos.

Os autores permanecem modestos em relação aos resultados de feedback de AGN, declarando explicitamente que a falta de impacto detectado é provavelmente devida ao pequeno volume da simulação, e não afirmam ter resolvido o mecanismo específico por trás da inesperada dependência de $\Omega_m$.