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

Resumen Técnico: Resolución de las Degeneraciones Núcleo-Envoltura en Planetas Gigantes y Sub-Neptunos

Planteamiento del Problema
La estructura interna de los planetas gigantes y los sub-neptunos representa un problema altamente degenerado y multidimensional en la ciencia planetaria. El reconocimiento de núcleos "difusos" o diluidos —donde los elementos pesados se distribuyen en gradientes en lugar de fronteras nítidas— y de mantos volátiles masivos ha complicado la interpretación de los radios ($R_p$) y masas ($M_p$) observados. Para cerrar la "brecha de radio" (radius valley), es necesario distinguir entre planetas con atmósferas espesas de H/He, océanos de agua supercrítica o composiciones mixtas. Los modelos evolutivos actuales (códigos de clase Henyey como MESA, APPLE, CEPAM) son computacionalmente costosos, lo que a menudo impide los barridos de parámetros densos necesarios para mapear todo el espacio de las degeneraciones estructurales estáticas, particularmente respecto al ancho de los gradientes de composición ($\sigma$). Además, la distribución espacial específica de los elementos pesados afecta significentemente el radio total y los momentos de gravedad ($J_n$), pero la topología funcional de estos gradientes permanece insuficientemente constreñida por la observación por sí sola.

Metodología
Los autores presentan fuzzycore, un integrador estructural hidrostático estático, multizona y de código abierto, diseñado para la exploración rápida de las degeneraciones de modelos directos (forward-modeling). A diferencia de los códigos evolutivos, fuzzycore desacopla el perfilado estructural del transporte de energía dependiente del tiempo, la mezcla convectiva y las curvas críticas de equilibrio de fase.

Componentes metodológicos clave:

• Equilibrio Hidrostático Multizona: El marco resuelve ecuaciones de 1d con simetría esférica para la conservación de la masa y el equilibrio hidrostático. Adapta dinámicamente la resolución de la rejilla para resolver gradientes de densidad pronunciados a través de las interfaces.
• Formulación de la Ecuación de Estado (EOS):
  • Núcleo Mixto: Utiliza una ley de volumen aditiva (ley de Amagat) para mezclar tablas ab initio de hierro y roca multifásica (separando MgSiO3 en MgO y SiO2 a altas presiones) para asegurar una estratificación de densidad monotónica.
  • Envoltura Fluida: Modela la envoltura como una mezcla ternaria (H, He, elementos pesados) utilizando la mezcla lineal ideal de volúmenes específicos. Los elementos pesados ($Z$) se aproximan mediante una EOS de H2O de amplio rango para capturar los efectos del peso molecular medio ($\mu$) y la compresibilidad a través de la transición desde una atmósfera prístina hasta las capas basales enriquecidas.
  • Manto Condensado: Soporta un manto de agua con fases separadas (Hielo VII/X o superiónico) para sub-neptunos ricos en volátiles, acoplado termodinámicamente a la base de la envoltura.
• Gradientes de Composición: En lugar de fronteras discretas, la fracción de masa de elementos pesados $Z(x)$ se modela mediante una función de decaimiento gaussiano continua. Para manejar la estabilidad numérica a medida que los gradientes se vuelven más agudos ($\sigma \to 0$), el marco emplea un promedio analítico de sub-rejilla para conservar el presupuesto de elementos pesados y prevenir picos de densidad no físicos.
• Integración Numérica: La envoltura se discretiza en una "escalera" de zonas convectivas totalmente mezcladas. Un método robusto de paso adiabático (adiabat stepper) mediante regresiones lineales de mínimos cuadrados locales sobre un KD-Tree del espacio de fase 3D ($P, T, S$) integra los estados termodinámicos.
• Disparo de Frontera (Boundary Shooting): Un método de disparo de frontera anidado (método de Brent) resuelve el problema de valor en la frontera, iterando sobre las presiones de interfaz para igualar la continuidad de masa, presión y temperatura entre la atmósfera, el manto y el núcleo. Soporta modos de fijación de masa y de fijación de gravedad superficial.
• Guardas de Estabilidad: Incluye una "guarda isotérmica" para evitar artefactos de enfriamiento numérico en las transiciones de fase de baja presión y un submuestreo adaptativo para envolturas ultra delgadas.

