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

Technische Zusammenfassung: Auswirkungen von 3D-Modellatmosphären mit hohen turbulenten Geschwindigkeiten auf die Bestimmung von $v \sin i$ bei O-Sternen

Problemstellung
Die Bestimmung der projizierten Rotationsgeschwindigkeit ($v \sin i$) massereicher O-Sterne ist grundlegend für das Verständnis ihrer Lebenszyklen und zur Einschränkung von Evolutionsmodellen. Diese Bestimmung steht jedoch vor einer primären Herausforderung: der Entartung zwischen rotationsbedingter Linienverbreiterung und anderen Verbreiterungsmechanismen, insbesondere der Makroturbulenz. Beobachtungen von Absorptionslinien bei O-Sternen zeigen Geschwindigkeitsdispersionen ($\sigma_V \sim 50$ km/s), die weit über thermischen Bewegungen liegen und nicht allein durch Rotation erklärt werden können. Während empirische Techniken wie die Fourier-Transformations-Methode (FT) und die Goodness-of-Fit-Methode (GOF) Standardverfahren zur Trennung dieser Komponenten sind, deuten frühere Arbeiten (z. B. Sundqvist et al. 2013b) darauf hin, dass diese Methoden unzuverlässig werden, wenn $v \sin i \lesssim 50$ km/s beträgt, aufgrund der Überlappung der Verbreiterungskomponenten. Jüngste 3D-Radiativ-Hydrodynamik-Simulationen (RHD) haben gezeigt, dass hohe turbulente Geschwindigkeiten diese „zusätzliche Verbreiterung" natürlich reproduzieren, doch diesen Modellen fehlte bisher die Rotation, sodass das Trennungsproblem in einem realistischen 3D-Kontext unbeantwortet blieb.

Methodik
Die Autoren führen Rotation in bestehende 3D-RHD-Modellatmosphären von O-Sternen ein, um synthetische Linienprofile zu erzeugen, die sowohl Rotation als auch intrinsische turbulente Geschwindigkeiten enthalten.

• Modellaufbau: Unter Verwendung eines 3D-Modells (Delbroek et al. 2025) mit einer festen inneren Randrotation ($v_{\text{rot,in}}$) conservieren die Autoren das azimutale Drehmoment, um Rotation zu simulieren. Der Fokus liegt auf der Absorptionslinie O III 5594 Å.
• Mock-Beobachtungen: Die Autoren generieren synthetische Linienprofile für verschiedene Neigungswinkel des Beobachters ($i$). Diese Profile dienen als „Mock-Beobachtungen", bei denen die wahre $v \sin i$ und die makroturbulente Geschwindigkeit ($v_{\text{macro}}$) aus der Simulation bekannt sind.
• Rückgewinnungsanalyse: Die Autoren wenden Standard-empirische 1D-Methoden auf diese 3D-Mock-Daten an, um deren Zuverlässigkeit zu testen:
  1. Fourier-Transformation (FT): Suche nach der ersten Nullstelle im Fourier-Leistungsspektrum zur Bestimmung von $v \sin i$.
  2. Goodness-of-Fit (GOF): Anpassung des Linienprofils mittels einer Faltung aus rotationsbedingter und isotroper gaußscher makroturbulenter Verbreiterung, um die besten Anpassungswerte zu finden.
• Einschränkungen: Die Studie geht davon aus, dass die 3D-Modelle kein Rauschen enthalten (im Gegensatz zu realen Beobachtungen) und erkennt an, dass 3D-Profile möglicherweise die Symmetrie von 1D-Modellen nicht aufweisen, was die Anpassung beeinflussen könnte.

Hauptergebnisse
Die Analyse zeigt signifikante Einschränkungen bei Standardmethoden, wenn turbulente Geschwindigkeiten hoch und unbekannt sind:

• Zuverlässigkeit der FT-Methode: Die FT-Methode liefert nur dann vernünftige Ergebnisse, wenn $v \sin i > v_{\text{macro}}$ gilt. Wenn $v \sin i < v_{\text{macro}}$ ist, wird die Methode höchst unzuverlässig und überschätzt $v \sin i$ häufig um Faktoren von 3 bis 10. Entscheidend ist, dass die FT-Methode aufgrund von Referenzminima, die durch hohe turbulente Geschwindigkeiten erzeugt werden, immer ein von Null verschiedenes Ergebnis liefert, selbst für polare Ansichten ($i=0^\circ$), bei denen Rotation nicht nachweisbar sein sollte.
• Entartung der GOF-Methode: Die GOF-Methode zeigt eine starke Entartung zwischen $v \sin i$ und $v_{\text{macro}}$. Wenn keiner der Parameter den anderen signifikant dominiert, sind die besten Anpassungswerte nicht eindeutig, sondern über einen großen Bereich des Parameterraums verteilt.
• Die „Howarth-Kurve"-Entartung: Der Bereich entarteter bester Anpassungswerte in der GOF-Analyse folgt eng der empirischen Formel, die von Howarth et al. (2007) abgeleitet wurde:
  $$C = \sqrt{\frac{v_{\text{macro}}^2}{2} + \frac{(v \sin i)^2}{4}}$$
  Der abgeleitete Wert ($C_i$) ist konsistent größer oder gleich dem erwarteten Wert ($C_e$). Dies impliziert, dass Standardanpassungen nur die Summe dieser Terme einschränken können, nicht jedoch die einzelnen Komponenten, es sei denn, eine dominiert klar.
• Parametereinschränkungen:
  • Wenn $v \sin i \gg v_{\text{macro}}$, ist $v \sin i$ gut eingeschränkt, aber $v_{\text{macro}}$ nicht.
  • Wenn $v_{\text{macro}} \gg v \sin i$, ist $v_{\text{macro}}$ gut eingeschränkt, aber $v \sin i$ wird überschätzt.
  • Wenn $v \sin i \approx v_{\text{macro}}$, kann keiner der Parameter einzeln zuverlässig bestimmt werden.

Bedeutung und Behauptungen
Die Studie kommt zu dem Schluss, dass Standard-spektroskopische Methoden (FT und GOF) im Allgemeinen $v \sin i$ und $v_{\text{macro}}$ nicht unabhängig voneinander bestimmen können, wenn hohe turbulente Geschwindigkeiten in 3D-O-Sternatmosphären vorhanden sind. Stattdessen approximieren diese Methoden primär die kombinierte Größe $\sqrt{v_{\text{macro}}^2/2 + (v \sin i)^2/4}$.

Die Autoren behaupten, dass diese Erkenntnis eine Neubearbeitung und Neuinterpretation zuvor etablierter empirischer Korrelationen zwischen $v \sin i$ und $v_{\text{macro}}$ (z. B. Simón-Díaz et al. 2017) sowie der statistischen Verteilungen beobachteter Rotationsraten in O-Sternpopulationen (z. B. Ramírez-Agudelo et al. 2013) erforderlich macht. Die Studie legt nahe, dass frühere Bestimmungen von Rotationsraten erheblich verändert sein könnten, wenn die makroturbulente Geschwindigkeit nicht unabhängig kalibriert wurde (z. B. via 3D-Simulationen) oder wenn sie die Verbreiterung nicht klar dominierte. Die Studie schlägt keine neuen Beobachtungstechniken vor, sondern unterstreicht die Notwendigkeit, bestehende Daten im Lichte dieser 3D-Atmosphäreneffekte neu zu interpretieren.