Quark and hybrid stars with renormalization group improvement of NNLO perturbative QCD

Énoncé du Problème
La description théorique de l'équation d'état (EoS) de la matière quarkique froide et dense est essentielle pour la modélisation des intérieurs des étoiles à neutrons (NS) et des étoiles hybrides. Bien que la chromodynamique quantique perturbative (pQCD) fournisse un cadre ab-initio pour les densités baryoniques élevées, son application aux densités intermédiaires pertinentes pour les cœurs d'étoiles à neutrons ($\sim 5-10 \rho_0$) est entravée par une sensibilité significative à l'échelle de renormalisation arbitraire ($\Lambda$). Spécifiquement, la pression de la pQCD au niveau suivant le suivant-suivant l'ordre de la perturbation (NNLO) présente une dépendance résiduelle à l'échelle importante, provenant principalement des contributions du secteur « dur ». Cette incertitude complique la détermination de l'existence ou non de la matière quarkique dans les cœurs des pulsars massifs observés (par exemple, PSR J0740+6620) et contraint la nature de la transition de phase hadron-quark. De plus, les calculs existants de la QCD sur réseau sont limités aux potentiels chimiques faibles en raison du problème de signe, et les modèles effectifs manquent souvent de liberté asymptotique.

Méthodologie
Les auteurs étendent le cadre de la Théorie de la Perturbation Optimisée par le Groupe de Renormalisation (RGOPT), précédemment appliquée à la matière symétrique, au cas de la matière quarkique non symétrique, en équilibre $\beta$ et de neutralité de charge avec $N_f = 2+1$ quarks massifs.

• Cadre RGOPT : Au lieu d'une expansion standard en faible couplage, la RGOPT introduit une expansion perturbative modifiée autour d'états de quasi-particules massives. Cela implique une masse de screening « dopée par le RG » ($\eta$) qui incorpore les propriétés de groupe de renormalisation (RG) à tous les ordres. La méthode utilise une substitution spécifique de paramètres (Éq. 2.3) et fixe le paramètre de masse arbitraire $\eta$ en résolvant l'équation de RG (Éq. 2.6) parallèlement à une condition de stationnarité, garantissant que la solution respecte la liberté asymptotique et reste réelle grâce à un changement de schéma de renormalisation (RSC) contrôlé.
• Consistance Thermodynamique : Pour appliquer la théorie aux étoiles compactes, les auteurs imposent l'équilibre $\beta$ et la neutralité de charge. Ils traitent le problème de l'inconsistance thermodynamique (où $dP/d\mu \neq \rho$) de la pression brute de la RGOPT en intégrant la différence entre les dérivées partielles et totales de la pression (Éq. 3.4).
• Formules Compactes : Reconnaissant la complexité de l'expression complète de la pression NNLO de la RGOPT, les auteurs dérivent des formules d'ajustement de type « pocket » pour la pression et la densité d'énergie thermodynamiquement cohérentes. Ces formules sont valables pour deux plages distinctes du paramètre d'échelle de renormalisation $X$ (où $\Lambda = X \mu_B/3$) : $2 \le X \le 3$ (pertinent pour les étoiles hybrides) et $3 \le X \le 6$ (pertinent pour les étoiles de quark pures).
• Construction d'Étoiles Hybrides : Pour modéliser les étoiles hybrides, les auteurs emploient une construction de Maxwell pour décrire une transition de phase du premier ordre entre les phases hadronique et quarkique. Ils combinent l'EoS de la RGOPT avec trois modèles hadroniques représentatifs basés sur la théorie du Champ Moyen Relativiste (RMF) : DD2 (rigide), SFHo (souple) et DIDY (souple, incluant les hyperons).

Contributions Clés

