FORGE'd in the Early Universe: The Effect of Protostellar Outflows on Pop III Accretion

Resumen Técnico: Simulación de una (proto)estrella de Población III con el marco FORGE'd in FIRE

Planteamiento del Problema
La formación y evolución de la primera generación de estrellas (Población III) permanece insuficientemente delimitada debido a la falta de observaciones directas y al inmenso rango dinámico requerido para conectar la formación de estructuras cosmológicas con la física detallada de los protoestrellas. Los modelos teóricos previos a menudo dependían de condiciones iniciales idealizadas, omitían física crítica (como campos magnéticos, radiación y chorros) o no lograban producir las Funciones Iniciales de Masa (IMF) realistas y las multiplicidades estelares observadas en entornos ricos en metales. Existe un debate específico sobre si las estrellas de la Pop. III pueden lanzar chorros (jets) protoestelares, un proceso que depende de la amplificación de campos magnéticos traza y de la presencia de un disco circunestelar rotatorio y magnetizado.

Metodología
Este estudio presenta una simulación de magnetohidrodinámica de radiación (RMHD) de "zoom-in" cosmológica utilizando el marco FORGE'd in FIRE. La simulación rastrea la formación y evolución de una única (proto)estrella de Pop. III libre de metales a un desplazamiento al rojo $z \sim 14$ dentro de un contexto cosmológico completo.

• Marco de Simulación: Se utiliza el código GIZMO, empleando el método de masa finita sin malla (MFM) para resolver las ecuaciones acopladas de gravedad, magnetohidrodinámica (MHD) y transporte de radiación.
• Física Incluida: La simulación incorpora el modelo de física FIRE-3, que incluye:
  • Radiación: Transporte multifrecuencia (fotoionizante, Lyman-Werner, etc.) acoplado de forma consistente con el gas y el polvo.
  • Termoquímica: Extendida a temperaturas de $1\text{--}10^{10}$ K y densidades de $10^{-8}\text{--}10^{24}$ cm$^{-3}$, rastreando las abundancias de H, He, HD y polvo.
  • Campos Magnéticos: Inicializados con campos semilla cosmológicos traza ($\sim 10^{-15}$ G) que se amplifican de forma consistente. Se omiten los efectos de la MHD no ideal (resistividad óhmica, efecto Hall, difusión ambipolar), justificado por comparaciones de escalas de tiempo en el post-procesamiento que muestran que son insignificantes en comparación con las escalas de tiempo dinámicas en el gas libre de metales.
  • Evolución Protoestelar: Las partículas sumidero (sink particles) representan (proto)estrellas, evolucionando a través de la acreción, las trayectorias de la secuencia principal y los mecanismos de retroalimentación (feedback).
• Modelado de Jets y Vientos: Los jets protoestelares se modelan siguiendo a Cunningham et al., desviando una fracción ($f_W = 0.3$) de la masa acrecionada hacia flujos bipolares colimados lanzados a una fracción ($f_K = 0.3$) de la velocidad Kepleriana. Se descartan los vientos estelares para las estrellas libres de metales, ya que su escala depende de la metalicidad.
• Resolución: La simulación utiliza "hiper-refinamiento", resolviendo celdas de gas hasta $\sim 10^{-3} M_\odot$ y escalas espaciales de $\sim 10^{-4}$ pc ($< 100$ au) cerca de la (proto)estrella, manteniendo al mismo tiempo el contexto cosmológico desde $z \sim 100$.

Resultados Clave

1. Crecimiento y Masa Estelar:

   • La (proto)estrella crece desde una tasa de acreción inicial de $\sim 5 \times 10^{-3} M_\odot \text{yr}^{-1}$ hasta una masa final de $\sim 27 M_\odot$ durante 31,000 años.
   • La tasa de crecimiento se estabiliza a medida que la estrella entra en la fase de "combustión de la capa de Shell D". La masa final es regulada por la retroalimentación (jets y radiación) en lugar por el agotamiento del reservorio circundante.
   • Al final de la simulación, la estrella tiene un radio de $\sim 10 R_\odot$, una temperatura efectiva de $\sim 3.5 \times 10^4$ K y una luminosidad de $\sim 1.3 \times 10^5 L_\odot$. El tiempo de escala de Kelvin-Helmholtz ($t_{KH}$) se vuelve más corto que el tiempo de escala de acreción ($t_{acc}$), lo que indica que la estrella alcanzará la secuencia principal antes de acrecer significativamente más masa.

