Pressure-regulated feedback-modulated star formation as a subgrid model for galaxy formation simulations

Resumo Técnico: Formação Estelar Modulada por Feedback e Regulada por Pressão como um Modelo de Subgrade

Definição do Problema
Simulações cosmológicas de grande volume da formação de galáxias (ex: IllustrisTNG, EAGLE, SIMBA) reproduzem com sucesso propriedades estatísticas de populações de galáxias, mas não conseguem resolver as escalas físicas das nuvens moleculares gigantes (GMCs) ou do meio interestelar (ISM) multifásico. Consequentemente, essas simulações dependem de prescrições de subgrade empiricamente calibradas para taxas de formação estelar (SFR) e para a equação de estado (EoS) efetiva. Modelos tradicionais tipicamente assumem que a SFR específica escala com a densidade do gás como $\dot{m}_*/m_g \propto \rho^{1/2}$ (tempo de queda livre inverso), com uma normalização definida por observações do universo local.

Os autores identificam duas limitações primárias nesta abordagem:

1. Falta de Conexão Física: A escala empírica carece de um vínculo direto com os processos físicos que regulam a formação estelar, especificamente o equilíbrio entre gravidade e feedback estelar.
2. Sensibilidade à Resolução e ao Ambiente: Essas prescrições frequentemente falham em regimes onde a densidade do gás e a pressão são significativamente mais altas do que no universo local (ex: galáxias de alto redshift). Além disso, modelos que dependem da densidade volumétrica local são sensíveis à resolução numérica; se a escala de altura do disco de gás não for resolvida, a densidade volumétrica medida torna-se imprecisa, levando a previsões de formação estelar inconsistentes entre diferentes resoluções.

Metodologia
O artigo introduz um novo framework de subgrade baseado na teoria da formação estelar Regulada por Pressão e Modulada por Feedback (PRFM) (Ostriker & Kim 2022). Neste framework, a SFR ajusta-se para manter o equilíbrio térmico e dinâmico vertical no ISM, onde a pressão total do plano médio ($P_{tot}$) equilibra o peso do ISM ($W$) e é proporcional à SFR por unidade de área ($\Sigma_{SFR}$) via um rendimento de feedback ($\Upsilon_{tot}$).

Os autores implementam duas versões complementares deste modelo no código de malha móvel Arepo, calibradas usando simulações magnetohidrodinâmicas de alta resolução da suíte TIGRESS:

1. PRFM-vol (Volumétrico): Projetado para simulações onde a espessura vertical do disco de gás é numericamente resolvida (tipicamente resolução de massa $\sim 10^5 M_\odot$).

   • Ele calcula a pressão efetiva ($P_{eff}$) em cada célula usando uma EoS calibrada pelo TIGRESS ($P_{tot} \propto n_H^{1.8}$).
   • O tempo dinâmico ($t_{dyn}$) é derivado das densidades volumétricas locais de gás, estrelas e matéria escura.
   • O tempo de esgotamento ($t_{dep}$) é definido com base na pressão local e um fator de renormalização ($R_f \approx 2$) para contabilizar a integração vertical da lei volumétrica.

2. PRFM-int (Integrado): Projetado para resoluções mais grosseiras onde a escala de altura do disco é não resolvida (tipicamente resolução de massa $\gtrsim 10^6 M_\odot$).

   • Em vez de depender da densidade volumétrica local, ele reconstrói a densidade e a pressão do plano médio através de considerações de equilíbrio vertical usando densidades superficiais ($\Sigma_g$, $\Sigma_*$) e o peso do ISM ($W$).
   • Ele resolve a escala de altura de gás de equilíbrio ($H_g$) e deriva o $t_{dep}$ diretamente dessas quantidades integradas, evitando o uso de densidades volumétricas sensíveis à resolução na prescrição da SFR.

Principais Contribuições

