Substellar Initial Mass Function of Trumpler 14

Resumo Técnico: Função de Massa Inicial Subestelar de Trumpler 14

Problema e Motivação
Aglomerados estelares jovens e massivos fornecem ambientes críticos para investigar a formação de anãs marrons (BDs) sob condições de alta densidade estelar e intensa radiação ultravioleta distante (FUV). Embora a Função de Massa Inicial (IMF) seja frequentemente considerada universal acima de 1 $M_\odot$, os regimes de baixa massa e subestelar exibem variações que podem depender de fatores ambientais, como densidade estelar, densidade gasosa e a presença de estrelas OB massivas. Trumpler 14 (Tr 14), um aglomerado de $\sim$1 Myr de idade na Nebulosa de Carina, rico em estrelas do tipo O e caracterizado por alta densidade estelar ($\sim$1000 estrelas/pc$^2$) e fluxo FUV extremo, serve como um laboratório ideal para investigar esses efeitos. Estudos anteriores de Tr 14 (ex: Ascenso et al. 2007; Hur et al. 2012) foram limitados por profundidades observacionais rasas, tipicamente alcançando apenas $\sim$1–2 $M_\odot$, ou careceram de cobertura robusta da população subestelar. Este estudo visa apresentar o levantamento mais profundo até o momento da IMF em Tr 14, estendendo-se até o regime subestelar para caracterizar a IMF de baixa massa e avaliar a eficiência de formação de BDs neste ambiente de alta radiação.

Metodologia
Os autores utilizaram imagens profundas no infravermelho próximo nas bandas $J$, $H$ e $K_s$ obtidas com o Gemini-South Adaptive Optics Imager (GSAOI), auxiliado pelo Gemini Multi-Conjugate Adaptive Optics System (GeMS). As observações cobrem a região central de $\sim$1.4' $\times$ 1.4' do aglomerado.

• Redução de Dados e Fotometria: As imagens foram processadas usando o pipeline THELI. A extração de fontes e a fotometria foram realizadas usando SourceExtractor e PSFEx com uma estratégia de detecção "quente" (hot) e "fria" (cold) para resolver tanto fontes brilhantes quanto tênues no campo congestionado. A calibração fotométrica foi ancorada ao catálogo de Ascenso et al. (2007), que por sua vez foi vinculado ao sistema 2MASS.
• Completude: Testes de estrelas artificiais foram conduzidos para determinar a completude de detecção. O levantamento atinge 50% de completude em $\sim$0,02 $M_\odot$ e 20% de completude em $\sim$0,01 $M_\odot$.
• Determinação de Membros: Para isolar os membros do aglomerado dos contaminantes de campo, dois métodos independentes de campo de controle (CF) foram empregados:
  1. Modelo de Galáxia de Besançon: Uma população sintética escalonada para a área do aglomerado.
  2. CF Observacional VISTA: Um campo real observado com o telescópio VISTA, cruzado com o Gaia DR3 para remover estrelas de primeiro plano via cortes de movimento próprio e paralaxe.
     Um procedimento de descontaminação estatística (Muži´c et al. 2019) foi aplicado a ambos os CFs, gerando probabilidades de associação para cada fonte.
• Derivação de Massa e Extinção: Massas e extinções individuais foram derivadas comparando a fotometria das fontes com uma isócrona de 1 Myr combinando modelos evolutivos PARSEC e BT-Settl. Uma abordagem de Monte Carlo (1000 execuções por fonte) contabilizou as incertezas fotométricas.
• Construção da IMF: A IMF foi derivada agrupando as distribuições de massa dos prováveis membros, corrigida pela completude. A análise focou no expoente da lei de potência ($\alpha$), onde $dN/dM \propto M^{-\alpha}$, e na razão estrela-BD.

Resultos Principais

• Inclinações da IMF:
  • Regime de Alta Massa (0,2–4,5 $M_\odot$): Ambos os métodos de descontaminação produzem uma inclinação consistente de $\alpha = 1,72 \pm 0,04$. Esta é mais rasa que a inclinação padrão de Salpeter ($\alpha = 2,35$), mas consistente com outros aglomerados jovens.
  • Regime de Baixa Massa/Subestelar: Os resultados dependem do método de descontaminação. Usando o CF VISTA (preferido pelos autores por ser mais realista), a inclinação para 0,01–0,2 $M_\odot$ é $\alpha = 0,14 \pm 0,19$. Usando o modelo de Besançon, obtém-se $\alpha = 0,03 \pm 0,24$.
  • Impacto da Completude: Os autores observam que a inclinação extremamente rasa no regime de baixa massa é amplamente impulsionada pelo bin de menor massa (0,01–0,02 $M_\odot$), que é o mais afetado pela incompleteness. Excluindo este bin, as inclinações derivadas entram em acordo com valores reportados para outros aglomerados jovens (ex: $\alpha \approx 0,26–0,33$ para 0,03–0,2 $M_\odot$).
• Razão Estrela-BD: Para a faixa de massa de 0,03–1 $M_\odot$, a razão mediana estrela-BD é 4,0 (intervalo de confiança de 95%: 2,8–5,8) usando o CF VISTA. Isso se enquadra na faixa típica de 3–6 observada em outros ambientes.
• Contexto Ambiental: Tr 14 exibe uma densidade de superfície de $\sim$1000 estrelas/pc$^2$ e um fluxo FUV de $\sim$10$^{5,54}$ $G_0$, excedendo o de outros aglomerados estudados como RCW 38 e NGC 2244.

Significância e Alegações
O artigo afirma fornecer a função de massa subestelar mais completa para Tr 14 até o momento. A conclusão primária é que, embora o Tr 14 possua uma IMF subestelar relativamente plana, essa planaridade é fortemente influenciada pelo bin de menor massa, que sofre com a maior incompleteness. Quando este bin é excluído, a IMF torna-se consistente com as de outros aglomerados jovens, sugerindo que anãs marrons com massas acima de $\sim$0,03 $M_\odot$ se formam com eficiência semelhante em diferentes ambientes.

No entanto, os autores destacam um potencial efeito ambiental em relação aos objetos de menor massa. A relativa falta de objetos abaixo do limiar de 0,03 $M_\odot$, combinada com a densidade extrema e a radiação FUV do aglomerado, sugere que tais ambientes podem suprimir a formação de BDs com massas $< 0,03 M_\odot$. Isso contrasta com teorias que sugerem que estrelas OB podem aumentar a formação de BDs (como proposto para NGC 2244), apontando, em vez disso, para um mecanismo de supressão nos ambientes mais extremos. O estudo ressalta a necessidade de considerar tanto a densidade estelar quanto o feedback radiativo ao modelar a extremidade de baixa massa da IMF.