Les sub-Neptunes sont la classe d'exoplanètes la plus répandue dans la Galaxie, mais leur structure interne reste un sujet de débat. On considérait traditionnellement que sous l'atmosphère d'hydrogène se trouvait une frontière nette avec le magma silicaté. Cependant, des calculs récents ont montré qu'aux pressions et températures caractéristiques des intérieurs des sub-Neptunes, le silicate et l'hydrogène peuvent se mélanger complètement. Cela pose le problème de la description des propriétés thermodynamiques d'un tel fluide supercritique, sans quoi il est impossible de construire des modèles corrects de la structure interne des planètes. Ce travail s'inscrit dans le cadre des recherches sur les équilibres de phases et les équations d'état pour les systèmes silicate-hydrogène, activement développées ces dernières années.
Les auteurs ont effectué des calculs de dynamique moléculaire ab initio dans le cadre de la théorie de la fonctionnelle de la densité, en utilisant le logiciel VASP et la fonctionnelle d'échange-corrélation PBEsol. Une cellule de 243 atomes (27 Mg, 27 Si, 81 O, 108 H) a été simulée, correspondant à la composition MgSiO3H4. D'abord, le MgSiO3 cristallin a été fondu, puis l'hydrogène a été ajouté. Les calculs ont été menés à des températures de 4000, 6000 et 8000 K et à des pressions allant jusqu'à ~600 GPa. Les données obtenues sur la pression et l'énergie interne ont été ajustées avec une énergie libre de Helmholtz, ce qui a permis de dériver la densité, l'entropie, la capacité thermique, la dilatation thermique, le module de compressibilité et le paramètre de Grüneisen à partir d'une surface thermodynamique unique.
Le fluide supercritique MgSiO3H4 s'est avéré moins dense que le silicate fondu sec, dans toutes les conditions étudiées. L'écart par rapport au mélange idéal des volumes de MgSiO3 et H2 atteint plusieurs pour cent et augmente avec la pression. L'analyse structurelle a montré que, avec l'augmentation de la pression, le nombre de coordination Si–H augmente, tandis que le nombre de liaisons H–H stables chute brusquement : à 4000 K, la fraction d'atomes d'hydrogène participant à des contacts H–H de longue durée tombe en dessous de 0,1 à des pressions d'environ 340 GPa. Cela indique une transition de l'hydrogène moléculaire vers une forme plus atomique, intégrée dans le réseau silicaté. L'application de l'équation d'état à des modèles de sub-Neptunes de masse 6 masses terrestres avec un stock total d'hydrogène identique de 4 % a montré que, selon la pression à la frontière atmosphère-intérieur (1 GPa contre 5 GPa), le rayon de la planète varie de 5,7 %.
Ce travail fournit pour la première fois une équation d'état auto-cohérente pour un fluide supercritique silicate-hydrogène, comblant ainsi une lacune importante dans la modélisation des sub-Neptunes. Contrairement aux travaux voisins qui évaluaient principalement les équilibres de phases ou le mélange au niveau de modèles thermodynamiques, ici des données ab initio directes sur la densité et les dérivées thermodynamiques ont été obtenues. Cela permet de prendre en compte quantitativement l'hydrogène dissous dans l'intérieur et son influence sur le rayon et l'évolution thermique des planètes.
La suite naturelle est l'extension de la composition avec le fer, qui peut fortement influencer la densité, les conditions d'oxydoréduction et les transitions de phase. Des vérifications expérimentales à hautes pressions et températures sont également nécessaires, mais elles sont encore très rares. Les efforts conjoints des théoriciens et des expérimentateurs permettront de construire des équations d'état pour toute la gamme de compositions, pressions et températures pertinentes pour les sub-Neptunes.
Les résultats influenceront directement les modèles de structure interne et d'évolution des sub-Neptunes, ainsi que l'interprétation des données sur les masses et les rayons des exoplanètes.
Inclure le fer dans les calculs et effectuer des mesures expérimentales de densité et de compressibilité des fluides supercritiques silicate-hydrogène.
Ce travail aborde un problème fondamental de la planétologie : comment déterminer la composition et la structure interne des exoplanètes à partir de données observationnelles limitées. Il montre qu'ignorer la non-idéalité du mélange silicate-hydrogène peut conduire à des erreurs systématiques dans l'estimation de la teneur en hydrogène et de l'état thermique des intérieurs.
Cette section s'adresse à l'auteur du travail et à ses collègues. Les travaux voisins ont été trouvés par une recherche sémantique dans notre archive, et non par correspondance de mots. Nous n'évaluons pas : nous montrons le champ autour.
Ce que fait ce travail
Les auteurs ont réalisé des calculs de dynamique moléculaire ab initio du mélange supercritique MgSiO3H4 à des pressions allant jusqu'à 600 GPa et des températures de 4000 à 8000 K. Ils ont obtenu l'équation d'état et montré que la densité du mélange est inférieure à celle du silicate sec, et que le comportement s'écarte du mélange idéal. En appliquant cette équation d'état à des modèles de sub-Neptunes de 6 masses terrestres, ils ont découvert que la distribution de l'hydrogène entre l'atmosphère et l'intérieur influence le rayon de la planète (une différence de 5,7 % pour des pressions limites de 1 et 5 GPa).
Où est sa force
La force de ce travail réside dans l'obtention d'une équation d'état auto-cohérente du fluide supercritique silicate-hydrogène à partir des premiers principes, ce qui comble une lacune dans la modélisation des sub-Neptunes et permet de prendre en compte la distribution réelle de l'hydrogène dans les profondeurs.
Vers où avancer
À essayer
Pourquoi cela peut compter
La prise en compte du mélange supercritique du silicate et de l'hydrogène pourrait modifier les estimations du rayon des sub-Neptunes de plusieurs pour cent, ce qui est important pour interpréter la relation observée rayon-masse et pour déterminer la structure interne des exoplanètes.
Travaux voisins dans notre archive
Le travail le plus proche de l'archive se situe à 0.65 par le sens — le niveau habituel de notre corpus.
À une pression d'environ 340 GPa et une température de 4000 K, seul un atome d'hydrogène sur dix conserve une liaison avec un autre hydrogène pendant plus de 100 femtosecondes ; les autres sont déjà intégrés dans le réseau silicaté.