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

Résumé technique : Les modèles d'évolution auto-cohérents montrent un mélange double-diffusif faible dans Jupiter et Saturne

Énoncé du problème
La structure interne des planètes géantes, en particulier l'existence de noyaux « flous » ou étendus dans Jupiter et Saturne, demeure un défi majeur pour les sciences planétaires. Bien que la convection double-diffusive ait été proposée comme un mécanisme capable de redistribuer les éléments lourds d'un noyau dense vers l'enveloppe sus-jacente, son efficacité dans des contextes évolutifs réalistes est incertaine. Les études précédentes se sont appuyées sur des structures compositionnelles idéalisées et n'ont pas traité le transport double-diffusif de manière auto-cohérente au sein des calculs d'évolution planétaire. De plus, il existe une divergence entre les simulations hydrodynamiques locales, qui prédisent souvent une fusion rapide des couches et un mélange efficace, et les exigences des modèles d'évolution globaux qui doivent reproduire les contraintes observationnelles actuelles.

Méthodologie
Les auteurs implémentent un traitement auto-cohérent du transport double-diffusif au sein du code d'évolution planétaire APPLE.

• Conditions initiales : L'étude utilise des modèles d'intérieur post-formation pour Jupiter (basés sur Stevenson et al. 2022) et Saturne (basés sur Bodenheimer et al. 2025). Ces modèles présentent des gradients compositionnels étendus (« noyaux flous ») résultant de la disruption partielle et de l'évaporation des planétésimaux accrétés lors de la formation, plutôt que des limites de noyau nettes.
• Physique du transport : Les auteurs incorporent la prescription de convection par couches dérivée de Spruit (2013). Cette approche modélise le transport à travers les escaliers de convection comme un processus de diffusion amélioré. Les diffusivités effectives compositionnelles ($D_{\text{eff}}$) et thermiques ($\kappa_{T,\text{eff}}$) sont dérivées de régimes asymptotiques où le transport est indépendant de l'épaisseur des couches, évoluant en fonction du rapport de stabilité ($R_\rho$) et des diffusivités microscopiques.
• Simulation évolutive : Les modèles sont évolués pendant 4,56 Gyr (l'âge du Système solaire). Le code résout simultanément l'évolution thermique et le transport compositionnel en utilisant une grille de masse lagrangienne. Les régions instables à la convection par renversement (critère de Ledoux) sont traitées avec la théorie de la longueur de mélange standard, tandis que les régions stables au sens de Ledoux ($1 < R_\rho < \tau^{-1/2}$) utilisent la prescription double-diffusive.
• Analyse de sensibilité : Pour tester la robustesse des résultats, les auteurs exécutent des simulations où les diffusivités effectives sont augmentées par des facteurs de 100 et 1000 par rapport aux valeurs fiduciaires de Spruit (2013).

Résultats clés

• Mélange inefficace : Les simulations révèlent que la convection double-diffusive produit un mélange très limité entre l'intérieur profond et l'enveloppe. Sur l'ensemble de l'histoire de refroidissement de 4,56 Gyr, moins d'une masse de $\sim 1 M_\oplus$ de matière lourde est redistribuée dans Jupiter et Saturne.
  • Dans le modèle fiduciel, seulement $\approx 0,11 M_\oplus$ (Jupiter) et $\approx 0,07 M_\oplus$ (Saturne) sont redistribués.
  • Même avec des diffusivités augmentées d'un facteur 1000, la masse redistribuée reste modeste ($\approx 2 M_\oplus$ pour Jupiter et $\approx 1,3 M_\oplus$ pour Saturne).
• Mécanisme d'auto-limitation : L'inefficacité provient d'une boucle de rétroaction négative. Le transport thermique amélioré modifie la stratification thermique, réduisant le rapport de stabilité ($R_\rho$) et supprimant ainsi les diffusivités effectives. Cela empêche un mélange incontrôlé, garantissant que même des coefficients de transport irréalistement élevés ne parviennent pas à homogénéiser les gradients compositionnels primordiaux.
• Contrainte énergétique : Les auteurs identifient une contrainte énergétique fondamentale. Dans l'approximation de Boussinesq utilisée dans de nombreuses simulations hydrodynamiques, le mélange est peu coûteux en énergie car l'épaisseur de la couche est petite par rapport à l'échelle de hauteur thermique. Cependant, dans les intérieurs de planètes géantes compressibles, la région stable s'étend sur une fraction substantielle du rayon planétaire ($H \sim 0,25\text{--}0,5 R_p$), rendant l'énergie requise pour mélanger le gradient compositionnel comparable au budget d'énergie thermique disponible. Par conséquent, le travail de flottabilité disponible pour piloter le transport compositionnel est sévèrement contraint.

Signification et affirmations
L'article affirme que la convection double-diffusive, lorsqu'elle est traitée de manière auto-cohérente dans un contexte d'évolution globale, ne peut pas éroder de manière significative les gradients compositionnels générés lors de la formation des planètes. Les distributions d'éléments lourds observées dans Jupiter et Saturne nécessitent probablement des mécanismes de transport additionnels ou des voies de formation produisant un mélange initial plus large que ce que prédisent les modèles d'accrétion standard.

Les auteurs suggèrent que les « noyaux flous » peuvent être un produit direct du processus de formation lui-même (par exemple, un mélange précoce plus étendu lors de l'ablation des planétésimaux) ou le résultat de collisions majeures, plutôt que le produit d'un mélange double-diffusif post-formation. L'étude souligne une limite critique dans l'application des résultats hydrodynamiques locaux (souvent basés sur l'approximation de Boussinesq) à l'évolution globale des planètes géantes, en soulignant que le budget énergétique de l'intérieur profond est le principal régulateur de l'efficacité du mélange.