Is the coexistence of strange quark stars and hadronic stars favored by astrophysical data? A Bayesian analysis

Résumé technique : Coexistence des étoiles de quarks étranges et des étoiles hadroniques

Énoncé du problème
La composition microscopique des étoiles compactes et l'équation d'état (EOS) de la matière fortement interagissante à haute densité demeurent des problèmes centraux en physique nucléaire et en astrophysique. Une tension spécifique existe entre les contraintes observationnelles exigeant une EOS « molle » aux densités intermédiaires (pour expliquer les petits rayons des objets de masse ordinaire et la déformabilité de marée de GW170817) et le besoin d'une EOS « raide » aux hautes densités (pour soutenir des pulsars dépassant $2 M_\odot$). De plus, l'inclusion d'hyperons et de résonances $\Delta$ dans la matière hadronique ramollit typiquement l'EOS, menant au « problème des hyperons » où les masses maximales tombent en dessous des valeurs observées.

Le scénario des « deux familles » (2F) propose une résolution : la coexistence de deux branches déconnectées d'étoiles compactes. Dans ce cadre, les étoiles hadroniques (HS) sont métastables et peuplent une branche avec des rayons plus petits, tandis que la matière de quarks étranges (SQM) absolument stable forme une branche auto-liée d'étoiles à quarks (QS) capable de supporter des masses plus élevées. Bien que ce scénario ait été exploré théoriquement, il n'avait pas encore été soumis à une analyse bayésienne rigoureuse comparant directement l'hypothèse standard de la « famille unique » (1F, où toutes les étoiles sont hadroniques) en utilisant un ensemble unifié de données astrophysiques et de laboratoire.

Méthodologie
Les auteurs réalisent la première analyse bayésienne détaillée pour contraindre les paramètres de l'EOS pour les scénarios 1F et 2F simultanément.

• Modèles :
  • Secteur Hadronique : Un modèle de champ moyen relativiste (RMF) non linéaire incluant des nucléons, des hyperons ($\Lambda, \Sigma, \Xi$) et des résonances $\Delta$. Le modèle utilise des paramètres empiriques nucléaires à la saturation et des couplages dépendants de la densité pour le méson $\rho$ afin de contrôler la pente de l'énergie de symétrie ($L_{sym}$).
  • Secteur de Quarks (2F uniquement) : Une EOS de modèle de sac de type confinement de saveur de couleur (CFL) décrivant la SQM absolument stable. Ce secteur inclut les paramètres du paramètre de sac ($B$), du gap de couplage CFL ($\Delta_{CFL}$) et des corrections de la QCD perturbative ($a_4$).
• Données et vraisemblances : L'analyse incorpore :
  • Contraintes Masse-Rayon (M-R) : Postérieurs conjoints de NICER (PSR J0030+0451, PSR J0740+6620, PSR J0614–3329) et des observations de rayons X de HESS J1731–347.
  • Contraintes de masse uniquement : PSR J0952–0607 et la composante secondaire de GW190814 (traitée comme un potentiel objet compact).
  • Ondes gravitationnelles : Le postérieur conjoint complet de GW170817 (masses et déformabilités de marée), traité par marginalisation sur les étiquettes de famille (HS-HS, HS-QS, QS-HS, QS-QS) dans le scénario 2F.
  • Contraintes de laboratoire/théoriques : Bandes de la théorie effective de champs chiraux ($\chi$EFT) pour la matière neutronique pure et données de collision d'ions lourds (HIC) pour la matière nucléaire symétrique.
• Cadre statistique :
  • L'étude utilise l'échantillonnage imbriqué dynamique (via le package dynesty) pour calculer les distributions postérieures et les évidences bayésiennes ($Z$).
  • Le facteur de Bayes ($B_{2F,1F} = Z_{2F}/Z_{1F}$) est utilisé pour comparer les modèles, en tenant compte des différentes dimensions de l'espace des paramètres (11 paramètres pour 1F contre 14 pour 2F).
  • Dans le scénario 2F, l'assignation de la famille de chaque source observée est traitée comme une variable latente et marginalisée, avec des priors définis sur la stabilité structurelle et les canaux de formation astrophysiques.

