Are NICER and GW170817 constraints suggesting a compactified scenario for Neutron stars?

Résumé Technique : Les contraintes de NICER et GW170817 suggèrent-elles un scénario de compactification pour les étoiles à neutrons ?

Énoncé du problème
L'équation d'état (EoS) de la matière à des densités dépassant la saturation nucléaire reste un problème fondamental non résolu en physique. Si la théorie du champ effectif chiral ($\chi$EFT) modélise de manière fiable les densités subnucléaires et que la QCD perturbative (pQCD) décrit les régimes de haute densité asymptotiques, la région de densité intermédiaire au sein des cœurs d'étoiles à neutrons (NS) est mal comprise. Des observations astrophysiques récentes provenant de l'instrument NICER (Neutron Star Interior Composition Explorer) et de l'événement d'ondes gravitationnelles GW170817 ont fourni des contraintes précises sur les masses, les rayons et les déformabilités de marée des étoiles à neutrons. Cependant, il reste incertain si ces contraintes favorisent une classe spécifique d'EoS (monotone, non monotone ou discontinue) ou si elles révèlent un motif structurel universel de la matière dense. Plus précisément, cet article examine si les données actuelles suggèrent un scénario de « compactification » où la gravité favorise naturellement les configurations qui maximisent la compacité stellaire plutôt que celles qui maximisent simplement la masse.

Méthodologie
Les auteurs utilisent un cadre d'inférence bayésienne complet pour analyser trois ensembles distincts et indépendants de classes d'EoS :

1. Monotone : Caractérisée par une vitesse du son ($c_s^2$) qui augmente de manière monotone avec la densité (typique des EoS hadroniques).
2. Non monotone : Caractérisée par un pic de $c_s^2$ suivi d'une diminution (associée à des transitions de phase lisses).
3. Discontinue : Caractérisée par un saut de densité ($\Delta n$) lors d'une transition de phase (associée à des transitions de phase de premier ordre fortes).

Construction de l'EoS :
L'EoS est construite selon une approche hybride. La matière à faible densité ($< 0,5 n_0$) est définie par l'EoS tabulée BPS. La région de $0,5 n_0$ à $1,1 n_0$ est modélisée par une polytrope dont l'indice adiabatique $\Gamma$ est contraint par les bandes d'incertitude de la $\chi$EFT. Pour les densités supérieures à $1,1 n_0$, le carré de la vitesse du son $c_s^2$ est paramétré comme une fonction piecewise-linéaire du potentiel chimique $\mu$ à travers quatre segments randomisés, s'étendant jusqu'au régime pQCD ($\mu \approx 2,6$ GeV).

Cadre d'inférence :
Contrairement aux méthodes antérieures de ré-échantillonnage par rejet qui génèrent des millions d'échantillons d'EoS pour les filtrer, cette étude utilise l'algorithme d'échantillonnage imbriqué UltraNest pour explorer directement la distribution postérieure. Le modèle contient 10 paramètres libres (les potentiels chimiques et les valeurs du carré de la vitesse du son pour les segments, plus $\Gamma_{CET}$ et un paramètre de saut de densité $\delta n$ pour la classe discontinue).
L'analyse intègre :

• Données NICER : Distributions postérieures conjointes pour la masse et le rayon de PSR J0030+0451, PSR J0740+6620, PSR J0437−4715, et le plus récent PSR J0614−3329.
• Données GW170817 : Échantillons postérieurs des masses des composantes et des déformabilités de marée de la fusion d'étoiles à neutrons binaires.
• Contraintes $\chi$EFT : Vraisemblances gaussiennes pour la pression à basses densités.

L'étude compare les distributions postérieures des paramètres de l'EoS et des propriétés stellaires dérivées (séquences Masse-Rayon) sous différentes combinaisons de contraintes observationnelles.

Résultats clés

1. Pouvoir de contrainte : L'analyse bayésienne confirme que les données observationnelles contraignent efficacement la région de faible densité de l'EoS (spécifiquement $\Gamma_{CET}$ et le premier segment de la vitesse du son $c_{s,1}^2$), qui reste robuste à travers toutes les classes d'EoS. En revanche, les paramètres de haute densité ($c_{s,2}^2, c_{s,3}^2$) présentent des distributions postérieures plus larges, indiquant une incertitude significative persistante.
2. Impact des nouvelles données : L'inclusion de PSR J0614−3329 et GW170817 déplace systématiquement les paramètres de potentiel chimique vers des valeurs plus élevées et réduit les limites de masse maximale. Cependant, l'évidence bayésienne ($\ln Z$) ne favorise pas fortement une classe d'EoS par rapport à une autre ; les différences d'évidence se situent dans la plage « inconclusive à faible ».
3. Masse maximale vs Compacité maximale : Une distinction cruciale est établie entre l'EoS qui produit la masse possible maximale ($M_{TOV}$) et l'EoS qui produit la compacité maximale ($C = GM/Rc^2$).
   • L'EoS de Masse Maximale est rigide dès les très faibles densités, fournissant un support de pression uniforme pour soutenir les masses les plus grandes. Ces configurations occupent le bord à grand rayon de la distribution postérieure M-R.
   • L'EoS de Compacité Maximale présente une EoS plus molle à basse/moyenne densité (minimisant le rayon) suivie d'un durcissement prononcé à haute densité (soutenant la masse).
4. Préférence Observationnelle : L'EoS la plus favorisée par les données combinées de NICER et GW170817 (le $\ln Z$ le plus élevé) ressemble structurellement à l'EoS de Compacité Maximale, et non à l'EoS de Masse Maximale. L'EoS favorisée par l'observation présente une relative mollesse aux densités intermédiaires, passant à une rigidité à haute densité.
5. Transitions de phase : Pour la classe discontinue, les données préfèrent de petits sauts de densité ($\delta n \approx 0,09-0,10$ fm$^{-3}$) et des transitions de phase tardives (survenant à des densités plus élevées). Bien qu'une transition de phase à faible densité puisse théoriquement produire une étoile hautement compacte avec un large cœur de quarks, les prédictions bayésiennes favorisent des configurations où la transition de phase survient plus tard, confinant la matière de quarks aux cœurs des étoiles suffisamment massives.

Signification et Revendications
L'article affirme que les contraintes astrophysiques actuelles sélectionnent des configurations d'EoS qui produisent des étoiles compactes plutôt que des étoiles maximalement massives. La convergence de l'EoS favorisée par l'observation avec la séquence de compacité maximale suggère une propriété fondamentale de la matière dense : l'EoS subit probablement un durcissement significatif aux densités intermédiaires (signalant potentiellement l'apparition de nouveaux degrés de liberté) avant de se rigidifier à nouveau pour soutenir les étoiles à neutrons observées de $\sim 2 M_\odot$.

Les auteurs concluent que cette structure « mou à rigide » est ce que la gravité favorise naturellement pour atteindre la compactification. Bien que les données ne rejettent pas définitivement des classes d'EoS spécifiques, elles pointent vers un motif universel où l'équation d'état passe de la mollesse à la rigidité, poussant les étoiles à neutrons vers leurs limites de compacité. L'article note avec modestie que de futures contraintes sur les pulsars de faible masse (pour mieux définir l'EoS de faible densité) et les pulsars de haute masse (pour contraindre l'EoS de haute densité) seront nécessaires pour affiner davantage ces conclusions et potentiellement rejeter certains scénarios, tels que les transitions de phase précoces.