Resolving core-envelope degeneracies in giant planets and sub-Neptunes: constraining the equivalence in the presence of dilute gradients

Sintesi Tecnica: Risoluzione delle Degenerazioni Core-Involucro nei Giganti Gassosi e nei Sub-Nettuno

Definizione del Problema
La struttura interna dei pianeti giganti e dei sub-Nettuno rappresenta un problema altamente degenerato e multidimensionale nella scienza planetaria. Il riconoscimento di nuclei "sfocati" o diluiti — dove gli elementi pesanti sono distribuiti in gradienti anziché con confini netti — e di massicci mantelli volatili ha complicato l'interpretazione dei raggi planetari osservati ($R_p$) e delle masse ($M_p$). Colmare il "gap del raggio" (radius valley) richiede di distinguere tra pianeti con spessi involucri di H/He, oceani di acqua supercritica o composizioni miste. Gli attuali modelli evolutivi (codici di classe Henyey come MESA, APPLE, CEPAM) sono computazionalmente costosi, il che spesso proibisce le densissime scansioni dei parametri necessarie per mappare l'intero spazio delle degenerazioni strutturali statiche, in particolare riguardo alla larghezza dei gradienti composizionali ($\sigma$). Inoltre, la specifica distribuzione spaziale degli elementi pesanti influenza significativamente il raggio totale e i momenti di gravità ($J_n$), sebbene la topologia funzionale di tali gradienti rimanga scarsamente vincolata dalla sola osservazione.

Metodologia
Gli autori presentano fuzzycore, un integratore strutturale idrostatico multi-zona, statico e open-source, progettato per l'esplorazione rapida dei modelli diretti delle degenerazioni interne. A differenza dei codici evolutivi, fuzzycore decoppia il profilo strutturale dal trasporto di energia tempo-dipendente, dal mixing convettivo e dalle curve critiche di equilibrio di fase.

Componenti metodologiche chiave includono:

• Equilibrio Idrostatico Multi-Zona: Il framework risolve equazioni 1D a simmetria sferica per la conservazione della massa e l'equilibrio idrostatico. Adatta dinamicamente la risoluzione della griglia per risolvere i ripidi gradienti di densità attraverso le interfacce.
• Formulazione dell'Equazione di Stato (EOS):
  • Nucleo Misto: Utilizza una legge del volume additivo (legge di Amagat) per miscelare tabelle ab initio di ferro e roccia multi-fase (separando MgSiO3 in MgO e SiO2 ad alte pressioni) per garantire una stratificazione della densità monotona.
  • Involucro Fluido: Modella l'involucro come una miscela ternaria (H, He, elementi pesanti) utilizzando la miscelazione lineare ideale dei volumi specifici. Gli elementi pesanti ($Z$) sono rappresentati da una EOS dell'H2O a vasto intervallo per catturare gli effetti della massa molecolare media ($\mu$) e della comprimibilità nel passaggio dall'atmosfera pristina agli strati basali arricchiti.
  • Mantello Condensato: Supporta un mantello d'acqua distinto e separato in fase (Ghiaccio VII/X o superionico) per i sub-Nettuno ricchi di volatili, accoppiato termodinamicamente alla base dell'involucro.
• Gradienti Composizionali: Invece di confini discreti, la frazione di massa degli elementi pesanti $Z(x)$ è modellata tramite una funzione di decadimento gaussiano continua. Per gestire la stabilità numerica quando i gradienti si fanno più ripidi ($\sigma \to 0$), il framework impiega una media analitica sub-griglia per conservare il budget degli elementi pesanti e prevenire picchi di densità non fisici.
• Integrazione Numerica: L'involucro è discretizzato in una serie di zone convettive completamente miscelate a "scala". Un robusto metodo di stepper adiabatico, basato su regressioni lineari dei minimi quadrati locali su un KD-Tree dello spazio di fase 3D ($P, T, S$), integra gli stati termodinamici.
• Boundary Shooting (Metodo del Tiro al Bersaglio sui Confini): Un metodo di boundary-shooting nidificato (metodo di Brent) risolve il problema del valore al contorno, iterando sulle pressioni di interfaccia per soddisfare la continuità di massa, pressione e temperatura tra l'atmosfera, il mantello e il nucleo. Supporta le modalità di puntamento della massa (mass-targeting) e di puntamento della gravità superficiale (surface-gravity-targeting).
• Guardie di Stabilità: Include una "guarda isotermica" per prevenire artefatti di raffreddamento numerico alle transizioni di fase a bassa pressione e un campionamento adattivo per involucri ultra-sottili.