Contributions clés

1. Première comparaison bayésienne : Ce travail fournit la première comparaison bayésienne rigoureuse entre les scénarios 1F et 2F en utilisant un ensemble complet de données incluant des contraintes nucléaires à basse densité et des observations astrophysiques de haute masse.
2. Cadre unifié : Il contraint simultanément l'EOS hadronique métastable et l'EOS de la SQM absolument stable au sein de l'hypothèse 2F, plutôt que de supposer la nature des sources observées a priori.
3. Résolution de la tension : L'analyse démontre comment le scénario 2F résout naturellement la tension entre la mollesse à densité intermédiaire et la rigidité à haute densité en séparant ces exigences en deux branches stellaires distinctes.

Résultats

• Préférence de modèle : Le scénario 2F est favorisé par rapport au scénario 1F pour toutes les combinaisons de jeux de données. Pour l'ensemble des données (incluant HIC et toutes les contraintes astrophysiques), le log-facteur de Bayes est $\ln B_{2F,1F} = 11,9$, indiquant un soutien « décisif » pour le scénario des deux familles.
• Caractéristiques de l'EOS :
  • Scénario 1F : Pour satisfaire à la fois les contraintes de petit rayon et la limite de $2 M_\odot$, l'EOS 1F nécessite un compromis finement équilibré. Ceci est réalisé par un début précoce des résonances $\Delta$ ($\sim 1,8 n_{sat}$) qui ramollit l'EOS aux densités intermédiaires, suivi d'un durcissement important aux hautes densités. La masse maximale est contrainte à $M_{max} \approx 2,05 M_\odot$.
  • Scénario 2F : La tension est levée. La branche hadronique reste molle (soutenue par les résonances $\Delta$ et les hyperons) avec une masse maximale de $M_{HS,max} \approx 1,54 M_\odot$ et un rayon de $R_{1,4} \approx 11,1$ km. Les pulsars massifs sont soutenus par la branche QS plus raide, atteignant $M_{QS,max} \approx 2,33 M_\odot$.
• Contraintes de paramètres : Le postérieur 2F préfère une masse effective de Dirac ($m^*/m$) plus grande et un ajustement moins extrême des couplages $\Delta$ par rapport au cas 1F. Les paramètres du secteur de quarks ($B, \Delta_{CFL}, a_4$) sont modérément contraints, avec une corrélation positive entre la constante de sac et le gap de couplage.
• Interprétation de GW170817 : Dans le scénario 2F, l'événement GW170817 est très probablement une fusion HS-HS, mais les canaux mixtes (HS-QS) contribuent de manière significative (~35 % de la vraisemblance), tandis que le canal QS-QS est défavorisé pour l'EOS la plus probable.

Signification et affirmations
L'article affirme que les données astrophysiques et de laboratoire actuelles sont plus naturellement organisées par deux familles d'étoiles compactes déconnectées que par une séquence purement hadronique unique. Le scénario 2F offre une explication naturelle à la coexistence d'objets de masse ordinaire très compacts (comme HESS J1731–347) et de pulsars très massifs, tout en résolvant efficacement le problème des hyperons.

Les auteurs soulignent que leur conclusion est robuste face aux variations des priors (potentiels des hyperons, plages de paramètres de la SQM) et des hypothèses de masse des sources, à condition d'utiliser l'ensemble des données. Cependant, ils notent que la préférence statistique est moins définitive si des ensembles de données restreints (excluant HESS J1731–347, PSR J0614–3329 et les données HIC) sont considérés.

L'étude exclut explicitement les étoiles hybrides (étoiles possédant des cœurs de quarks) de la comparaison actuelle, notant qu'une sélection de modèle complète devrait inclure les étoiles hybrides comme une troisième hypothèse dans des travaux futurs. Les auteurs concluent que les détections futures d'objets plus massifs, d'objets de masse ordinaire très compacts, ou de mesures précises de deux masses distinctes avec le même rayon, fourniraient des indications encore plus fortes en faveur du scénario des deux familles.