Investigating Twin Star Equation of States in Light of Recent Astrophysical Observations

Résumé technique : Investigation des équations d'état à étoiles jumelles à la lumière des récentes observations astrophysiques

Énoncé du problème
La composition interne des étoiles à neutrons (EN) et la nature de la transition de phase (TP) de la matière hadronique (MH) vers la matière de quarks (MQ) demeurent non résolues. Bien que la chromodynamique quantique (QCD) prédise une transition à des densités élevées, les conditions spécifiques — telles que la densité d'apparition, la pression de transition et l'intensité de la discontinuité — sont inconnues. Une transition de phase du premier ordre forte pourrait conduire à l'existence d'« étoiles jumelles » (EJ) : un troisième rameau stable d'objets compacts dans le diagramme masse-rayon ($M-R$), séparé du rameau standard des étoiles à neutrons par une région instable. Ces étoiles jumelles posséderaient une masse égale mais des rayons plus petits que leurs homologues de neutrons. Le défi consiste à identifier les équations d'état (EoS) hybrides qui satisfont à la fois l'équilibre hydrostatique de la relativité générale et les contraintes de plus en plus précises de l'astrophysique multi-messagers, incluant les mesures de masse-rayon de NICER et les observations d'ondes gravitationnelles (GW170817).

Méthodologie
Les auteurs adoptent une approche indépendante du modèle pour construire un ensemble complet d'EoS hybrides capables de produire des solutions d'étoiles jumelles. La méthodologie repose sur les équations de Tolman-Oppenheimer-Volkoff (TOV) pour résoudre l'équilibre hydrostatique. La construction de l'EoS comprend quatre régions distinctes :

1. Région de la croûte : Une EoS polytropique par morceaux fixe en dessous de $\sim 0,4\rho_{sat}$.
2. Région hadronique : Modélisée par une polytrope avec un indice adiabatique $\gamma$ libre, déterminé par les conditions limites à la transition.
3. Région de transition de phase : Modélisée via une construction de Maxwell, caractérisée par une pression de transition constante ($P_{tr}$) et un saut discontinu de densité d'énergie ($\Delta e$).
4. Région de quarks : Modélisée en utilisant une paramétrisation à vitesse du son constante (CSS) ($c^2_{s,QM}$).

L'espace des paramètres est défini par quatre variables indépendantes : la densité d'énergie de transition ($e_{tr}$), la pression de transition ($P_{tr}$), l'intensité de la discontinuité de densité d'énergie ($\Delta e$), et la vitesse du son dans la matière de quarks ($c^2_{s,QM}$). L'étude fait varier systématiquement ces paramètres dans les limites de causalité pour identifier les conditions qui génèrent un troisième rameau stable dans la séquence $M-R$.

Principales contributions

• Définition des courbes en « chapeau de sorcière » (Witch-Hat) : Le papier introduit une représentation géométrique dans l'espace des paramètres $(P_{tr}, \Delta e)$ pour des valeurs fixes de $e_{tr}$ et $c^2_{s,QM}$. Ces courbes délimitent les frontières des intensités de discontinuité autorisées permettant la formation d'étoiles jumelles. Les auteurs identifient trois régions distinctes dans cet espace :
  • Région sans solution : Aucune étoile jumelle ne se forme, quel que soit $\Delta e$.
  • Région bornée : Des étoiles jumelles se forment pour une plage spécifique de $\Delta e$ en dessous d'une pression de transition maximale.
  • Région non bornée : Des étoiles jumelles se forment pour $\Delta e$ au-dessus d'une limite minimale, sans limite supérieure sur $\Delta e$ (dans le cadre numérique).
• Analyse de la causalité et de la complétude : L'étude révèle que les courbes en « chapeau de sorcière » peuvent être « ponctuées » ou incomplètes. Cela se produit lorsque la pression de transition maximale autorisée ($P_{tr}^{max}$) requise pour former une étoile jumelle excède la limite de causalité ($c_s = 1$) ou la pression centrale de la configuration de masse maximale hadronique.
• Intégration des récentes contraintes observationnelles : Les auteurs appliquent les intervalles de confiance à 2$\sigma$ (ainsi que 1$\sigma$/3$\sigma$) des mesures récentes de NICER (PSR J0614–3329, PSR J0437–4715, PSR J0740+6620, PSR J0030+0451), du pulsar massif PSR J0348+0432, et de la déformabilité de marée de GW170817 pour contraindre l'espace des paramètres.

