Diffuse gamma-ray emissions around the stellar cluster Berkeley 59

Sintesi Tecnica: Emissioni gamma diffuse attorno all'ammasso stellare Berkeley 59

Problema e Motivazione
L'origine dei raggi cosmici (CR) galattici rimane una questione aperta fondamentale in astrofisica. Sebbene i resti di supernova (SNR) siano stati a lungo considerati i principali siti di accelerazione, crescenti evidenze suggeriscono che gli ammassi stellari, in particolare quelli che ospitano stelle OB massive, fungano da fonti energetiche significative per i raggi cosmici ad altissima energia (VHE) e ultralalta energia (UHE). Questi ammassi generano potenti venti stellari ed esplosioni di supernova che creano superbolle, potenzialmente accelerando le particelle ad alte energie. Tuttavia, distinguere tra origini leptoniche e adroniche per le emissioni ad alta energia in questi ambienti è impegnativo, e la distribuzione spaziale dei CR all'interno degli ambienti degli ammassi è spesso semplificata eccessivamente nei modelli di diffusione. Questo studio investiga il giovane ammasso stellare Berkeley 59 (Ber 59) per determinare se agisca come sorgente di emissione gamma diffusa e per vincolare i meccanismi fisici di accelerazione e propagazione dei CR all'interno di tale ammasso.

Metodologia
Gli autori hanno condotto un'analisi completa utilizzando 16 anni di dati Fermi-LAT (Pass 8, raccolti da agosto 2008 a ottobre 2024) all'interno di una regione di interesse di $15^\circ \times 15^\circ$ centrata sulla sorgente 4FGL J0002.1+6721, che è spostata di soli 4,8 minuti d'arco da Ber 59.

1. Analisi dei raggi gamma: L'analisi ha utilizzato i pacchetti Fermipy e Fermi Science Tools. È stata eseguita un'analisi di verosimiglianza binnata nell'intervallo di energia da 300 MeV a 300 GeV. Il team ha testato vari modelli spaziali (sorgente puntuale, Gaussiana radiale e disco radiale) contro i modelli di fondo (emissione diffusa galattica gll_iem_v07 ed emissione extragalattica isotropa P8_R3_V3). Hanno impiegato la Statistica del Test (TS) e il Criterio di Informazione di Akaike (AIC) per valutare la significatività delle estensioni spaziali e la preferenza per specifici modelli spettrali (Legge di Potenza, LogParabola e PLSuperExpCutoff).
2. Caratterizzazione del contenuto gassoso: Per indagare l'origine adronica, gli autori hanno mappato la distribuzione dell'idrogeno molecolare ($H_2$), dell'idrogeno atomico neutro ($HI$) e dell'idrogeno ionizzato ($HII$) attorno a Ber 59.
   • La densità di colonna di $H_2$ è stata derivata dai dati CO (Dame et al. 2001) utilizzando un fattore di conversione $X_{CO} = 2.0 \times 10^{20} \text{ cm}^{-2} \text{ K}^{-1} \text{ km}^{-1} \text{ s}$.
   • I dati $HI$ sono stati ottenuti dal sondaggio all-sky HI4PI.
   • La densità di colonna di $HII$ è stata derivata dalla mappa free-free di Planck (353 GHz).
3. Modellazione dell'iniezione di CR: Gli autori hanno modellato la distribuzione spettrale di energia (SED) utilizzando il pacchetto NAIMA. Hanno testato sia scenari adronici (interazioni protone-protone che portano al decadimento di pioni neutri) sia scenari leptonici (scattering Compton inverso e bremsstrahlung). Inoltre, hanno derivato la relazione tra l'efficienza di accelerazione dei raggi cosmici ($\eta_{CR}$) e il coefficiente di diffusione ($D$) basandosi sul tasso di iniezione di energia cinetica dai venti stellari.

Risultati Chiave

