Maximal mass of neutron stars constrained by neutron star observations

En employant un cadre de pondération bayésien qui intègre des observations multimessagers (incluant GW170817, NICER et des objets compacts candidats) avec des équations d'état hybrides, cette étude détermine que la masse maximale des étoiles à neutrons est robustement contrainte à environ 2,2–2,3 masses solaires, tandis que le rayon correspondant dépend plus fortement du modèle hadronique sous-jacent, tombant typiquement près de 12 km.

L'essentiel

Imaginez que l'univers soit rempli d'un matériau mystérieux et super-dense appelé « matière d'étoile à neutrons ». Il est si lourd qu'une seule cuillère à café pèserait autant qu'une montagne. Depuis longtemps, les physiciens tentent de déterminer les « règles du jeu » de ce matériau — spécifiquement, sa rigidité ou sa compressibilité. Cet ensemble de règles est appelé l'Équation d'État (EOS).

La grande question que cet article tente de répondre est la suivante : Quelle est la masse absolue maximale qu'une étoile à neutrons peut atteindre avant de s'effondrer en un trou noir ?

Voici l'histoire de la manière dont les auteurs ont résolu ce puzzle, expliquée en termes simples :

1. Les Deux Points de Départ (Les Recettes)

Pour déterminer les règles, les scientifiques ont commencé avec deux « recettes » différentes décrivant le comportement de cette matière dense à des densités plus faibles. Imaginez-les comme deux théories différentes sur la façon dont les ingrédients se mélangent :

• Recette A (SFHo) : Une recette « plus molle », ce qui signifie que la matière est un peu plus facile à comprimer.
• Recette B (DD2) : Une recette « plus rigide », ce qui signifie que la matière résiste davantage à la compression.

Ils savaient que ces recettes fonctionnaient bien au « début » de l'échelle de densité, mais ils ignoraient ce qui se produisait aux densités extrêmes, ultra-élevées, trouvées au centre d'une étoile à neutrons. Pour combler ce vide, ils ont utilisé un « pont » mathématique pour relier leurs recettes à ce que nous savons de la physique des particules aux énergies les plus élevées possibles.

2. Le Travail d'Enquête (Utilisation de Vraies Indices)

Au lieu de simplement deviner, les auteurs ont agi comme des détectives. Ils ont pris leurs deux recettes et les ont testées contre de véritables indices recueillis par des télescopes et des détecteurs d'ondes gravitationnelles. Ils ont utilisé une méthode statistique spéciale (appelée pondération bayésienne) pour voir quelles versions de leurs recettes résistaient à l'épreuve.

Voici les indices qu'ils ont utilisés :

• Le « Grand Crash » (GW170817) : Lorsque deux étoiles à neutrons se sont percutées, elles ont envoyé des ondulations dans l'espace. La manière dont ces ondulations se sont comportées a révélé aux scientifiques à quel point les étoiles étaient « compressibles ».
• La « Lampe de Poche » (NICER) : Un télescope spatial a pris des photos de points chauds sur des étoiles à neutrons en rotation. En mesurant la taille apparente des étoiles et leur masse, ils ont obtenu un rapport direct taille-poids.
• Le Candidat « Léger » (HESS J1731–347) : Un objet très petit et léger qui pourrait être une étoile à neutrons.
• Le Candidat « Poids Lourd » (GW190814) : Un objet mystérieux plus lourd que la plupart des étoiles à neutrons mais plus léger que la plupart des trous noirs. Les scientifiques se sont demandé : Cet objet pourrait-il en réalité être une étoile à neutrons ultra-lourde ?

3. Les Résultats : Ce que les Indices Leur Ont Dit

Les scientifiques ont soumis leurs deux recettes à ces indices et ont examiné les résultats.

La Limite de Poids (Masse Maximale) :

• La Surprise : Peu importait quelle recette de départ (Molle ou Rigide) ils utilisaient. Les indices réels étaient si puissants qu'ils ont forcé les deux recettes à s'accorder sur la même réponse.
• Le Verdict : Lorsqu'ils ont utilisé les indices les plus fiables (le « Grand Crash » et la « Lampe de Poche »), la masse maximale qu'une étoile à neutrons peut supporter est d'environ 2,2 à 2,3 fois la masse de notre Soleil.
• La « Twist » du Poids Lourd : S'ils supposent que cet objet lourd mystérieux (GW190814) est une étoile à neutrons, la limite grimpe à environ 2,6 à 2,7 fois la masse du Soleil. Cependant, cela crée un conflit avec les indices de « compressibilité » du Grand Crash, rendant la situation délicate.

La Limite de Taille (Rayon) :

• La Différence : Contrairement au poids, la taille de l'étoile dépendait de la recette de départ utilisée.
• Le Verdict : La recette « Molle » prédisait un rayon d'environ 11,8 km, tandis que la recette « Rigide » prédisait environ 12,4 km.
• Le Point Doux : Lorsque tous les meilleurs indices sont combinés, la taille la plus probable pour ces étoiles est d'environ 12 kilomètres (plus ou moins 1 km).

