Cent mille étoiles dans le noyau dense d’un amas globulaire se déplacent comme des billes d’acier dans un flipper géant. Chaque « choc » — une rencontre gravitationnelle entre un couple binaire et une étoile isolée — déclenche une cascade d’événements imprévisibles. Le duo peut se disloquer, l’intrus peut éjecter l’un des partenaires et prendre sa place, et parfois les trois se figent dans une danse hiérarchique complexe. Ces rendez-vous durent quelques jours, mais leurs conséquences scellent le destin de millions d’astres et donnent naissance à des objets exotiques — des trou noir aux sources d’ondes gravitationnelles.

Le classique problème à trois corps est chaotique : un infime écart dans la vitesse initiale, et le dénouement change du tout au tout. Henri Poincaré fut le premier à le comprendre dès la fin du XIXe siècle, en démontrant qu’il n’existe pas de solution analytique générale. Mais un groupe d’astrophysiciens a prouvé que le chaos n’est pas total. En réalisant près d’un million de simulations numériques avec le code REBOUND, ils ont entraîné un modèle d’apprentissage automatique XGBoost à deviner l’issue. Les données d’entrée étaient des grandeurs physiques sans dimension — la rigidité du couple binaire, le coefficient de focalisation, le nombre de Safronov. La rigidité, rapport entre l’énergie de liaison du couple et l’énergie cinétique de l’étoile incidente, s’est avérée être le chef d’orchestre : si le couple est « dur », il survivra presque à coup sûr au choc ou procédera à un échange de partenaires.

Les erreurs ne sont pas réparties uniformément. Elles se concentrent sur une frontière étroite entre l’échange et les systèmes triples hiérarchiques. Là, la redistribution d’énergie et de moment cinétique est si sensible que les captures temporaires d’étoiles provoquent une résonance chaotique. L’issue dépend de fluctuations microscopiques — comme une bille qui hésiterait sur la crête entre deux cuvettes. L’analyse de l’entropie des bassins a confirmé que cette frontière possède une structure fractale. C’est précisément là que le modèle vacille et affiche une faible confiance. Non un défaut, mais le reflet de l’imprévisibilité fondamentale découverte par Poincaré.

Le gain pratique est colossal. Avant, simuler une seule de ces rencontres prenait des minutes de temps de calcul. Désormais, le substitut s’exécute en millisecondes. Une accélération d’un facteur cent mille permet d’intégrer l’algorithme dans des simulations globales d’amas, où chaque événement était auparavant remplacé par une approximation grossière. Désormais, la dynamique des résonances et la formation des systèmes triples seront reproduites avec une fidélité physique. La prochaine étape consistera à ajouter les forces de marée, les tailles finies des étoiles et les corrections relativistes pour prédire les fusions d’objets compacts avec un niveau de détail inédit. Cela influencera directement les études de population des sources d’ondes gravitationnelles.

Le gigantesque flipper stellaire trouve ses règles. Le chaos, qui semblait absolu, obéit à un ordre caché, accessible à un algorithme entraîné. Là où les mathématiques capitulent, la physique laisse des indices — et le modèle apprend à les déchiffrer. Ainsi, le fracas des boules qui s’entrechoquent dans l’obscurité des amas ne sera plus un bruit aléatoire.