Turbulent infall onto class 0 disks as cause of CAI brief condensation episode in the solar system

Resumo Técnico: Queda Turbulenta sobre Discos de Classe 0 como Causa do Episódio Breve de Condensação de Inclusões Ricas em Cálcio e Alumínio no Sistema Solar

Declaração do Problema
As inclusões ricas em cálcio e alumínio (CAIs) encontradas em meteoritos condritos carbonáceos representam os sólidos mais antigos do sistema solar, caracterizados por um episódio de condensação notavelmente breve, com duração aproximada de $10^5$ anos. Um desafio central na teoria da formação do sistema solar é conciliar essa curta escala de tempo de formação com a poluição de longa duração (Myr) das regiões de formação estelar por radionuclídeos de vida curta (SLRs) provenientes de fontes externas, como supernovas ou estrelas do ramo gigante assintótico. Embora modelos anteriores, como o modelo de vento X, tenham tentado explicar a formação de CAIs por meio de ventos magnetocentrífugos, a datação isotópica precisa restringiu a formação de CAIs a uma janela estreita que desafia modelos baseados em irradiação solar prolongada. Além disso, embora vários mecanismos de aquecimento (por exemplo, erupções de FU Orionis, supernovas próximas) tenham sido propostos, eles frequentemente lutam para explicar simultaneamente a extensão espacial (várias UA), a natureza transitória e a intensidade do aquecimento necessários para sublimar e recondensar grãos refratários. Os autores investigam se a queda dinâmica de gás sobre discos protostelares de Classe 0 pode gerar as condições térmicas necessárias para esse episódio breve e de alta temperatura.

Metodologia
O estudo emprega simulações hidrodinâmicas de alta resolução utilizando o código GIZMO no modo de Massa Finita sem Malha (MFM). Este esquema Lagrangiano é selecionado por sua captura superior de choques, dissipação numérica reduzida e excelente conservação de momento angular, que são críticos para modelar discos torcidos e queda turbulenta.

• Condições Iniciais: As simulações utilizam condições iniciais derivadas de um estudo de síntese populacional da formação de aglomerados estelares (Bate 2012, 2019). Os autores traçam a evolução de estrelas de massa solar de volta ao estágio de Classe 0, capturando histórias de queda realistas, dependentes do tempo e anisotrópicas provenientes de fluxos de acreção filamentares.
• Configuração da Simulação: Duas simulações principais são conduzidas:
  • Execução 1: Apresenta uma fase de queda rápida e caótica com fluxos de momento angular significativamente desalinhados impactando o disco.
  • Execução 2: Um caso de controle com configuração inicial similar, mas sem queda vertical significativa de gás a partir de fluxos, servindo como linha de base para a evolução isolada do disco.
• Modelos Físicos:
  • Transferência Radiativa: Em vez de uma equação de estado barotrópica simples, os autores empregam um método de resfriamento radiativo aproximado (Young et al. 2024) que hibridiza estimadores de densidade de coluna (Stamatellos et al. 2007; Lombardi Jr et al. 2015) para lidar com geometrias esféricas e semelhantes a discos. Isso permite uma evolução térmica realista onde o aquecimento viscoso domina sobre a irradiação no plano médio opticamente espesso.
  • Resolução: Um estudo de resolução (Apêndice A) confirma que as características físicas principais, especificamente ondas de densidade espiral e aquecimento por choque associado, convergem em contagens de partículas de $\sim 10$ milhões. As simulações finais utilizam 10 milhões de partículas.

Principais Resultados
As simulações revelam uma evolução distinta em três fases na Execução 1 (o cenário de queda), contrastando fortemente com a evolução estável da Execução 2.

1. Evolução Global e Reorientação do Disco:

   • Execução 1: O sistema passa por uma "fase principal de aquecimento" ($t \approx 6000-9000$ anos) impulsionada por uma queda vigorosa e desalinhada. A interação entre os fluxos de queda e o disco rotativo causa uma reorientação global do disco em quase $30^\circ$. Isso leva à formação de torções e ondas de densidade espiral globais.
   • Execução 2: Sem queda externa significativa, o disco permanece geometricamente fino e estável, com apenas inclinações menores ($10^\circ-20^\circ$) e nenhuma reorientação em grande escala.

2. Evolução Térmica:

   • Execução 1: A fase principal de aquecimento é caracterizada por um rápido aumento de temperatura. A temperatura máxima do disco atinge $\sim 2000$ K, suficiente para sublimar grãos de CAI. Esse aquecimento é sustentado por vários milhares de anos. O mecanismo de aquecimento é identificado como dissipação de choque na interface disco-envelope e dentro do plano médio do disco devido a choques espirais. O disco denso e verticalmente espessado durante o pico da queda reduz a eficiência do resfriamento radiativo, aprisionando o calor.
   • Execução 2: A temperatura do disco permanece significativamente mais baixa, atingindo um pico apenas em torno de $400-600$ K. A ausência de queda externa impede a geração do intenso aquecimento por choque observado na Execução 1.

3. Mecanismos de Aquecimento:

   • Choques de Superfície: Fatias verticais mostram fluxos de queda convergentes colidindo com a superfície do disco, criando choques fortes e zonas de alta entropia.
   • Choques Espirais: A análise do plano médio revela que ondas de densidade espiral (impulsionadas por instabilidade gravitacional com Toomre $Q \approx 1$) geram saltos localizados de entropia e aquecimento. No disco interno ($r \lesssim 3$ UA), as temperaturas atingem $\sim 2000$ K, enquanto choques espirais na faixa de 2–3 UA criam pontos quentes superiores a 1500 K.
   • Momento Angular: O desalinhamento dos vetores de momento angular em queda impulsiona a reorientação e a turbulência, que por sua vez aumentam a dissipação viscosa e a formação de choques.

Significado e Alegações
O artigo alega que a queda turbulenta sobre discos de Classe 0 fornece um mecanismo viável para o episódio breve e de alta temperatura de condensação necessário para formar CAIs. Especificamente:

• Aquecimento Transitório: O aquecimento impulsionado pela queda é transitório, durando $\sim 10^3-10^4$ anos, o que se alinha com as restrições de idade radiogênica das CAIs.
• Extensão Espacial: O aquecimento se estende a $\sim 2-3$ UA, cobrindo a região onde as CAIs são pensadas para se formar e serem preservadas.
• Recondensação: O modelo sugere um ciclo onde grãos carregados com CAIs "pré-solares" em nuvens em queda são sublimados pelo aquecimento por choque e subsequentemente recondensam à medida que o sistema esfria durante a fase de subsídio da queda.
• Dicotomia de Sistemas: A comparação entre a Execução 1 e a Execução 2 implica que a formação de um reservatório de CAIs de população única não é um resultado inevitável de toda formação de disco. Depende de condições de contorno específicas, como a presença de fluxos de queda caóticos e desalinhados. Isso pode explicar por que as CAIs são enriquecidas em condritos carbonáceos (CC), mas escassas em condritos não carbonáceos (NC), sugerindo que apenas uma fração dos sistemas planetários experimenta as dinâmicas de queda específicas necessárias para redefinir globalmente o relógio das CAIs.

Os autores concluem que, embora suas simulações cubram o estágio muito inicial da formação do disco (Classe 0, massa central $\sim 0,2 M_\odot$), o mecanismo de aquecimento impulsionado pela queda oferece uma explicação robusta para as condições térmicas necessárias para a formação de CAIs, ligando a história dinâmica do envelope protostelar diretamente ao registro meteorítico.