The contribution to Galactic Centre {\gamma}-ray excess from cluster-born millisecond pulsars. Constraints from direct N-body simulations

Utilizzando simulazioni N-body dirette ad alta risoluzione, lo studio dimostra che un'associazione di pulsar millisecondo nate in ammassi globulari, sia sopravvissuti che distrutti, può spiegare l'eccesso di raggi gamma nel centro galattico, favorendo un'origine astrofisica rispetto a quella della materia oscura.

L'essenziale

Ecco una spiegazione semplice e creativa di questo studio scientifico, pensata per chiunque, anche senza un background in astronomia.

🌌 Il Mistero del "Bagliore Fantasma" al Centro della Galassia

Immaginate di guardare il cielo notturno e di vedere la nostra galassia, la Via Lattea, come una grande spirale di luci. Al centro di questa spirale c'è un "cuore" molto denso e attivo, dove si trova un buco nero gigante (Sagittarius A*).

Da anni, i telescopi spaziali come il Fermi-LAT vedono qualcosa di strano in quel cuore: c'è un eccesso di raggi gamma (una forma di luce molto energetica e invisibile all'occhio umano) che è più intenso di quanto ci si aspetterebbe. È come se qualcuno avesse acceso una lampadina potente al centro della galassia, ma non sappiamo chi l'abbia accesa.

Finora, c'erano due teorie principali su chi fosse il "colpevole":

1. La teoria "Fantasma": Forse è la materia oscura che si sta annichilando (sparando via energia) in quel punto.
2. La teoria "Folla": Forse non è un'unica fonte misteriosa, ma una folla enorme di piccole stelle morte e velocissime (dette pulsar a millisecondi) che brillano tutte insieme, creando l'illusione di un unico bagliore gigante.

🚚 Il Viaggio dei "Camionisti" Stellari

Questo nuovo studio, condotto da un team di ricercatori internazionali, decide di dare un'occhiata alla seconda teoria con un approccio molto pratico: la dinamica.

Immaginate che la nostra galassia sia una grande città e i Ammassi Globulari (grandi gruppi di stelle che girano intorno al centro) siano dei camionisti che viaggiano verso il centro città.

• Questi camionisti trasportano un carico speciale: stelle di neutroni (resti di stelle esplose, molto dense).
• Durante il viaggio, a causa delle strade sconnesse (le forze di marea della galassia) e delle frenate brusche, i camionisti perdono parte del carico.
• Il carico che cade a terra finisce proprio nel centro della città (il Centro Galattico).

Gli scienziati hanno usato dei supercomputer per fare una simulazione incredibilmente dettagliata di questo viaggio. Hanno creato un "mondo virtuale" dove hanno fatto viaggiare:

1. I camionisti attuali: Gli ammassi che vediamo ancora oggi e che stanno ancora viaggiando.
2. I camionisti fantasma: Ammassi che esistevano miliardi di anni fa, si sono distrutti completamente durante il viaggio, ma il loro carico (le stelle di neutroni) è arrivato a destinazione.

🔍 Cosa hanno scoperto?

Ecco i punti chiave, tradotti in immagini semplici:

• Il carico è arrivato: La simulazione ha mostrato che questi "camionisti" hanno depositato un numero enorme di stelle di neutroni nel centro della galassia.
• Da stelle morte a lampadine viventi: Non tutte le stelle di neutroni brillano. Alcune devono essere "riciclate" (come se qualcuno le avesse rimesso in moto) per diventare pulsar che emettono raggi gamma. Gli scienziati hanno usato un trucco statistico: hanno guardato quanti pulsar ci sono negli ammassi che vediamo oggi e hanno detto: "Ok, se ne vediamo X qui, probabilmente ce ne sono altrettanti nel carico perso nel viaggio".
• Il risultato è un "Sì" alla teoria della folla: Quando hanno sommato la luce di tutte queste pulsar immaginarie (sia quelle degli ammassi vivi che di quelli morti), il risultato è stato sorprendente. La luce totale generata corrisponde quasi perfettamente al "Bagliore Fantasma" (GCE) che vediamo nei dati reali.
• La forma è giusta: Non solo la quantità di luce è giusta, ma anche la forma. La simulazione mostra che queste stelle formano una nuvola luminosa che è più densa nel piano della galassia (come un disco) e ha anche una struttura verticale. Questo assomiglia molto a ciò che vediamo nei telescopi.

🧩 Perché è importante?

