Nane bianche: superanto drasticamente limite di Chandrasekhar che resisteva da decenni.

Nel vasto panorama dell’astrofisica moderna, la massa limite di Chandrasekhar ($M_{Ch} \approx 1.44 M_{\odot}$) ha rappresentato per quasi un secolo una pietra miliare inamovibile. Formulata dal fisico Subrahmanyan Chandrasekhar negli anni ’30, questa soglia definisce la massa massima che una nana bianca può sostenere prima che la pressione di degenerazione degli elettroni ceda alla gravità, portando la stella a collassare o a esplodere. Recentemente, un team di ricercatori guidato da Zenia Zuraiq, Banibrata Mukhopadhyay e Christopher A. Tout ha pubblicato un lavoro scientifico rivoluzionario basato sul codice di evoluzione stellare STARS di Cambridge. Questo studio dimostra che le stelle dotate di forti campi magnetici possono superare tale limite fisicamente stabilito, dando origine a nane bianche super-Chandrasekhar o B-WDs con masse fino a $2.4 M_{\odot}$.

La portata di questa scoperta risuona direttamente nello studio delle supernovae di Tipo Ia, le celebri esplosioni cosmiche impiegate come candele standard per misurare le distanze nell’Universo e per comprendere l’espansione accelerata guidata dall’energia oscura. L’esistenza dimostrata di nane bianche magnetiche altamente massicce offre una spiegazione teorica rigorosa a eventi straordinariamente luminosi registrati dai telescopi, come le supernovae SN 2003fg e SN 2009dc, la cui origine era finora rimasta un enigma.

Stelle a fine vita ed equilibrio della materia degenerata

Per comprendere la portata del risultato, occorre esplorare cosa sia una nana bianca. Quando una stella con una massa simile a quella del nostro Sole o fino a circa otto masse solari esaurisce il suo combustibile nucleare al centro, espelle i suoi strati esterni nello spazio, formando una nebulosa planetaria. Ciò che rimane è un nucleo denso e caldissimo, composto prevalentemente da carbonio e ossigeno (oppure elio nelle stelle meno massicce).

In una stella comune come il Sole, la forza di gravità che schiaccia la materia verso il centro è bilanciata dalla pressione verso l’esterno generata dalle reazioni di fusione nucleare. In una nana bianca, dove le reazioni nucleari si sono spente, la gravità viene contrastata dalla pressione di degenerazione elettronica, un effetto derivante dalla meccanica quantistica e dal principio di esclusione di Pauli. Questo principio impedisce agli elettroni di occupare lo stesso stato quantico, creando una resistenza estrema alla compressione.

Tuttavia, il modello matematico tradizionale di Chandrasekhar ipotizzava che la stella fosse non rotante e priva di campo magnetico significativo. Con queste condizioni, se una nana bianca accresce materia fino a circa $1.38 M_{\odot}$, la densità centrale diventa tale da innescare la fusione termonucleare incontrollata del carbonio, distruggendo la stella in un’esplosione di tipo Ia.

Forza magnetica e struttura interna delle stelle B-WDs

L’elemento chiave introdotto dalle nuove simulazioni del codice STARS riguarda il ruolo attivo del campo magnetico centrale e superficiale. Quando un campo magnetico dipendente dalla densità viene incorporato nel bilancio idrostatico, agisce fornendo un supporto ulteriore contro la gravità.

La pressione totale all’interno della stella diventa la somma della pressione della materia degenerata e della pressione magnetica ($P_B = B^2 / 2\mu_0$), oltre all’effetto della tensione magnetica. Questa ulteriore spinta permette alla struttura stellare di sostenere una massa nettamente superiore rispetto a quella prevista per oggetti non magnetici, spingendo la curva massa-raggio verso regioni mai esplorate prima.

I ricercatori hanno modellato la geometria magnetica interna combinando una componente toroidale dominante nella zona profonda del nucleo e una componente poloidale limitata alla parte più esterna (l’ultimo 10%-30% del raggio). Questa configurazione mista garantisce la stabilità magnetoidrodinamica dell’astro nel corso dell’evoluzione, evitando che la stella subisca deformazioni asimmetriche distruttive o deviazioni marcate dalla forma sferica.

Evoluzione stellare e accrescimento nei sistemi binari

