Non-linear Evolution of Dark Plasma Subhalos

Resumen Técnico: Evolución No Lineal de Subhalos de Plasma Oscuro

Planteamiento del Problema
Este artículo investiga la evolución dinámica de los subhalos de materia oscura (DM) que orbitan dentro del potencial de la Vía Láctea (MW) bajo la suposición de que la DM constituye un "plasma oscuro". En modelos donde la DM interactúa consigo misma mediante fuerzas de largo alcance mediadas por un fotón oculto (un sector $U(1)$ oscuro), la dinámica se asemeja a la de los plasmas del Modelo Estándar. Específicamente, cuando un subhalo de DM se mueve a través del plasma oscuro ambiental de la MW, la velocidad de deriva relativa entre las dos poblaciones puede exceder sus dispersiones de velocidad internas. Esta configuración es susceptible a inestabilidades de plasma (por ejemplo, inestabilidades de dos corrientes y de haz-plasma), las cuales generan campos electrostáticos y electromagnéticos oscuros. Estos campos dispersan las partículas, modifican la distribución en el espacio de fases y conducen al sistema hacia el equilibrio térmico. El problema central abordado es la evolución no lineal de estas inestabilidades y su impacto en la pérdida de masa de los subhalos, comparando específicamente el calentamiento inducido por plasma frente al despojo de marea gravitacional estándar.

Metodología
Los autores emplean un enfoque híbrido semianalítico y numérico para modelar la evolución de los subhalos:

1. Marco Semianalítico:

   • Condiciones Iniciales: Los subhalos se modelan como plasmas esféricos de carga neta neutra con perfiles de densidad NFW y distribuciones de velocidad de Maxwell-Boltzmann. Se consideran dos masas de referencia: $10^7 M_\odot$ y $10^9 M_\odot$.
   • Dinámica Orbital: Las órbitas de los subhalos se integran utilizando el paquete galpy dentro de un potencial de la MW (MWPotential2014). Las óridades se inicializan en el radio virial con diversas excentricidades ($\epsilon = 0.1$ a $0.9$).
   • Acoplamiento de la Pérdida de Masa: El perfil de densidad del subhalo se discretiza en capas radiales concéntricas. La evolución de cada capa se rastrea en pasos de tiempo ($\delta t = 7$ Myr) hasta $1.2 \times t_{\text{peri}}$ (tiempo del primer paso por el pericentro).
   • Mecanismos: La pérdida de masa se calcula a través de dos canales distintos:
     • Despojo de Marea (Tidal Stripping): Se modela utilizando un parámetro de eficiencia estándar ($\alpha = 0.55$) aplicado a las capas situadas fuera del radio de marea instantáneo.
     • Calentamiento por Plasma: La fracción de partículas ligadas ($f_{\text{bound}}$) en cada capa se determina mediante la interpolación de resultados de una suite de simulaciones de Partícula en Celda (PIC). Las partículas calentadas por encima de la velocidad de escape local ($v_{\text{esc}} \approx 4\sigma_{\text{SH}}$) son eliminadas, y la masa restante ligada se re-virializa.

2. Simulaciones de Partícula en Celda (PIC):

   • Código: El estudio utiliza el código TRISTAN-MP, adaptado para plasmas oscuros mediante la invariancia de reescalado de coordenadas de las ecuaciones de Vlasov-Poisson.
   • Configuración: Se realizan simulaciones electrostáticas unidimensionales (1D) con condiciones de contorno periódicas. La configuración modela haces de contravía que representan el subhalo y el fondo de la MW.
   • Parámetros: Las simulaciones varían la relación de densidad ($\rho_{\text{SH}}/\rho_{\text{MW}} \in [10^{-3}, 10^2]$) y la velocidad de deriva relativa ($v_{\text{SH}}/\sigma_{\text{MW}}$). La relación carga-masa de referencia es $q_\chi/m_\chi = 10^{-14} e/m_p$.
   • Dimensionalidad: Aunque los resultados principales son 1D, los autores realizan una simulación limitada 2D-2V (dos dimensiones espaciales, dos de velocidad) para verificar que la eficiencia de calentamiento aumenta en dimensiones superiores, lo que sugiere que los resultados 1D son conservadores.

Contribuciones Clave

