Formation of Isotopically Heterogeneous Molecular Cloud Cores in Filamentary Molecular Clouds

Resumo Técnico: Formação de Núcleos de Nuvens Moleculares Isotopicamente Heterogêneos em Nuvens Moleculares Filamentares

Definição do Problema
A análise de meteoritos revela que os sólidos do sistema solar primitivo, como as inclusões ricas em cálcio-alumínio (CAIs), exibem heterogeneidade isotópica (por exemplo, variações nas razões $^{54}\text{Cr}/^{52}\text{Cr}$ e $^{50}\text{Ti}/^{47}\text{Ti}$). Embora o processamento térmico em discos protoplanetários e a herança de núcleos de nuvens moleculares progenitoras heterogêneas sejam explicações propostas, permanece incerto o quanto a heterogeneidade isotópica sobrevive ao colapso de nuvens moleculares filamentares em núcleos densos. Especificamente, é necessário quantificar quanto da variação isotópica inicial presente em uma nuvem filamentar é preservada ou suavizada durante a fragmentação gravitacional e a formação de núcleos e como essa heterogeneidade se traduz no disco circunestelar resultante.

Metodologia
Os autores utilizaram simulações de Hidrodinâmica de Partículas Suavizadas (SPH), empregando especificamente o método magnetohidrodinâmico de partículas suavizadas de Godunov (GSPM), para modelar a fragmentação gravitacional de nuvens moleculares filamentares.

• Configuração da Simulação: As simulações modelaram filamentos com uma massa linear crítica ($M_{\text{line,crit}} \approx 18 M_\odot \text{pc}^{-1}$) e uma largura de $\sim 0,1$ pc, consistentes com observações. As condições iniciais incluíram um perfil de densidade em equilíbrio hidrostático, turbulência de Kolmogorov ($\sigma = 2c_s$) e um campo magnético inicial de $10 \mu\text{G}$ perpendicular ao eixo do filamento. Casos puramente hidrodinâmicos também foram executados para comparação.
• Rastreamento Isotópico: As espécies isotópicas foram tratadas como traçadores passivos. Os autores construíram seis modelos distintos (Tabela 1) para impor gradientes de razão isotópica inicial ($C_0$) ao longo de diferentes eixos: o eixo longitudinal ($z$), os eixos menores ($x, y$) e combinações destes. A razão isotópica foi normalizada entre 0 e 1 sobre o domínio computacional.
• Análise: A evolução da razão isotópica foi rastreada seguindo as partículas SPH. Os autores definiram cascas esféricas ao redor do pico de densidade dos núcleos em formação para calcular a razão isotópica média ($\bar{C}$) e a diferença em relação ao centro do núcleo ($|\Delta C|$). Eles também empregaram uma abordagem semi-analítica para estimar o deslocamento inicial do centro de massa das cascas ($\delta z_{\text{ini}}$) baseado em campos de velocidade turbulentos e espectros de potência de densidade.

Principais Resultados

1. Sobrevivência da Heterogeneidade: As simulações demonstram que as heterogeneidades isotópicas não são completamente suavizadas durante a fragmentação do filamento. Os núcleos retêm aproximadamente 0,1% a 1% da heterogeneidade isotópica inicial presente ao longo do comprimento inicial do filamento de 1 pc.
2. Dependência Geométrica: A magnitude da heterogeneidade retida depende da orientação do gradiente isotópico inicial em relação à geometria do filamento.
   • Gradientes ao longo do eixo longitudinal (direção $z$, Modelo $M_z$) resultam nas maiores diferenças isotópicas entre as cascas ($|\Delta C| \sim 0,01$).
   • Gradientes ao longo dos eixos menores ($x$ ou $y$) resultam em diferenças significativamente menores. Isso é atribuído à geometria do filamento, que restringe o deslocamento dos elementos de fluido ao longo dos eixos menores, limitando assim a mistura de material com diferentes razões isotópicas.
3. Mecanismo de Inhomogeneidade: O principal motor da diferença isotópica residual entre as cascas é a diferença inicial no centro de massa das cascas, causada pelo campo de velocidade turbulenta. As estimativas semi-analíticas, que consideram o deslocamento inicial dos elementos de fluido, concordam com os resultados da simulação dentro de um fator de 2, confirmando que a configuração turbulenta inicial dita o perfil isotópico final.
4. Influência do Campo Magnético: A presença de um campo magnético ($10 \mu\text{G}$) não altera significativamente o perfil interno da razão isotópica em comparação com o caso puramente hidrodinâmico.
5. Implicações para o Disco Circunestelar: Ao estimar o raio centrífugo do material que accreta, os autores descobriram que a heterogeneidade isotópica persiste no disco circunestelar. Embora o disco herde a heterogeneidade do núcleo, a distribuição anisotrópica do momento angular significa que o gás accreta de várias cascas esféricas. Os resultados sugerem que aproximadamente 1% da heterogeneidade filamentar inicial permanece no disco.

Significância e Alegações
O artigo afirma fornecer uma ligação quantitativa entre as anomalias isotópicas observadas em meteoritos e as condições astrofísicas da nuvem molecular progenitora. Combinando seus achados com o fator de redução relatado por Jacquet et al. (2019) (que sugere uma redução de um fator de $\sim 3$ do núcleo para o disco), os autores inferem que a heterogeneidade residual no disco circunestelar final é de aproximadamente 0,03% a 0,3% da variação filamentar inicial.

Os autores afirmam que este nível de heterogeneidade residual é suficiente para explicar as anomalias isotópicas observadas nos materiais do Sistema Solar (por exemplo, variações de $^{54}\text{Cr}/^{52}\text{Cr}$ e $^{50}\text{Ti}/^{47}\text{Ti}$ de 1 a 10 partes em 10.000) sem exigir processamento térmico extremo ou cenários exóticos de injeção nucleossintética pós-colapso. O estudo conclui que a inhomogeneidade isotópica dos blocos de construção do sistema solar pode sobreviver naturalmente a todo o processo, desde a fragmentação do filamento até a formação do disco, desde que a nuvem progenitora possuísse um nível apropriado de inhomogeneidade inicial. Os autores observam que simulações de alta resolução são necessárias para trabalhos futuros para resolver plenamente as propriedades realistas das razões isotópicas no disco circunestelar, particularmente no que diz respeito à queda (infall) e fluxos (outflows) anisotrópicos.