Contribuciones Clave y Resultados

1. Robustez de la Forma Funcional: La comparación con un perfil de Júpiter derivado de la evolución (Vazan et al. 2018) demuestra que, cuando se igualan las metalicidades de la frontera de la envoltura y los presupuestos integrados de elementos pesados, el radio planetario macroscópico es robusto a la topología exacta del gradiente. Seis familias paramétricas distintas (decaimiento gaussiano, sigmoide, exponencial, bilineal, de paso único y escaleras hidrodinámicas digitalizadas) produjeron una dispersión de radio de solo un 2.3%. Esto valida el uso de parametrizaciones analíticas suaves (como el decaimiento gaussiano) para mapear las degeneraciones.
2. Atlas de Degeneración de Sub-Neptunos: El marco se aplicó a un planeta fijo de $8 M_\oplus$ para generar un "atlas de degeneración de mundos de agua" a través de una rejilla de $15 \times 15 \times 12$ de masa de agua, masa de envoltura y ancho de dilución ($\sigma$).
   • $\mu$-Compresión: El estudio cuantifica cómo la redistribución de elementos pesados desde un núcleo compacto hacia una envoltura diluida aumenta el peso molecular medio, compactando físicamente la atmósfera. Esta "$\mu$-compresión" puede suprimir la expansión térmica, reduciendo el radio planetario hasta en un 45% para presupuestos de elementos pesados fijos.
   • Manto Térmico: Las envolturas ricas en metales actúan como frenos térmicos, aumentando la diferencia de temperatura entre la superficie y el núcleo en aproximadamente 5000 K, aunque este es un hallazgo estructural estático, no una historia de enfriamiento.
   • Límites Observacionales: Para un $(M_p, R_p)$ observado fijo, el ancho de dilución $\sigma$ es esencialmente no restringido. El atlas revela una familia continua de arquitecturas estáticas (variando la masa de agua, la masa de la envoltura y $\sigma$) que reproducen el mismo radio, resaltando que $(M_p, R_p)$ por sí solo no puede resolver la estructura interna.
3. Aplicación a Objetivos Específicos:
   • GJ 1214 b (Rico en volátiles): Modelado como un cuerpo de agua-roca compacto con una envoltura delgada y altamente enriquecida, convergiendo a un radio dentro del 1% de las observaciones.
   • Kepler-11 e (Rico en gas/Super-puff): Modelado como un núcleo masivo con una envoltura diluida y fuertemente contaminada. El modelo demuestra que se requiere un gradiente de elementos pesados profundamente suspendido para "lastrar" la envoltura masiva hacia su radio compacto observado, contrarrestando la inflación esperada de una envoltura de H/He prístina.
4. Acoplamiento Atmósfera-Interior: El marco está diseñado para servir como trasfondo para las recuperaciones atmosféricas (ej. acoplamiento con HADES). Puede vincular los priors de metalicidad atmosférica de JWST/Ariel con arquitecturas de interiores profundos mediante la resolución de la gravedad superficial y las condiciones de frontera autoconsecuentes.

Significancia y Alcance
El artículo posiciona a fuzzycore no como un reemplazo para los códigos evolutivos totalmente implícitos, sino como una herramienta complementaria necesaria. Al eliminar la carga computacional del mezclado dependiente del tiempo y el transporte de energía, fuzzycore permite el mapeo sistemático y de alta resolución de las degeneraciones estructurales estáticas que de otro modo serían demasiado costosas de explorar con resolvedores de clase Henyey.

Los autores enfatizan que fuzzycore identifica configuraciones hidrostáticas mecánicamente permitidas. No determina si estas configuraciones son termodinámicamente consistentes con la física de desmezcla o hidrodinámicamente estables contra la erosión turbulenta a lo largo de escalas de tiempo geológicas; esas preguntas siguen siendo el dominio del modelado evolutivo y los estudios de equilibrio de fase ab initio. Sin embargo, el marco proporciona el "atlas de degeneración" esencial requerido para interpretar los datos observacionales actuales y futuros, demostrando que la distribución espacial de los elementos pesados es una dimensión crítica, aunque observacionalmente degenerada, en la caracterización planetaria.