• Primer Estudio No Lineal: Este trabajo presenta el primer estudio de la evolución no lineal de las inestabilidades electrostáticas en subhalos de DM utilizando simulaciones cinéticas PIC completas.
• Cuantificación del Calentamiento por Plasma: Los autores cuantifican la eficiencia de la pérdida de masa inducida por plasma, demostrando que puede exceder significativamente el despojo de marea, particularmente para subhalos de baja masa en órbitas excéntricas.
• Modelo Semianalítico: Se desarrolla un marco novedoso que acopla el crecimiento de la inestabilidad de plasma local (derivado de los datos PIC) con la evolución orbital global para predecir la evolución de la masa y el perfil de densidad del subhalo.

Resultos

• Eficiencia de Pérdida de Masa: Para subhalos en órbitas altamente excéntricas ($\epsilon = 0.9$), el calentamiento por plasma es el mecanismo dominante de pérdida de masa. Para el momento del primer pericentro ($t = 1.2 t_{\text{peri}}$):
  • Un subhalo de $10^7 M_\odot$ pierde $\sim 97\%$ de su masa inicial.
  • Un subhalo de $10^9 M_\odot$ pierde $\sim 84\%$ de su masa inicial.
  • En contraste, el despojo de marea por sí solo elimina solo $\sim 30\%$ de la masa para ambos casos.
• Dependencia Radial ("Despojo de Afuera hacia Adentro"): La pérdida de masa procede desde el exterior hacia el interior. Las inestabilidades crecen más eficientemente donde el contraste de densidad entre el subhalo y la MW es modesto ($\rho_{\text{SH}}/\rho_{\text{MW}} \sim 1$) y la velocidad relativa es alta. Consecuentemente, las capas externas son calentadas y despojadas antes que el núcleo. Para órbitas altamente excéntricas, el núcleo eventualmente se ve afectado cerca del pericentro; para órbitas circulares, la pérdida de masa se confina a las regiones exteriores.
• Implicaciones en la Función de Masa: El despojo preferencial de subhalos de baja masa (debido a sus menores dispersiones de velocidad interna y velocidades de escape) sugiere una supresión del extremo de baja masa de la Función de Masa de Subhalos (SHMF). Los autores estiman que, para los parámetros de referencia, la SHMF podría exhibir una pendiente más plana o un corte abrupto por debajo de $M_{\text{SH}} \sim 10^7 M_\odot$.
• Excentricidad Orbital: El efecto es altamente sensible a la excentricidad. Para órbitas circulares ($\epsilon = 0.1$), el calentamiento por plasma elimina solo $\sim 50\%$ ($10^7 M_\odot$) y $\sim 30\%$ ($10^9 M_\odot$) de la masa, valores comparables o ligeramente superiores a los efectos de marea.

Significado y Reivindicaciones
El artículo sostiene que la física del plasma oscuro ofrece un mecanismo distinto para la evolución de los subhalos que puede dejar firmas observables en la población de satélites de la MW. Específicamente, los autores argumentan que el calentamiento y el despojo inducidos por plasma pueden:

1. Remodelar la SHMF: Al destruir preferencialmente los subhalos de baja masa y "fríos", los modelos de plasma oscuro pueden resolver o alterar el problema de los "satélites faltantes" o el problema de "demasiado grande para fallar" (too big to fail), lo que potencialmente conduce a una pendiente de baja masa suprimida en la población observada de subhalos.
2. Proporcionar Restricciones Observacionales: La supresión de subhalos de baja masa y la consecuente modificación de la SHMF proporcionan una nueva vía para restringir los parámetros del sector oscuro (específicamente la relación carga-masa $q_\chi/m_\chi$) utilizando datos de sondeos como ELVES, SPARC, SAGA y Merian.
3. Destacar la Dinámica No Lineal: El estudio enfatiza que los efectos colectivos de plasma, impulsados por inestabilidades de corriente, pueden operar en escalas de tiempo comparables a los tiempos dinámicos galácticos, lo que conduce a una evaporación rápida que las simulaciones de CDM colisionales estándar no capturan.

Los autores mantienen un tono modesto respecto al alcance de sus hallazgos, señalando que sus resultados se basan en simulaciones 1D (aunque las pruebas 2D sugieren que el modelo 1D es conservador) y que una predicción completa a nivel de población requiere un trabajo futuro que incorpore la gravedad autoconsciente y los efectos de plasma en dimensiones superiores.