Perché alcune galassie non hanno la barra a spirale?

Le galassie a spirale spesso mostrano una “barra” centrale, ma non tutte la hanno. Gli astronomi hanno scoperto che la presenza di questa barra dipende da due fattori chiave nei loro modelli al computer: il parametro di ammorbidimento (ε) e la frazione di massa del disco (md). Il parametro ε controlla quanto le forze gravitazionali tra le particelle simulano la realtà; valori troppo piccoli provocano un “riscaldamento numerico”, cioè particelle che accelerano troppo al centro del disco. Questo crea orbite caotiche e impedisce la formazione della barra. D’altro canto, aumentare md di solito ritarda la comparsa della barra, ma solo se ε è moderato (≥10 pc). Se ε è molto piccolo, anche un disco più massiccio non riesce a far nascere la barra, perché il riscaldamento numerico domina. Gli esperimenti con milioni di particelle (N≈10⁷) mostrano che, con ε ≤5 pc, la barra è quasi sempre soppressa, indipendentemente da md. Un altro indizio utile è il profilo di accelerazione verticale del disco: quando questo profilo si “infiamma”, è segno che il riscaldamento numerico sta interferendo. In sintesi, la barra galattica è il risultato di un delicato equilibrio tra la massa del disco e la scelta del parametro ε nei simulazioni, e capire questo aiuta a interpretare le diverse forme osservate nelle galassie reali.

Le simulazioni N‑body sono il principale strumento per studiare la dinamica galattica. In queste simulazioni, ogni particella rappresenta una porzione di materia (stelle, gas o materia oscura) e interagisce gravitazionalmente con le altre. Il parametro di ammorbidimento ε è introdotto per evitare divergenze numeriche quando due particelle si avvicinano troppo: la forza gravitazionale viene “smussata” su una scala ε, rendendo il potenziale più regolare. Tuttavia, se ε è troppo piccolo rispetto alla risoluzione spaziale della simulazione, le interazioni a corto raggio diventano troppo intense, generando accelerazioni artificiali (numerical heating). Questo fenomeno è particolarmente evidente nei modelli a risoluzione elevata (N≈10⁷), dove le particelle hanno masse molto piccole e quindi sono più sensibili a variazioni di ε. Il risultato è una popolazione di orbite quasi‑circolari al centro del disco, che non supportano la crescita di una struttura bar‑like.

La frazione di massa del disco md influisce sulla stabilità dinamica: un disco più massiccio aumenta la gravità auto‑coerente, favorendo l’instabilità bar. La dipendenza teorica è spesso descritta da una relazione inversa, bar strength ∝ md⁻¹, ma le simulazioni mostrano che questa dipendenza è modulata da ε. Quando ε è moderato (≥10 pc), l’aumento di md riduce il drift dell’accelerazione verticale senza alterare significativamente la crescita della barra. Al contrario, per ε ≤5 pc, l’effetto di riscaldamento numerico sovrasta l’influenza di md, sopprimendo la barra anche in dischi molto massicci.

Un indicatore diagnostico efficace è il profilo di accelerazione verticale a_z(z) del disco. In assenza di riscaldamento numerico, a_z(z) rimane stabile e segue la distribuzione di massa prevista dal modello di equilibrio. Un picco improvviso o una crescita esponenziale di a_z indica che le particelle stanno ricevendo impulsi energetici non fisici, tipici del riscaldamento numerico. Monitorare questo profilo permette di identificare in tempo reale quando ε è troppo piccolo e di intervenire, ad esempio aumentando ε o riducendo la risoluzione temporale, per mantenere la simulazione fisicamente affidabile.

Infine, la scelta di ε deve bilanciare due esigenze: minimizzare il riscaldamento numerico e preservare la capacità della simulazione di risolvere le scale di interesse, come la larghezza della barra (spesso dell’ordine di 1‑2 kpc). Tecniche avanzate, come l’uso di softening adattivo (ε variabile in funzione della densità locale), stanno emergendo per ottimizzare questo compromesso, consentendo di studiare con maggiore precisione la formazione e l’evoluzione delle barre galattiche in un contesto cosmologico più ampio.