Résultats
L'application des données observationnelles restreint considérablement l'espace des paramètres viable pour les étoiles jumelles :

• Contraintes sur la densité de transition : Les mesures récentes de PSR J0614–3329 et PSR J0437–4715 suggèrent une EoS plus molle autour de la masse canonique, ce qui contraint étroitement la densité d'énergie de transition. L'étude établit une limite supérieure sur la densité d'énergie de transition de $e_{tr}^{max} \approx 1,40 \times 10^{15}$ g/cm$^3$ ($\approx 4,03\rho_{sat}$) pour l'EoS de quarks la plus raide ($c^2_{s,QM}=1,0$).
• Contraintes sur l'intensité de la discontinuité : La discontinuité maximale autorisée de la densité de masse au repos est trouvée à $\Delta\rho_{max} \approx 7,76\rho_{sat}$ ($\Delta e_{max} \approx 2,94 \times 10^{15}$ g/cm$^3$). Cette limite est régie par les EoS produisant des étoiles jumelles de Catégorie II (où le second rameau supporte $>2,0 M_\odot$ mais pas le troisième).
• Masse maximale des étoiles jumelles : La masse maximale atteignable par les étoiles jumelles est calculée à $M_{TS}^{max} \approx 2,05 M_\odot$ (pour $c^2_{s,QM}=1,0$ et 2$\sigma$ de confiance). Cette limite est fixée par les étoiles jumelles de Catégorie I (où les deux branches supportent $>2,0 M_\odot$).
• Impact des niveaux de confiance : Les résultats sont sensibles au niveau de confiance des inférences observationnelles. À un niveau de confiance strict de 1$\sigma$, les étoiles jumelles de Catégorie I et II sont entièrement exclues, ne laissant que des solutions de Catégorie III et IV avec des masses maximales nettement plus faibles ($\approx 1,37 M_\odot$). Inversement, à 3$\sigma$, les contraintes se relâchent, permettant une gamme plus large de paramètres, incluant des cas où la pression de transition se situe dans la région « non bornée ».
• Dépendance aux catégories : L'étude catégorise les étoiles jumelles en fonction des masses maximales de leurs seconds ($M_{TOV,2}$) et troisièmes ($M_{TOV,3}$) rameaux. Elle conclut que les étoiles jumelles de Catégorie I (massives sur les deux branches) ne sont permises qu'à des densités de transition extrêmement élevées ($\approx 4\rho_{sat}$), tandis que des densités de transition plus douces favorisent les catégories II, III et IV.

Signification
L'article affirme que l'observation d'étoiles jumelles, particulièrement massives, constituerait une signature de type « smoking gun » d'une forte transition de phase du premier ordre hadron-quark dans les cœurs des étoiles à neutrons. En cartographiant les courbes en « chapeau de sorcière » et en appliquant des contraintes observationnelles rigoureuses, l'étude fournit des limites thermodynamiques spécifiques sur la densité de transition ($e_{tr}$), l'intensité de la discontinuité ($\Delta e$) et la masse maximale des étoiles jumelles. Ces limites offrent un cadre critique pour interpréter les futures données astrophysiques et restreindre la nature possible de la matière dense. Les auteurs notent que bien que leur analyse soit basée sur l'équilibre hydrostatique, l'existence des étoiles jumelles de Catégorie II est particulièrement fragile ; des études dynamiques antérieures suggèrent qu'elles pourraient être instables face aux perturbations, migrant vers les branches hadroniques ou s'effondrant en trous noirs. Ainsi, l'identification d'EoS d'étoiles jumelles viables est une étape nécessaire, bien que potentiellement insuffisante, vers la confirmation de l'existence de la matière de quarks dans les étoiles à neutrons.