An Increase in the Galactic Planet Host Fraction Fails to Reproduce the Galactic Height Trend in Planet Occurrence

Résumé Technique : Une augmentation de la fraction d'étoiles hôtes galactiques ne parvient pas à reproduire la tendance de la hauteur galactique dans l'occurrence des planètes

Énoncé du problème
Des résultats observationnels récents, spécifiquement Zink et al. (2023) (ci-après Z23), ont identifié une forte tendance où le taux d'occurrence des petites planètes (Super-Terres et Sous-Neptunes) autour des naines FGK diminue significativement avec l'augmentation de l'amplitude d'oscillation verticale maximale ($Z_{max}$) par rapport au plan moyen galactique. Bien que la métallicité stellaire soit connue pour être corrélée à la formation planétaire, le seul gradient de métallicité galactique est insuffisant pour expliquer la raideur de ce déclin observé. La question centrale abordée par cette étude est de savoir si une augmentation dépendante du temps de la fraction d'étoiles hôtes de planètes ($f$) au cours de l'histoire de la Voie Lactée pourrait reproduire la tendance occurrence-$Z_{max}$ observée. Les auteurs examinent si un « boost » de l'efficacité de la formation planétaire à un certain moment du passé de la Galaxie, combiné au chauffage cinétique des orbites stellaires au fil du temps, peut rendre compte du déficit de planètes aux hauteurs galactiques élevées.

Méthodologie
Les auteurs ont développé un nouveau code de synthèse de population de systèmes planétaires, psps, pour simuler des populations planétaires synthétiques autour d'un échantillon de naines FGK. La méthodologie comprenait les étapes suivantes :

1. Échantillons Stellaires : Deux échantillons stellaires distincts ont été construits pour tester la robustesse des résultats face aux incertitudes sur l'âge :

   • Échantillon Réel : Un échantillon conjoint Kepler et K2 (22 184 étoiles) tiré du catalogue de Hardegree-Ullman et al. (2025), mis en correspondance avec les âges d'isochrones de Berger et al. (2020b, 2025). Les valeurs de $Z_{max}$ ont été calculées à l'aide de la cinématique Gaia DR3 et d'un modèle de potentiel galactique via le package Gala.
   • Échantillon Synthétique : Un échantillon de synthèse de population TRILEGAL (49 109 étoiles) ajusté pour ressembler aux champs Kepler et K2. Cet échantillon a fourni un scénario « optimiste » avec des âges stellaires arbitrairement précis et des modèles de chauffage cinétique connus, permettant aux auteurs d'isoler les effets de la précision de l'âge des modèles physiques.

2. Modèles d'Occurrence Dépendants du Temps : L'étude a testé deux familles de modèles pour l'évolution de la fraction d'étoiles hôtes de planètes ($f$) en fonction de l'âge stellaire (temps de regard en arrière) :

   • Fonctions Échelonnées (Step Functions) : Un saut soudain d'une fraction initiale $f_1$ vers une fraction actuelle $f_2$ à un temps de regard en arrière $t$ spécifique.
   • Fonctions par Morceaux (Piecewise Functions) : Une augmentation linéaire graduelle de $f_1$ vers $f_2$ commençant au temps $t$.
   • Tous les modèles ont été contraints pour donner une fraction actuelle d'étoiles hôtes de planètes de $\sim 30-35\,\%$, cohérente avec les contraintes observationnelles.

3. Synthèse de Population et Observation : En utilisant psps, les auteurs ont généré des systèmes planétaires basés sur la fraction $f$ prescrite pour chaque étoile. Les systèmes ont reçu des architectures (systèmes multi-planétaires compacts vs planètes uniques) et des paramètres orbitaux (excentricité, inclinaison) basés sur des contraintes démographiques existantes (ex: Lam & Ballard 2024). Les populations synthétiques ont ensuite été « observées » en utilisant des méthodes d'efficacité de détection inverse (IDEM) qui tenaient compte de la probabilité de transit géométrique et des fonctions de sensibilité spécifiques des missions Kepler et K2.