Questo studio è come un detective che risolve un caso.
Invece di dire "Forse è materia oscura", gli scienziati dicono: "Guardate, abbiamo ricostruito il viaggio di milioni di stelle morte. Se mettiamo insieme la loro luce, otteniamo esattamente il mistero che stiamo cercando di risolvere."

Per far combaciare perfettamente i numeri, hanno dovuto ipotizzare che ci fossero un po' più di "camionisti fantasma" (ammassi distrutti) di quanto pensassimo prima, o che le loro lampadine fossero un po' più luminose. Ma anche con queste piccole aggiustature, la teoria delle pulsar sembra molto più solida e naturale della teoria della materia oscura.

🚀 Cosa succede dopo?

Ora che sappiamo dove guardare, i telescopi radio futuri (come il potente SKA o FAST) potranno cercare direttamente queste stelle nel centro della galassia. Se troveranno la "folla" di pulsar che la simulazione prevede, il caso del "Bagliore Fantasma" sarà chiuso definitivamente: non è un fantasma, è solo una folla di stelle molto veloci e luminose che ci hanno nascosto la loro identità per troppo tempo!

In sintesi: La galassia non è illuminata da una magia misteriosa, ma dal lavoro di milioni di "camionisti" stellari che hanno consegnato il loro carico di stelle morte proprio nel cuore della Via Lattea.

Sommario tecnico

Ecco una sintesi tecnica dettagliata del documento scientifico in italiano, strutturata secondo le sezioni richieste.

Titolo: Il contributo dell'eccesso di raggi $\gamma$ del Centro Galattico da pulsar millisecondo nate negli ammassi: Vincoli da simulazioni N-body dirette

1. Il Problema: L'Eccesso di Raggi $\gamma$ del Centro Galattico (GCE)

Il Centro Galattico (GalC) è una regione di intensa attività astrofisica. Le osservazioni del telescopio spaziale Fermi-LAT hanno rivelato un eccesso anomalo di emissione di raggi $\gamma$ nell'intervallo di energia GeV, noto come Galactic Centre GeV Excess (GCE). Questo segnale, sfericamente simmetrico attorno a Sagittarius A* e con un picco spettrale tra 1 e 4 GeV, supera le previsioni dei modelli astrofisici standard basati sulle interazioni dei raggi cosmici con il mezzo interstellare.
Le interpretazioni principali del GCE sono due:

1. Fisica Esotica: Annichilazione di Materia Oscura (WIMP) concentrata nel centro galattico.
2. Origine Astrofisica: Una popolazione non risolta di Pulsar Millisecondo (MSP) nel bulge galattico.
   L'obiettivo di questo studio è riesaminare l'ipotesi delle MSP, seguendo un quadro puramente dinamico per determinare se le pulsar prodotte negli ammassi globulari (GC) e successivamente rilasciate nel centro galattico possano spiegare l'intero segnale osservato.

2. Metodologia

Gli autori hanno condotto simulazioni N-body dirette ad alta risoluzione utilizzando il codice parallelo $\phi$-GPU, che integra uno schema di integrazione di Hermite del quarto ordine e include l'evoluzione stellare aggiornata (perdita di massa, supernove, calci natali).

La metodologia si basa su due canali complementari:

• Ammassi Globulari Attuali (MW GCs): Simulazione di 12 ammassi globulari reali della Via Lattea (es. Terzan 5, NGC 6642) su orbite reali, tracciando la loro evoluzione dinamica in un potenziale galattico variabile nel tempo (derivato dalla simulazione IllustrisTNG-100).
• Ammassi Distrutti (Early GCs): Simulazione di un campione teorico di 54 ammassi (50 di massa $6 \times 10^4 M_\odot$e 4 di massa$1.8 \times 10^5 M_\odot$) che si sono completamente distrutti nei primi 5 Gyr, depositando i loro resti nel centro galattico.

Processo di modellazione:

1. Dinamica delle Stelle di Neutroni (NS): Le simulazioni tracciano la fase-spaziale completa (posizione e velocità) di tutte le stelle di neutroni (NS) formatesi negli ammassi. Si tiene conto dei "natal kicks" (velocità di espulsione dalle supernove), che causano la perdita della maggior parte delle NS dall'ammasso.
2. Conversione NS $\to$ MSP: Poiché le simulazioni non includono l'evoluzione binaria necessaria per "riciclare" le NS in MSP, gli autori hanno adottato un approccio statistico. Hanno calibrato il rapporto tra il numero di MSP osservate/predette e il numero di NS legate agli ammassi reali, ottenendo un coefficiente mediano empirico:
   $$k = \text{median}(N_{MSP}/N_{bound\_NS}) \approx 0.0114$$
   Questo coefficiente è stato applicato uniformemente alle NS depositate nel centro galattico (sia da ammassi vivi che distrutti) per stimare il numero totale di MSP.
