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

Resumen Técnico: Efectos en las Determinaciones de $v \sin i$ de Estrellas O a partir de Modelos Atmosféricos 3D con Altas Velocidades Turbulentas

Planteamiento del Problema
Determinar la velocidad de rotación proyectada ($v \sin i$) de estrellas O masivas es fundamental para comprender sus ciclos de vida y restringir los modelos evolutivos. Sin embargo, esta determinación enfrenta un desafío principal: la degeneración entre el ensanchamiento rotacional de las líneas y otros mecanismos de ensanchamiento, específicamente la macroturbulencia. Las observaciones de líneas de absorción de estrellas O revelan dispersiones de velocidad ($\sigma_V \sim 50$ km/s) que superan con creces los movimientos térmicos, las cuales no pueden explicarse únicamente por la rotación. Si bien técnicas empíricas como el método de Transformada de Fourier (FT) y el método de Bondad de Ajuste (GOF) son estándar para separar estos componentes, trabajos anteriores (por ejemplo, Sundqvist et al. 2013b) sugieren que estos métodos se vuelven poco fiables cuando $v \sin i \lesssim 50$ km/s debido a la superposición de los componentes de ensanchamiento. Simulaciones recientes de hidrodinámica de radiación 3D (RHD) han demostrado que altas velocidades turbulentas reproducen naturalmente este "ensanchamiento extra", pero estos modelos carecían previamente de rotación, dejando el problema de la separación sin abordar en un contexto 3D realista.

Metodología
Los autores introducen rotación en modelos atmosféricos 3D RHD existentes de estrellas O para generar perfiles de línea sintéticos que incluyen tanto la rotación como las velocidades turbulentas intrínsecas.

• Configuración del Modelo: Utilizando un modelo 3D (Delbroek et al. 2025) con una rotación fija en el límite interior ($v_{\text{rot,in}}$), los autores conservan el momento angular azimutal para simular la rotación. Se centran en la línea de absorción O III 5594 Å.
• Observaciones Simuladas: Los autores generan perfiles de línea sintéticos para varios ángulos de inclinación del observador ($i$). Estos perfiles sirven como "observaciones simuladas" donde la verdadera $v \sin i$ y la velocidad macroturbulenta ($v_{\text{macro}}$) son conocidas a partir de la simulación.
• Análisis de Recuperación: Los autores aplican métodos empíricos 1D estándar a estos datos simulados 3D para probar su fiabilidad:
  1. Transformada de Fourier (FT): Búsqueda del primer cero en el espectro de potencia de Fourier para determinar $v \sin i$.
  2. Bondad de Ajuste (GOF): Ajuste del perfil de línea utilizando una convolución de ensanchamiento macroturbulento rotacional y gaussiano isotrópico para encontrar los parámetros de mejor ajuste.
• Restricciones: El estudio asume que los modelos 3D no tienen ruido (a diferencia de las observaciones reales) y reconoce que los perfiles 3D pueden carecer de la simetría de los modelos 1D, lo que podría afectar el ajuste.

Resultados Clave
El análisis revela limitaciones significativas en los métodos estándar cuando las velocidades turbulentas son altas y desconocidas:

• Fiabilidad del Método FT: El método FT produce resultados razonables solo cuando $v \sin i > v_{\text{macro}}$. Cuando $v \sin i < v_{\text{macro}}$, el método se vuelve altamente poco fiable, a menudo sobreestimando $v \sin i$ por factores de 3 a 10. Crucialmente, el método FT siempre produce un resultado no nulo debido a mínimos fiduciales creados por altas velocidades turbulentas, incluso para vistas polares ($i=0^\circ$) donde la rotación debería ser indetectable.
• Degeneración del Método GOF: El método GOF exhibe una fuerte degeneración entre $v \sin i$ y $v_{\text{macro}}$. Cuando ninguno de los parámetros domina significativamente al otro, los valores de mejor ajuste no son únicos, sino que se distribuyen a través de una gran región del espacio de parámetros.
• La Degeneración de la "Curva de Howarth": La región de valores de mejor ajuste degenerados en el análisis GOF sigue estrechamente la fórmula empírica derivada por Howarth et al. (2007):
  $$C = \sqrt{\frac{v_{\text{macro}}^2}{2} + \frac{(v \sin i)^2}{4}}$$
  El valor inferido ($C_i$) es consistentemente mayor o igual al valor esperado ($C_e$). Esto implica que el ajuste estándar solo puede restringir la suma de estos términos, no los componentes individuales, a menos que uno domine claramente.
• Restricciones de Parámetros:
  • Si $v \sin i \gg v_{\text{macro}}$, $v \sin i$ está bien restringido, pero $v_{\text{macro}}$ no lo está.
  • Si $v_{\text{macro}} \gg v \sin i$, $v_{\text{macro}}$ está bien restringido, pero $v \sin i$ se sobreestima.
  • Si $v \sin i \approx v_{\text{macro}}$, ninguno de los parámetros puede determinarse individualmente de manera fiable.

Significado y Afirmaciones
El artículo concluye que los métodos espectroscópicos estándar (FT y GOF) generalmente no pueden determinar independientemente $v \sin i$ y $v_{\text{macro}}$ en presencia de altas velocidades turbulentas encontradas en atmósferas 3D de estrellas O. En cambio, estos métodos aproximan principalmente la cantidad combinada $\sqrt{v_{\text{macro}}^2/2 + (v \sin i)^2/4}$.

Los autores afirman que este hallazgo requiere un reanálisis y una reinterpretación de las correlaciones empíricas previamente establecidas entre $v \sin i$ y $v_{\text{macro}}$ (por ejemplo, Simón-Díaz et al. 2017) y de las distribuciones estadísticas de las tasas de rotación observadas en poblaciones de estrellas O (por ejemplo, Ramírez-Agudelo et al. 2013). El estudio sugiere que las determinaciones anteriores de las tasas de rotación podrían alterarse significativamente si la velocidad macroturbulenta no se calibró independientemente (por ejemplo, mediante simulaciones 3D) o si no dominó claramente el ensanchamiento. El artículo no propone nuevas técnicas observacionales, sino que destaca la necesidad de reinterpretar los datos existentes a la luz de estos efectos atmosféricos 3D.