La formazione di una nana bianca magnetica super-Chandrasekhar avviene principalmente attraverso processi di evoluzione in sistemi binari con trasferimento di massa. Il flusso sequenziale che porta alla nascita di questi oggetti cosmici si articola nelle seguenti fasi:

1. Fase di sequenza principale: Una stella di massa compresa tra $3 M_{\odot}$ e $8 M_{\odot}$ possiede un campo magnetico iniziale dormiente, con valori centrali compresi tra $10^6$ G e $10^7$ G e valori superficiali moderati.

2. Evoluzione lungo il ramo gigante e AGB: Durante la fase di ramo asintotico delle giganti (AGB), la stella perde i suoi strati esterni tramite intensi venti stellari o spogliazione causata dall’interazione gravitazionale con una stella compagna.

3. Formazione della nana bianca con campo magnetico dormiente: Si crea un nucleo degenerato di carbonio-ossigeno (CO) o elio-carbonio-ossigeno (HeCO) con una massa sub-Chandrasekhar (ad esempio $0.45 M_{\odot}$ o $1.02 M_{\odot}$). Il campo magnetico rimane idrodinamicamente dormiente ($\le 5 \times 10^{13}$ G) e non altera l’equilibrio iniziale.

4. Fase di accrescimento di materia: La nana bianca inizia ad aspirare gas idrogeno o elio/carbonio dalla stella compagna a un tasso costante ($\dot{M} \approx 10^{-9} M_{\odot} \text{ yr}^{-1}$).

5. Attivazione del campo magnetico per congelamento del flusso: Man mano che la stella accumula massa e si contrae, la densità centrale aumenta. Per la conservazione del flusso magnetico (flux freezing), il campo magnetico centrale cresce fino a superare $10^{14}$ G.

6. Raggiungimento della massa super-Chandrasekhar: Grazie al sostegno magnetico aggiunto, la stella supera la soglia di $1.44 M_{\odot}$ senza innescare precocemente la fusione del carbonio, stabilizzandosi su valori fino a $2.4 M_{\odot}$ prima del punto di saturazione.

Anomalie dei raggi osservati nelle nane bianche fredde

Oltre a spiegare l’esistenza di nane bianche super-massicce, il modello teorico sviluppato risolve un dilemma astrofisico legato alle nane bianche sub-Chandrasekhar che mostrano un raggio anomalo rispetto alla loro massa.

L’analisi dei dati di osservazione su stelle come SDSS J123204.19+522548.2 (con una temperatura efficace $T_e = 8750\text{ K}$) aveva rivelato una discrepanza evidente rispetto alla relazione massa-raggio standard. Considerare soltanto gli effetti della temperatura finita non era sufficiente per spiegare l’espansione del raggio, lasciando un margine di errore vicino al 7%.

L’introduzione congiunta della temperatura finita e della pressione magnetica toroidale permette di riprodurre perfettamente le dimensioni misurate. Il campo magnetico agisce dilatando la struttura esterna dell’astro. Nel caso specifico di SDSS J123204.19+522548.2, l’integrazione del fattore magnetico riduce lo scarto teorico a un valore inferiore allo 0.7%, fornendo la prima spiegazione fisica completa coerente con le osservazioni.

Composizione chimica ed effetti dei bruciamenti di elio

La capacità di una nana bianca magnetica di raggiungere il massimo limite di massa dipende fortemente dalla composizione chimica del nucleo e del materiale accresciuto. Le simulazioni eseguite con STARS hanno messo a confronto due scenari distinti:

• Nane bianche a carbonio-ossigeno (CO): Formate a partire da stelle di massa iniziale attorno a $8 M_{\odot}$, mantengono strati di elio estremamente ridotti. L’accrescimento di materiale ricco di carbonio e ossigeno procede senza che si inneschino reazioni termonucleari instabili in superficie. In questo caso, il supporto magnetico consente di raggiungere una massa finale di circa $2.4 M_{\odot}$ prima che la densità centrale inneschi la fusione del carbonio.

