Resolving Star Cluster Formation in Galaxy Simulations with Cosmic Ray Feedback

Resumo Técnico: Resolução da Formação de Aglomerados Estelares em Simulações Galácticas com Feedback de Raios Cósmicos

Definição do Problema
Os aglomerados estelares são os principais locais de formação de estrelas massivas, que impulsionam processos de feedback que moldam o meio interestelar (ISM) e regulam a formação estelar futura. Uma compreensão abrangente do feedback galáctico requer a resolução da formação de aglomerados em escala de parsec e do ISM multifásico. Sabe-se que os raios cósmicos (CRs) contribuem para a pressão e o aquecimento do ISM, mas seu impacto específico nas taxas de formação de aglomerados (SFRs), nas funções de massa de aglomerados e na ligação gravitacional permanece incerto. Estudos anteriores exploraram o feedback de CRs no contexto de fluxos galácticos ou utilizaram modelos de subgrade para aglomerados, mas poucos resolveram simultaneamente partículas estelares individuais, CRs acoplados dinamicamente e um ISM multifásico regulado por feedback. Este trabalho aborda a lacuna no entendimento de como os modelos de transporte de CR influenciam a formação e as propriedades dos aglomerados estelares.

Metodologia
Os autores analisam simulações de alta resolução do tipo "tallbox" realizadas utilizando o código de malha móvel Arepo. As simulações empregam o modelo de termoquímica não-equilibrada e feedback Crisp, que rastreia processos de resfriamento/aquecimento para hidrogênio e metais. O estudo utiliza três casos de física distintos para isolar os efeitos do transporte de CR:

1. MHD: Magnetohidrodinâmica com CRs tratados como um termo de aquecimento/ionização uniforme.
2. CR-NL: Magnetohidrodinâmica de dois momentos para CR (CRMHD) com amortecimento de Landau não linear (NLLD) como o único mecanismo de amortecimento, resultando em um acoplamento uniforme entre CR e gás.
3. CR-NL-IN (Fiducial): CRMHD com ambos, NLLD e amortecimento íon-neutro (IND). Neste caso, os CRs se desacoplam do gás neutro frio, reduzindo seu impacto local na formação estelar.

As simulações resolvem o gás até uma resolução de massa de $10 M_\odot$ e uma escala espacial de $\sim 0,5$ pc. A formação estelar é implementada via uma prescrição do tipo Schmidt, onde células de gás que excedem um limiar de densidade ($n_H = 10^3 \text{ cm}^{-3}$) são convertidas em partículas estelares colisionais com uma eficiência local de 100% por tempo de free-fall. O feedback estelar é modelado exclusivamente através de supernovas (SNe); o feedback radiativo e ventos estelares são omitidos para isolar o impacto de CRs e SNe. As SNe são injetadas estocasticamente para estrelas com idades entre 3–38 Myr, com uma adição de 5% da energia de SN injetada como CR nos casos CRMHD.

Os aglomerados estelares são identificados usando um algoritmo Friends-of-Friends (FoF) modificado que incorpora uma coordenada temporal, garantindo que os grupos identificados consistam de estrelas coevales dentro de um comprimento de ligação espacial de 2 pc. A análise foca em aglomerados "jovens" (idades $\sim 10$ Myr) para caracterizar suas propriedades iniciais antes de uma disrupção tidal significativa.

Resultos Principais

1. Funções de Massa de Aglomerados: Todos os três casos de simulação produzem funções de massa de aglomerados que seguem uma lei de potência ($dN/dM \propto M^{\alpha_{PL}}$) com inclinações entre $-2,0$e$-2,4$ para massas $> 10^3 M_\odot$. Isso coincide com tendências observacionais onde a inclinação correlaciona-se positivamente com a densidade superficial da taxa de formação estelar ($\Sigma_{SFR}$). A inclusão de CRs desloca os pontos de dados ao longo dessas relações observadas, mas não produz resultados não físicos.
2. Eficiência de Formação de Aglomerados (CFE): A CFE ($\Gamma$) mostra uma correlação positiva com $\Sigma_{SFR}$ quando calculada para aglomerados massivos ($M > 5000 M_\odot$), consistente com previsões teóricas (ex: Kruijssen 2012) e observações. No entanto, para limiares de massa mais baixos ($M > 1000 M_\odot$), os casos CRMHD mostram uma CFE alta e constante, enquanto o caso MHD mostra valores menores. Os autores sugerem que a alta CFE nos casos CRMHD pode ser um artefato do alto limiar de densidade e da eficiência utilizada na prescrição de formação estelar.
3. Agrupamento de Supernovas: Aproximadamente 50% das SNe ocorrem em gás difuso, pré-condicionado ($n_H < 10^{-2} \text{ cm}^{-3}$) em todos os três casos. Essa fração é amplamente independente do modelo de CR, indicando que o feedback de CR não altera significativamente o ambiente local onde as primeiras SNe explodem, apesar de alterar a fração de volume global do ISM difuso.
4. Propriedades dos Aglomerados e Parâmetros Viriais:
   • Raios e Dispersões de Velocidade: Os aglomerados no caso MHD exibem raios de massa média maiores e dispersões de velocidade 1D mais altas comparados aos casos CRMHD.
   • Ligação Gravitacional: Os aglomerados MHD são, em média, desvinculados ($\alpha_{vir} > 1$), enquanto os aglomerados CRMHD são, em média, vinculados ($\alpha_{vir} < 1$).
   • Mecanismo: Os autores atribuem essas diferenças ao orçamento de energia turbulenta do ISM. O caso MHD possui a maior SFR, levando a uma maior taxa de injeção de energia de SN e um ISM mais turbulento. Esta turbulência ambiente mais alta é herdada pelas nuvens moleculares em formação, resultando em aglomerados estelares com dispersões de velocidade mais altas que são menos vinculados gravitacionalmente. O feedback de CR reduz a SFR, diminuindo o orçamento de energia turbulenta e produzindo aglomerados mais fortemente vinculados.

Significância e Alegações
O artigo alega apresentar as primeiras simulações de formação de aglomerados estelares que incluem simultaneamente CRs acoplados dinamicamente, um ISM multifásico e partículas estelares individuais representativas de estrelas massivas. A principal significância reside em demonstrar que:

• O modelo de feedback Crisp pode reproduzir de forma autossuficiente a formação de aglomerados estelares com inclinações de função de massa e dependências ambientais observadas sem ajuste fino, mesmo na ausência de feedback radiativo pré-SN.
• O impacto primário do feedback de CR neste contexto é a redução da SFR global. Esta redução diminui o orçamento de energia turbulenta do ISM, o que, por sua vez, produz aglomerados estelares que são mais fortemente vinculados (menores dispersões de velocidade) comparados a simulações sem feedback de CR eficaz.
• Os CRs não parecem interromper dinamicamente o colapso de nuvens moleculares densas individuais ou alterar significativamente o pré-condicionamento local de SNe; seu efeito é mediado através da regulação global do estado do ISM.

Os autores concluem que, embora seu setup idealizado omita o feedback radiativo e ventos estelares, as diferenças relativas entre os casos quanto ao impacto de CR na ligação de aglomerados e na SFR são robustas. Os resultados fornecem novas restrições para modelos de transporte de CR e sugerem que o feedback de CR atua para mover os resultados das simulações ao longo das relações observadas, em vez de criar novos regimes não observados.