Leçon
La théorie — Probabilité de planètes semblables à la Terre
Il s’agit d’un effectif attendu, et non d’une probabilité : on part d’un nombre d’étoiles et on le multiplie par une chaîne de fractions, dont chacune répond à la question “parmi celles-ci, quelle proportion également… ?”. Le résultat correspond au nombre d’objets auxquels on s’attendrait si chaque fraction était exacte — et il s’agit d’un calcul de la même forme que l’équation de Drake, avec la même force et la même faiblesse.
Six barres empilées, chacune deux fois moins large que celle du dessus, des étoiles observées jusqu’aux détections attendues.
- Les étoiles observées — la barre entière, et le seul nombre ici qui soit compté plutôt que supposé.
- Chaque fraction prend une part de la barre au-dessus d’elle, jamais de la barre du haut. C’est ce qui rend la multiplication légitime, et c’est l’étape que l’on saute.
- Après quatre fractions la barre vaut un seizième du départ :
156256objets attendus. - La détection s’applique après la vérité, comme une étape distincte —
78128. Ce qui existe et ce qu’un relevé verrait sont deux barres différentes.
Ce que signifie chaque symbole
N*- le nombre d’étoiles observées,
150,000avec les valeurs par défaut. fractions- la chaîne —
f_pavec des planètes,f_HZdans la zone habitable,f⊕de la taille de la Terre,f_rockyrocheuses. Chacune est une proportion conditionnelle de la ligne précédente, et non du total. E- la population attendue — combien existent réellement :
3,029. p_det- l’efficacité de détection : la proportion que l’on détecterait réellement. Elle transforme une réalité en une observation.
E_det- les détections attendues,
E × p_det = 242.
D’où vient la formule
- La multiplication n’est la bonne opération que parce que chaque fraction est conditionnée par la précédente.
f_HZest la proportion d’étoiles abritant des planètes qui possèdent une planète en zone habitable, et non la proportion de toutes les étoiles. - En les enchaînant :
E = N* · f_p · f_HZ · f⊕ · f_rocky. Avec les valeurs par défaut, cela donne150,000 × 0.85 × 0.18 × 0.22 × 0.6 = 3,029. - La détection est une étape distincte, appliquée après la réalité sous-jacente :
E_det = E · p_det = 3,029 × 0.08 = 242. Tout l’intérêt de les séparer est de distinguer ce qui existe de ce qu’un relevé observerait.
- Suppose
- Que chaque fraction est indépendante des autres sachant l’étape précédente, et qu’une valeur unique peut représenter chacune d’elles. Aucune de ces hypothèses n’est sûre — l’occurrence des planètes est corrélée au type spectral et à la métallicité de l’étoile, si bien que ces fractions ne peuvent pas vraiment être réglées indépendamment les unes des autres.
- Ne tient plus quand
- Considérez la précision avec scepticisme. Le résultat est donné à l’unité près —
242— à partir de cinq fractions réglées chacune par un curseur à une décimale, et multiplier cinq estimations multiplie leurs erreurs au lieu de les lisser. Divisez deux des fractions par deux et le résultat est divisé par quatre. La fourchette affichée est étroite parce qu’elle explore la taille du relevé et non les fractions, alors que c’est dans ces dernières que réside la quasi-totalité de l’incertitude.
Exercice
Vérifiez-vous
Prédisez d’abord la réponse, puis utilisez les commandes ci-dessus pour vérifier. N’affichez la solution qu’après vous être engagé sur une hypothèse : c’est ce qui en fait un exercice.
-
Mettez toutes les fractions à 0.5 et notez les deux comptes. Montez maintenant une seule de fp, fHZ, fEarth ou frocky à 0.6 — puis remettez-la et montez pdetect à 0.6 à la place. Qu’y a-t-il de différent avec la cinquième ?
Afficher la réponse
Les quatre fractions de population sont interchangeables. Quelle que soit celle que vous montez, le nombre détecté passe de 78,128 à 93,754 et la vérité de population de 156,256 à 187,508, car la multiplication se moque de l’ordre : une amélioration de 20% vaut 20% où que vous l’appliquiez. Monter pdetect donne les mêmes 93,754 détections mais laisse la vérité à 156,256 — l’efficacité de détection est un fait sur votre télescope, pas sur la galaxie. C’est exactement à quoi servent les deux lignes : quatre facteurs disent combien de planètes existent, et le cinquième combien vous en verriez. C’est aussi pourquoi ce genre d’argument tourne toujours autour du facteur le plus *incertain* et non du plus petit. Le levier est identique ; seules les barres d’erreur les distinguent. -
La fourchette affiche 155,481 – 157,031 autour d’un centre à 156,256 — soit plus ou moins 775. D’où vient ce 775, et qu’est-ce qu’il omet ?
Afficher la réponse
√156,256 = 395.3, et 1.96 × 395.3 = 775 : un intervalle à 95% pour le bruit de comptage poissonien. C’est la dispersion que vous verriez si les cinq fractions étaient connues exactement et que seul le hasard décidait combien de planètes tombent dans votre échantillon. Il omet toute incertitude sur les fractions elles-mêmes — et celle-ci domine complètement. fHZ et fEarth se discutent dans la littérature à des facteurs deux près, pas à un demi pour cent près. Faites passer fHZ de 0.5 à 0.6 et l’estimation bondit de 31,252, environ quarante fois la largeur de la fourchette que l’outil venait d’afficher. La fourchette est honnête sur la seule chose qu’elle mesure, et follement optimiste sur le nombre dans son ensemble.
Problème entièrement résolu
-
Planètes de type terrestre à détecter dans un relevé de 150 000 étoiles 6 étapes
Observez 150 000 étoiles. En partant des fractions retenues par cet outil, calculez combien de planètes similaires à la Terre existent et combien vous en détecteriez réellement — puis déterminez dans quelle mesure on peut croire à la bande d'incertitude affichée.
-
Le modèle est un simple produit : on part du nombre d'étoiles et on le multiplie successivement par chaque fraction. Rien y concernant n'est conditionné par un autre facteur, ce qui constitue en soi une hypothèse.
-
Parcourez le calcul de gauche à droite. Chaque facteur réduit l'effectif restant, et les trois derniers éliminent à eux seuls 96 % des étoiles qui possédaient des planètes.
-
La détection constitue un sixième facteur, appliqué à la population de planètes plutôt qu'aux étoiles. Environ une planète similaire à la Terre sur douze est détectée.
-
Passons à la bande d'incertitude. Traiter le décompte selon une loi de Poisson attribue comme écart-type la racine carrée de la moyenne, et l'intervalle affiché par l'outil correspond à deux écarts-types de part et d'autre.
-
Comparez cette efficacité de détection avec la seule géométrie. Un transit exige que l'orbite soit vue par la tranche à R☉/1 au près, soit une chance sur 215 — ainsi, 0,08 n'est pas une probabilité géométrique, mais toute une stratégie d'observation condensée en un seul chiffre, et le résultat purement géométrique des transits donnerait 14 planètes au lieu de 242.
-
Enfin, faites varier une seule fraction dans un intervalle qu'aucun astronome ne jugerait irraisonnable, et observez la réponse varier d'un facteur quatre.
Réponse
3 029 présents, 242 détectés, et la marge de ± 108 ne vaut presque rien. Cette marge correspond au bruit de comptage poissonien, c’est-à-dire à la dispersion que vous observeriez si les cinq fractions étaient parfaitement exactes et que seul le hasard des dés intervenait. Or ces fractions comportent elles-mêmes des incertitudes. Il suffit que l’une d’elles, la fraction de la zone habitable, puisse être divisée ou multipliée par 2 pour que l’estimation s’étende de 1 515 à 6 059, soit une plage vingt fois plus large que celle affichée. C’est le défaut classique de tout modèle fondé sur un produit de fractions, équation de Drake comprise : le calcul produit une précision à quatre chiffres à partir de données qui n’en justifient qu’un seul. Cette précision n’est qu’un artefact de la multiplication.
-
Références (1)
- The measured occurrence rate of rocky habitable-zone planets, the one factor here that is not a guess: S. Bryson, M. Kunimoto, R. K. Kopparapu, J. L. Coughlin, W. J. Borucki & D. Koch, "The Occurrence of Rocky Habitable-zone Planets around Solar-like Stars from Kepler Data." The Astronomical Journal 161(1), 36, 2020.