A comprehensive study of $δ$ Scuti-type pulsators in eclipsing binaries: oscillating eclipsing Algols

Résumé Technique : Une Étude Approfondie des Pulsateurs de Type $\delta$ Scuti dans les Binaires Éclipsantes

Énoncé du Problème
Les binaires éclipsantes (BE) qui sont également des binaires spectroscopiques à deux composantes (SB2) représentent un laboratoire critique pour la détermination des paramètres stellaires absolus, tels que la masse et le rayon, qui sont essentiels pour tester les modèles d'évolution stellaire. Cependant, les paramètres fondamentaux seuls échouent souvent à révéler la structure interne ou la composition chimique d'une étoile. Bien que les pulsateurs de type $\delta$ Scuti offrent une voie pour sonder l'intérieur des étoiles via l'astérosismologie, leur caractérisation dans les systèmes binaires est entravée par le manque de données de vitesse radiale (VR) de haute qualité. Sans mesures de VR, les masses absolues ne peuvent être déterminées, laissant l'histoire évolutive de ces systèmes incomplète. De plus, un sous-ensemble spécifique de ces systèmes — les Algols éclipsants oscillants (oEA) — implique des configurations semi-détachées avec des rapports de masse extrêmes, indiquant une histoire de transfert de masse qui modifie la durée de vie et la position sur le diagramme de Hertzsprung–Russell (HR). Une compréhension globale de la manière dont le transfert de masse affecte les propriétés de pulsation et l'évolution de ces étoiles de masse intermédiaire nécessite une analyse simultanée de la photométrie, de la spectroscopie et de la modélisation évolutive.

Méthodologie
Les auteurs présentent une analyse multi-facettes de cinq systèmes oEA identifiés à partir du relevé CRÉME (Comprehensive Research with Échelles on the Most Interesting Eclipsing Binaries) : TIC 64437380 (HD 139774), TIC 82474821 (SW Pup), TIC 126945917 (HD 202042), TIC 10756751 (GP Cet), et TIC 37817410 (TZ Eri).

1. Acquisition des Données :

   • Photométrie : Des courbes de lumière (LC) de haute précision à cadence courte (2 minutes) ont été récupérées auprès du satellite Transiting Exoplanet Survey Satellite (TESS). Les données ont été traitées à l'aide des packages lightkurve et wotan pour normaliser le flux et éliminer les tendances à long terme.
   • Spectroscopie : Des spectres multi-époques à haute résolution ont été obtenus à l'aide de divers spectrographes, principalement FEROS (MPG-2.2m) et CHIRON (SMARTS 1.5m), avec des données supplémentaires pour GP Cet provenant de HIDES et IRCS.

2. Dérivation des Paramètres Orbitaux et Stellaires :

   • Vitesses Radiales : Les VR ont été calculées en utilisant la technique de corrélation croisée bidimensionnelle (todcor). Les solutions orbitales ont été dérivées à l'aide de v2fit, en ajustant la période, la vitesse systémique et les demi-amplitudes de vitesse.
   • Analyse Spectrale : Les paramètres atmosphériques ($T_{\text{eff}}$, $\log g$, $v \sin i$) pour les étoiles primaires ont été déterminés à l'aide du cadre iSpec et du module GSSP (Grid Search in Stellar Parameters). Le désenchevêtrement spectral a été réalisé à l'aide du code disentangling_shift_and_add pour isoler les spectres de chaque composante.
   • Modélisation de la Courbe de Lumière : Le code phoebe2 a été employé pour modéliser les LC, en incorporant la géométrie de Roche, l'assombrissement au limb et l'assombrissement gravitationnel. La modélisation a utilisé une contrainte de semi-détaché où l'étoile secondaire remplit son lobe de Roche. Les paramètres tels que l'inclinaison, les rapports de rayon et les rapports de température ont été optimisés à l'aide des algorithmes Nelder-Mead et d'Évolution Différentielle, les incertitudes étant estimées via un échantillonnage par Monte Carlo par chaînes de Markov (MCMC).