4. La Grande Image

L'article conclut qu'en examinant les « points d'extrémité » (les étoiles les plus lourdes et les plus grandes possibles) et en utilisant un mélange de données astronomiques réelles, nous pouvons affiner les règles de la matière la plus dense de l'univers.

• Le Poids : L'univers semble avoir une « limite de vitesse » pour la masse maximale d'une étoile à neutrons, se situant confortablement autour de 2,2 à 2,3 masses solaires. Cela correspond à l'étoile à neutrons la plus lourde que nous ayons réellement observée jusqu'à présent.
• La Taille : Elles ont à peu près la taille d'une petite ville, soit environ 12 km de diamètre.
• L'Essentiel : Les observations réelles (les indices) sont bien plus puissantes que les hypothèses théoriques de départ. Peu importe la théorie avec laquelle vous commencez, les données provenant des étoiles elles-mêmes forcent la réponse à converger vers les mêmes chiffres.

En résumé, l'univers nous a donné une réponse très claire : les étoiles à neutrons peuvent devenir incroyablement lourdes, mais il existe un plafond strict, et elles sont étonnamment petites compte tenu de leur masse.

Résumé technique

1. Énoncé du problème

Le défi central en physique nucléaire et des ions lourds est de déterminer l'Équation d'État (EOS) de la matière fortement interactive aux densités et potentiels chimiques intermédiaires. Alors que la matière à faible densité est contrainte par les expériences nucléaires et la QCD sur réseau, et la matière à haute densité par la QCD perturbative (pQCD), la région pertinente pour les cœurs d'étoiles à neutrons (plusieurs fois la densité de saturation nucléaire, $\rho_0$) reste mal comprise.

Le problème spécifique abordé est la détermination de la masse maximale ($M_{TOV}$) et du rayon correspondant ($R_{TOV}$) des étoiles à neutrons. Ces points extrêmes des séquences masse-rayon ($M(R)$) sont cruciaux pour comprendre la rigidité de l'EOS aux densités supranucléaires. Les auteurs visent à quantifier comment les observations multimessagers actuelles contraignent ces points extrêmes et dans quelle mesure les résultats dépendent du choix du modèle hadronique de base sous-jacent.

2. Méthodologie

A. Construction d'EOS hybrides

Les auteurs construisent des familles d'EOS hybrides causales en reliant la matière hadronique à faible densité à la matière quark à haute densité :

• Régime de faible densité : Deux modèles hadroniques représentatifs sont utilisés comme bases :
  • SFHo : Un modèle de champ moyen relativiste (EOS plus souple, incompressibilité $K=245$ MeV).
  • DD2 : Une EOS hadronique plus rigide.
• Régime de haute densité : Un modèle Sigma linéaire étendu (eLSM) basé sur la symétrie chirale $U(3) \times U(3)$, incorporant des quarks constituants et des nonets de mésons. Ce modèle est contraint pour correspondre à la thermodynamique de la QCD sur réseau à faible potentiel chimique et aux résultats de la pQCD à haut potentiel chimique.
• Interpolation : Pour relier les deux régimes, une interpolation polynomiale du cinquième ordre est appliquée à la densité d'énergie en fonction de la densité baryonique. Cela assure la cohérence thermodynamique et la continuité de la pression et de la vitesse du son ($c_s$) à travers la région de transition, évitant les transitions de phase du premier ordre importantes.
• Contrainte asymptotique : L'EOS quark-méson est requise pour se connecter de manière fluide aux résultats de la pQCD aux densités asymptotiquement élevées ($\mu \approx 2.6$ GeV).

B. Cadre bayésien

L'étude emploie un cadre de pondération bayésienne pour analyser les distributions de probabilité des points extrêmes de la séquence $M(R)$.

1. A priori : Un large ensemble d'EOS hybrides causales est généré, couvrant une gamme de rigidité aux densités supranucléaires.
2. Pondération par vraisemblance : L'ensemble est pondéré par rapport aux données observationnelles pour générer des distributions de probabilité a posteriori pour $M_{TOV}$ et $R_{TOV}$.

C. Contraintes observationnelles

Les contraintes suivantes sont appliquées séquentiellement et en combinaison :