2. Disco Circunestelar:

   • Se forma un disco soportado por rotación, que se extiende hasta $\sim 0.01$ pc ($\sim 2000$ au).
   • El disco permanece gravitacionalmente estable durante toda la simulación, con un parámetro de Toomre $Q$ consistentemente $>> 1$.
   • La estabilidad se mantiene mediante una combinación de soporte térmico (en las regiones internas, $T \sim 10^4$ K) y soporte turbulento (en las regiones externas, $R > 1000$ au).
   • El perfil de densidad superficial transiciona de $\Sigma \propto R^{-1}$ en el disco exterior a un perfil de disco-$\alpha$ en las regiones internas, con un aumento abrupto de la densidad a $\sim 20$ au.

3. Jets Protoestelares y Flujos (Outflows):

   • Lanzamiento de Jets: Los jets se forman porque la relación masa-flujo magnética ($\mu$) cae dentro del régimen ($\mu \sim 10^{-1} \text{--} 10^1$) que permite el lanzamiento de jets, según lo definido por Machida y Doi. Esto ocurre a pesar del uso de campos magnéticos semilla cosmológicos traza.
   • Dinámica: Los jets se propagan hasta $\sim 1$ pc, colisionando con el gas entrante para formar choques de terminación.
   • Regulación: Los jets impulsan un momento sustancial, barriendo el gas y reduciendo la tasa de acreción neta hacia adentro en aproximadamente un orden de magnitud en radios grandes. Esta retroalimentación impulsada por el momento regula la masa estelar final.
   • La tasa de flujo alcanza su máximo en $\sim 10^{-3} M_\odot \text{yr}^{-1}$ dentro de los primeros 1,000 años antes de declinar, mientras que la tasa de flujo hacia adentro se estabiliza alrededor de $3\text{--}5 \times 10^{-4} M_\odot \text{yr}^{-1}$.

4. Estructura Química y Térmica:

   • En el disco interno, las altas temperaturas ($T \sim 10^4$ K) suprimen el hidrógeno molecular ($H_2$) en favor del hidrógeno ionizado, desplazando los mecanismos de enfriamiento de $H_2$ hacia el enfriamiento de Lyman.
   • El disco permanece ópticamente delgado a su radiación de enfriamiento dominante, permitiendo que los fotones de IR escapen eficientemente.

Significancia y Reivindicaciones
El artículo afirma presentar la primera aplicación del marco FORGE'd in FIRE a una (proto)estrella de Pop. III resuelta individualmente dentro de una configuración de "zoom-in" cosmológica completa. La significancia principal radica en demostrar que:

• Los jets son viables en la formación de la Pop. III: Incluso con campos magnéticos semilla cosmológicos traza, se cumplen las condiciones para un fuerte lanzamiento de jets, desafiando las suposiciones previas de que los jets podrían estar ausentes en entornos libres de metales.
• La retroalimentación regula la masa: Los jets protoestelares juegan un papel crítico en la configuración de la masa final de las estrellas de la Pop. III al regular la acreción mediante flujos impulsados por el momento, en lugar de que la estrella simplemente consuma todo su reservorio.
• Estabilidad del Disco: Es posible la formación de un disco circunestelar estable y no fragmentado, mantenido por soporte térmico y turbulento.

Estos resultados proporcionan una imagen físicamente más completa de las regiones de formación estelar más tempranas, ofreciendo perspectivas sobre la Función Inicial de Masa (IMF) de la Pop. III y los procesos de retroalimentación. Los autores señalan que el trabajo futuro expandirá esto a un grupo de estrellas de la Pop. III para caracterizar completamente la IMF y la multiplicidad estelar predicha.