Formation of Isotopically Heterogeneous Molecular Cloud Cores in Filamentary Molecular Clouds

Resumen Técnico: Formación de núcleos de nubes moleculares isotópicamente heterogéneos en nubes moleculares filamentosas

Planteamiento del problema
El análisis de meteoritos revela que los sólidos del sistema solar temprano, como las inclusiones ricas en calcio y aluminio (CAIs), exhiben inhomogeneidad isotópica (por ejemplo, variaciones en las razones $^{54}\text{Cr}/^{52}\text{Cr}$ y $^{50}\text{Ti}/^{47}\text{Ti}$). Se han propuesto tanto el procesamiento térmico en discos protoplanetarios como la herencia de núcleos de nubes moleculares progenitoras inhomogéneas como explicaciones, pero la medida en que la heterogeneidad isotópica sobrevive al colapso de las nubes moleculares filamentosas hacia núcleos densos sigue sin estar clara. Específicamente, es necesario cuantificar cuánto de la variación isotópica inicial presente en una nube filamentosa se preserva o se suaviza durante la fragmentación gravitacional y la formación de núcleos, y cómo esta heterogeneidad se traduce al disco circunestelar resultante.

Metodología
Los autores utilizaron simulaciones de Hidrodinámica de Partículas Suavizadas (SPH), empleando específicamente el método de magnetohidrodinámica de partículas suavizadas de Godunov (GSPM), para modelar la fragmentación gravitacional de nubes moleculares filamentosas.

• Configuración de la simulación: Las simulaciones modelaron filamentos con una masa lineal crítica ($M_{\text{line,crit}} \approx 18 M_\odot \text{pc}^{-1}$) y un ancho de $\sim 0.1$ pc, consistentes con las observaciones. Las condiciones iniciales incluyeron un perfil de densidad en equilibrio hidrostático, turbulencia de Kolmogorov ($\sigma = 2c_s$), y un campo magnético inicial de $10 \mu\text{G}$ perpendicular al eje del filamento. También se realizaron casos puramente hidrodinámicos para comparación.
• Rastreo isotópico: Las especies isotópicas fueron tratadas como trazadores pasivos. Los autores construyeron seis modelos distintos (Tabla 1) para imponer gradientes de razón isotópica iniciales ($C_0$) a lo largo de diferentes ejes: el eje longitudinal ($z$), los ejes menores ($x, y$), y combinaciones de estos. La razón isotópica se normalizó entre 0 y 1 sobre el dominio computacional.
• Análisis: Se rastreó la evolución de la razón isotópica siguiendo las partículas SPH. Los autores definieron capas esféricas alrededor del pico de densidad de los núcleos en formación para calcular la razón isotópica promedio ($\bar{C}$) y la diferencia respecto al centro del núcleo ($|\Delta C|$). También emplearon un enfoque semianalítico para estimar el desplazamiento inicial del centro de masa de las capas ($\delta z_{\text{ini}}$) basado en campos de velocidad turbulentos y espectros de potencia de densidad.

Result mayoría de resultados clave

1. Supervivencia de la heterogeneidad: Las simulaciones demuestran que las inhomogeneidades isotópicas no se suavizan completamente durante la fragmentación del filamento. Los núcleos retienen aproximadamente del 0.1% al 1% de la inhomogeneidad isotópica inicial presente en la longitud inicial del filamento de 1 pc.
2. Dependencia geométrica: La magnitud de la heterogeneidad retenida depende de la orientación del gradiente isotópico inicial respecto a la geometría del filamento.
   • Los gradientes a lo largo del eje longitudinal (dirección $z$, Modelo $M_z$) resultan en las mayores diferencias isotópicas entre capas ($|\Delta C| \sim 0.01$).
   • Los gradientes a lo largo de los ejes menores ($x$ o $y$) resultan en diferencias significativamente menores. Esto se atribuye a la geometría del filamento, que restringe el desplazamiento de los elementos de fluido a lo largo de los ejes menores, limitando así la mezcla de material con diferentes razones isotópicas.
3. Mecanismo de la inhomogeneidad: El principal motor de la diferencia isotópica residual entre las capas es la diferencia inicial en el centro de masa de las capas, causada por el campo de velocidad turbulenta. Las estimaciones semianalíticas, que consideran el desplazamiento inicial de los elementos de fluido, concuerdan con los resultados de la simulación dentro de un factor de 2, confirmando que la configuración turbulenta inicial dicta el perfil isotópico final.
4. Influencia del campo magnético: La presencia de un campo magnético ($10 \mu\text{G}$) no altera significativamente el perfil interno de la razón isotópica en comparación con el caso puramente hidrodinámico.
5. Implicaciones para el disco circunestelar: Al estimar el radio centrífugo del material que acrece, los autores encontraron que la heterogeneidad isotópica persiste en el disco circunestelar. Aunque el disco hereda la heterogeneidad del núcleo, la distribución anisotrópica del momento angular implica que el gas acrece desde varias capas esféricas. Los resultados sugieren que aproximadamente el 1% de la heterogeneidad inicial del filamento permanece en el disco.

Significancia y afirmaciones
El artículo afirma proporcionar un vínculo cuantitativo entre las anomalías isotópicas observadas en meteoritos y las condiciones astrofísicas de la nube molecular progenitora. Al combinar sus hallazgos con el factor de reducción reportado por Jacquet et al. (2019) (que sugiere una reducción de un factor de $\sim 3$ del núcleo al disco), los autores infieren que la heterogeneidad residual en el disco circunestelar final es aproximadamente del 0.03% al 0.3% de la variación inicial del filamento.

Los autores aseveran que este nivel de heterogeneidad residual es suficiente para explicar las anomalías isotópicas observadas en los materiales del Sistema Solar (por ejemplo, variaciones de $^{54}\text{Cr}/^{52}\text{Cr}$ y $^{50}\text{Ti}/^{47}\text{Ti}$ de 1 a 10 partes en 10,000) sin requerir un procesamiento térmico extremo o escenarios de inyección nucleosintética exótica post-colapso. El estudio concluye que la inhomogeneidad isotópica de los componentes de construcción del sistema solar podría sobrevivir naturalmente a todo el proceso, desde la fragmentación del filamento hasta la formación del disco, siempre que la nube progenitora poseyera un nivel apropiado de inhomogeneidad inicial. Los autores señalan que se requieren simulaciones de alta resolución para trabajos futuros para resolver completamente las propiedades realistas de las razones isotópicas en el disco circunestelar, particularmente en lo que respecta a la caída (infall) y salidas (outflows) anisotrópicas.