Effects on vsini determinations of O stars from 3D model atmospheres with high turbulent velocities

Résumé technique : Effets sur les déterminations de $v \sin i$ des étoiles O à partir d'atmosphères modèles 3D avec de fortes vitesses turbulentes

Énoncé du problème
La détermination de la vitesse de rotation projetée ($v \sin i$) des étoiles O massives est fondamentale pour comprendre leurs cycles de vie et contraindre les modèles d'évolution. Cependant, cette détermination fait face à un défi majeur : la dégénérescence entre l'élargissement rotationnel des raies et d'autres mécanismes d'élargissement, spécifiquement la macroturbulence. Les observations des raies d'absorption des étoiles O révèlent des dispersions de vitesses ($\sigma_V \sim 50$ km/s) bien supérieures aux mouvements thermiques, ce qui ne peut s'expliquer par la rotation seule. Bien que des techniques empiriques comme la méthode de Transformée de Fourier (TF) et la méthode d'Adéquation aux Données (GOF) soient standards pour séparer ces composantes, des travaux antérieurs (par exemple, Sundqvist et al. 2013b) suggèrent que ces méthodes deviennent peu fiables lorsque $v \sin i \lesssim 50$ km/s en raison du chevauchement des composantes d'élargissement. Des simulations récentes d'hydrodynamique radiative 3D (RHD) ont montré que de fortes vitesses turbulentes reproduisent naturellement cet « élargissement supplémentaire », mais ces modèles manquaient auparavant de rotation, laissant le problème de la séparation non résolu dans un contexte 3D réaliste.

Méthodologie
Les auteurs introduisent la rotation dans des modèles d'atmosphères 3D RHD existants d'étoiles O pour générer des profils de raies synthétiques incluant à la fois la rotation et les vitesses turbulentes intrinsèques.

• Configuration du modèle : En utilisant un modèle 3D (Delbroek et al. 2025) avec une rotation fixe à la limite intérieure ($v_{\text{rot,in}}$), les auteurs conservent le moment angulaire azimutal pour simuler la rotation. Ils se concentrent sur la raie d'absorption O III 5594 Å.
• Observations simulées : Les auteurs génèrent des profils de raies synthétiques pour divers angles d'inclinaison de l'observateur ($i$). Ces profils servent d'« observations simulées » où la vraie $v \sin i$ et la vitesse macroturbulente ($v_{\text{macro}}$) sont connues grâce à la simulation.
• Analyse de récupération : Les auteurs appliquent des méthodes empiriques 1D standard à ces données simulées 3D pour tester leur fiabilité :
  1. Transformée de Fourier (TF) : Recherche du premier zéro dans le spectre de puissance de Fourier pour déterminer $v \sin i$.
  2. Adéquation aux Données (GOF) : Ajustement du profil de raie en utilisant une convolution d'élargissement rotationnel et d'élargissement macroturbulent isotrope gaussien pour trouver les paramètres d'ajustement optimal.
• Contraintes : L'étude suppose que les modèles 3D ne comportent pas de bruit (contrairement aux observations réelles) et reconnaît que les profils 3D peuvent manquer de la symétrie des modèles 1D, affectant potentiellement l'ajustement.

Résultats clés
L'analyse révèle des limitations significatives des méthodes standard lorsque les vitesses turbulentes sont élevées et inconnues :

• Fiabilité de la méthode TF : La méthode TF ne donne des résultats raisonnables que lorsque $v \sin i > v_{\text{macro}}$. Lorsque $v \sin i < v_{\text{macro}}$, la méthode devient hautement peu fiable, surestimant souvent $v \sin i$ par des facteurs de 3 à 10. Crucialement, la méthode TF produit toujours un résultat non nul en raison de minima de référence créés par les fortes vitesses turbulentes, même pour des vues polaires ($i=0^\circ$) où la rotation devrait être indétectable.
• Dégénérescence de la méthode GOF : La méthode GOF présente une forte dégénérescence entre $v \sin i$ et $v_{\text{macro}}$. Lorsque ni l'un ni l'autre paramètre ne domine significativement l'autre, les valeurs d'ajustement optimal ne sont pas uniques mais sont distribuées sur une vaste région de l'espace des paramètres.
• La dégénérescence de la « courbe de Howarth » : La région des valeurs d'ajustement optimal dégénérées dans l'analyse GOF suit étroitement la formule empirique dérivée par Howarth et al. (2007) :
  $$C = \sqrt{\frac{v_{\text{macro}}^2}{2} + \frac{(v \sin i)^2}{4}}$$
  La valeur déduite ($C_i$) est systématiquement supérieure ou égale à la valeur attendue ($C_e$). Cela implique que l'ajustement standard ne peut contraindre que la somme de ces termes, et non les composantes individuelles, à moins que l'une ne domine clairement.
• Contraintes de paramètres :
  • Si $v \sin i \gg v_{\text{macro}}$, $v \sin i$ est bien contraint, mais $v_{\text{macro}}$ ne l'est pas.
  • Si $v_{\text{macro}} \gg v \sin i$, $v_{\text{macro}}$ est bien contraint, mais $v \sin i$ est surestimé.
  • Si $v \sin i \approx v_{\text{macro}}$, aucun des deux paramètres ne peut être déterminé de manière fiable individuellement.

Signification et affirmations
L'article conclut que les méthodes spectroscopiques standard (TF et GOF) ne peuvent généralement pas déterminer indépendamment $v \sin i$ et $v_{\text{macro}}$ en présence des fortes vitesses turbulentes trouvées dans les atmosphères 3D des étoiles O. Au lieu de cela, ces méthodes approximent principalement la quantité combinée $\sqrt{v_{\text{macro}}^2/2 + (v \sin i)^2/4}$.

Les auteurs affirment que cette découverte nécessite une réanalyse et une réinterprétation des corrélations empiriques précédemment établies entre $v \sin i$ et $v_{\text{macro}}$ (par exemple, Simón-Díaz et al. 2017) et des distributions statistiques des taux de rotation observés dans les populations d'étoiles O (par exemple, Ramírez-Agudelo et al. 2013). L'étude suggère que les déterminations antérieures des taux de rotation pourraient être significativement modifiées si la vitesse macroturbulente n'a pas été calibrée indépendamment (par exemple, via des simulations 3D) ou si elle ne dominait pas clairement l'élargissement. L'article ne propose pas de nouvelles techniques observationnelles mais met en évidence la nécessité de réinterpréter les données existantes à la lumière de ces effets d'atmosphère 3D.