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

Resumo Técnico: Modelagem de Coagulação de Poeira, Arrasto Dinâmico e Mistura Turbulenta durante a Formação de Estrelas e Discos

Definição do Problema
A poeira desempenha um papel crítico tanto na formação de estrelas quanto de planetas, influenciando o resfriamento do gás, os níveis de ionização e o acoplamento dos campos magnéticos ao gás. Embora a coagulação da poeira em grãos maiores (seixos) seja essencial para a formação de planetesimais e o crescimento de núcleos, simulações hidrodinâmicas tridimensionais anteriores de formação estelar negligenciaram amplamente a evolução dinâmica da poeira em relação ao gás. Métodos anteriores, como os de Bate (2022), tratavam as populações de poeira como evoluindo independentemente dentro de cada elemento de fluido (partícula SPH), assumindo que a poeira estava perfeitamente acoplada ao gás. Essa aproximação era válida para as fases iniciais de colapso, de baixa densidade, mas falhou em considerar a migração, o assentamento e a mistura da poeira ao longo de escalas de tempo mais longas, particularmente em discos protoplanetários onde as densidades são mais altas e o crescimento de grãos é rápido. Além disso, modelos existentes frequentemente careciam de um tratamento autoconsistente de como a turbulência impulsiona tanto as velocidades relativas necessárias para o crescimento de grãos quanto a difusão turbulenta que mistura a poeira verticalmente.

Metodologia
Os autores apresentam uma nova implementação no código de hidrodinâmica de partículas suavizadas (SPH) SPHNG, que combina três processos físicos distintos:

1. Coagulação de Poeira: Utilizando o método da equação de Smoluchowski de campo médio discretizado (Hirashita & Yan 2009; Bate 2022), o código rastreia a evolução de uma população de grãos de poeira através de $N$ classes de tamanho espaçadas logaritmicamente para cada partícula de gás SPH. O método considera colisões binárias com uma probabilidade de adesão unitária abaixo de uma velocidade de limiar de coagulação específica.
2. Arrasto Dinâmico (MULTIGRAIN): Os autores integram o algoritmo MULTIGRAIN (Hutchison et al. 2018), que resolve as equações de fluido contínuo para uma mistura de gás e múltiplas espécies de poeira na aproximação de velocidade terminal. Isso permite o cálculo de velocidades diferenciais globais entre os tamanhos de poeira e o gás devido a gradientes de pressão e gravidade.
3. Integração Implícita e Difusão Turbulenta: Uma inovação fundamental é o uso de um método de integração implícita recentemente desenvolvido (Elsender & Bate 2024) para resolver as equações de arrasto da poeira. Isso evita as restrições severas de passo de tempo (time-step) de métodos explícitos para grãos com tempos de parada longos e evita frações de poeira negativas não físicas. Adicionalmente, os autores introduzem um novo método implícito para modelar a difusão turbulenta. Eles tratam a agitação turbulenta da poeira como um processo de difusão ($D_i = \alpha_{SS} c_s H / Sc$), onde o número de Schmidt ($Sc$) contabiliza a eficiência de acoplamento de diferentes tamanhos de grãos. Este termo de difusão é resolvido usando uma abordagem de separação de operadores de Gauss-Seidel, permitindo que a poeira se misture entre elementos de fluido.

O código rastreia a evolução da poeira desde o colapso de um núcleo de nuvem molecular, passando pela fase do primeiro núcleo hidrostático, até a formação e evolução de um disco protoplanetário, utilizando partículas de sumidouro (sink particles) para acompanhar o sistema além da formação do núcleo estelar.

Principais Contribuições

• Primeira Simulação Acoplada 3D: Este trabalho representa a primeira simulação hidrodinâmica tridimensional a modelar simultaneamente a coagulação de poeira, o arrasto dinâmico do gás e a mistura turbulenta durante os estágios iniciais de formação de estrelas e seus discos.
• Esquema de Difusão Implícita: O desenvolvimento de um solver implícito robusto para a difusão de poeira turbulenta permite simulações estáveis da mistura de poeira sem as severas limitações de passo de tempo associadas a esquemas de difusão explícitos.
• Tratamento Unificado: O código unifica o tratamento das velocidades relativas (que impulsionam o crescimento) e da difusão turbulenta (que impulsiona a mistura) dentro de um único arcabouço, utilizando o mesmo parâmetro de turbulência ($\alpha_{SS}$) para ambos os processos nos casos de teste.

