Numerical calculations of neutron star mountains supported by crustal lattice pressure

Questo studio calcola in modo coerente le asimmetrie di temperatura e di pressione reticolare in stelle di neutroni accrescenti, utilizzando equazioni di stato realistiche, e conclude che le "montagne" risultanti sono troppo piccole per determinare l'equilibrio di rotazione delle binarie X a bassa massa o per generare un'emissione significativa di onde gravitazionali.

L'essenziale

Immagina una stella di neutroni come una palla da biliardo cosmica, ma fatta di materia così densa che un cucchiaino peserebbe quanto una montagna. Queste stelle ruotano su se stesse a velocità folli, come trottole che non vogliono fermarsi.

Gli scienziati sperano che alcune di queste stelle emettano un "canto" continuo sotto forma di onde gravitazionali (increspature nello spazio-tempo), che potremmo un giorno ascoltare con i nostri rivelatori come LIGO. Per farlo, però, la stella non deve essere una sfera perfetta: deve avere una "gobba", una piccola irregolarità chiamata "montagna".

Il problema è: come si forma questa montagna su una stella così dura e calda?

La vecchia teoria: Il "terremoto" nucleare

Fino a poco tempo fa, la teoria più accreditata (il modello UCB) diceva che la montagna si formava perché la stella, mangiando materia dalla compagna, subiva dei "cambi di stato" nucleari. Immagina di avere una torta a strati: se alcune parti della torta si cuociono prima di altre, gli strati si spostano e la torta diventa storta.
In questa vecchia teoria, si pensava che gli elettroni venissero "catturati" dai nuclei atomici in punti specifici della crosta della stella, creando degli strati ondulati. Tuttavia, calcoli più recenti hanno mostrato che questi "strati nucleari" sono molto meno frequenti e meno potenti di quanto si pensasse. È come scoprire che la torta non ha abbastanza strati per diventare storta.

La nuova teoria: La "crosta termica"

In questo nuovo studio, gli autori (Hutchins e Jones) hanno provato un approccio diverso. Invece di guardare i cambiamenti nucleari, hanno guardato la pressione del reticolo cristallino.

Immagina la crosta della stella di neutroni non come una roccia, ma come un gigantesco cristallo di ghiaccio o una griglia metallica perfetta.

1. Il calore non è uniforme: A causa dei forti campi magnetici interni della stella, il calore non si diffonde in modo uniforme. È come se avessi un radiatore che scalda molto di più un lato della stanza rispetto all'altro.
2. La pressione cambia: Quando una parte di questo "cristallo" è più calda, si espande leggermente e cambia la pressione interna, proprio come l'aria calda in un palloncino che spinge contro le pareti.
3. La deformazione: Questa differenza di pressione spinge la crosta solida a deformarsi, creando quella piccola "gobba" o montagna.

Gli autori hanno usato modelli matematici molto precisi (le equazioni di stato BSk) per calcolare esattamente quanto questa "gobba" sarebbe grande.

Il risultato: Una montagna... minuscola

La sorpresa? Le montagne create da questo meccanismo sono piccolissime.
Pensaci così: se la superficie della stella fosse liscia come un tavolo da biliardo, questa "montagna" sarebbe alta meno di un atomo. È così piccola che, anche se la stella ruota velocissima, l'onda gravitazionale che produce è debole come un sussurro in mezzo a un uragano.

Cosa significa per noi?