• Estensione Spaziale: L'analisi ha rivelato un'emissione GeV estesa significativa attorno a Berkeley 59. L'emissione è meglio modellata da un disco radiale con un raggio di $1.02^\circ \pm 0.17^\circ$, producendo una significatività di estensione di $10.6\sigma$ (TS estensione = 97 per il modello a disco). Questa estensione è robusta anche tenendo conto di sorgenti vicine come 4FGL J0005.6+6746c.
• Proprietà Spettrali: Lo spettro dei raggi gamma è ben descritto da una semplice Legge di Potenza (PL) con un indice fotoni di $\Gamma = 2.88 \pm 0.06$ e una luminosità di $L_\gamma = 1.15 \times 10^{33} \text{ erg/s}$. La mancanza di una curvatura spettrale significativa ($TS_{LP} = 0$, $TS_{PLSC} = 1.6$) favorisce un modello PL rispetto al modello LogParabola utilizzato nel catalogo 4FGL per la sorgente vicina.
• Massa e Densità del Gas: La massa totale del gas all'interno della regione estesa (derivata da $H_2$ e $HII$) è di circa $2.89 \times 10^2 M_\odot$, con una densità numerica media di protoni di $99.41 \text{ cm}^{-3}$. La massa stellare totale dell'ammasso è stimata in $\sim 5.4 \times 10^4 M_\odot$.
• Efficienza di Accelerazione e Diffusione: Assumendo un'origine adronica, l'indice dei protoni migliore è $\alpha_p = 2.91 \pm 0.05$, con un'energia totale dei protoni superiore a 1 GeV di $W_p = 1.51 \times 10^{49} \text{ erg}$. Gli autori hanno derivato un coefficiente di diffusione $D$ nell'intervallo da $8.3 \times 10^{26}$ a $3.3 \times 10^{27} \text{ cm}^2 \text{ s}^{-1}$, corrispondente a un'efficienza di accelerazione $\eta_{CR}$ da 0,1 a 0,4. Questo coefficiente di diffusione è inferiore alla media galattica, suggerendo un confinamento potenziato dei CR all'interno dell'ambiente dell'ammasso.
• Contributo Leptonico: Sebbene uno scenario leptonico (scattering Compton inverso e bremsstrahlung) sia stato testato e abbia prodotto un indice elettronico migliore di $\alpha_e = 2.19 \pm 0.04$, gli autori notano che questa origine non può essere definitivamente esclusa, sebbene lo scenario adronico sia coerente con la densità del gas osservata e l'indice spettrale.

Significato e Affermazioni
Il lavoro afferma che Berkeley 59 è una sorgente di emissione gamma estesa originata da raggi cosmici accelerati dai venti dell'ammasso che interagiscono con il gas circostante. Lo studio fornisce un vincolo dettagliato sull'efficienza di accelerazione dei CR e sul coefficiente di diffusione in un ambiente di ammasso stellare giovane. Il coefficiente di diffusione derivato, essendo inferiore alla media galattica, supporta l'ipotesi che i campi magnetici turbolenti all'interno dell'ammasso (guidati da venti stellari e supernove) migliorino il confinamento dei CR, facilitando un'accelerazione efficiente delle particelle.

Gli autori concludono che l'energia cinetica totale iniettata dai venti stellari durante la vita di $\sim 2$ Myr dell'ammasso ($\sim 6.07 \times 10^{49} \text{ erg}$) è sufficiente a spiegare le emissioni gamma osservate. Questo lavoro contribuisce al crescente corpo di prove che gli ammassi stellari sono contributori significativi al bilancio dei raggi cosmici galattici, offrendo un laboratorio unico per studiare l'interazione tra formazione stellare, struttura del mezzo interstellare e fisica delle particelle ad alta energia. Gli autori riconoscono che le incertezze nella densità efficace del gas e nell'efficienza di accelerazione richiedono un ulteriore affinamento attraverso una mappatura del gas più dettagliata nelle osservazioni future.