Resultados
Os autores apresentam resultados de tanto cálculos de teste (assentamento de poeira em uma coluna vertical) quanto de simulações completas de formação estelar:

• Cálculos de Teste: O código recupera com sucesso soluções analíticas para o equilíbrio entre o assentamento de poeira e a difusão turbulenta. Os testes confirmam que a inclusão da difusão turbulenta evita a formação de camadas de poeira irrealisticamente finas e mantém uma escala de altura de poeira finita.
• Simulações de Formação Estelar:
  • Primeiro Núcleo Hidrostático: Consistente com trabalhos anteriores (Bate 2022), ocorre um crescimento rápido de poeira dentro do primeiro núcleo hidrostático, produzindo populações localmente monodispersas com tamanhos que excedem 100 $\mu$m no interior do núcleo. A inclusão de arrasto e difusão teve impacto negligenciável nesta fase inicial devido às curtas escalas de tempo evolutivas.
  • Evolução do Disco Protoplanetário: Uma vez formado o disco, a inclusão da difusão turbulenta altera significativamente a evolução da poeira.
    • Crescimento Aumentado: A difusão turbulenta aumenta substancialmente a taxa de crescimento dos grãos de poeira. Ao fornecer continuamente grãos pequenos e de tamanho intermediário para regiões de alta densidade (onde grãos grandes estão crescendo), a difusão evita a depleção local de "alimento" para a coagulação.
    • Dependência da Turbulência: O nível de turbulência ($\alpha_{SS}$) é crítico. Uma turbulência alta ($\alpha_{SS} = 10^{-3}$) leva ao crescimento rápido para populações monodispersas de $\sim30-100$ $\mu$m em todo o disco. Uma turbulência menor ($\alpha_{SS} = 10^{-4}$) resulta em um crescimento mais lento e em uma população de poeira mais estratificada verticalmente, com grãos maiores confinados mais próximos ao plano médio.
    • Assentamento de Poeira: Em simulações com turbulência fraca, o assentamento da poeira leva a um aumento significativo na razão poeira-gás no plano médio (em mais de 50%), enquanto grãos pequenos permanecem empobrecidos perto do plano médio e confinados às camadas superficiais.
  • Arrasto vs. Difusão: Comparações de simulações com e sem processos específicos revelam que o arrasto dinâmico poeira-gás tem pouco efeito nas distribuições de poeira nestas escalas de tempo, pois os grãos permanecem bem acoplados. Em contraste, a inclusão da difusão turbulenta tem um efeito enorme, suavizando as distribuições de poeira e acelerando a coagulação.

Significância e Alegações
O artigo afirma que esta nova metodologia permite, pela primeira vez, a modelagem autoconsistente do crescimento, mistura e migração de poeira durante a formação de estrelas e seus discos. Os autores enfatizam que, embora a inclusão de arrasto e difusão não altere as descobertas fundamentais sobre o crescimento de poeira no primeiro núcleo hidrostático, isso é crucial para a evolução subsequente do disco protoplanetário. Especificamente, eles afirmam que a difusão de poeira turbulenta proporciona um aumento substancial na taxa de crescimento dos grãos de poeira ao atuar como uma fonte de grãos pequenos para as regiões onde grãos grandes estão se formando.

Os autores mantêm-se modestos quanto ao escopo de suas descobertas, observando que:

• Os resultados dependem do modelo de turbulência adotado (ex: disco-$\alpha$ vs. turbulência MHD).
• As simulações consideram atualmente a coagulação, mas não a fragmentação; diferentes limiares de coagulação poderiam levar à fragmentação e repopulação de grãos pequenos.
• Os cálculos de transferência radiativa ainda assumem opacidades padrão do meio interestelar em vez de usar as distribuições de tamanho de grão locais em evolução, um passo identificado como uma extensão futura necessária.

Em última análise, o trabalho fornece uma ferramenta para comparar simulações hidrodinâmicas 3D com a riqueza de novos dados observacionais de instalações como ALMA, VLT-SPHERE e JWST, particularmente no que diz respeito ao assentamento de poeira e agitação turbulenta em discos vistos de perfil (edge-on).