Modelling dust coagulation, dynamical drag and turbulent mixing during star and disc formation

Resumen Técnico: Modelado de la Coagulación de Polvo, Arrastre Dinámico y Mezcla Turbulenta durante la Formación de Estrellas y Discos

Planteamiento del Problema
El polvo desempeña un papel crítico tanto en la formación de estrellas como de planetas, influyendo en el enfriamiento del gas, los niveles de ionización y el acoplamiento de los campos magnéticos al gas. Si bien la coagulación del polvo en granos más grandes (guijarros) es esencial para la formación de planetesimales y el crecimiento de núcleos, las simulaciones hidrodinámicas tridimensionales previas de la formación estelar han omitido en gran medida la evolución dinámica del polvo con respecto al gas. Los métodos anteriores, como los de Bate (2022), trataban las poblaciones de polvo como entidades que evolucionaban de forma independiente dentro de cada elemento de fluido (partícula SPH), asumiendo que el polvo estaba perfectamente acoplado al gas. Esta aproximación fue válida para las fases tempranas y de baja densidad del colapso, pero no logró dar cuenta de la migración, el asentamiento y la mezcla del polvo en escalas de tiempo más largas, particularmente dentro de los discos protoplanetarios donde las densidades son más altas y el crecimiento de los granos es rápido. Además, los modelos existentes a menudo carecían de un tratamiento consistente de cómo la turbulencia impulsa tanto las velocidades relativas necesarias para el crecimiento de los granos como la difusión turbulenta que mezcla el polvo verticalmente.

Metodología
Los autores presentan una nueva implementación dentro del código de hidrodinámica de partículas suavizadas (SPH) SPHNG que combina tres procesos físicos distintos:

1. Coagulación de Polvo: Utilizando el método de la ecuación de Smoluchowski de campo medio discretizada (Hirashita & Yan 2009; Bate 2022), el código rastrea la evolución de una población de granos de polvo a través de $N$ contenedores de tamaño espaciados logarítmicamente para cada partícula de gas SPH. El método considera colisiones binarias con una probabilidad de adhesión de unidad por debajo de un umbral de velocidad de coagulación específico.
2. Arrastre Dinámico (MULTIGRAIN): Los autores integran el algoritmo MULTIGRAIN (Hutchison et al. 2018), que resuelve las ecuaciones de fluido continuo para una mezcla de gas y múltiples especies de polvo bajo la aproximación de velocidad terminal. Esto permite calcular las velocidades diferenciales globales entre los diferentes tamaños de polvo y el gas debido a gradientes de presión y gravedad.
3. Integración Implícita y Difusión Turbulenta: Una innovación clave es el uso de un método de integración implícita recientemente desarrollado (Elsender & Bate 2024) para resolver las ecuaciones de arrastre del polvo. Esto evita las restricciones de paso de tiempo restrictivas de los métodos explícitos para granos con tiempos de parada largos y previene fracciones de polvo negativas no físicas. Además, los autores introducen un nuevo método implícito para modelar la difusión turbulenta. Tratan el agitación turbulenta del polvo como un proceso de difusión ($D_i = \alpha_{SS} c_s H / Sc$), donde el número de Schmidt ($Sc$) contabiliza la eficiencia de acoplamiento de diferentes tamaños de grano. Este término de difusión se resuelve utilizando un enfoque de división de operadores de Gauss-Seidel, lo que permite que el polvo se mezcle entre elementos de fluido.

El código rastrea la evolución del polvo desde el colapso de un núcleo de nube molecular, pasando por la fase del primer núcleo hidrostático, hasta la formación y evolución de un disco protoplanetario, utilizando partículas sumidero (sink particles) para seguir el sistema más allá de la formación del núcleo estelar.

Contribuciones Clave

• Primera Simulación Acoplada 3D: Este trabajo representa la primera simulación hidrodinámica tridimensional que modela simultáneamente la coagulación de polvo, el arrastre dinámico del gas y la mezcla turbulenta durante las etapas tempranas de la formación de estrellas y sus discos.
• Esquema de Difusión Implícita: El desarrollo de un robusto resolvedor implícito para la difusión turbulenta del polvo permite simulaciones estables de la mezcla de polvo sin las severas limitaciones de paso de tiempo asociadas con los esquemas de difusión explícitos.
• Tratamiento Unificado: El código unifica el tratamiento de las velocidades relativas (que impulsan el crecimiento) y la difusión turbulenta (que impulsa la mezcla) dentro de un solo marco, utilizando el mismo parámetro de turbulencia ($\alpha_{SS}$) para ambos procesos en los casos de prueba.

