Astrophysical constraints on the cold equation of state of the strongly interacting matter

Cette étude utilise des observations astrophysiques, notamment des mesures de pulsars massifs, des données NICER et des contraintes de déformabilité tidale de GW170817, ainsi que des calculs de QCD perturbative, pour restreindre considérablement l'espace paramétrique admissible pour l'équation d'état de la matière froide, dense et fortement interactive.

L'essentiel

Imaginez l'univers comme une gigantesque cuisine avec une recette très spécifique, impossible : une matière froide et dense.

Sur Terre, nous ne pouvons pas cuisiner ce plat. Nos collisionneurs de particules les plus puissants (comme le LHC) sont comme des fours à haute température ; ils peuvent fracasser des atomes les uns contre les autres, mais ils créent une soupe chaude et chaotique qui ne nous apprend pas grand-chose sur ce qui se passe lorsque la matière est comprimée à froid et serrée.

L'unique endroit dans l'univers où cette recette « froide et dense » existe réellement se trouve à l'intérieur des étoiles à neutrons. Ce sont les restes cosmiques d'étoiles massives qui se sont effondrées. Elles sont si lourdes qu'une cuillère à café de leur matière pèserait un milliard de tonnes sur Terre. Parce qu'elles sont si denses, elles agissent comme le seul laboratoire naturel de l'univers pour étudier comment la matière se comporte sous une pression extrême.

Le Mystère : « L'Équation d'État »

Les physiciens veulent connaître l'« Équation d'État » (EoS). Considérez l'EoS comme le manuel d'instructions de cette matière dense. Il nous dit : Si vous serrez cette matière plus fort, combien résiste-t-elle ? Devient-elle molle, ou se transforme-t-elle en quelque chose de plus dur que le diamant ?

Le problème est que nous n'avons pas le manuel. Nous devons deviner les règles en observant les étoiles à neutrons et en voyant comment elles se comportent.

Le Travail de Détective : Utiliser des Indices pour Affiner les Hypothèses

Les auteurs de cet article ont agi comme des détectives tentant de résoudre un mystère. Ils ont commencé par une immense bibliothèque de 10 000 manuels d'instructions possibles (des théories sur le comportement de la matière). La plupart de ces manuels n'étaient que des hypothèses basées sur des principes mathématiques et physiques.

Ensuite, ils ont utilisé des indices réels provenant de l'espace pour rayer les manuels qui ne correspondaient pas. Voici les indices qu'ils ont utilisés :

1. L'Indice « Poids Lourd » (Masse) :
   Nous savons qu'il existe une étoile à neutrons appelée pulsar « Veuve Noire » qui est incroyablement lourde (environ 2,22 fois la masse de notre Soleil).

   • L'Analogie : Imaginez que vous avez une pile de 10 000 ponts différents. Vous savez avec certitude qu'un camion pesant 2,22 tonnes a traversé l'un d'eux sans qu'il ne s'effondre. Tout design de pont qui se serait effondré sous ce poids est immédiatement jeté à la poubelle.
   • Résultat : Cet unique indice a éliminé environ 80 % des manuels possibles.

2. L'Indice « Limites de Vitesse » (pQCD) :
   Au cœur même d'une étoile à neutrons, la matière est si dense que les lois de la physique changent, et nous pouvons utiliser un type spécifique de mathématiques (la QCD perturbative) pour prédire ce qui se passe.

   • L'Analogie : C'est comme savoir que peu importe comment vous concevez une voiture, elle ne peut pas légalement rouler plus vite que la vitesse de la lumière. Si un design de pont implique que la voiture dépasserait la vitesse de la lumière, il est invalide.
   • Résultat : Cela a éliminé quelques manuels supplémentaires qui étaient physiquement impossibles.

