pop-cosmos: Galaxy size evolution across structural and star-formation classifications in COSMOS-Web

Resumo Técnico: pop-cosmos: Evolução do tamanho das galáxias através de classificações estruturais e de formação estelar em COSMOS-Web

Problema e Motivação
Sabe-se que o tamanho das galáxias correlaciona-se com a massa estelar e o redshift, caracterizado por relações de escala de tamanho. No entanto, as formas inferidas dessas relações são sensíveis à maneira como as galáxias são classificadas. Tradicionalmente, estudos separam galáxias por atividade de formação estelar (ex: taxa específica de formação estelar, sSFR) ou por marcadores morfológicos (ex: razão bojo-total, índice de Sérsic), tratando frequentemente essas classificações como equivalentes. Uma questão fundamental não resolvida é se esses esquemas identificam as mesmas populações estruturais ou se capturam dimensões diferentes da evolução galáctica. Além disso, trabalhos anteriores basearam-se em seleções de cor no referencial de repouso (ex: $UVJ$), que podem sofrer contaminação por galáxias formadoras de estrelas com poeira, e em imagens baseadas em solo que podem borrar detalhes estruturais. Este artigo investiga se as classificações baseadas em sSFR e em morfologia produzem relações de tamanho–massa e tamanho–redshift quantitativamente diferentes, e se essas diferenças são recuperáveis apenas com resolução baseada no espaço.

Metodologia
Os autores combinam dois conjuntos de dados distintos para o campo COSMOS para construir uma amostra de 99.369 galáxias:

1. Parâmetros Físicos: Massa estelar ($M_*$), redshift ($z$) e estimativas de sSFR são derivados do modelo generativo pop-cosmos (Thorp et al. 2025). Este modelo é treinado em fotometria de 26 bandas do COSMOS2020 (Weaver et al. 2022) e utiliza um modelo de difusão baseado em pontuação (score-based) para aprender a distribuição conjunta de 16 parâmetros de energia espectral (SED), evitando formas funcionais restritivas.
2. Tamanho e Morfologia: Raios de meia luz ($R_{50}$), índices de Sérsic ($n$) e razões bojo-total ($B/T$) são obtidos do catálogo COSMOS2025 (Shuntov et al. 2025), baseado em imagens de alta resolução do JWST/NIRCam do levantamento COSMOS-Web.

A amostra de análise é cruzada com uma abertura de 0,5" e sujeita a cortes de qualidade para remover estrelas, fontes fundidas e medições de tamanho não confiáveis. A amostra é restrita a $0,0 < z < 3,0$.
Os autores ajustam as relações tamanho–massa e tamanho–redshift usando inferência Bayesiana com verossimilhanças Gaussianas que contabilizam incertezas de medição e dispersão intrínseca ($\sigma_{int}$).

• Populações formadoras de estrelas/dominadas por disco/tipo-latente são ajustadas com uma lei de potência única: $R_{50} = R_0 (M_*/5\times10^{10} M_\odot)^B$.
• Populações quiescentes/dominadas por bojo/tipo-precoce são ajustadas com uma lei de potência dupla para capturar potenciais quebras: $R_{50} = \gamma (M_*/M_\odot)^\alpha (1 + M_*/10^\delta M_\odot)^{\beta-\alpha}$.
• Relações tamanho–redshift são modeladas como $R_0 = B_z (1+z_m)^{-\beta_z}$.

As classificações são aplicadas usando:

• sSFR: Quiescente ($\log_{10}(\text{sSFR}/\text{yr}^{-1}) \le -11$) vs. Formadora de estrelas.
• B/T: Dominada por disco ($B/T \le 0,2$), Intermediária, Dominada por bojo ($B/T \ge 0,6$).
• Índice de Sérsic: Tipo-latente ($n \le 2,5$) vs. Tipo-precoce ($n > 2,5$).

Principais Resultados

