Que révèlent les isotopes de calcium et d’étain sur la matière invisible des étoiles à neutrons ?

Une densité que le laboratoire ne saurait atteindre

Une étoile à neutrons concentre davantage de masse que le Soleil dans une sphère d’à peine quelques dizaines de kilomètres de diamètre. Cette densité échappe aux capacités de tout laboratoire. Les astronomes doivent donc la déduire des étoiles elles-mêmes. Le pulsar J0740+6620, le plus massif connu à ce jour, pèse désormais 2,08 ± 0,07 masses solaires dans cette lecture. Pour qu’une étoile de ce poids ne s’effondre pas en trou noir, la matière la plus dense à son cœur doit être suffisamment rigide pour supporter la charge. D’autres mesures, notamment quatre déterminations du rayon par NICER et un signal de fusion, empêchent cette rigidité de devenir arbitraire.

La calibration bayésienne transforme ces observations en distributions de probabilité plutôt qu’en une réponse unique sur la matière nucléaire. La masse du pulsar fixe une limite supérieure à la souplesse possible de la partie haute densité, tandis que les rayons et la déformation lors de la fusion servent de garde-fous quant au volume qu’une telle étoile peut occuper. Ensemble, ils resserrent l’ensemble des modèles acceptables sans le réduire à une seule courbe. Chaque modèle survivant conserve une incertitude quantifiée. Le résultat est donc un éventail de comportements possibles pour les étoiles denses, assorti de probabilités.

Cette forme probabiliste est essentielle car le même calcul doit rester cohérent avec la matière nucléaire ordinaire près de la densité de saturation. Une estimation libre à haute densité pourrait soutenir des étoiles massives tout en contredisant ce que nous savons des noyaux au laboratoire. Le protocole transmet donc l’information basse densité vers l’avant, afin que le cœur dense et les couches externes restent partie d’une même description physique. Ainsi, le pulsar le plus massif n’ajoute pas simplement un nombre : il ancre une région de matière autrement inaccessible à des preuves que nous pouvons effectivement mesurer.

Des noyaux de calcium et d’étain à coquille ouverte contraindront le modèle nucléaire

La calibration nucléaire probabiliste intègre désormais des noyaux de calcium et d’étain à coquille ouverte afin de tester sa capacité à traiter les particules non appariées. Ces noyaux sont déterminants : leurs masses et tailles mesurées réagissent à l’appariement ainsi qu’à la variation de l’énergie nucléaire avec la densité. Quatre isotopes du calcium (42, 44, 46 et 50) et trois de l’étain (112, 120 et 124) ont été ajoutés à la calibration. Leurs énergies de liaison et rayons de charge fournissent au modèle un ensemble plus large de points d’ancrage à basse densité que les noyaux à coquille fermée seuls. Un signal d’appariement des neutrons issu de l’étain-120 permet également de fixer la force d’appariement introduite dans le calcul. Ce terme est crucial car les noyaux à coquille ouverte peuvent présenter des neutrons non appariés, et leurs masses révèlent l’intensité avec laquelle ces derniers se lient. L’inclusion de ces isotopes a fait baisser la valeur de l’énergie de liaison à densité nucléaire ordinaire. Les paramètres de l’énergie d’asymétrie ont, eux, légèrement augmenté. Rappelons que cette énergie représente le coût supplémentaire lié à un excès de neutrons par rapport aux protons dans un noyau. Sa pente décrit la rapidité avec laquelle cette pénalité croît lorsque la matière se densifie. La pente déduite ici reste basse comparée à d’autres estimations. Toutefois, dès que les données des coquilles ouvertes entrent dans l’ajustement, les valeurs supérieures à 50 MeV ne sont plus fortement exclues. Le résultat est donc une distribution de probabilité sur les paramètres nucléaires possibles, et non une réponse unique pour la matière dense. Cette dispersion a son importance : les énergies de liaison tirent vers une région du paramètre tandis que les observables riches en neutrons tirent différemment depuis une autre. Les mesures du calcium et de l’étain resserrent ainsi certaines parties du modèle sans effacer la tension entre elles.

La croûte s’alourdit sans perdre de rigidité

