Self-Consistent Evolution Models Show Weak Double-Diffusive Mixing in Jupiter and Saturn

Resumen Técnico: Los modelos de evolución autoconsistentes muestran una mezcla doble-difusiva débil en Júpiter y Saturno

Planteamiento del problema
La estructura interna de los planetas gigantes, particularmente la existencia de núcleos "difusos" o extendidos en Júpiter y Saturno, sigue siendo un desafío significativo para las ciencias planetarias. Aunque se ha propuesto la convección doble-difusiva como un mecanismo para redistribuir elementos pesados desde un núcleo denso hacia el sobrecapa, su eficiencia en contextos evolutivos realistas es incierta. Estudios previos han dependido de estructuras composicionales idealizadas y no han tratado el transporte doble-difusivo de manera autoconsistente dentro de los cálculos de evolución planetaria. Además, existe una discrepancia entre las simulaciones hidrodinámicas locales, que a menudo predicen una fusión rápida de capas y una mezcla eficiente, y los requisitos de los modelos evolutivos globales que deben reproducir las restricciones observacionales actuales.

Metodología
Los autores implementan un tratamiento autoconsistente del transporte doble-difusivo dentro del código de evolución planetaria APPLE.

• Condiciones iniciales: El estudio utiliza modelos de interior post-formación para Júpiter (basados en Stevenson et al. 2022) y Saturno (basados en Bodenheimer et al. 2025). Estos modelos presentan gradientes composicionales extendidos ("núcleos difusos") resultantes de la disrupción parcial y evaporación de planetesimales acrecionados durante la formación, en lugar de límites de núcleo nítidos.
• Física del transporte: Los autores incorporan la prescripción de convección por capas derivada por Spruit (2013). Este enfoque modela el transporte a través de escaleras de convección como un proceso de difusión mejorado. Las difusividades efectivas de composición ($D_{\text{eff}}$) y térmica ($\kappa_{T,\text{eff}}$) se derivan de regímenes asintóticos donde el transporte es independiente del espesor de la capa, escalando con la razón de estabilidad ($R_\rho$) y las difusividades microscópicas.
• Simulación evolutiva: Los modelos se evolucionan durante 4.56 Gyr (la edad del Sistema Solar). El código resuelve simultáneamente la evolución térmica y el transporte composicional utilizando una rejilla de masa lagrangiana. Las regiones inestables a la convección por inversión (criterio de Ledoux) se tratan con la teoría de longitud de mezcla estándar, mientras que las regiones estables según Ledoux con $1 < R_\rho < \tau^{-1/2}$ utilizan la prescripción doble-difusiva.
• Análisis de sensibilidad: Para probar la robustez de los resultados, los autores ejecutan simulaciones donde las difusividades efectivas se incrementan por factores de 100 y 1000 respecto a los valores fiduciarios de Spruit (2013).

Resultados clave

• Mezcla ineficiente: Las simulaciones revelan que la convección doble-difusiva produce una mezcla muy limitada entre el interior profundo y la envoltura. A lo largo de la historia de enfriamiento de 4.56 Gyr, se redistribuye menos de $\sim 1 M_\oplus$ de material pesado tanto en Júpiter como en Saturno.
  • En el modelo fiduciario, solo se redistribuyen aproximadamente $\approx 0.11 M_\oplus$ (Júpiter) y $\approx 0.07 M_\oplus$ (Saturno).
  • Incluso con difusividades incrementadas por un factor de 1000, la masa redistribuida sigue siendo modesta ($\approx 2 M_\oplus$ para Júpiter y $\approx 1.3 M_\oplus$ para Saturno).
• Mecanismo de autolimitación: La ineficiencia surge de un bucle de retroalimentación negativa. El transporte térmico mejorado modifica la estratificación térmica, reduciendo la razón de estabilidad ($R_\rho$) y, por ende, suprimiendo las difusividades efectivas. Esto evita una mezcla descontrolada, asegurando que incluso coeficientes de transporte irrealmente grandes no logren homogeneizar los gradientes composicionales primordiales.
• Restricción energética: Los autores identifican una restricción energética fundamental. En la aproximación de Boussinesq utilizada en muchas simulaciones hidrodinámicas, la mezcla es energéticamente barata porque el espesor de la capa es pequeño comparado con la escala térmica. Sin embargo, en los interiores compresibles de los planetas gigantes, la región estable abarca una fracción sustancial del radio planetario ($H \sim 0.25\text{--}0.5 R_p$), lo que hace que la energía requerida para mezclar el gradiente composicional sea comparable al presupuesto de energía térmica disponible. Consecuentemente, el trabajo de flotabilidad disponible para impulsar el transporte composicional está severamente restringido.

Significado y afirmaciones
El artículo afirma que la convección doble-difusiva, cuando se trata de forma autoconsistente en un contexto evolutivo global, no puede erosionar significativamente los gradientes composicionales generados durante la formación planetaria. Las distribuciones de elementos pesados observadas en Júpiter y Saturno probablemente requieren mecanismos de transporte adicionales o vías de formación que produzcan una mezcla inicial más amplia de lo que predicen los modelos de acreción estándar.

Los autores sugieren que los núcleos "difusos" pueden ser un producto directo del proceso de formación (por ejemplo, una mezcla temprana más extensa durante la ablación de planetesimales) o el resultado de grandes eventos colisionales, en lugar de ser un producto de la mezcla doble-difusiva post-formación. El estudio destaca una limitación crítica al aplicar los resultados de la hidrodinámica local (a menudo basados en la aproximación de Boussinesq) a la evolución global de los planetas gigantes, enfatizando que el presupuesto energético del interior profundo es el regulador principal de la eficiencia de la mezcla.