Euclid Quick Data Release (Q1). The average far-infrared properties of Euclid-selected star-forming galaxies

Resumo Técnico: As Propriedades Médias de Infravermelho Distante das Galáxias Formadoras de Estrelas Selecionadas pelo Euclid

Problema e Motivação
O primeiro Lançamento Rápido de Dados (Q1) da missão Euclid fornece um conjunto de dados de óptico e infravermelho próximo (NIR) sem precedentes, contendo milhões de galáxias com redshifts fotométricos robustos, massas estelares e taxas de formação estelar (SFRs) derivadas do ajuste de densidade espectral de energia (SED). No entanto, essas SFRs derivadas de óptico/NIR dependem de suposições sobre a atenuação de poeira e a população estelar subjacente. Para restringir independentemente as propriedades médias de infravermelho distante (far-IR) dessas galáxias e verificar a consistência das estimativas de SFR, este estudo investiga a emissão térmica de grãos de poeira. Embora análises de empilhamento (stacking) anteriores de dados de far-IR (por exemplo, nos campos UDS, COSMOS ou GAMA) tenham sido limitadas por pequenas áreas de céu ou intervalos de redshift, o catálogo Euclid Q1 permite o empilhamento de milhões de galáxias através de dezenas de graus quadrados, permitindo uma caracterização estatisticamente robusta da sequência principal (MS) de galáxias formadoras de estrelas até $z \approx 3$.

Metodologia
Os autores utilizaram o catálogo de fusão Euclid Q1 (especificamente a amostra descrita por Euclid Collaboration: Enia et al. 2025), selecionando galáxias formadoras de estrelas dentro do intervalo de redshift $0.2 < z < 3.0$ e intervalo de massa estelar $8.3 < \log_{10}(M_*/M_\odot) < 11.5$. A amostra foi cruzada com imagens de far-IR existentes do Observatório Espacial Herschel (PACS em 100/160 $\mu$m; SPIRE em 250/350/500 $\mu$m) e do instrumento SCUBA-2 (450/850 $\mu$m). A análise cobriu uma área de sobreposição total de 17,6 deg$^2$ para o SPIRE e 2,4 deg$^2$ para o SCUBA-2.

Para abordar a resolução angular grosseira dos instrumentos de far-IR (10–30 arcsegundos) em comparação com a resolução do Euclid (~0,2 arcsegundos), os autores empregaram uma análise de empilhamento usando o algoritmo SimStack. Este método resolve simultaneamente as densidades de fluxo médias de múltiplos bins definidos por massa estelar e redshift. Características metodológicas principais incluíram:

• Implementação do SimStack: Resolução de um sistema linear $y = S_1 X_1 + \dots + S_N X_N$ para determinar as amplitudes de empilhamento ($S_b$) para $N$ bins, contabilizando o agrupamento (clustering) de galáxias e correlações entre bins sem assumir uma distribuição de Poisson.
• Mascaramento: Mascaramento rigoroso de regiões contaminadas por estrelas brilhantes, efeitos de borda e fontes galácticas (ex: NGC 6543) para evitar viés no fluxo empilhado.
• Ajuste de SED: As densidades de fluxo empilhadas foram ajustadas com um modelo de corpo negro modificado transitando para uma lei de potência em altas frequências. O índice de emissividade de poeira ($\beta$) e a inclinação da lei de potência ($\alpha$) foram fixados em valores derivados de galáxias locais selecionadas em far-IR ($\beta=1.96, \alpha=2.3$) devido às limitações de restrição fotométrica.
• Derivação de Parâmetros: A partir dos melhores ajustes de SED, os autores derivaram temperaturas de poeira médias ($T_d$), massas de poeira ($M_d$) e SFRs derivadas de far-IR.

Resultos Principais

1. Consistência de SFR: As SFRs derivadas de far-IR mostram boa concordância com as SFRs preditas pela sequência principal (MS) de formação estelar derivada de fotometria de óptico/NIR, com uma razão média ponderada de 1,0.
2. Evolução da Temperatura de Poeira: A temperatura média de poeira ($T_d$) diminui com o tempo cósmico (aumento do redshift). A evolução é bem descrita pela função $T_d(t) = T_2 + (T_1 - T_2)e^{-t/\tau}$, onde $T_1 = 79,7 \pm 7,4$ K (extrapolado para $t=0$), $T_2 = 23,2 \pm 0,1$ K (valor assintótico), e $\tau = 1,6 \pm 0,1$ Gyr.
   • Crucialmente, $T_d$ é encontrado como independente da massa estelar em todos os bins estudados.
   • A convergência para um valor não nulo ($T_2 \approx 23$ K) abaixo de $z=1$ sugere que a poeira em galáxias da MS é cada vez mais aquecida pela população estelar existente, mais fria e antiga, em vez de apenas por estrelas jovens e quentes em regiões de formação estelar.
3. Massa de Poeira e Razões: As massas de poeira ($M_d$) e as razões de poeira para massa estelar ($M_d/M_*$) aumentam de $z=0,2$ até $z \approx 1$, após o qual estabilizam ou declinam. A razão poeira-massa estelar diminui com o aumento da massa estelar em todos os redshifts, atribuída a menores sSFRs (taxas específicas de formação estelar) em galáxias mais massivas.
4. Fundo Cósmico de Infravermelho (CIB): O catálogo Euclid Q1 resolve mais de 60% do CIB em 250, 350 e 500 $\mu$m. Em 100 $\mu$m, a fração resolvida é estimada em 30–80%, enquanto em 850 $\mu$m é de 30–40%.
5. Comparação com Simulações: Quando comparado ao catálogo mock MAMBO (baseado na Simulação Millennium e no Gerador Empírico de Galáxias), as massas de poeira e SFRs simuladas concordam razoavelmente bem. No entanto, a simulação prevê uma diminuição monotônica da temperatura de poeira com o tempo, falhando em reproduzir o platô observado em baixos redshifts ($z < 1$).

Significância e Alegações
O artigo afirma que este estudo representa a maior análise de empilhamento deste tipo até o momento, aproveitando a profundidade e área únicas do lançamento Q1 do Euclid. A principal significância reside na descoberta de que a temperatura média de poeira das galáxias de formação estelar converge para um valor constante em baixos redshifts, independente da massa estelar. Esta descoberta desafia modelos onde a temperatura da poeira é estritamente proporcional à SFR ou sSFR, sugerindo uma transição no mecanismo de aquecimento da poeira de estrelas jovens e massivas para o campo de radiação estelar geral de populações estelares mais antigas conforme o universo envelhece.

Além disso, os autores demonstram que a correlação entre a temperatura de poeira e a SFR é mediada pela massa estelar, visto que a própria $T_d$ é independente da massa. Os resultados também fornecem um parâmetro de referência para futuros lançamentos de dados do Euclid, indicando que será necessária uma imagem mais profunda para resolver totalmente o CIB e sondar as propriedades de far-IR de galáxias em redshifts mais altos ($z > 1,5$) e massas estelares menores. O estudo conclui que, embora os dados atuais do Q1 não resolvam todo o CIB, eles explicam a maior parte do fundo em comprimentos de onda do SPIRE, validando o potencial do Euclid para mapear a história térmica da população formadora de estrelas do universo.