The FUor Mass Distribution Matches the Solar Neighborhood IMF: Evidence for a Universal Eruptive Phase

1. El Problema

La formación estelar implica una fase de acreción de masa, pero la línea temporal de este proceso no está bien comprendida. Existe una tensión observacional entre las tasas de acreción medidas en objetos estelares jóvenes (YSOs) de clase I y II (típicamente $10^{-8} - 10^{-7} M_{\odot} \text{yr}^{-1}$) y las masas finales esperadas en la secuencia principal. Si la acreción fuera constante, los discos durarían demasiado tiempo o no acumularían suficiente masa.

Los objetos FU Ori (FUors) representan la solución a esta tensión mediante eventos de acreción eruptiva extrema, donde las tasas aumentan a $10^{-5} - 10^{-4} M_{\odot} \text{yr}^{-1}$ durante décadas. Sin embargo, durante estas erupciones, el disco interno es ~100 veces más brillante que la estrella central, haciendo imposible observar la estrella directamente. Esto ha dificultado la caracterización de la población subyacente de FUors y la determinación de si todas las estrellas jóvenes pasan por esta fase eruptiva.

El desafío principal es inferir la masa estelar ($M_*$) sin ver la estrella, utilizando únicamente las líneas de absorción del disco, las cuales están afectadas por la rotación kepleriana, la inclinación del disco ($i$) y el radio de truncación interno ($R_{inner}$).

2. Metodología

Los autores emplearon un enfoque estadístico basado en la rotación diferencial del disco para superar las limitaciones de las mediciones individuales:

• Muestra de Datos: Se obtuvieron espectros de alta resolución ($R \approx 15,000 - 25,000$) en el infrarrojo cercano (NIR, 1.0–2.4 µm) utilizando el instrumento NIRSPEC en el telescopio Keck para 23 FUors. Se complementó con datos archivados del espectrógrafo IGRINS ($R \approx 45,000$) para 8 objetos adicionales, totalizando una muestra de 28 FUors.
• Análisis de Rotación Kepleriana:
  • Se asume que el gas del disco interno rota keplerianamente ($v \propto r^{-1/2}$).
  • Se utilizó la técnica de función de correlación cruzada (CCF) en la banda H (1.5–1.8 µm), considerada óptima por su mezcla de características atómicas y moleculares y menor contaminación por vientos estelares o polvo caliente en comparación con otras bandas.
  • Se compararon los espectros observados con una cuadrícula de modelos de disco de acreción (basados en FU Ori) con diferentes anchos de línea rotacional ($v_{max}$).
  • Se midió el ensanchamiento diferencial de las líneas espectrales entre las bandas Y/J (más azules, origen en radios más pequeños y calientes) y la banda K (más roja, origen en radios más grandes y fríos) para confirmar la rotación kepleriana.
• Simulación de Distribuciones de Masa:
  • Se simularon distribuciones de $v_{max}$ asumiendo diferentes Funciones de Masa Inicial (IMF): una "plana" (uniforme) y la IMF del vecindario solar (K24, basada en Kirkpatrick et al. 2024), tanto con como sin enanas marrones (BD).
  • Se asignaron radios estelares basados en isocronas (PARSEC y Baraffe) e inclinaciones aleatorias para generar distribuciones simuladas de $v_{max}$ y compararlas con la distribución observada.

3. Contribuciones Clave

• Validación de la Rotación Kepleriana: Demostración observacional inequívoca de la rotación kepleriana en los discos internos de 20 FUors mediante el ensanchamiento diferencial de líneas en el NIR, confirmando que el modelo de disco de acreción caliente es válido a nivel poblacional.
• Técnica de Medición de $v_{max}$: Desarrollo y aplicación de un método robusto basado en CCF en la banda H para medir la velocidad máxima proyectada ($v_{max} = \sqrt{GM_*/R_{inner}} \sin i$) en una población de objetos, minimizando la degeneración entre masa, radio e inclinación.
• Análisis Poblacional: El primer estudio que no intenta derivar una masa única para cada FUor, sino que compara la distribución estadística de las velocidades medidas contra distribuciones de masa teóricas.

4. Resultados

• Distribución Bimodal de $v_{max}$: La distribución observada de velocidades proyectadas muestra dos picos distintivos: uno alrededor de $\sim30 \text{ km s}^{-1}$ y otro alrededor de $\sim120 \text{ km s}^{-1}$.
• Coincidencia con la IMF Solar:
  • La distribución observada coincide notablemente bien con la simulación basada en la IMF del vecindario solar (K24) restringida a masas estelares de 0.1 a 2 $M_{\odot}$.
  • Una distribución de masa "plana" (uniforme) fue rechazada, ya que sobreestimaba los objetos de alta velocidad.
  • La inclusión de enanas marrones (BD) en la simulación (K24 completa) generó un exceso de objetos de muy baja velocidad ($<30 \text{ km s}^{-1}$) que no se observa en la muestra actual, sugiriendo un sesgo de detección hacia objetos más luminosos o la ausencia real de FUors de masa de enana marrón en la muestra actual.
• Radio Estelar: Las simulaciones indicaron que asumir que las estrellas FUor están significativamente "hinchadas" (radios un 20% mayores) no altera sustancialmente la distribución de $v_{max}$ predicha, lo que sugiere que la relación masa-radio estándar es suficiente para explicar los datos.

5. Significado

• Fase Eruptiva Universal: El hallazgo de que la distribución de masas de los FUors coincide con la IMF del vecindario solar sugiere fuertemente que todas las estrellas jóvenes (o al menos la gran mayoría en el rango de masas estelares) pasan por una fase eruptiva tipo FU Ori durante su evolución pre-secuencia principal. Esto resuelve la tensión entre las tasas de acreción observadas en YSOs tranquilos y la masa final de las estrellas.
• Sesgo de Detección: La falta de FUors de masa de enana marrón en la muestra podría deberse a que sus erupciones son demasiado tenues para ser detectadas con los sondeos actuales, no a que no existan.
• Herramienta Futura: La técnica desarrollada permite medir $v_{max}$ con espectros de banda H de moderada relación señal-ruido, lo que facilitará la caracterización de la población completa de FUors (~40 objetos conocidos) y mejorará los ajustes de la distribución de energía espectral (SED) para obtener parámetros físicos precisos de objetos individuales.

En conclusión, el estudio proporciona la primera evidencia observacional sólida de que los eventos eruptivos FU Ori son una fase evolutiva común y necesaria en la formación estelar, con una distribución de masas que refleja la población estelar general.