Oubliez les images habituelles des noyaux planétaires — la frontière nette entre roche et gaz est révolue. La sonde Juno, en plongeant dans le champ gravitationnel de Jupiter, et Cassini, dans les bras de Saturne, ont découvert que les éléments lourds ne sont pas concentrés au centre mais étalés sur une vaste région, comme si quelqu'un avait mal remué un cocktail. Pendant longtemps, on a cru que cet étalement était l'œuvre de la convection à double diffusion. Ce processus, observé dans l'océan où des couches d'eau de salinité et de température différentes forment une structure en escalier, semblait idéal pour mélanger les intérieurs. Mais les calculs disent le contraire.
Les auteurs se sont armés de la recette des échelles convectives de Spruit (2013) et l'ont intégrée dans le code d'évolution APPLE, qui modélise toute la vie de la planète — de l'enfance en fusion jusqu'à aujourd'hui, 4,56 milliards d'années. Les conditions initiales incluaient déjà des noyaux dilués : les planétésimaux, en bombardant la protoplanète, s'évaporaient partiellement, laissant derrière eux de longs panaches de eau et de carbone, sortes de fantômes de comètes en chute. Ensuite, le code calculait comment ces couches évoluaient sous l'effet de la convection ordinaire — celle qui obéit au critère de Schwarzschild —, de la diffusion moléculaire et, bien sûr, de la double diffusion. Résultat ? Jupiter n'a remué que 0,11 masse terrestre, Saturne encore moins. Même si un expérimentateur de laboratoire amplifiait mentalement les coefficients de transport d'un facteur mille, les chiffres atteignent à peine deux masses — et ce, sur toute l'histoire du système solaire.
Cette inefficacité a une justification physique rigoureuse. Le paramètre de stabilité Rρ = (αμ/αT) (∇μ/(∇ – ∇ad)) détermine si des couches se forment. Pour Rρ à peine supérieur à 1, la stratification s'active, mais dès que la double diffusion tente d'augmenter le transfert de chaleur, le gradient superadiabatique chute et Rρ diminue — le système s'auto-atténue, comme le fait de remuer des braises refroidissantes avec un tisonnier éteint la chaleur. La deuxième équation, $E_{mix}$/(ρ₀ $c_P$ ΔT H) ~ H/$H_T$, met carrément un terme aux ambitions : elle montre que dans un intérieur compressible, l'épaisseur de la couche convective H est comparable à l'échelle de hauteur thermique $H_T$, ce qui signifie que le coût du mélange frôle le budget thermique total. Les planètes géantes, composées principalement de hydrogène et d'hélium, poussent ce principe à l'extrême, et la production d'entropie due au mélange est négligeable. Chandrasekhar, dans ses travaux sur les étoiles, avertissait : dans un fluide compressible, même de petites hétérogénéités de composition tiennent bon.
Qu'est-ce que cela signifie pour notre compréhension des planètes ? Les noyaux dilués ne sont pas le produit d'une lente combustion, mais soit une caractéristique innée, imprimée à l'époque de l'accrétion, soit une cicatrice d'un impact colossal — par exemple, une collision frontale avec un corps d'une dizaine de masses terrestres. Les modèles qui reposaient sur un mélange post-formation efficace sont relégués aux oubliettes. Nos conceptions des exoplanètes lointaines devront également être révisées : si les couches sont si stables, la spectroscopie de leurs atmosphères lira non pas des réactions chimiques actuelles, mais une recette ancienne, conservée depuis la naissance. Dans les entrailles des planètes, comme dans les fouilles archéologiques, chaque couche est une page de l'histoire.
À l'avenir, des simulations hydrodynamiques sans compromis sur l'approximation de Boussinesq pourront capturer toute la rigidité de la convection compressible. Des expériences de laboratoire avec stratification multicouche reproduisent pas à pas des conditions où chaque centimètre de hauteur dévore des joules. La prochaine génération de modèles théoriques de transport turbulent pourrait apprendre à prédire comment l'évolution thermique distribue les éléments non seulement dans Jupiter et Saturne, mais aussi dans les naines brunes et les mondes gonflés autour d'étoiles lointaines. Une chose est sûre : le mixeur cosmique s'est avéré désespérément faible, et chaque corps céleste conserve son histoire dans un emballage stratifié intact.
L'idée de l'« échelle » à double diffusion vient de l'océanographie : là-bas, des profils en escalier de température et de salinité sont observés dans les mers et les océans. Dans les entrailles de Jupiter, de telles couches, si elles existaient, pourraient s'étendre sur des centaines de kilomètres. Cependant, l'énergie nécessaire pour les maintenir serait prodigieuse — et c'est précisément ce qui les tue.