• Calibração com Fundamentação Física: O modelo deriva a SFR e a EoS a partir da condição de equilíbrio PRFM ($W = P_{tot} = \Upsilon_{tot}\Sigma_{SFR}$), calibrada contra simulações TIGRESS em vez de ajustes empíricos a galáxias locais.
• Robustez à Resolução: O desenvolvimento da abordagem PRFM-int aborda especificamente a falha dos modelos volumétricos em baixa resolução ao utilizar densidades superficiais e restrições de equilíbrio, que são menos sensíveis à incapacidade de resolver a estrutura vertical do disco.
• Eficiência Variável: Ao contrário dos modelos tradicionais que assumem uma eficiência constante por tempo de queda livre, o modelo PRFM prevê que a eficiência de formação estelar por tempo dinâmico ($\epsilon_{dyn}$) aumenta com a pressão ($\epsilon_{dyn} \propto P_{tot}^{0.43}$).

Resultados
Os autores testaram ambas as implementações em simulações de disco isolado do tipo Via Láctea, através de resoluções de massa de $10^5$, $10^6$ e $10^7 M_\odot$, comparando-as com a prescrição padrão da IllustrisTNG.

• Concordância de Alta Resolução: Em $10^5 M_\odot$, tanto o PRFM-vol quanto o PRFM-int concordam estreitamente entre si e com os resultados do TIGRESS para SFRs e tempos de esgotamento. Ambos recuperam com sucesso a relação Ostriker-Kim calibrada pelo TIGRESS ($\Sigma_{SFR} \propto P_{tot}^{1.2}$).
• Sensibilidade à Resolução de Modelos Volumétricos: O PRFM-vol e os modelos TNG mostram sensibilidade significativa à resolução. À medida que a resolução torna-se mais grosseira para $10^6$ e $10^7 M_\odot$, a incapacidade de resolver a escala de altura do disco leva a medições imprecisas da densidade volumétrica. Isso faz com que o modelo PRFM-vol superestime os tempos de esgotamento e falhe em corresponder ao histórico de SFR de alta resolução.
• Robustez do PRFM-int: A implementação PRFM-int permanece robusta em todas as resoluções testadas. Ela rastreia com sucesso o histórico global de formação estelar e as distribuições de densidade superficial da simulação PRFM-vol de alta resolução ($10^5 M_\odot$), mesmo em $10^7 M_\odot$.
• Comparação com TNG: Comparado ao modelo TNG, os modelos PRFM preveem tempos de esgotamento de gás mais curtos, particularmente em regiões de alta pressão e alta densidade. Os modelos PRFM exibem uma gama dinâmica muito mais ampla de tempos de esgotamento e uma dependência mais íngreme da pressão, enquanto o modelo TNG mantém uma eficiência constante por tempo de queda livre.

Significância
O artigo afirma que as prescrições de subgrade derivadas do PRFM fornecem um framework fisicamente fundamentado e numericamente estável para a formação estelar através da gama dinâmica de simulações de formação de galáxias.

• Para Simulações de Alta Resolução: O modelo PRFM-vol oferece um tratamento mais realista da física do ISM do que as leis empíricas tradicionais, ligando a SFR diretamente ao equilíbrio dinâmico da galáxia.
• Para Simulações Cosmológicas de Grande Escala: O modelo PRFM-int resolve o problema crítico da dependência de resolução. Ele permite que simulações cosmológicas de grandes volumes (que tipicamente não conseguem resolver as escalas de altura de disco) utilizem leis de formação estelar motivadas fisicamente que permanecem consistentes com a física de alta resolução.
• Implicações para Alto Redshift: Os autores observam que a eficiência dependente da pressão do modelo PRFM sugere que a formação estelar em galáxias de alto redshift (onde as pressões são mais altas) pode ser significativamente mais eficiente do que o previsto pelos atuais modelos do tipo TNG, potencialmente impactando a interpretação de galáxias luminosas observadas pelo JWST.

O trabalho é apresentado como um passo em direção ao objetivo da colaboração "Learning the Universe" de desenvolver modelos de subgrade preditivos, em vez de empiricamente ajustados, que possam ser atualizados conforme as simulações de ISM de alta resolução (como o TIGRESS) expandem para cobrir espaços de parâmetros mais amplos (ex: metalicidade variável).