As above, so below: assessing extremeness of the neutron-star equation of state based on the unstable branch

Résumé technique : Évaluation de l'extrémisme de l'équation d'état des étoiles à neutrons basée sur la branche instable

Énoncé du problème
Les modèles microphysiques de la matière des étoiles à neutrons (EN) sont essentiels pour interpréter les observations astrophysiques, mais ils sont typiquement construits et validés uniquement jusqu'aux densités maximales atteintes par les étoiles à neutrons stables. Le comportement de l'équation d'état (EoS) aux densités supérieures à ces limites a historiquement été traité comme non contraint ou supposé suivre la trajectoire originale du modèle. Cependant, l'EoS doit ultimement se connecter au régime de la chromodynamique quantique perturbative (pQCD) à des densités asymptotiquement élevées ($n \sim 20\text{--}40 n_{\text{sat}}$). L'exigence que cette connexion soit causale, stable et thermodynamiquement cohérente impose des contraintes sur l'EoS à des densités intermédiaires, même dans les régions où ni les modèles phénoménologiques ni les calculs pQCD ne sont directement applicables. Cet article examine si les modèles hadroniques standards, lorsqu'ils sont étendus aux densités centrales maximales des EN stables, impliquent des continuations physiquement plausibles vers le régime pQCD ou s'ils nécessitent des comportements « extrêmes » et improbables.

Méthodologie
Les auteurs utilisent une approche non paramétrique et agnostique du modèle pour échantillonner l'espace des extensions valides de l'EoS. La méthodologie implique :

1. Cadre : Utilisation d'un pont par processus gaussien (GPB) pour générer un large ensemble d'extensions d'EoS reliant la densité de terminaison d'un modèle ($n_{\text{term}}$) au régime pQCD à $n_{\text{pQCD}} = 30 n_{\text{sat}}$.
2. Contraintes : L'échantillonnage impose la stabilité mécanique, la causalité (vitesse du son $c_s \leq 1$) et la cohérence thermodynamique. Les extensions sont diffusées pour imposer des longueurs de corrélation, permettant des caractéristiques nettes à la transition entre le modèle et l'extension tout en assurant un comportement lisse vers la limite pQCD.
3. Indice de tension pQCD ($I_{\text{pQCD}}$) : Pour quantifier l'« extrémisme » d'une connexion, les auteurs utilisent l'indice défini par Komoltsev & Kurkela (2022) et Komoltsev et al. (2024). Cet indice mesure à quel point une extension d'EoS valide est proche des constructions de différence de pression minimale ou maximale requises pour relier les points de basse et haute densité. $I_{\text{pQCD}} \approx 1$ indique une connexion nécessitant une rigidité maximale ou de fortes transitions de phase.
4. Sélection de modèles : Un ensemble représentatif d'EoSs froides en équilibre $\beta$ a été sélectionné à partir de la base de données CompOSE, classé en :
   • Purement hadroniques : Calculs microscopiques (APR, BL(chiral)), fonctionnelles de densité non relativistes (Sly4, BSk24) et fonctionnelles de densité relativistes (DD2, SFHo, FSU2R).
   • Modèles avec degrés de liberté additionnels : Modèles incluant de la matière de quarks (QHC18, V-QCD), des hyperons (DDHdeltaY4, CMF-7) ou de la matière quarkyonique.
5. Analyse : L'étude examine les extensions à partir de deux points de terminaison : la densité centrale d'une étoile canonique de $1.4 M_\odot$ et la densité centrale maximale ($n_{\text{TOV}}$) de l'étoile stable la plus massive supportée par chaque modèle.

Résultats clés