• Nane bianche ibride ad elio-carbonio-ossigeno (HeCO): Derivate da stelle di sequenza principale meno massicce (circa $3 M_{\odot}$), mantengono una frazione residua di elio non trascurabile ($\sim 0.1 M_{\odot}$). Durante la fase di accrescimento, l’aumento di pressione e temperatura nella zona esterna provoca l’accensione instabile dell’elio. Questo fenomeno destabilizza la stella ben prima di raggiungere il limite teorico di massa, provocando eruzioni di tipo nova o detonazioni superficiali che interrompono la crescita.

Le simulazioni confermano che, disattivando artificialmente le reazioni nucleari dell’elio nel codice, la nana bianca HeCO riprende a crescere lungo la curva di equilibrio magnetico, dimostrando in modo inequivocabile che l’arresto prematuro dell’accrescimento è causato dalla reattività del combustibile nucleare e non da un’instabilità magnetica o gravitazionale.

Implicazioni cosmologiche per le candele standard e l’energia oscura

L’impatto di questa scoperta si estende al campo della cosmologia osservativa. Le supernovae di Tipo Ia sono state lo strumento fondamentale che ha consentito di scoprire l’espansione accelerata dell’Universo e di postulari l’esistenza dell’energia oscura. La loro utilità come misuratori cosmici si basa sul presupposto che la stella progenitrice esploda sempre al raggiungimento della medesima massa critica ($M_{Ch} \approx 1.44 M_{\odot}$), generando una luminosità di picco uniforme e calibrabile.

L’accertata presenza di progenitori magnetici super-Chandrasekhar altera questo paradigma:

• Le nane bianche altamente magnetiche possono esplodere a masse vicine a $2.4 M_{\odot}$, liberando una quantità decisamente maggiore di nichel-56 ($^{56}\text{Ni}$) e producendo eventi straordinariamente più luminosi della media.

• Le supernovae di Tipo Ia non sono tutte identiche, ma presentano una variabilità intrinseca legata al campo magnetico della stella progenitrice.

• I parametri cosmologici derivati, compresa la stima precisa della costante di Hubble e della tensione ad essa associata, potrebbero richiedere una ricalibrazione teorica che tenga conto della popolazione di B-WDs.

Approfondimenti tecnici e formalismo matematico dal codice STARS

Per gli specialisti della materia, vengono riportati i dettagli tecnici, i coefficienti e le equazioni differenziali fondamentali implementate nel codice STARS per simulare la struttura magnetica ed evolutiva delle nane bianche.

Equazioni di struttura stellare idrostatica e termica

Il sistema di equazioni differenziali alla base della struttura unidimensionale modificate per tenere conto degli effetti magnetici e delle perdite di energia è espresso come:

Dove:

• $m$ rappresenta la massa racchiusa all’interno del raggio $r$ al tempo $t$.

• $P$ è la pressione totale, che combina la pressione della materia degenerata $P_m$ e il contributo della pressione magnetica $P_B$.

• $T$ è la temperatura locale, $\rho$ la densità di massa, $L$ la luminosità integrata.

• $\nabla$ indica il gradiente termico logaritmico con la pressione, definito in base al regime locale (radiativo o convettivo).

• $\epsilon$ descrive il tasso di generazione di energia tramite la rete di 21 reazioni nucleari integrate nel codice.

• $\epsilon_{\nu}$ rappresenta la perdita di energia causata dall’emissione di neutrini.

• $u$ è l’energia interna specifica per unità di massa.

Profilo di campo magnetico dipendente dalla densità

Il campo magnetico interno $B(\rho)$ viene modellato tramite una funzione analitica legata alla densità locale del plasma, atta a simulare l’effetto di congelamento del flusso durante la contrazione stellare:

I parametri del modello adottati nelle simulazioni principali sono:

• $B_0 = 10^{14}\text{ G}$, rappresenta l’intensità di scala del campo magnetico centrale.

• $\rho_0 = 10^9\text{ g cm}^{-3}$, rappresenta la densità di riferimento di saturazione.

• $B_s$ indica il campo magnetico superficiale, mantenuto coerente con la conservazione del flusso magnetico ideale ($B_s R^2 = \text{costante}$). Per una stella di sequenza principale con raggio pari a quello solare ($R = R_{\odot}$) si imposta $B_s \sim 10\text{ G} – 1000\text{ G}$, portando a valori superficiali di nana bianca compresi tra $10^4\text{ G}$ e $10^9\text{ G}$.