Resultados
Los autores presentan resultados tanto de cálculos de prueba (asentamiento de polvo en una columna vertical) como de simulaciones completas de formación estelar:

• Cálculos de Prueba: El código recupera con éxito las soluciones analíticas para el equilibrio entre el asentamiento de polvo y la difusión turbulenta. Las pruebas confirman que la inclusión de la difusión turbulenta evita la formación de capas de polvo irrealmente delgadas y mantiene una escala de altura de polvo finita.
• Simulaciones de Formación Estelar:
  • Primer Núcleo Hidrostático: Consistente con trabajos previos (Bate 2022), ocurre un crecimiento rápido de polvo dentro del primer núcleo hidrostático, produciendo poblaciones localmente monodispersas con tamaños que superan los 100 $\mu$m en el interior del núcleo. La inclusión de arrastre y difusión tuvo un impacto insignificante en esta fase temprana debido a los cortos tiempos evolutivos.
  • Evolución del Disco Protoplanetario: Una vez que se forma el disco, la inclusión de la difusión turbulenta altera significativamente la evolución del polvo.
    • Crecimiento Mejorado: La difusión turbulenta potencia sustancialmente la tasa de crecimiento de los granos de polvo. Al suministrar continuamente granos pequeños e intermedios a las regiones de alta densidad (donde los granos grandes están creciendo), la difusión evita el agotamiento local del "alimento" para la coagulación.
    • Dependencia de la Turbulencia: El nivel de turbulencia ($\alpha_{SS}$) es crítico. Una turbulencia alta ($\alpha_{SS} = 10^{-3}$) conduce a un crecimiento rápido hacia poblaciones monodispersas de $\sim30-100$ $\mu$m en todo el disco. Una turbulencia menor ($\alpha_{SS} = 10^{-4}$) resulta en un crecimiento más lento y una población de polvo más estratificada verticalmente, con granos más grandes confinados cerca del plano medio.
    • Asentamiento de Polvo: En simulaciones con turbulencia débil, el asentamiento de polvo conduce a un aumento significativo en la relación polvo-gas en el plano medio (en más del 50%), mientras que los granos pequeños permanecen empobrecidos cerca del plano medio y confinados en las capas superficiales.
  • Arrastre vs. Difusión: Las comparaciones de simulaciones con y sin procesos específicos revelan que el arrastre dinámico polvo-gas tiene poco efecto en las distribuciones de polvo en estas escalas de tiempo, ya que los granos permanecen bien acoplados. En contraste, la inclusión de la difusión turbulenta tiene un efecto enorme, suavizando las distribuciones de polvo y acelerando la coagulación.

Significancia y Reivindicaciones
El artículo afirma que esta nueva metodología permite, por primera vez, el modelado autosistente del crecimiento, la mezcla y la migración del polvo durante la formación de estrellas y sus discos. Los autores enfatizan que, si bien la inclusión de arrastre y difusión no altera los hallazgos fundamentales respecto al crecimiento de polvo en el primer núcleo hidrostático, es crucial para la evolución posterior del disco protoplanetario. Específicamente, afirman que la difusión de polvo turbulenta proporciona un aumento sustancial a la tasa de crecimiento de los granos de polvo al actuar como una fuente de granos pequeños para las regiones donde se están formando granos grandes.

Los autores mantienen la modestia respecto al alcance de sus hallazgos, señalando que:

• Los resultados dependen del modelo de turbulencia adoptado (por ejemplo, disco-$\alpha$ vs. turbulencia MHD).
• Las simulaciones consideran actualmente la coagulación pero no la fragmentación; diferentes umbrales de coagulación podrían conducir a la fragmentación y repoblación de granos pequeños.
• Los cálculos de transferencia radiativa todavía asumen opacidades interestelares estándar en lugar de utilizar las distribuciones locales de tamaño de grano en evolución, un paso identificado como una extensión futura necesaria.

En última instancia, este trabajo proporciona una herramienta para comparar las simulaciones hidrodinámicas 3D con la riqueza de nuevos datos observacionales de instalaciones como ALMA, VLT-SPHERE y JWST, particularmente en lo que respecta al asentamiento de polvo y la agitación turbulenta en discos vistos de canto.