1. Diferenças Quantitativas nas Relações de Escala: As divisões por sSFR e morfologia produzem inclinações, interceptos e comportamentos de dispersão intrínseca quantitativamente diferentes.

   • Inclinações: Sistemas formadores de estrelas/dominados por disco seguem inclinações mais rasas ($B \sim 0,12\text{--}0,27$), enquanto sistemas quiescentes/dominados por bojo exibem inclinações mais íngremes em altas massas ($B \sim 0,46\text{--}0,62$).
   • Dispersão Intrínseca: Crucialmente, $\sigma_{int}$ depende da morfologia estrutural, mas não da sSFR. Sistemas dominados por disco mostram menor dispersão ($\sigma_{int} \sim 0,13\text{--}0,14$ dex) comparados a sistemas dominados por bojo ($\sigma_{int} \sim 0,18\text{--}0,20$ dex). Em contraste, populações formadoras de estrelas e quiescentes mostem dispersão semelhante ($\sim 0,15$ dex).
   • Quebras: Populações quiescentes e dominadas por bojo mostram quebras de lei de potência dupla, mas em massas de pivô diferentes. A massa de pivô quiescente é $\delta_Q \sim 10^{10,7} M_\odot$, enquanto a do bojo é $\delta_B \sim 10^{11,1} M_\odot$.

2. Desacoplamento do Encerramento (Quenching) e Transformação Estrutural: O deslocamento entre a massa de pivô quiescente ($\delta_Q$) e a massa de pivô dominada por bojo ($\delta_B$) sugere que o encerramento da formação estelar (quenching) e a transformação morfológica ocorrem em escalas de tempo diferentes. O quenching precede a transformação morfológica.

3. Tendências Dependentes da Morfologia Requerem Imagens Espaciais: As tendências estruturais identificadas (especificamente a correlação entre $B/T$, massa e redshift) são suavizadas em imagens baseadas em solo (HSC-I), mas são claramente recuperadas nos dados do JWST/COSMOS-Web. Medições baseadas em solo falham em resolver a diversidade estrutural dentro das classes de sSFR, particularmente em baixo redshift, onde uma fração significativa de galáxias quiescentes é dominada por disco.

4. Conexão com o Feedback de AGN:

   • A massa de pivô quiescente ($\sim 10^{10,7} M_\odot$) coincide com a escala de massa onde a fração de luminosidade bolométrica do toro de infravermelho do AGN atinge o pico em galáxias em transição e onde o feedback de AGN torna-se eficaz em halos ($M_{200} \sim 10^{12} M_\odot$).
   • A massa de pivô dominada por bojo ($\sim 10^{11,1} M_\odot$) coincide com a massa do halo ($M_{200} \sim 10^{12,7} M_\odot$) onde a redistribuição bariônica impulsionada por AGN atinge o pico, criando condições para o aquecimento de discos e transformação estrutural.
   • Em $z < 0,5$, o pivô da lei de potência dupla torna-se não restringido para as galáxias mais massivas, sugerindo que nos halos de maior massa ($M_{200} \gtrsim 10^{13,5} M_\odot$), o feedback de AGN pode falhar em prevenir a reacumulação de gás, levando a uma quebra na transição aguda de quenching.

Significância e Alegações
O artigo alega que as relações de escala de tamanho das galáxias codificam uma sequência de limiares fisicamente distintos, em vez de uma única transição evolutiva. A principal significância reside em demonstrar que:

• Esquemas de Classificação Não São Intercambiáveis: sSFR, $B/T$ e o índice de Sérsic capturam diferentes estágios da sequência "quenching-então-transformação". Tratá-los como equivalentes obscurece o desacoplamento dependente da massa desses processos.
• Dispersão Intrínseca como Sonda de Halo: A dependência da dispersão intrínseca na morfologia (mas não na sSFR) implica que o tamanho da galáxia está conectado ao spin do halo primariamente através da configuração estrutural (ex: conservação de momento angular em discos vs. processos seculares em buljos), e não através da atividade de formação estelar.
• Necessidade de Resolução Baseada no Espaço: Recuperar as tendências dependentes da morfologia e a verdadeira diversidade estrutural dentro das classes de sSFR requer a resolução fornecida por imagens espaciais (JWST), uma vez que os dados baseados em solo borram essas distinções.
• Evolução em Dois Estágios: Os resultados apoiam um quadro de dois estágios onde o feedback de AGN desencadeia o quenching em uma escala de massa específica, seguido por uma transformação estrutural atrasada em sistemas dominados por bojo, impulsionada por evolução secular e aquecimento de disco, com o ritmo modulado pelo ambiente.

Os autores concluem que o uso de um modelo direto de nível populacional (pop-cosmos) combinado com morfologia de alta resolução baseada no espaço fornece um arcabouço coerente para desmembrar esses caminhos evolutivos.