The Sensitivity of Substructure Lensing to SIDM Core-collapse Model Variation

Resumen Técnico: La Sensibilidad del Lenteado de Subestructura a la Variación del Modelo de Colapso del Núcleo de SIDM

Enunciado del Problema
El lenteado gravitacional fuerte ha surgido como una sonda poderosa para la subestructura de materia oscura, ofreciendo una prueba específica para los modelos de Materia Oscura Auto-Interactuante (SIDM). En escenarios SIDM con colisiones elásticas, la transferencia de calor conduce a la formación de núcleos isotérmicos, seguidos por una "catástrofe gravotermal" (colapso del núcleo) que produce núcleos centrales densos. Se predice que estos halos compactos dejan huellas distintivas en las observables de lenteado. Sin embargo, la predicción precisa de estas señales se ve obstaculizada por la dificultad de modelar la línea temporal del colapso del núcleo y el estado final de los subhalos. Los enfoques actuales a menudo dependen de suposiciones simplificadas, como tratar el colapso del núcleo como una transición instantánea o utilizar asignaciones probabilísticas para los estados de colapso. Los autores identifican una desconexión entre las pipelines existentes de inferencia de lenteado y los modelos físicos subyacentes, señalando que pequeñas variaciones en cómo se modelan los subhalos en colapso pueden conducir a diferencias significativas en la potencia de lenteado. Sin cuantificar estas incertidumbres de modelado, la capacidad de las observaciones de lenteado para restringir la naturaleza de la materia oscura es limitada.

Metodología
Los autores cuantifican este desafío de modelado probando varias modificaciones a los perfiles de densidad de halos en colapso del núcleo SIDM dentro del contexto del lenteado de subestructura. Su metodología implica:

1. Simulación y Calibración: Utilizan datos de simulación de $N$-cuerpos (generados con el código Arepo y el módulo SIDM de la Ref. [20]) para calibrar un tratamiento aproximado de los efectos de marea sobre el colapso gravotermal del núcleo. Derivan un "factor de aceleración de colapso" para tener en cuenta la diferencia en las escalas de tiempo de colapso entre halos aislados y subhalos que evolucionan en un potencial huésped.
2. Realizaciones de Lenteado: Utilizando el paquete de código abierto pyHalo, generan 200 realizaciones aleatorias de sistemas de lentes a escala galáctica (en $z_{lente}=0.5$) que contienen subhalos del deflector principal y halos a lo largo de la línea de visión (LOS). Reemplazan los perfiles estándar de Materia Oscura Fría (CDM) con perfiles SIDM basados en una solución de fluido gravotermal calibrada con simulaciones de $N$-cuerpos.
3. Variaciones del Modelo: El estudio prueba varias opciones de modelado:
   • Continuo vs. Instantáneo: Comparando perfiles de densidad que evolucionan suavemente contra modelos binarios donde los halos están completamente sin colapsar o completamente colapsados.
   • Umbral de Colapso: Variando el momento en que el perfil de densidad se "congela" (detiene) o se "impulsa" al estado final, probando diferentes estados finales del colapso.
   • Probabilístico vs. Determinista: Comparando el seguimiento de la evolución individual de halos contra la asignación aleatoria de estados de colapso basados en bins de masa.
   • Perfiles Iniciales: Probando halos sin colapsar modelados ya sea como TNFW (NFW truncado) estándar o como habiendo alcanzado la fase de núcleo maximamente desarrollado.
4. Análisis Estadístico: Calculan el potencial de lenteado gravitacional efectivo de múltiples planos, calculando la divergencia ($\kappa_{div}$) y el rotacional ($\kappa_{curl}$) del campo de deflexión efectivo. La estadística principal es la función de correlación de dos puntos de estos campos, descompuesta en momentos multipolares (monopolo y cuadrupolo) utilizando bases ortonormales para desentrañar los efectos de la LOS del lente principal.

Contribuciones Clave