1. Extension de la RGOPT à l'équilibre $\beta$ : L'article fournit la première EoS de la RGOPT au niveau NNLO pour la matière froide et dense avec $N_f=2+1$ quarks massifs en équilibre $\beta$ et neutralité de charge, une étape nécessaire pour une modélisation stellaire réaliste.
2. Réduction de l'Incertitude d'Échelle : L'étude démontre que la pression de la RGOPT est nettement moins sensible aux variations du paramètre d'échelle de renormalisation $X$ par rapport à la pQCD NNLO standard.
3. Outils Analytiques : Les auteurs fournissent des formules d'ajustement compactes et explicites (Éqs. 3.9 et 3.13) pour la pression et la densité baryonique de la RGOPT, permettant une application immédiate dans les simulations astrophysiques sans recalculer les intégrales perturbatives complexes.
4. Comparaison Systématique : Une comparaison exhaustive est faite entre les étoiles de quark pures et les étoiles hybrides en utilisant à la fois la RGOPT et la pQCD standard, contraintes par les observations astrophysiques récentes.

Résultats

• Étoiles de Quark Pures (QS) :

  • Pour reproduire les masses maximales de $2,0, 2,3$ et $2,6 M_\odot$, la RGOPT nécessite des paramètres d'échelle dans la plage $X \approx 3,08 - 3,58$. En revanche, la pQCD standard nécessite une plage plus étroite et plus basse ($X \approx 2,95 - 3,26$) et présente une dépendance plus forte de la masse maximale vis-à-vis de $X$.
  • Si l'objet de faible masse de l'événement GW190814 est interprété comme une étoile à neutrons, la RGOPT nécessite $X \approx 4,10$, tandis que la pQCD nécessite $X \approx 3,56$.
  • Les prédictions de la RGOPT pour les étoiles de quark pures sont plus stables contre les variations d'échelle, avec une masse maximale changeant de $\sim 0,59 M_\odot$ par unité de $X$, contre $\sim 0,98 M_\odot$ pour la pQCD.
  • Les deux approches satisfont l'hypothèse de Bodmer-Witten pour la matière étrange stable.

• Étoiles Hybrides :

  • Les étoiles hybrides ont été construites en utilisant la construction de Maxwell avec les trois modèles hadroniques.
  • Pour les étoiles hybrides compatibles avec le pulsar massif PSR J0740+6620, le paramètre d'échelle de la RGOPT varie d'environ $1,99$à$2,98$ selon le modèle hadronique.
  • Les valeurs de $X$ les plus grandes produisent des étoiles hybrides avec de grands cœurs de quark (rayons de $5$à$8$ km). Les valeurs de $X$ les plus petites produisent des étoiles avec de très petits cœurs de quark.
  • La densité de transition et la force de la transition de phase dépendent du modèle hadronique. Pour le modèle rigide DD2, la transition se produit à des densités plus basses ($\sim 2 \rho_0$) comparé aux modèles plus souples SFHo et DIDY ($\sim 4 \rho_0$).
  • Les étoiles hybrides avec de grands cœurs de quark (compatibles avec PSR J0740+6620) ont généralement des enveloppes hadroniques de $3-5$ km et sont incompatibles avec l'objet de faible masse HESS J1731-347. Inversement, les étoiles hybrides avec de petits cœurs de quark (petit $X$) sont compatibles avec HESS J1731-347.
  • La vitesse du son dans les EoS hybrides dépasse la limite conforme ($\sqrt{1/3}$) dans la branche hadronique mais chute en dessous de celle-ci dans la phase quark après la transition.

Signification et Revendications
L'article affirme que le cadre RGOPT offre un outil théorique plus robuste pour décrire le secteur quarkique des étoiles compactes que la pQCD standard au même ordre de perturbation. En incorporant les propriétés de RG directement dans l'expansion perturbative, la RGOPT atténue la grande dépendance résiduelle à l'échelle qui affecte les prédictions de la pQCD standard aux densités intermédiaires. Cette réduction de l'incertitude théorique permet de contraindre plus étroitement l'échelle de renormalisation lors de la mise en correspondance avec les données astrophysiques. Les auteurs affirment que leurs formules compactes dérivées fournissent un moyen pratique et précis d'incorporer ces résultats QCD améliorés dans les calculs de structure stellaire, facilitant l'investigation de l'existence de cœurs de matière quarkique dans les étoiles à neutrons observées. Le travail souligne que, bien que la RGOPT et la pQCD puissent toutes deux accommoder les contraintes observationnelles actuelles, la RGOPT offre une base théorique plus stable pour explorer l'équation d'état à haute densité.