1. Non è la causa principale: Questi calcoli suggeriscono che le "gobbe" termiche non sono abbastanza grandi da spiegare perché le stelle di neutroni smettono di accelerare (il cosiddetto equilibrio di coppia). C'è probabilmente qualcos'altro che le frena.
2. Servono orecchie più sensibili: Se vogliamo "sentire" queste montagne, i nostri attuali rivelatori di onde gravitazionali non sono abbastanza sensibili. Servirà la prossima generazione di strumenti (come l'Einstein Telescope) che saranno molto più precisi.
3. Un nuovo campo di caccia: Anche se le stelle normali (in sistemi binari) hanno montagne troppo piccole, gli autori suggeriscono di guardare altrove. Forse le stelle di neutroni in sistemi ultra-luminosi (ULX), che hanno campi magnetici enormi e mangiano materia a velocità incredibili, potrebbero avere montagne più grandi e quindi essere più facili da ascoltare.

In sintesi: Gli scienziati hanno scoperto che le "gobbe" sulle stelle di neutroni causate dal calore e dal magnetismo sono molto più piccole di quanto speravamo. È come cercare di sentire il battito di un'ala di mosca da un chilometro di distanza: è possibile in teoria, ma ci serviranno orecchie (e strumenti) molto, molto migliori per riuscirci.

Sommario tecnico

Titolo: Calcoli numerici di montagne di stelle di neutroni sostenute dalla pressione del reticolo cristallino della crosta

1. Il Problema

Le onde gravitazionali (GW) continue potrebbero essere la chiave per spiegare le frequenze di rotazione delle stelle di neutroni (NS) in sistemi binari a raggi X di bassa massa (LMXB). Un meccanismo proposto per generare tali onde è la formazione di una "montagna" di massa (un'asimmetria quadrupolare) sostenuta da deformazioni elastiche nella crosta solida della stella.
Il problema centrale affrontato da questo studio è la mancanza di modelli auto-consistenti che colleghino l'origine delle asimmetrie termiche alla generazione effettiva della deformazione elastica.

• Limiti dei modelli precedenti: La maggior parte degli studi si è concentrata separatamente sull'origine delle tensioni elastiche o sull'asimmetria di temperatura. In particolare, il modello classico (UCB, Ushomirsky et al. 2000) si basa su uno spostamento degli strati di cattura elettronica dovuto a gradienti termici. Tuttavia, recenti calcoli dell'Equazione di Stato (EOS) della materia nucleare suggeriscono che tali strati di cattura siano molto meno numerosi e abbiano energie di soglia più basse rispetto a quanto ipotizzato, rendendo il modello UCB potenzialmente sovrastimato.
• L'ipotesi alternativa: In un lavoro precedente, gli autori avevano proposto che la pressione termica del reticolo ionico della crosta (anziché gli strati di cattura) potesse essere una fonte competitiva di asimmetria. Tuttavia, mancava un calcolo dettagliato e auto-consistente di questo scenario.

2. Metodologia

Gli autori hanno sviluppato un quadro di calcolo auto-consistente che integra tre fasi principali:

1. Struttura Idrostatica e EOS:

   • Utilizzo di equazioni di stato realistiche e unificate (BSk19, BSk20, BSk21) basate sui funzionali di densità nucleare di Bruxelles-Montreal. Queste descrivono coerentemente tutte le regioni della stella (crosta esterna, crosta interna e nucleo) in modo termodinamicamente consistente.
   • Per evitare problemi di convergenza numerica dovuti ai salti di densità nelle transizioni nucleari, l'EOS tabulata è stata approssimata mediante una funzione analitica complessa (fit a 23 parametri).
   • La struttura idrostatica è risolta con l'equazione di Tolman-Oppenheimer-Volkoff (TOV), ma le fasi termiche ed elastiche sono trattate in un contesto newtoniano (con una mappatura delle coordinate per preservare la struttura relativistica di base).

2. Struttura Termica e Asimmetrie:

   • Calcolo del profilo termico di fondo per stelle di neutroni in accrescimento stazionario, includendo il riscaldamento profondo (reazioni nucleari non in equilibrio) e il riscaldamento superficiale (shallow heating).
   • Introduzione di asimmetrie termiche su larga scala causate da conduttività termica anisotropa dovuta a campi magnetici interni (poloidali o toroidali). I campi magnetici deflettono il flusso di calore degli elettroni relativistici, creando gradienti di temperatura latitudinali.

3. Deformazione Elastica (Montagne):

   • Calcolo della risposta elastica della crosta alle perturbazioni di temperatura.
   • Meccanismo chiave: Invece di usare lo spostamento degli strati di cattura, gli autori calcolano le perturbazioni di pressione generate direttamente dalla componente termica della pressione del reticolo ionico ($P_{th}$). Questa pressione dipende dalla temperatura e dalla struttura cristallina della crosta.
   • Risoluzione delle equazioni di perturbazione elastica per ottenere il campo di spostamento lagrangiano e, di conseguenza, il momento di quadrupolo di massa ($Q_{22}$) e l'ellitticità ($\varepsilon$).

3. Contributi Chiave

• Primo calcolo auto-consistente completo: Per la prima volta, asimmetrie termiche, perturbazioni di pressione del reticolo e masse asimmetriche sono calcolate all'interno di un unico modello, utilizzando EOS unificate realistiche (BSk).
• Superamento del modello UCB: Il lavoro valida e quantifica l'ipotesi che la pressione del reticolo termico possa generare montagne, offrendo un'alternativa al modello degli strati di cattura che è stato messo in discussione dalle nuove EOS.
• Analisi comparativa: Confronto diretto tra le montagne generate dalla pressione del reticolo e quelle stimate dal vecchio modello degli strati di cattura (aggiornato con le nuove EOS), dimostrando come le stime tradizionali siano state sovrastimate a causa della mancanza di strati di cattura ad alta energia nelle regioni interne.
• Parametrizzazione EOS: Sviluppo di un fit analitico preciso per le EOS accresciute (BSk19-21) per permettere integrazioni numeriche stabili attraverso le discontinuità di densità.

4. Risultati

• Dimensioni delle montagne: Le montagne generate dal meccanismo di pressione del reticolo termico sono estremamente piccole. L'ellitticità massima calcolata è dell'ordine di $\varepsilon \lesssim 10^{-11}$.
• Confronto con il limite di equilibrio: Questo valore è di diversi ordini di grandezza inferiore al limite di equilibrio del torque gravitazionale ($\varepsilon_{TB} \sim 10^{-9} - 10^{-7}$) necessario per spiegare la frenata della rotazione delle LMXB tramite onde gravitazionali.
• Dipendenza dai parametri: L'ellitticità aumenta con il tasso di accrescimento e l'entità del riscaldamento superficiale, ma raggiunge un plateau quando il tasso di accrescimento è così alto da fondere la crosta esterna ($\dot{M} \gtrsim 5 \times 10^{-9} M_\odot/\text{anno}$), limitando la regione solida disponibile per la deformazione.
• Confronto con il modello UCB: Anche aggiornando il modello UCB con le nuove EOS (che riducono drasticamente il numero di strati di cattura), le montagne risultanti rimangono molto più grandi di quelle generate dalla pressione del reticolo termico (fino a 3 ordini di grandezza in più per il modello UCB aggiornato, ma comunque spesso insufficienti per il torque balance). Tuttavia, il modello UCB aggiornato è comunque inferiore ai limiti osservativi attuali.
• Implicazioni per i campi magnetici: Le asimmetrie termiche sono più elevate per campi magnetici nella crosta ($B \sim 10^{12}$ G) rispetto a quelli nel nucleo ($B \sim 10^8$ G), ma in entrambi i casi le montagne risultano troppo piccole per essere rilevate con la sensibilità attuale.

5. Significato e Conclusioni

• Rilevanza per le Onde Gravitazionali: Le montagne termiche sostenute dalla pressione del reticolo, sebbene fisicamente reali e calcolate in modo rigoroso, sono troppo piccole per essere la causa principale dell'equilibrio di rotazione delle LMXB o per essere rilevate dai rivelatori attuali (LIGO/Virgo/KAGRA). La loro emissione di onde gravitazionali è trascurabile rispetto ai limiti osservativi attuali ($\varepsilon \sim 10^{-7}$).
• Impatto sui modelli teorici: Lo studio conferma che i modelli basati sugli strati di cattura elettronica (UCB) sono probabilmente sovrastimati se si utilizzano EOS moderne, e che anche il meccanismo alternativo della pressione del reticolo non è sufficiente a spiegare le osservazioni di spin-down.
• Prospettive future:
  • La ricerca di montagne termiche potrebbe richiedere rivelatori di prossima generazione (come Einstein Telescope o Cosmic Explorer) con sensibilità molto superiori.
  • Il meccanismo potrebbe essere più rilevante in sistemi con campi magnetici estremamente forti e tassi di accrescimento super-Eddington, come le sorgenti X ultraluminose (ULX), sebbene la loro bassa frequenza di rotazione possa compensare il guadagno in massa quadrupolare.
  • Sono necessari miglioramenti nel modello, inclusa l'estensione alla Relatività Generale completa e il rilassamento dell'approssimazione di Cowling, per affinare ulteriormente le stime.

In sintesi, il paper chiude un capitolo importante nella modellizzazione delle montagne di stelle di neutroni, dimostrando che i meccanismi termici "realistici" basati sulle attuali conoscenze della fisica nucleare producono deformazioni troppo deboli per essere la fonte dominante delle onde gravitazionali continue attese dalle LMXB.