Substellar Initial Mass Function of Trumpler 14

Sintesi Tecnica: Funzione di Massa Iniziale Sottostellare di Trumpler 14

Problema e Motivazione
I giovani ammassi stellari massicci forniscono ambienti critici per investigare la formazione di nane brune (BD) in condizioni di alta densità stellare e intensa radiazione ultravioletta lontano (FUV). Sebbene la Funzione di Massa Iniziale (IMF) sia spesso considerata universale sopra 1 $M_\odot$, i regimi di bassa massa e sottostellare esibiscono variazioni che possono dipendere da fattori ambientali come la densità stellare, la densità del gas e la presenza di stelle massicce di tipo OB. Trumpler 14 (Tr 14), un ammasso di circa 1 Myr nella Nebulosa della Carena, ricco di stelle di tipo O e caratterizzato da un'alta densità stellare ($\sim$1000 stelle/pc$^2$) ed estremo flusso FUV, funge da laboratorio ideale per indagare questi effetti. Studi precedenti su Tr 14 (ad es. Ascenso et al. 2007; Hur et al. 2012) erano limitati da profondità osservative superficiali, raggiungendo tipicamente solo $\sim$1–2 $M_\odot$, o mancavano di una copertura robusta della popolazione sottostellare. Questo studio mira a presentare il sondaggio più profondo ad oggi dell'IMF in Tr 14, estendendosi fino al regime sottostellare per caratterizzare l'IMF a bassa massa e valutare l'efficienza di formazione delle BD in questo ambiente ad alta radiazione.

Metodologia
Gli autori hanno utilizzato immagini nel vicino infrarosso profondo nelle bande $J$, $H$ e $K_s$ ottenute con il Gemini-South Adaptive Optics Imager (GSAOI$, assistito dal sistema Gemini Multi-Conjugate Adaptive Optics (GeMS). Le osservazioni coprono la regione centrale di $\sim$1.4' $\times$ 1.4' dell'ammasso.

• Riduzione dei Dati e Fotometria: Le immagini sono state elaborate utilizzando la pipeline THELI. L'estrazione delle sorgenti e la fotometria sono state eseguite utilizzando SourceExtractor e PSFEx con una strategia di rilevamento "hot" e "cold" per risolvere sia le sorgenti brillanti che quelle deboli in un campo affollato. La calibrazione fotometrica è stata ancorata al catalogo di Ascenso et al. (2007), il quale era a sua volta legato al sistema 2MASS.
• Completezza: Test di stelle artificiali sono stati condotti per determinare la completezza del rilevamento. Il sondaggio raggiunge una completezza del 50% a $\sim$0.02 $M_\odot$ e del 20% a $\sim$0.01 $M_\odot$.
• Determinazione dell'Appartenenza: Per isolare i membri dell'ammasso dai contaminanti del campo, sono stati impiegati due approcci indipendenti di campo di controllo (CF):
  1. Modello di Galassia di Besançon: Una popolazione sintetica scalata sull'area dell'ammasso.
  2. CF Osservativo VISTA: Un campo reale osservato con il telescopio VISTA, incrociato con Gaia DR3 per rimuovere le stelle in primo piano tramite tagli di moto proprio e parallasse.
     Una procedura di decontaminazione statistica (Mužić et al. 2019) è stata applicata a entrambi i CF, generando probabilità di appartenenza per ogni sorgente.
• Derivazione di Massa ed Estinzione: Masse ed estinzioni individuali sono state derivate confrontando la fotometria delle sorgenti con un isocrona di 1 Myr che combina i modelli evolutivi PARSEC e BT-Settl. Un approccio Monte Carlo (1000 iterazioni per sorgente) ha tenuto conto delle incertezze fotometriche.
• Costruzione dell'IMF: L'IMF è stata derivata raggruppando le distribuzioni di massa dei probabili membri, corretta per la completezza. L'analisi si è concentrata sulla pendenza della legge di potenza ($\alpha$) dove $dN/dM \propto M^{-\alpha}$, e sul rapporto stella-BD.

Risultati Chiave

• Pendenze dell'IMF:
  • Regime di Alta Massa (0.2–4.5 $M_\odot$): Entrambi i metodi di decontaminazione forniscono una pendenza coerente di $\alpha = 1.72 \pm 0.04$. Questa è più piatta rispetto alla standard pendenza di Salpeter ($\alpha = 2.35$) ma coerente con altri ammassi giovani.
  • Regime di Bassa Massa/Sottostellare: I risultati dipendono dal metodo di decontaminazione. Utilizzando il CF VISTA (preferito dagli autori in quanto più realistico), la pendenza per 0.01–0.2 $M_\odot$ è $\alpha = 0.14 \pm 0.19$. Utilizzando il modello di Besançon, si ottiene $\alpha = 0.03 \pm 0.24$.
  • Impatto della Completezza: Gli autori notano che la pendenza estremamente piatta nel regime di bassa massa è largamente guidata dal bin di massa più basso (0.01–0.02 $M_\odot$), che è il più influenzato dall'incompletezza. Escludendo questo bin, le pendenze derivate entrano in accordo con i valori riportati per altri ammassi giovani (ad es. $\alpha \approx 0.26–0.33$ per 0.03–0.2 $M_\odot$).
• Rapporto Stella-BD: Per l'intervallo di massa 0.03–1 $M_\odot$, il rapporto mediano stella-BD è 4.0 (intervallo di confidenza al 95%: 2.8–5.8) usando il CF VISTA. Questo rientra nel tipico intervallo 3–6 osservato in altri ambienti.
• Contesto Ambientale: Tr 14 presenta una densità superficiale di $\sim$1000 stelle/pc$^2$ e un flusso FUV di $\sim$10$^{5.54}$ $G_0$, superando quello di altri ammassi studiati come RCW 38 e NGC 2244.

Significatività e Rivendicazioni
L'articolo rivendica di fornire la funzione di massa sottostellare più completa per Tr 14 ad oggi. La conclusione primaria è che, sebbene Tr 14 ospiti un'IMF sottostellare relativamente piatta, questa piattezza è fortemente influenzata dal bin di massa più bassa, che soffre della maggiore incompletezza. Quando questo bin viene escluso, l'IMF diventa coerente con quelle di altri ammassi giovani, suggerendo che le nane brune con masse sopra $\sim$0.03 $M_\odot$ si formino con un'efficienza simile in diversi ambienti.

Tuttavia, gli autori evidenziano un potenziale effetto ambientale riguardante gli oggetti di massa più bassa. La relativa mancanza di oggetti sotto la soglia di 0.03 $M_\odot$, combinata con l'estrema densità e la radiazione FUV dell'ammasso, suggerisce che tali ambienti possano sopprimere la formazione di BD a masse $< 0.03 M_\odot$. Ciò contrasta con le teorie che suggeriscono che le stelle OB possano potenziare la formazione di BD (come proposto per NGC 2244), puntando invece a un potenziale meccanismo di soppressione negli ambienti più estremi. Lo studio sottolinea la necessità di considerare sia la densità stellare che il feedback radiativo quando si modellano l'estremità di bassa massa dell'IMF.