Controlling artificial surface heating in neutron star simulations: Application to hybrid equations of state

Le stelle di neutroni sono gli oggetti più densi dell'universo, i nuclei collassati di stelle massicce che hanno esaurito il proprio combustibile. Immaginate una sfera grande quanto una città contenente più massa del nostro intero Sole, compattata così strettamente che un singolo cucchiaino del suo materiale peserebbe miliardi di tonnellate. Questi laboratori cosmici permettono agli scienziati di studiare la materia in condizioni di pressione e densità che non possono essere replicate sulla Terra. Per comprendere come queste stelle si comportino, come si muovano e come possano scontrarsi tra loro, i ricercatori si affidano a potenti simulazioni al computer. Questi modelli digitali agiscono come laboratori virtuali, risolvendo equazioni complesse per prevedere l'evoluzione delle stelle e le onde gravitazionali che emettono quando collidono. Tuttavia, affinché queste previsioni siano utili, i modelli informatici devono essere incredibilmente precisi, specialmente al limite estremo della stella, dove la materia densa incontra il vuoto dello spazio circostante.

Il confine tra una stella di neutroni e il vuoto circostante è un luogo di estrema difficoltà per le simulazioni al computer. Nel mondo reale, questa transizione è netta, ma in un computer la stella è rappresentata da una griglia di punti, come una maglia digitale. Mentre la stella si muove o cambia forma all'interno di questa griglia, il computer fatica a gestire il salto improvviso dalla materia densa al nulla. Questa lotta crea un errore sottile ma persistente: il computer aggiunge accidentalmente calore alla superficie della stella dove non dovrebbe esistere alcun calore fisico. Questo riscaldamento artificiale non è causato dall'attrito o da reazioni nucleari; è un errore matematico. Con il tempo, questo calore falso si diffonde verso l'interno, facendo espandere leggermente la stella e alterando la sua energia interna. Se lasciato senza controllo, questo errore numerico può distorcere la simulazione, portando gli scienziati a trarre conclusioni errate sulla struttura della stella o sui segnali che essa invia all'universo.

In uno studio recente, un ricercatore della Goethe University Frankfurt ha indagato su come fermare questo riscaldamento artificiale. Il lavoro si è concentrato su un tipo specifico di codice informatico utilizzato per simulare le stelle di neutroni, noto come codice BAM. Questo codice è ampiamente utilizzato per modellare le violente collisioni di stelle di neutroni binarie, dove due tali stelle spiraleggiano l'una verso l'altra e infine si fondono. Il ricercatore ha testato una nuova strategia matematica chiamata limitazione del flusso basata sull'entropia (entropy-based flux-limiting). In termini semplici, questo metodo agisce come un filtro intelligente per i calcoli del computer. Monitora il flusso di energia e materia alla superficie della stella e regola automaticamente il metodo di calcolo quando rileva il tipo di cambiamenti bruschi che di solito causano errori. In questo modo, impedisce al computer di generare quel calore spurio fin dall'inizio.

Per vedere se questo metodo funzionasse, il ricercatore ha eseguito due serie di simulazioni. La prima serie coinvolgeva stelle di neutroni isolate, che si trovano sole nello spazio, mentre la seconda simulava coppie di stelle che spiraleggiano l'una verso l'altra. In entrambi gli scenari, al computer è stata data una scelta: utilizzare il metodo di calcolo standard o utilizzare il nuovo filtro basato sull'entropia. I risultati sono stati chiari e coerenti. Quando veniva utilizzato il metodo standard, l'energia interna delle stelle aumentava costantemente nel tempo e le stelle stesse diventavano più grandi, come se si stessero gonfiando a causa del calore falso. Questa espansione non era una realtà fisica, ma un sintomo dell'errore numerico. Al contrario, quando veniva applicato il nuovo filtro, le stelle rimanevano stabili. La loro energia interna rimaneva vicina al valore iniziale e le loro dimensioni non cambiavano artificialmente. Il nuovo metodo ha mantenuto con successo le stelle dal riscaldarsi ed espandersi, rappresentando molto più fedelmente il loro stato fisico.

Lo studio ha esaminato una vasta gamma di modelli di stelle di neutroni, ciascuno basato su diverse teorie riguardanti il comportamento della materia a densità così estreme. Alcuni modelli assumevano che le stelle fossero composte da un tipo di particella, mentre altri includevano miscele più complesse. Indipendentemente dal modello utilizzato, il nuovo metodo di filtraggio ha performato ugualmente bene. Ha ridotto il riscaldamento artificiale di un margine significativo, mantenendo la crescita dell'energia interna a una frazione minuscola rispetto a quella osservata con il vecchio metodo. Nelle simulazioni di stelle binarie, dove la forza di gravità e le forze mareali sono intense, il nuovo metodo ha continuato a sopprimere il riscaldamento durante l'intera danza dell'inseguimento a spirale, fino al momento della collisione. Le stelle simulate con il nuovo metodo rimanevano più piccole e più fredde rispetto a quelle simulate con il vecchio metodo, suggerendo che gli artefatti digitali venivano eliminati efficacemente.

Questo lavoro non pretende di aver risolto ogni problema nelle simulazioni delle stelle di neutroni, né suggerisce che il nuovo metodo sia una soluzione magica per tutte le sfide numeriche. Al contrario, fornisce una prova forte che una specifica fonte di errore — il riscaldamento artificiale alla superficie — può essere controllata. Il ricercatore osserva che, sebbene questi test siano stati esegtti utilizzando un modo semplificato di modellare il calore della stella, il successo del metodo suggerisce che funzionerà altrettanto bene quando applicato a modelli più complessi e realistici che tracciano la temperatura con maggiore dettaglio. Le scoperte offrono una strada da seguire per creare simulazioni più accurate delle collisioni tra stelle di neutroni. Rimuovendo queste distorsioni digitali, gli scienziati possono fidarsi più profondamente dei loro modelli, assicurando che i segnali che prevedono per i rilevatori di onde gravitazionali siano basati sulla vera fisica dell'universo piuttosto che sui limiti del codice informatico.