Probing the inner disk edge of rotating neutron stars with dipole excitations in finite nuclei

Résumé technique : Sondage du bord interne du disque de l'étoile à neutrons en rotation par les excitations dipolaires dans les noyaux finis

Énoncé du problème
La structure des étoiles à neutrons en rotation rapide et la dynamique orbitale au sein de leur espace-temps sont régies par l'équation d'état (EOS) nucléaire. Une incertitude critique dans les modèles actuels concerne le point de transition où l'orbite circulaire stable la plus interne (ISCO) atteint la surface stellaire. Cette « transition surface–ISCO » détermine si le bord interne d'un disque d'accrétion est défini par la surface stellaire ou par une instabilité orbitale (un ISCO externe). Cette frontière dicte la fréquence Keplerienne maximale et l'énergie de liaison orbitale disponible pour être libérée. Cependant, les prédictions concernant la masse et le rayon auxquels cette transition se produit sont actuellement limitées par les incertitudes liées à l'EOS, particulièrement la dépendance en densité de l'énergie de symétrie. Bien que les observations astrophysiques contraignent le comportement à haute densité, il est nécessaire de quantifier comment les mesures nucléaires proches de la densité de saturation peuvent affiner ces prédictions pour les configurations en rotation.

Méthodologie
Les auteurs modélisent des étoiles à neutrons en rotation uniforme dans le cadre de la relativité générale en utilisant le code RNS. Ils emploient une famille d'équations d'état dérivées de fonctionnelles de densité d'énergie relativistes avec des interactions de couplage ponctuel dépendant de la densité (DD–PC). Cette famille de fonctionnelles fait varier systématiquement l'énergie de symétrie à la densité de saturation ($J$) de 29 à 36 MeV, avec des paramètres de pente ($L$) correspondants allant de 29 à 94,1 MeV. Toutes les fonctionnelles sont calibrées sur des données de l'état fondamental nucléaire et supportent des masses maximales non rotatives d'au moins $2 M_\odot$.

Pour contraindre cette famille, les auteurs utilisent des mesures de la polarisabilité dipolaire électrique ($\alpha_D$) dans des noyaux riches en neutrons : $^{48}$Ca, $^{68}$Ni, $^{120}$Sn et $^{208}$Pb (collectivement l'ensemble CNSP-4). Des travaux antérieurs (Koliogiannis et al. 2026a) ont déduit une énergie de symétrie de $J = 31,34 \pm 0,79$ MeV à partir de ces données. L'étude compare la famille complète DD–PC par rapport au sous-ensemble cohérent avec les contraintes CNSP-4. Pour chaque configuration, les auteurs localisent l'ISCO prograde par rapport à la surface stellaire afin de définir la masse de transition $M_{tr}(f)$ et le rayon équatorial $R_{eq}$. Ils calculent ensuite la fréquence orbitale ($\nu_{in}$) et l'efficacité de liaison ($\eta_{in}$) au bord interne du disque.

Résultats clés

1. Réduction de l'incertitude : L'application de la contrainte de polarisabilité dipolaire réduit considérablement l'écart des propriétés de transition prédites. À travers les fréquences de spin considérées, la contrainte réduit l'incertitude sur la masse de transition d'un facteur d'environ sept à dix (passant d'écarts de 0,12–0,2 $M_\odot$ à moins de 0,025 $M_\odot$) et sur le rayon de transition par des facteurs similaires (de 0,5–1 km à moins de 0,15 km). Cela représente une réduction de 85 à 90 % de l'écart dépendant de l'EOS.
2. Contraintes sur la fréquence orbitale : Le long de la transition surface–ISCO, la contrainte CNSP-4 confine la fréquence orbitale du bord interne dans une bande plus étroite que 20 Hz à chaque fréquence de spin. En revanche, la famille complète DD–PC produit des largeurs d'environ 38–134 Hz.
3. Cas spécifique (716 Hz) : Pour l'étoile à neutrons la plus rapide connue (PSR J1748−2446ad, $f=716$ Hz), la plage de masse de transition se contracte de $1,586–1,731 M_\odot$ (famille complète) à $1,620–1,635 M_\odot$ (contrainte CNSP-4). La plage du rayon équatorial se contracte de $12,58–13,46$ km à $12,70–12,80$ km.
4. Application à 4U 0614+09 : En combinant la mesure de spin ($f \approx 414$ Hz) et une limite inférieure sur la fréquence la plus élevée des QPO kHz ($\nu_{in} \ge 1267$ Hz) pour le système binaire X de faible masse 4U 0614+09, les auteurs dérivent les intervalles de masse et de rayon autorisés. Sans la contrainte nucléaire, la limite de masse inférieure varie de 0,81 à 1,23 $M_\odot$ ; avec la contrainte CNSP-4, elle se resserre à 0,91–0,98 $M_\odot$ (une réduction de l'écart de 82 %). La limite de masse supérieure reste presque indépendante de l'EOS à $\approx 1,97 M_\odot$. La contrainte nucléaire resserre les limites du rayon équatorial d'environ 82–95 %.

Signification et affirmations
L'article démontre que la réponse dipolaire électrique des noyaux finis constitue une contrainte puissante et indépendante sur la structure des étoiles à neutrons en rotation et sur les propriétés orbitales de leurs disques d'accrétion internes. En reliant les observables nucléaires à la dynamique du disque d'accrétion interne, les auteurs montrent que les mesures de polarisabilité dipolaire permettent de :

• Définir précisément la masse et le rayon auxquels se produit la transition surface–ISCO.
• Confiner la fréquence orbitale prédite au bord du disque interne dans une fenêtre étroite, facilitant l'interprétation des oscillations quasi-périodiques kilohertz (QPO kHz).
• Permettre une estimation de masse indépendante des étoiles à neutrons accrétantes lorsqu'elles sont combinées avec les mesures de fréquence de spin et de données QPO.

Les auteurs concluent que ces contraintes nucléaires informent à la fois la structure stellaire et la mécanique orbitale du disque interne, offrant une voie pour tester la contrainte de l'énergie de symétrie par des observations astrophysiques d'étoiles à neutrons en rotation. Ils notent que des mesures plus précises des réponses dipolaires et des QPO à plus haute fréquence (par exemple, via de futures missions de rayons X) pourraient affiner davantage ces prédictions.