1. Densité canonique ($1.4 M_\odot$) : Lorsque les modèles ne sont considérés valides que jusqu'à la densité centrale d'une étoile de $1.4 M_\odot$, les extensions requises vers le régime pQCD sont relativement peu contraintes. Une large gamme de comportements, incluant divers profils de vitesse du son et d'anomalies de trace, est compatible avec les contraintes pQCD pour toutes les échelles de renormalisation ($X$).
2. Densité maximale ($n_{\text{TOV}}$) : La situation change radicalement lorsque les modèles sont supposés valides jusqu'à la densité de masse maximale.
   • Modèles purement hadroniques : La plupart des modèles purement nucléoniques (Sly4, BSk24, DD2, SFHo, BL) sont incompatibles avec les valeurs plus élevées de l'échelle de renormalisation pQCD ($X_{\text{high}}$). Pour les valeurs de $X$ plus faibles autorisées, les extensions nécessitent un adoucissement abrupt et dramatique immédiatement au-dessus de la densité TOV. Cela se manifeste par un saut de densité fort, de type transition de phase du premier ordre ( $\Delta n \approx 15 n_{\text{sat}}$), et un décrochement (kink) dans la relation masse-rayon en entrant dans la branche instable.
   • L'exception FSU2R : Le modèle FSU2R permet des extensions sur toute la gamme de $X$ car il présente un comportement plus mou et une anomalie de trace positive jusqu'à $n_{\text{TOV}}$. Cependant, les auteurs notent que ce modèle prédit des rayons $> 12$ km, ce qui entre en conflit avec les observations de NICER.
   • Modèles avec degrés de liberté additionnels : Les modèles incorporant des quarks, des hyperons ou de la matière quarkyonique (V-QCD, Quarkyonic) s'adoucissent généralement avant d'atteindre la densité TOV. Par conséquent, ils ne nécessitent pas la transition de phase spectaculaire et coïncidente immédiatement au-dessus de $n_{\text{TOV}}$ dont les modèles purement hadroniques ont besoin. Leurs extensions vers la pQCD sont moins restrictives.
3. Tendance quantitative : L'analyse de l'ensemble de la base de données CompOSE via l'indice $I_{\text{pQCD}}$ révèle que les EoSs purement hadroniques, considérées valides jusqu'à $n_{\text{TOV}}$, se regroupent près de $I_{\text{pQCD}} = 1$, indiquant un besoin générique de comportement extrême, de type premier ordre. Les modèles avec des degrés de liberté d'adoucissement se regroupent à des valeurs de $I_{\text{pQCD}}$ plus basses.

Signification et affirmations
L'article soutient que l'exigence de connecter la matière des étoiles à neutrons stables au régime pQCD de manière causale et thermodynamiquement cohérente agit comme un filtre puissant pour les modèles microphysiques. Les auteurs affirment que :

• Les modèles purement nucléoniques, si l'on suppose qu'ils sont valides jusqu'à la densité de masse maximale, impliquent un comportement très spécifique et « artificiel » : une forte transition de phase du premier ordre se produisant de manière fortuite juste au moment où l'étoile atteint sa limite de masse maximale. Les auteurs jugent cette coïncidence hautement improbable.
• Par conséquent, l'étude défavorise les descriptions purement nucléoniques pour toutes les étoiles à neutrons massives stables.
• Les résultats soutiennent l'apparition de degrés de liberté additionnels (tels que les quarks ou les hyperons) qui provoquent l'adoucissement de l'EoS avant la densité TOV, évitant ainsi le besoin d'extensions extrêmes et coïncidentes vers le régime pQCD.
• Ce travail fournit un outil (via le cadre GPB et l'indice $I_{\text{pQCD}}$) pour construire et visualiser ces contraintes, offrant des orientations pour la construction future de modèles d'EoS respectant les contraintes fondamentales de la QCD.

Les auteurs ne prétendent pas identifier un type spécifique de nouvelle matière, mais affirment que quelque mécanisme d'adoucissement sous $M_{\text{TOV}}$ est nécessaire pour éviter les contraintes extrêmes imposées par la cohérence pQCD sur les modèles purement hadroniques.