Contributi Chiave e Risultati

1. Robustezza della Forma Funzionale: Il benchmarking rispetto a un profilo di Giove derivato dall'evoluzione (Vazan et al. 2018) dimostra che, quando vengono corrisponduti le metallicità del confine dell'involucro e i budget integrati degli elementi pesanti, il raggio macroscopico è robusto rispetto alla precisa topologia funzionale del gradiente. Sei diverse famiglie parametriche (gaussiana, sigmoide, esponenziale, bilineare, a singolo step e scale idrodinamiche digitalizzate) hanno prodotto una dispersione del raggio di solo il 2,3%. Ciò valida l'uso di parametrizzazioni analitiche lisce (come il decadimento gaussiano) per mappare le degenerazioni.
2. Atlante delle Degenerazioni dei Sub-Nettuno: Il framework è stato applicato a un pianeta fisso di $8 M_\oplus$ per generare un "atlante delle degenerazioni dell'acqua" attraverso una griglia $15 \times 15 \times 12$ di massa d'acqua, massa dell'involucro e larghezza di diluizione ($\sigma$).
   • $\mu$-Compattazione: Lo studio quantifica come la ridistribuzione degli elementi pesanti da un nucleo compatto a un involucro diluito aumenti la massa molecolare media, compattando fisicamente l'atmosfera. Questa "$\mu$-compattazione" può sopprimere l'espansione termica, riducendo il raggio planetario fino al 45% per budget di elementi pesanti fissi.
   • Copertura Termica (Thermal Blanketing): Gli involucri ricchi di metalli agiscono come freni termici, aumentando la differenza di temperatura tra la superficie e il nucleo di circa $\sim 5000$ K, sebbene si tratti di un reperto strutturale statico, non di una storia di raffreddamento.
   • Limiti Osservativi: Per un dato $(M_p, R_p)$ osservato, la larghezza di diluizione $\sigma$ è essenzialmente non vincolata. L'atlante rivela una famiglia continua di architetture statiche (variando massa d'acqua, massa dell'involucro e $\sigma$) che riproducono lo stesso raggio, evidenziando che $(M_p, R_p)$ da soli non possono risolvere la struttura interna.
3. Applicazione a Target Specifici:
   • GJ 1214 b (Ricco di volatili): Modellato come un corpo acqua-roccia compatto con un involucro sottile e altamente arricchito, convergendo a un raggio entro l'1% delle osservazioni.
   • Kepler-11 e (Ricco di gas/Super-puff): Modellato come un nucleo massiccio con un involucro diluito e pesantemente inquinato. Il modello dimostra che è richiesto un gradiente di elementi pesanti sospeso profondamente per "zavorrare" l'involucro massiccio verso il raggio compatto osservato, contrastando l'inflazione attesa da un involucro di H/He pristino.
4. Accoppiamento Atmosfera-Interno: Il framework è progettato per servire come backend per i recuperi atmosferici (es. accoppiamento con HADES). Può collegare i prior di metallicità atmosferica di JWST/Ariel alle architetture interne profonde, risolvendo per la gravità superficiale e le condizioni al contorno auto-coerenti.

Significato e Ambito
L'articolo posiziona fuzzycore non come un sostituto dei codici evolutivi completamente impliciti, ma come uno strumento complementare necessario. Eliminando l'overhead computazionale del mixing tempo-dipendente e del trasporto di energia, fuzzycore consente la mappatura sistematica e ad alta risoluzione delle degenerazioni strutturali statiche che sarebbero altrimenti troppo costose da esplorare con i solver di classe Henyey.

Gli autori sottolineano che fuzzycore identifica configurazioni idrostatiche meccanicamente ammesse. Non determina se tali configurazioni siano termodinamicamente coerenti con la fisica del demixing o idrodinamicamente stabili contro l'erosione turbolenta su scale temporali geologiche; queste domande rimangono dominio della modellazione evolutiva e degli studi di equilibrio di fase ab initio. Tuttavia, il framework fornisce l'essenziale "atlante delle degenerazioni" necessario per interpretare i dati osservativi attuali e futuri, dimostrando che la distribuzione spaziale degli elementi pesanti è una dimensione critica, sebbene osservativamente degenerata, nella caratterizzazione planetaria.