The Mass Function of Neutron Stars from Core-Collapse Supernova Simulations

Au cœur des étoiles massives, un drame final et violent se joue. Lorsqu'une étoile bien plus lourde que notre Soleil vient à manquer de carburant, son noyau s'effondre sous son propre poids, déclenchant une explosion de supernova capable d'éclipser une galaxie entière. Dans le sillage de cet événement, le noyau survit souvent sous la forme d'une étoile à neutrons, un objet si dense qu'une seule cuillère à café de sa matière pèserait l'équivalent d'une montagne. Pendant des décades, les astronomes ont tenté de comprendre la naissance de ces vestiges cosmiques. Ils ont mesuré les masses des étoiles à neutrons trouvées dans des systèmes binaires, où deux étoiles orbitent l'une autour de l'autre, mais les raisons sous-jacentes à la distribution spécifique de ces poids sont restées insaisissables. Existe-t-il un poids standard pour une étoile à neutrons nouveau-née, ou naissent-elles avec une grande variété de masses ? Et comment le processus chaotique de la mort d'une étoile détermine-t-il si elle devient une étoile à neutrons ou s'effondre complètement en un trou noir ?

Une nouvelle étude menée par une équipe de chercheurs réunit les simulations informatiques les plus avancées de ces explosions stellaires avec les dernières observations d'étoiles à neutrons pour répondre à ces questions. En reliant la structure interne de l'étoile mourante à la masse du vestige qu'elle laisse derrière elle, l'équipe a construit une image théorique de la façon dont les étoiles à neutrons naissent. Leurs travaux suggèrent que la variété des masses que nous observons aujourd'hui n'est pas seulement le résultat d'un gain de poids plus tard dans la vie des étoiles, mais qu'elle est inscrite dès l'instant même de leur création. Les conclusions offrent une explication cohérente sur la raison pour laquelle la plupart des étoiles à neutrons se regroupent autour d'un poids spécifique, tandis qu'un groupe plus restreint s'étire vers des masses beaucoup plus lourdes, et elles fournissent une nouvelle estimation de la fréquence à laquelle les étoiles massives ne parviennent pas à devenir des étoiles à neutrons, s'effondrant plutôt en trous noirs.

Les chercheurs ont commencé par examiner la relation complexe entre la taille initiale d'une étoile massive et la structure de son noyau juste avant qu'elle n'explose. Ils ont utilisé des modèles informatiques tridimensionnels sophistiqués qui simulent la physique complexe d'une supernova, suivant le comportement du noyau de l'étoile dans les derniers instants de sa vie. Ces simulations ont révélé un schéma clair : la « compacité » interne du noyau de l'étoile — une mesure de la densité de la matière à un point spécifique à l'intérieur — détermine directement la masse de l'étoile à neutrons qui se forme. Un noyau moins compact a tendance à laisser derrière lui une étoile à neutrons plus légère, tandis qu'un noyau plus compact produit une étoile plus lourde. Cette connexion a permis à l'équipe de prendre la distribution connue des étoiles massives dans l'univers et de prédire la distribution des masses des étoiles à neutrons au moment de leur naissance.

Lorsqu'ils ont comparé cette distribution de naissance théorique avec les masses réelles des étoiles à neutrons observées dans le ciel, une divergence est apparue. Les simulations prédisaient que les étoiles à neutrons nouveau-nées les plus communes devraient être légèrement plus légères que la masse moyenne observée, avec un groupe secondaire d'étoiles plus lourdes. Cependant, les données observées montraient un pic à un poids légèrement supérieur. Les chercheurs ont réalisé que cette différence pouvait s'expliquer par deux facteurs principaux : le gain de masse des étoiles à neutrons après leur naissance, et la manière dont leur naissance affecte leur maintien dans un système binaire.

Les étoiles à neutrons dans les systèmes binaires aspirent souvent de la matière de leurs étoiles compagnes, un processus appelé accrétion, ce qui ajoute du poids à l'étoile à neutrons au fil du temps. L'équipe a appliqué des corrections aux données observées pour estimer quelles auraient été les masses à la naissance, en retirant le poids gagné grâce à leurs compagnes. Cet ajustement a déplacé le pic observé vers le bas pour l'aligner beaucoup plus étroitement avec la prédiction théorique. De plus, l'étude a mis en évidence le rôle des « coups de pied » (kicks) reçus par les étoiles à neutrons à leur naissance. Tout comme une fusée est propulsée vers l'avant par ses gaz d'échappement, une étoile à neutrons reçoit une poussée soudaine dans la direction opposée à l'explosion. Les simulations suggèrent que les étoiles à neutrons plus lourdes reçoivent en moyenne des coups de pied plus puissants. Ces coups de pied puissants peuvent briser le lien gravitationnel qui maintient le système binaire ensemble, éjectant les étoiles à neutrons les plus lourdes dans le vide où elles sont plus difficiles à trouver. Ce retrait sélectif signifie que les systèmes binaires que nous observons sont biaisés en faveur des étoiles à neutrons plus légères qui ont reçu des coups de pied plus doux et sont restées en orbite.

En combinant les effets de la formation de trous noirs, de l'accrétion de masse et de ces coups de pied de naissance, l'équipe a créé un modèle raffiné qui correspond remarquablement bien aux données observées. Leur modèle suggère que l'univers produit un large spectre de masses d'étoiles à neutrons à la naissance, allant d'environ 1,2 à 2,0 fois la masse de notre Soleil. La masse de naissance la plus courante est d'environ 1,35 fois la masse du Soleil, ce qui explique le pic proéminent observé dans les observations. La queue plus lourde de la distribution, s'étendant jusqu'à 2,0 masses solaires, est un résultat naturel des structures variables des étoiles progénitrices, plutôt que du seul résultat d'une accrétion ultérieure.

L'étude a également fourni une nouvelle estimation de la fraction des étoiles massives qui terminent leur vie sous forme de trous noirs au lieu d'étoiles à neutrons. En analysant quels types d'étoiles dans leurs simulations n'ont pas réussi à produire une étoile à neutrons, les chercheurs ont calculé qu'environ 21 % des étoiles massives dans la plage de masse pertinente s'effondrent directement en trous noirs. Ce chiffre n'est pas une règle simple basée sur la taille initiale de l'étoile ; il dépend de la structure complexe et non linéaire du noyau de l'étoile, ce qui signifie que des étoiles de tailles similaires peuvent avoir des destins très différents.

Ce travail représente une étape importante pour relier la théorie de la mort des étoiles à la réalité de ce que nous observons. Il dépasse l'ancienne idée selon laquelle toutes les étoiles à neutrons naissent avec presque la même masse, en la remplaçant par une image dynamique où la structure initiale de l'étoile dicte une large gamme de résultats. L'accord entre les simulations complexes et les données observées suggère que notre compréhension des explosions de supernovae arrive à maturité. Bien qu'il subsiste des incertitudes, telles que l'influence précise de la composition chimique de l'étoile ou les détails des interactions binaires, l'étude démontre que la théorie moderne des supernovas peut désormais expliquer la population d'objets compacts que nous voyons dans le ciel. Elle ouvre une nouvelle ère où la naissance de ces objets extrêmes peut être retracée jusqu'aux conditions spécifiques des étoiles qui les ont créés, transformant une collection de mesures en un récit cohérent de mort et de renaissance stellaires.