• Cuantificación de la Incertidumbre de Modelado: El artículo proporciona una cuantificación sistemática de cómo suposiciones simplificadas específicas en el modelado del colapso del núcleo SIDM afectan las estadísticas de lenteado.
• Calibración de Marea: Introduce un método probabilístico para tener en cuenta la aceleración/desaceleración por marea del colapso del núcleo en subhalos, calibrado contra simulaciones de $N$-cuerpos.
• Comparación Binaria vs. Continua: Compara rigurosamente la aproximación de "colapso binario" (común en la literatura) contra un modelo gravotermal continuo y evolutivo en el tiempo, destacando dónde falla la aproximación.
• Sensibilidad del Perfil: El estudio demuestra que la elección del perfil de densidad inicial (TNFW vs. con núcleo) y el método de asignación de estados de colapso alteran significativamente la señal de lenteado predicha.

Resultados

• Sensibilidad a la Densidad vs. Tiempo: La función de correlación de dos puntos (específicamente los momentos multipolares) es altamente sensible a la densidad central de los halos colapsados pero relativamente insensible a la escala de tiempo exacta del colapso. Las variaciones en el tiempo de colapso (por ejemplo, deteniendo la evolución al 20% del tiempo total) no producen resultados significativamente diferentes en comparación con la evolución completa, siempre que el perfil de densidad final sea preciso.
• Supresión por Modelos Binarios: Los modelos de colapso del núcleo binario suprimen consistentemente la potencia a pequeña escala en las estadísticas de lenteado en comparación con el modelo de evolución temporal completa. Esta supresión es más severa cuando:
  • Los halos sin colapsar se modelan con perfiles iniciales con núcleo en lugar de perfiles TNFW.
  • Los estados de colapso se asignan aleatoriamente en lugar de basarse en la evolución individual del halo.
• Dependencia de la Escala: Las diferencias entre modelos son más distinguibles a escalas de longitud pequeñas (por debajo de 0.1 segundos de arco). El momento cuadrupolar del campo de divergencia ($\xi_{div,2}$) permanece sensible a las variaciones del modelo a escalas mayores que el monopolo o los momentos de rotacional.
• Correlación entre Colapso y Potencia de Lenteado: El estudio encuentra que los halos con un efecto de lenteado más pronunciado (mayor masa/concentración) tienden a colapsar más rápido. La aleatorización de los estados de colapso suprime la mejora SIDM porque desacopla el estado de colapso de las propiedades del halo que impulsan la señal de lenteado.
• Validación con N-cuerpos: El modelo de aceleración por marea aproximado de los autores produce estadísticas de lenteado consistentes con las simulaciones de $N$-cuerpos dentro de una desviación estándar, validando su uso como una aproximación de primer orden para este estudio.

Significado y Afirmaciones
El artículo afirma que la función de correlación de dos puntos es una sonda sensible de las simplificaciones específicas realizadas en el modelado del colapso del núcleo SIDM. Los autores afirman que:

1. Robustez del Estado Final: Las incertidumbres sistemáticas modestas en el tiempo del colapso del núcleo pueden no afectar significativamente las predicciones de lenteado, siempre que el perfil de densidad final del halo colapsado se modele con precisión.
2. Limitaciones de las Aproximaciones Binarias: Aunque los modelos de colapso binario son computacionalmente eficientes, pueden introducir sesgos significativos al suprimir la potencia a pequeña escala, particularmente si asumen condiciones iniciales con núcleo o aleatorizan los estados de colapso.
3. Necesidad de Motivación Física: Los futuros análisis de lenteado destinados a probar SIDM deben realizar elecciones intencionales y motivadas físicamente con respecto a los perfiles de subestructura y halos a lo largo de la línea de visión. Los autores recomiendan que tales trabajos incluyan pruebas análogas de sus elecciones de perfiles y comparen sus perfiles de ángulo de deflexión contra modelos de evolución completa motivados físicamente.
4. Dependencia Observacional: La detectabilidad de estas diferencias de modelado depende en gran medida de la sonda observacional y de las escalas que se analizan. El artículo no afirma que los datos actuales puedan descartar definitivamente modelos SIDM específicos, sino que la interpretación de dichos datos depende en gran medida de la robustez de los modelos de subestructura subyacentes.

Los autores concluyen que refinar los modelos del colapso del núcleo SIDM en entornos realistas es esencial, ya que las predicciones comprobables pueden depender fuertemente del estado final de los halos con colapso del núcleo. Sugieren que su marco puede utilizarse para implementar estos modelos de evolución completa en predicciones para anomalías de relación de flujo y otras observables de lenteado en trabajos futuros.