4. Analyse : Les taux d'occurrence simulés ont été regroupés par classes de $Z_{max}$ (100–1000 pc) et ajustés avec une loi de puissance ($y \propto Z_{max}^\tau$) pour déterminer la pente ($\tau$) et la normalisation ($\eta$). Ces résultats ont été comparés à la tendance observée de Z23 ($\tau \approx -0,28$).

Résultats Clés

• Insuffisance des Augmentations de la Fraction d'Hôtes : L'étude conclut qu'augmenter la fraction d'étoiles hôtes de planètes ($f$) au fil du temps, que ce soit de manière instantanée ou graduelle, est généralement insuffisant pour reproduire la raideur de la tendance observée occurrence-$Z_{max}$.
• Contraintes des Modèles :
  • Pour l'échantillon réel Kepler-K2, la grande majorité des modèles produisaient des pentes ($\tau$) trop faibles (proches de 0) par rapport au $\tau \approx -0,28$ observé. Les modèles qui se rapprochaient le plus de la tendance nécessitaient des paramètres extrêmes et physiquement implausibles (ex: $f_1 \approx 1\,\%$ ou $f_2 \approx 100\,\%$) et des temps de saut ($t$) restreints à une fenêtre étroite entre 2,2 et 3,2 milliards d'années.
  • Pour l'échantillon synthétique TRILEGAL (avec des âges précis), les résultats étaient similaires. Bien qu'une augmentation tardive et instantanée de $f$ puisse correspondre à la pente de Z23, elle nécessitait à nouveau un ajustement vers des extrêmes théoriques ($f_1=1\,\%$ ou $f_2=100\,\%$). Cela suggère que même avec une précision d'âge parfaite, une simple augmentation de $f$ ne peut pas, à elle seule, piloter la tendance observée.
• Analyse Kepler-Uniquement : Lors de l'analyse du seul champ Kepler (qui possède une distribution d'âge différente et manque de données K2), la pente occurrence-$Z_{max}$ dérivée était similaire à celle de l'échantillon conjoint. Une augmentation par saut de $f$ correspondait légèrement mieux à cette tendance Kepler-uniquement, mais la difficulté fondamentale de reproduire la raideur persistait.
• Rôle de la Précision de l'Âge : La comparaison entre l'échantillon réel et l'échantillon synthétique indique que l'incapacité des modèles à correspondre à la tendance n'est pas principalement due à une précision insuffisante des mesures d'âge stellaire. Même avec des âges idéalisés, les modèles d'« augmentation graduelle » ont été écartés, et les modèles de « saut » nécessitaient des paramètres irréalistes.

Signification et Revendications
L'article conclut que la forte corrélation négative observée entre l'occurrence des petites planètes et la hauteur galactique ne peut être expliquée uniquement par une augmentation dépendante du temps de la fraction d'étoiles hébergeant des planètes. Les auteurs affirment que, bien qu'un boost tardif de $f$ (spécifiquement une fonction en escalier survenant il y a 2 à 3 milliards d'années) puisse marginalement reproduire la tendance, cela nécessite un ajustement vers les limites physiques du modèle (ex: presque 100 % des étoiles hébergeant des planètes aujourd'hui).

Les auteurs soutiennent que cette divergence implique que d'autres facteurs doivent être en jeu. Ils suggèrent que les futurs modèles pourraient obtenir un meilleur ajustement en :

1. Distinguant les Super-Terres des Sous-Neptunes (car Z23 a trouvé que la tendance est principalement pilotée par les Sous-Neptunes).
2. Permettant à la fraction de systèmes multi-planétaires compacts « intacts » d'évoluer au fil du temps, plutôt que de la maintenir constante.
3. Investiguant des mécanismes physiques qui soit augmentent la formation planétaire récente, soit perturbent progressivement les systèmes plus anciens (ex: via des passages stellaires, des marées galactiques ou des changements dans la turbulence du milieu interstellaire).

L'étude souligne que l'image canonique de l'enrichissement métallique local comme seul lien entre le contexte galactique et la démographie planétaire est insuffisante à l'ère de Gaia. Les résultats mettent en évidence la nécessité de connecter l'histoire dynamique à l'échelle galactique (fusions, chutes de satellites) avec les processus locaux affectant les disques protoplanétaires et la stabilité des systèmes pour expliquer pleinement la démographie des exoplanètes.