1. Contrainte pQCD : Assure que l'EOS approche le point de référence N3LO de la pQCD ($\mu=2.6$ GeV, $\rho=6.47$ fm$^{-3}$, $p=3823$ MeV/fm$^3$) tout en maintenant la causalité ($c_s < 1$).
2. GW170817 : Contrainte de déformabilité tidale ($\tilde{\Lambda} < 720$) et une contrainte combinée plus stricte ("$\Lambda_{190}$" : $70 < \Lambda_{1.4} < 580$ et $9.1 < R_{1.4} < 12.8$ km).
3. NICER : Mesures masse-rayon pour cinq pulsars (J0030+0451, J0740+6620, J0614–3329, J1231–1411, J0437–4715) avec une précision d'environ $10\%$.
4. HESS J1731–347 : Un candidat pour une étoile à neutrons de très faible masse ($M \approx 0.77 M_\odot$), traité avec prudence en raison des incertitudes d'interprétation.
5. GW190814 : L'objet du "trou de masse" ($M \approx 2.59 M_\odot$). Les auteurs testent l'hypothèse que cet objet est une étoile à neutrons en appliquant une contrainte gaussienne centrée sur $2.59 M_\odot$.

3. Contributions clés

• Diagnostic unifié : Démontre que les distributions de probabilité des points extrêmes de la séquence $M(R)$ servent de diagnostic sensible et complémentaire pour contraindre le comportement de l'EOS à haute densité dans un cadre multimessager.
• Indépendance vs dépendance de la base : Dissocie systématiquement l'influence des données observationnelles du choix théorique de la base hadronique (SFHo vs DD2).
• Identification de tensions : Identifie une tension spécifique entre l'EOS rigide DD2, l'hypothèse de l'étoile à neutrons du "trou de masse" et les contraintes strictes de déformabilité tidale ($\Lambda_{190}$).

4. Résultats clés

Distributions de la masse maximale ($M_{TOV}$)

• Dominance des observations : Le pic de la distribution $M_{TOV}$ est principalement piloté par les contraintes observationnelles plutôt que par le choix de l'EOS hadronique de base (SFHo ou DD2).
• Progression des contraintes :
  • pQCD uniquement : Produit un pic légèrement supérieur à $1.5 M_\odot$.
  • GW170817/HESS : Aucun décalage significatif de la position du pic.
  • NICER : Déplace le pic vers 2.1–2.2 $M_\odot$ et affine la distribution.
  • Trou de masse (GW190814) : Déplace significativement le pic vers 2.6–2.7 $M_\odot$.
• Conclusion robuste : Lors de l'application des contraintes les plus robustes (pQCD + GW170817 + NICER), les deux modèles EOS convergent vers une masse maximale préférée de 2.2–2.3 $M_\odot$. Cela est cohérent avec la masse observée de PSR J0952–0607 ($2.35 \pm 0.17 M_\odot$).

Distributions du rayon ($R_{TOV}$)

• Dépendance de la base : Contrairement à la masse, la distribution du rayon montre une forte dépendance à l'EOS hadronique sous-jacente.
  • SFHo+eLSM : Pic près de 11.8 km (avec des contraintes robustes).
  • DD2+eLSM : Pic près de 12.4 km (avec des contraintes robustes).
• Préférence générale : Les contraintes combinées favorisent une plage de rayons de $12 \pm 1$ km.
• Impact de la déformabilité tidale : L'application de la contrainte $\Lambda_{190}$ plus stricte rétrécit les distributions de rayons et abaisse les valeurs préférées, car une déformabilité tidale plus faible est corrélée à des rayons plus petits.

La tension du "trou de masse"

• Si l'objet GW190814 est interprété comme une étoile à neutrons, l'EOS basée sur DD2 échoue à fournir une description cohérente lorsqu'elle est combinée avec la contrainte tidale $\Lambda_{190}$. La rigidité requise pour soutenir une étoile de $2.6 M_\odot$ dans le modèle DD2 entraîne des déformabilités tidales trop grandes pour satisfaire les limites de GW170817.
• L'EOS basée sur SFHo (plus souple) gère cette combinaison légèrement mieux, mais fait toujours face à des défis.

5. Importance

Ce travail fournit une évaluation statistique rigoureuse de la masse maximale des étoiles à neutrons, allant au-delà des estimations ponctuelles uniques pour atteindre des distributions de probabilité.

1. Validation des modèles actuels : Les résultats soutiennent l'existence d'étoiles à neutrons jusqu'à $\sim 2.2-2.3 M_\odot$, validant les modèles de physique nucléaire actuels face aux données astrophysiques.
2. Contrainte sur les transitions de phase : La procédure d'appariement fluide et les distributions résultantes suggèrent que de grandes transitions de phase du premier ordre sont peu susceptibles de se produire dans la plage de densité sondée par les observations actuelles, ou sont du moins fortement contraintes.
3. Perspectives futures : L'étude souligne que la précision future des mesures de rayon (par exemple, de NICER ou de futures missions X) et la nature définitive de l'objet GW190814 (étoile à neutrons vs trou noir) sont cruciales pour résoudre les incertitudes restantes dans l'EOS à haute densité. La tension entre l'hypothèse du trou de masse et les contraintes de déformabilité tidale suggère que si GW190814 est effectivement une étoile à neutrons, l'EOS doit être plus souple que ce que le modèle DD2 prédit, ou que les contraintes tidales nécessitent une réévaluation.