3. L'Indice « Mollesse » (Déformabilité de Marée) :
   Lorsque deux étoiles à neutrons entrent en collision (comme l'événement GW170817), elles s'étirent mutuellement comme du sucre filé avant de fusionner. Cette « élasticité » est appelée déformabilité de marée.

   • L'Analogie : Imaginez deux guimauves entrant en collision. Si elles sont très rigides, elles changent à peine de forme. Si elles sont molles, elles s'écrasent beaucoup. Les ondes gravitationnelles de la collision nous disent exactement combien elles se sont écrasées.
   • Résultat : Ce fut le filtre le plus important. Il s'est avéré que la plupart des manuels restants prédisaient des étoiles à neutrons soit trop rigides, soit trop molles par rapport à ce que nous avons observé lors de la collision. Cet indice seul a réduit la liste des manuels valides à moins de 2 %.

4. L'Indice « Taille » (NICER) :
   Le télescope NICER sur la Station spatiale internationale prend des images en rayons X des étoiles à neutrons pour mesurer leur taille (rayon).

   • L'Analogie : C'est comme mesurer la circonférence de la guimauve.
   • Résultat : Bien qu'utiles, les mesures de NICER comportent encore un certain « flou » (incertitude). Elles ont aidé à réduire la liste, mais elles n'étaient pas aussi strictes que l'indice de « mollesse ».

Qu'Ont-ils Découvert ?

Après avoir appliqué tous ces filtres, les auteurs ont découvert que le « manuel d'instructions » de la matière dense est beaucoup plus spécifique que nous ne le pensions.

• Le « Point Doux » : La matière à l'intérieur de ces étoiles semble subir une transition. Elle commence comme une matière atomique normale (hadrons) puis se transforme en une soupe de quarks (les composants de base des protons et des neutrons).
• La Transition : Ce changement ne se produit pas instantanément comme un interrupteur (un saut net) ; il se produit progressivement, comme une transition douce. Les auteurs ont constaté que cette transition se produit probablement à une densité environ 4,8 fois supérieure à la densité d'un noyau atomique normal.
• La Taille : Les manuels valides suggèrent que les étoiles à neutrons sont généralement assez grandes (environ 12 à 13 km de rayon) et pas aussi petites que le suggéraient certaines autres théories.

Les Scénarios « Et Si »

Les auteurs ont également testé deux cartes blanches :

1. L'Étoile « Minuscule » : Il existe un objet candidat qui pourrait être une étoile à neutrons très légère. Si cela est réel, cela obligerait les règles à changer encore davantage. Cependant, les auteurs notent que cet objet est controversé et pourrait même ne pas être une étoile à neutrons.
2. L'Étoile « Écart » : Un objet mystérieux a été détecté lors d'une collision (GW190814) qui est plus lourd que n'importe quelle étoile à neutrons connue mais plus léger qu'un trou noir. Si cet objet est une étoile à neutrons, cela constituerait une contrainte majeure, obligeant le « manuel d'instructions » à être très rigide pour supporter ce poids.

La Conclusion

L'article conclut que les observations d'étoiles à neutrons sont le filtre ultime. Bien que nous ayons de nombreuses théories sur le fonctionnement de la matière, l'univers est très exigeant. La combinaison des étoiles les plus lourdes connues et de la « mollesse » observée lors des collisions d'étoiles a considérablement réduit les possibilités.

Actuellement, les indices les plus restrictifs sont la masse des étoiles les plus lourdes et la déformabilité de marée provenant des collisions. Les mesures de « taille » provenant des télescopes sont utiles mais encore un peu trop floues pour être le facteur décisif. Les auteurs se retrouvent avec un ensemble spécifique de règles que la matière doit suivre, mais ils admettent qu'il reste encore du travail à faire pour comprendre exactement pourquoi la matière se comporte ainsi.

Résumé technique

1. Énoncé du problème

Les propriétés fondamentales de la matière fortement interactive, froide et dense, demeurent l'un des problèmes non résolus centraux de la physique nucléaire. Bien que la Chromodynamique Quantique (QCD) soit bien comprise aux faibles densités (via la physique nucléaire) et aux densités asymptotiquement élevées (via la QCD perturbative), le régime intermédiaire — où l'on s'attend à ce que réside le point critique final — est inaccessible aux expériences terrestres. Les collisions d'ions lourds sondent cette région mais à des températures élevées ($\sim 100$ MeV), rendant difficile l'isolement de l'équation d'état (EOS) froide.

Les étoiles à neutrons (EN) servent de seuls laboratoires naturels pour étudier la matière à des densités plusieurs fois supérieures à la densité de saturation nucléaire ($\rho_0$) à des températures proches de zéro. Cependant, l'EOS de la matière dense n'est pas uniquement déterminée par les observations actuelles. L'article vise à quantifier comment diverses observables astrophysiques contraignent l'EOS, en investiguant spécifiquement la transition de la matière hadronique à la matière quark et les contraintes résultantes sur la masse, le rayon et la déformabilité tidale des EN.

2. Méthodologie

2.1 Paramétrisation de l'EOS

Les auteurs construisent une EOS unifiée en reliant trois régimes de densité distincts :

1. Faible densité (Croûte) : Modélisée à l'aide de l'EOS de croûte unifiée de Douchin et Haensel.
2. Phase hadronique ($\rho \lesssim \rho_{BL}$) : Modélisée à l'aide de l'EOS relativiste de champ moyen SFHo (basée sur les nucléons et les mésons $\omega, \rho, \sigma$). Les auteurs testent également la robustesse en utilisant l'EOS DD2.
3. Phase quark ($\rho \gtrsim \rho_{BU}$) : Modélisée à l'aide d'un modèle quark-méson basé sur la symétrie chirale $U(3) \times U(3)$ (eLSM), qui inclut les quarks constituants et les nonets de mésons.
4. Haute densité asymptotique : Contrainte par les calculs de QCD perturbative (pQCD) à $\mu \approx 2.6$ GeV ($\rho \approx 40\rho_0$).

Interpolation : La transition entre les phases hadronique et quark est gérée via une interpolation polynomiale du cinquième ordre de la densité d'énergie $\varepsilon(\rho_B)$ dans une région de densité intermédiaire définie par une densité centrale $\bar{\rho}$ et une largeur $\Gamma$. Cela assure la continuité de la pression et de la vitesse du son ($c_s$).

Paramètres libres : L'analyse fait varier trois paramètres clés :

• $g_V$ : Constante de couplage vectoriel entre les mésons vectoriels et les quarks constituants ($0 \le g_V \le 10$).
• $\bar{\rho}$ : Densité centrale de la région de transition ($2\rho_0 \le \bar{\rho} \le 5\rho_0$).
• $\Gamma$ : Largeur de la région de transition ($\rho_0 \le \Gamma \le 4\rho_0$).

2.2 Contraintes et données

L'étude emploie une inférence bayésienne pour contraindre l'espace des paramètres en utilisant un ensemble de données d'environ 10 000 EOS échantillonnées aléatoirement. Les contraintes sont appliquées soit comme des conditions strictes (coupes dures), soit comme des vraisemblances probabilistes :

• Conditions strictes :

  • Causalité et stabilité : $c_s < 1$ et stabilité thermodynamique.
  • Connexion pQCD : L'EOS doit se connecter de manière causale au point de référence pQCD à $\mu = 2.6$ GeV.
  • Masse maximale : L'EOS doit supporter une étoile non rotative (TOV) avec $M_{TOV} \ge 2.22 M_\odot$ (basé sur une analyse de population incluant PSR J0952–0607).

• Contraintes probabilistes (Vraisemblances) :

  • NICER : Mesures masse-rayon pour cinq pulsars (PSR J0030+0451, J0740+6620, J0614-3329, J1231-1411, J0437-4715).
  • GW170817 : Contraintes de déformabilité tidale. L'étude utilise une borne conservative $\tilde{\Lambda} < 720$ et une plage plus stricte $70 < \Lambda_{1.4} < 580$.
  • Scénarios hypothétiques :
    • HESS J1731–347 : Un objet compact très léger ($M \approx 0.77 M_\odot$).
    • GW190814 : Un objet de « trou de masse » ($M \approx 2.59 M_\odot$) interprété comme une étoile à neutrons.

3. Contributions clés

1. Analyse bayésienne systématique : L'article fournit un cadre statistique rigoureux pour cartographier la variété EOS autorisée en faisant varier la densité et la largeur de transition entre la matière hadronique et la matière quark, plutôt que d'assumer un modèle de transition fixe.
2. Quantification de la force des contraintes : Il classe explicitement le pouvoir contraignant des différentes observations astrophysiques, démontrant que la déformabilité tidale et la masse maximale sont les facteurs les plus restrictifs, tandis que les mesures de rayon actuelles de NICER (en raison de grandes incertitudes) ont un impact moindre sur le rétrécissement de la variété EOS.
3. Robustesse du modèle : Les auteurs valident leurs résultats en comparant le modèle SFHo-eLSM à des constructions alternatives (DD2-eLSM, SFHo-MHT, et un modèle de transition de phase du premier ordre BBKF), montrant que la préférence pour des paramètres de transition spécifiques est robuste à travers différents modèles hadroniques.
4. Caractéristiques de la transition : L'étude identifie que les données favorisent une transition de recouvrement large plutôt qu'une transition de phase du premier ordre abrupte.

4. Résultats clés

• Réduction de la variété EOS : À partir d'environ 9 261 EOS initiales, l'application de la causalité et de la stabilité strictes réduit ce nombre à environ 1 993. La contrainte de masse minimale ($M_{TOV} \ge 2.22 M_\odot$) est le filtre unique le plus puissant, réduisant l'ensemble viable à environ 600.
• Déformabilité tidale : La contrainte de déformabilité tidale de GW170817 est la deuxième plus restrictive. L'application de la borne conservative $\tilde{\Lambda} < 720$ réduit les EOS autorisées à moins de 2 % de l'ensemble original. La borne plus stricte ($70 < \Lambda_{1.4} < 580$) réduit encore davantage l'espace viable.
• Paramètres privilégiés : L'ensemble de paramètres le plus probable trouvé est :
  • Densité de transition : $\bar{\rho} \approx 4.6 - 5.0 \rho_0$
  • Largeur de transition : $\Gamma \approx 2.8 - 3.0 \rho_0$
  • Couplage vectoriel : $g_V \approx 4.0 - 6.0$
• Prédictions de rayon : Le modèle prédit des rayons pour une étoile à neutrons de $1.4 M_\odot$ dans la plage de 12,5–12,8 km. Cela crée une tension avec la mesure NICER pour PSR J0614–3329, qui suggère un rayon plus petit ($\sim 10.3$ km), bien que les grandes incertitudes dans les données NICER empêchent cela d'être une exclusion définitive.
• Implications du trou de masse : Si l'objet GW190814 ($2.59 M_\odot$) est une étoile à neutrons, les contraintes EOS deviennent extrêmement sévères, éliminant efficacement l'espace de paramètres privilégié du modèle actuel.
• Objet HESS : L'existence d'une EN très légère ($0.77 M_\odot$) dans HESS J1731–347 a un impact négligeable sur les contraintes EOS dans ce cadre.

5. Signification

Ce travail démontre que les observations astrophysiques sont désormais le moteur principal pour contraindre l'EOS froide de la matière fortement interactive. L'étude met en évidence une tension critique : les modèles actuels qui satisfont la contrainte de haute masse ($>2.2 M_\odot$) et les limites pQCD tendent à prédire des rayons plus grands et des déformabilités tidales plus élevées que ce que suggèrent certaines mesures spécifiques de NICER.

La découverte qu'une transition de recouvrement large ($\Gamma \approx 3\rho_0$) est favorisée par rapport à une transition de phase du premier ordre abrupte ou à un recouvrement étroit suggère que la transition de phase de la matière hadronique à la matière quark dans les cœurs d'étoiles à neutrons est probablement un recouvrement lisse se produisant à des densités élevées ($\sim 5\rho_0$). Cela fournit une référence cruciale pour les développements théoriques futurs en QCD et guide l'interprétation des futures données multi-messagers provenant des détecteurs d'ondes gravitationnelles et des télescopes à rayons X.