3. Calcolo del Flusso: Sono state generate mappe del cielo e profili di flusso cumulativo assumendo diverse luminosità medie per MSP ($\langle L \rangle \sim 1, 4, 8 \times 10^{33}$ erg s$^{-1}$) e confrontandole con i dati osservativi di Fermi-LAT.

3. Contributi Chiave

• Approccio Dinamico Completo: A differenza di studi precedenti che assumevano frazioni fisse di NS trattenute, questo lavoro calcola dinamicamente la frazione di NS rimaste nel raggio di 1 kpc per ogni singolo ammasso, tenendo conto delle orbite specifiche e della storia di stripping mareale.
• Inclusione degli Ammassi Distrutti: Il modello quantifica esplicitamente il contributo degli ammassi globulari che si sono già dissolti, i quali costituiscono una fonte significativa di NS nel bulge.
• Calibrazione Empirica: Uso di dati osservativi recenti (inclusi nuovi MSP scoperti da FAST) per calibrare il rapporto MSP/NS, riducendo le incertezze teoriche sulla popolazione di pulsar.

4. Risultati

• Distribuzione delle NS: Le NS provenienti dagli ammassi distrutti mostrano una morfologia assimmetrica con una sovraddensità pronunciata nel piano galattico e una struttura meno evidente perpendicolarmente ad esso. Questo è coerente con recenti report sulla morfologia del GCE.
• Contributo al Flusso $\gamma$:
  • Gli ammassi attuali (vivi) forniscono già un segnale $\gamma$ sostanziale.
  • Aggiungendo il contributo degli ammassi distrutti, l'ampiezza e la concentrazione centrale del segnale aumentano significativamente.
  • Utilizzando una luminosità media di $\langle L \rangle \approx 8 \times 10^{33}$ erg s$^{-1}$ e il coefficiente mediano $k$, il modello riproduce bene l'ampiezza e la concentrazione del GCE osservato.
  • Per ottenere un accordo quasi perfetto con i dati osservativi (dopo aver sottratto il contributo del buco nero supermassiccio), è necessario aumentare il numero di NS provenienti dagli ammassi distrutti di un fattore circa 2 rispetto al campione simulato, oppure adottare i valori massimi del rapporto MSP/NS ($k_{theor2}$).
• Confronto con la Materia Oscura: Il modello MSP riesce a riprodurre naturalmente sia l'ampiezza che la concentrazione spaziale del segnale senza richiedere parametri fisici esotici.

5. Significato e Conclusioni

Lo studio conclude che l'origine astrofisica (MSP) dell'eccesso di raggi $\gamma$ del Centro Galattico è favorita rispetto all'ipotesi di annichilazione della Materia Oscura.

• Robustezza: La combinazione di MSP consegnate sia da ammassi sopravvissuti che da ammassi distrutti spiega il segnale osservato sotto ipotesi ragionevoli.
• Implicazioni Future: I risultati suggeriscono che future osservazioni radio profonde (con telescopi come FAST, MeerKAT, SKA) dovrebbero cercare pulsar nelle regioni di sovraddensità previste nel piano galattico. Inoltre, la necessità di una luminosità media relativamente alta ($\sim 8 \times 10^{33}$ erg s$^{-1}$) o di un numero maggiore di ammassi distrutti indica che potrebbero esserci meccanismi di formazione in situ nel bulge o una storia di distruzione degli ammassi più intensa di quanto simulato.
• Limiti e Sviluppi: Il modello attuale si basa su un coefficiente $k$ globale e non include l'evoluzione binaria dettagliata. Futuri lavori dovrebbero integrare l'evoluzione binaria diretta nelle simulazioni N-body e considerare la geometria del fascio delle pulsar per ridurre le sistematiche.

In sintesi, il lavoro fornisce una prova dinamica solida a sostegno dell'ipotesi che le pulsar millisecondo, rilasciate dalla distruzione e dall'evoluzione degli ammassi globulari, siano la fonte primaria del misterioso eccesso di raggi $\gamma$ al centro della nostra galassia.