3. Modélisation Évolutive :

   • Une grille de 5 000 modèles d'évolution binaire a été construite à l'aide du code MESA (Modules for Experiments in Stellar Astrophysics).
   • La grille a fait varier les périodes orbitales initiales (1,2–5,0 jours), les masses initiales (0,5–2,5 $M_\odot$), la métallicité, l'overshooting convectif et l'efficacité du transfert de masse ($\beta$).
   • Les modèles ont été comparés aux masses et rayons observés en utilisant la distance de Mahalanobis pour identifier les trajectes évolutives de meilleur ajustement, en distinguant entre les scénarios de Cas A (transfert de masse pendant la combustion de l'hydrogène central) et de Cas B (transfert de masse après l'épuisement de l'hydrogène central).

4. Analyse de la Pulsation :

   • L'analyse des fréquences a été effectuée sur les résidus des ajustements du modèle binaire.
   • Les résidus ont été nettoyés des harmoniques orbitales et des tendances à long terme.
   • Le logiciel Period04 a été utilisé pour rechercher les fréquences d'oscillation (0–80 d$^{-1}$) à l'aide d'une procédure de blanchiment (pre-whitening). Les fréquences avec un rapport signal sur bruit (SNR) $\ge$ 5 ont été considérées comme significatives.

Résultats Clés

• Paramètres du Système : Des paramètres orbitaux et stellaires précis ont été dérivés pour les cinq systèmes. Les systèmes présentent de faibles rapports de masse ($q \approx 0,1$–$0,3$), confirmant leur nature d'Algol semi-détaché. Les composantes primaires sont des pulsateurs de type $\delta$ Scuti (masses $\sim$1,7–2,4 $M_\odot$), tandis que les composantes secondaires sont plus froides, sous-lumineuses et gonflées, remplissant leurs lobes de Roche.
• Histoire Évolutive : Les simulations MESA suggèrent que tous les systèmes ont subi un transfert de masse significatif.
  • TIC 82474821 : Les meilleurs ajustements suggèrent une histoire de transfert de masse de Cas B non conservatrice, avec un âge du système estimé entre 0,5 et 3 Gyr.
  • TIC 10756751 : Les modèles indiquent un transfert de masse de Cas A presque conservateur, avec une plage d'âge de 1,86–8,13 Gyr.
  • TIC 37817410 : Les ajustements pointent vers un scénario de Cas B avec des métallicités sur-solaires et des cycles de combustion CNO instables dans l'enveloppe de la donneuse.
  • TIC 64437380 & TIC 126945917 : Bien que les ajustements soient moins précis (écartant jusqu'à 50 %), les résultats soutiennent toujours une classification de Cas B, bien que les auteurs notent que l'affinement des plages de paramètres initiaux (spécifiquement la période orbitale initiale et la masse de la secondaire) pourrait améliorer la qualité de l'ajustement.
• Pulsation : Des fréquences d'oscillation significatives ont été détectées dans les résidus de tous les systèmes, confirmant leur statut de binaires pulsantes. L'analyse a réussi à isoler ces fréquences du signal orbital binaire.

Signification et Revendications
L'article affirme contribuer à l'enrichissement du catalogue des systèmes oEA bien caractérisés, fournissant des points de référence essentiels pour tester les modèles théoriques de structure et d'évolution stellaires. En combinant les contraintes géométriques précises des binaires éclipsantes avec les données astérosismiques, l'étude démontre la faisabilité d'une contrainte simultanée de la géométrie stellaire et de la structure interne.

Les auteurs concluent modestement que, bien que leurs modèles évolutifs reproduisent généralement les systèmes observés, des divergences subsistent. Ils attribuent cela aux incertitudes dans les traitements physiques tels que l'overshooting convectif, l'efficacité du transfert de masse, la perte de moment angulaire et la rotation stellaire. Par conséquent, le papier se positionne comme une étape fondamentale plutôt que comme une solution définitive, notant que des travaux futurs impliquant la photométrie multi-couleurs pour l'identification des modes et la modélisation séismique sont nécessaires pour contraindre davantage les scénarios de formation et comprendre les changements de structure interne induits par le transfert de masse. L'étude souligne le rôle critique des interactions binaires dans le façonnement des trajectes évolutives des étoiles de masse intermédiaire.