L’inférence statistique mise à jour, fondée sur les mesures des noyaux atomiques, décale la plage de probabilité de la charge ionique dans la croûte vers le haut. Ce glissement traduit une valeur plus élevée du terme d’énergie de symétrie à basse densité, qui quantifie le coût énergétique d’un déséquilibre entre neutrons et protons. Les ions de la croûte adoptent ainsi des charges plus lourdes, non pas parce que le réseau cristallin est reconstruit, mais parce que la fraction de protons autorisée augmente. En revanche, les propriétés mécaniques globales de la croûte, notamment sa pression et son module de cisaillement, restent quasi inchangées.

La taille des cellules qui détermine l’agencement de ces ions dans le solide demeure proche de sa valeur antérieure. Densité centrale, densité du gaz et rayon ne présentent aucune variation significative lorsque les nouvelles données nucléaires sont intégrées. Comme les ions demeurent ancrés dans un réseau dont l’écartement est quasi identique, la résistance au cisaillement de la croûte n’augmente pas avec la charge plus lourde. Même les corrections de taille finie appliquées aux ions laissent ce module de cisaillement essentiellement inchangé.

L’effet est donc local, limité à la composition, et non global sur la structure. L’équation d’état et la relation masse-rayon restent proches des analyses antérieures des étoiles à neutrons et demeurent compatibles avec les mesures de rayon de NICER. Cette dissociation est cruciale : une croûte peut porter des ions différents sans que la taille totale ou le profil de pression de l’étoile ne soient modifiés.

Une énergie de symétrie plus souple à basse densité, mais un comportement plus rigide aux hautes densités

La calibration bayésienne d’un modèle de Skyrme aboutit à une distribution de probabilité pour la pente de l’énergie de symétrie, qui reste relativement faible aux abords de la densité nucléaire ordinaire. Ce type de modèle constitue une recette mathématique décrivant l’énergie du noyau. Ici, l’ajustement s’appuie sur des données variées : les noyaux de calcium et d’étain à coquilles ouvertes, d’autres grandeurs mesurables sur des noyaux finis, des calculs directs de la matière neutronique pure, les mesures de rayons par NICER ainsi que les contraintes issues des ondes gravitationnelles. L’intérêt majeur réside dans le fait que les corrélations entre paramètres sont conservées d’une étape à l’autre, au lieu d’être réinitialisées. Un résumé statistique compact permet ainsi aux autres équipes de reprendre directement ces incertitudes liées, sans avoir à supposer que chaque paramètre nucléaire varie indépendamment des autres.

Cette inférence transforme un problème complexe à nombreux paramètres en un point de départ commun. Les isotopes de calcium et d’étain à coquilles ouvertes apportent des énergies de liaison, des rayons de charge et une contrainte sur le gap d’appariement. Le résultat ne repose donc pas sur une seule famille nucléaire ni sur une masse mal décrite. L’étape suivante relie ces évolutions nucléaires à la croûte et au cœur de l’étoile, tout en maintenant la courbe globale masse-rayon compatible avec les mesures de rayons de NICER.

La limite du procédé est que les données des noyaux finis contraignent le modèle très fortement près de la densité nucléaire ordinaire, tandis que les densités du cœur sont fournies par une extension flexible à haute densité. Dans ce cadre, le comportement autorisé à haute densité tend à se rigidifier au-delà de la densité nucléaire ordinaire si l’on veut rester compatible avec les mesures de rayons de NICER et avec une masse de pulsar d’environ 2,08 masses solaires utilisée comme contrainte de masse maximale. Les contraintes marées issues des ondes gravitationnelles limitent ensuite le degré de rigidité que cette région à haute densité peut atteindre. Le résumé gaussien compact omet par ailleurs l’asymétrie et la multimodalité, notamment pour un paramètre de surface dont la distribution de probabilité est large et irrégulière. Les équipes qui utilisent ce résumé doivent donc garder à l’esprit qu’il ne représente pas la forme complète de l’incertitude.

Comme l’inférence reste probabiliste, les futures mesures de rayons et de marées testeront jusqu’à quel point la matière nucléaire à haute densité doit se rigidifier tout en maintenant une énergie de symétrie relativement souple à basse densité.