• $\eta$ e $\gamma$ sono parametri adimensionali che regolano la pendenza e la forma del profilo di campo magnetico. Nelle simulazioni esaminate si sono impiegati valori di $\eta = 0.7, 1.1, 1.3$ e $\gamma = 0.8, 0.9$.

Valutazione dell’ellitticità e della deviazione dalla simmetria sferica

Per garantire la validità dell’approssimazione unidimensionale del codice STARS, l’ellitticità indotta dalla perturbazione di densità magnetica $\delta\rho$ è stata calcolata tramite l’integrale:

Dove:

• $\delta\rho = \frac{B^2}{2 \mu_0}$ è la perturbazione di densità causata dal campo magnetico.

• $I_0$ è il momento di inerzia imperturbato della stella.

I risultati numerici confermano che $\epsilon \le 10^{-7}$ per le nane bianche appena formate e $\epsilon \le 10^{-5}$ per le nane bianche super-massicce da $2.4 M_{\odot}$. Tali valori estremamente bassi confermano che la deviazione dalla simmetria sferica è trascurabile, giustificando pienamente l’uso di codici idrodinamici 1D.

Scale temporali di decadimento del campo magnetico

La stabilità temporale del campo magnetico è stata verificata integrando gli effetti accoppiati della dissipazione Ohmica e dell’effetto Hall:

Le due scale temporali caratteristiche sono definite da:

Dove:

• $\rho_{c,n} = \rho_c / 10^n\text{ g cm}^{-3}$, con $\rho_c$ densità centrale.

• $\rho_{\text{av}}$ è la densità media dell’astro.

• $R_4 = R / 10^4\text{ km}$ descrive il raggio scalato.

• $T_{c,7} = T_c / 10^7\text{ K}$ rappresenta la temperatura centrale.

• $B_{0,14} = B_0 / 10^{14}\text{ G}$ è l’intensità del campo centrale.

• $l_8 = l / 10^8\text{ cm}$ indica la scala di lunghezza caratteristica dei cappi di flusso magnetico.

I calcoli indicano che il decadimento significativo del campo magnetico richiede oltre 5 miliardi di anni ($5\text{ Gyr}$), mentre la fase di accrescimento che porta alla formazione delle nane bianche super-Chandrasekhar dura al massimo 1 miliardo di anni ($1\text{ Gyr}$). Questo divario temporale garantisce che il campo magnetico rimanga integro e protetto per l’intera durata della crescita della stella.

Traduzione abstract. La disponibilità dello studio completo, in accesso libero o riservato, è verificabile tramite il link del DOI nel box sotto se disponibile.

Abstract

Esploriamo l’evoluzione temporale delle stelle di sequenza principale (MS) magnetizzate in nane bianche (WD) magnetizzate non rotanti con una serie di modelli che descrivono un diagramma di Hertzsprung-Russell modificato. Utilizziamo il Cambridge Stellar Evolution Code, STARS, con modifiche che introducono effetti magnetici e raffreddamento. Indaghiamo la possibile esistenza di nane bianche super-Chandrasekhar, massicce, stabili e magnetizzate, approssimativamente sferiche, al di fuori della curva massa-raggio di Chandrasekhar. Queste vengono definite B-WD. La loro esistenza è dedotta dalle osservazioni di peculiari supernove di tipo Ia sovraluminose, ad esempio SN 2003fg. Dimostriamo che gli effetti classici dei campi magnetici influenzano significativamente le MS e la loro evoluzione temporale in WD. Inoltre, il campo magnetico potrebbe essere decisivo per spiegare il raggio insolitamente grande osservato in alcune nane bianche sub-Chandrasekhar, ad esempio SDSS J123204.19+522548.2, insieme all’effetto della temperatura finita. Abbiamo ulteriormente studiato la stabilità delle nane bianche magnetizzate, sia al di sotto che al di sopra del limite di Chandrasekhar, e confermato la possibilità di nane bianche di tipo B super-Chandrasekhar con nuovi limiti di massa che dipendono dalla geometria del campo magnetico, dall’accrescimento di massa da una compagna e dai tassi di raffreddamento. Pertanto, sebbene il “limite di Chandrasekhar” possa essere sacrosanto, il suo valore dipende fortemente dalla fisica sottostante, in questo caso il campo magnetico.

by