Collapse versus Disruption: The Fate of Compact Stellar Systems in Ultralight Dark Matter Halos

Resumen Técnico: Colapso frente a Disrupción en Halos de Materia Oscura Ultraligera

Planteamiento del Problema
Los modelos de materia oscura ultraligera (ULDM, por sus siglas en inglés) predicen que la interferencia de las ondas de ULDM genera gránulos de densidad estocásticos, lo que genera fluctuaciones en el potencial gravitatorio que varían en el tiempo. Estas fluctuaciones inducen un calentamiento dinámico en los sistemas estelares incrustados, provocando su expansión. Si bien este mecanismo se ha utilizado para establecer límites inferiores estrictos sobre la masa de la partícula de ULDM ($m_a$) utilizando sistemas pequeños como cúmulos estelares nucleares y galaxias enanas ultra débiles (UFDs), el comportamiento evolutivo de sistemas significativamente más pequeños que la longitud de onda de de Broglie ($\lambda_{dB}$) de la ULDM sigue siendo poco comprendido. Específicamente, no está claro cómo la interacción entre el calentamiento inducido por la ULDM y la relajación estelar interna (impulsada por la dispersión gravitatoria de dos cuerpos) afecta el destino de los sistemas compactos. Estudios previos sugirieron que el calentamiento podría ser ineficaz para sistemas muy pequeños, mientras que otros señalaron que los sistemas compactos tienen tiempos de relajación más cortos que la edad del Universo, lo que potencialmente permite que el colapso del núcleo domine.

Metodología
Los autores emplean un marco de simulación numérica de dos etapas para resolver las escalas espaciales y temporales dispares involucradas:

1. Etapa 1 (Evolución del Halo): El halo de ULDM se modela como un todo, evolucionando dentro de los campos de marea de la Vía Láctea (MW) y la Gran Nube de Magallanes (LMC). El sistema estelar se trata como una masa puntual ($M_\star = 10^3 M_\odot$) para registrar su trayectoria y el tensor de marea dependiente del tiempo inducido por las inhomogeneidades de la ULDM. La simulación adopta una historia de ensamblaje realista donde el halo evoluciona fuera de la MW durante $\sim 3$ Gyr antes de entrar en una trayectoria similar a la de la UFD Segue 1.
2. Etapa 2 (Evolución Interna): Utilizando el código de $N$-cuerpos directos PeTar, la evolución interna del sistema estelar se simula bajo el tensor de marea dependiente del tiempo derivado de la Etapa 1. El sistema estelar se inicializa con un perfil Plummer en equilibrio dinámico aproximado.
   • Parámetros: Se simulan seis sistemas representativos (S1–S6) con variaciones en la masa estelar total ($M_\star$), el número de partículas ($N$) y los radios de semimasa ($R_h$).
   • Comparación: Cada sistema se evoluciona tanto con el calentamiento inducido por la ULDM como en aislamiento (sin ULDM) para aislar el efecto de calentamiento.
   • Física: La aceleración de la ULDM se aproxima mediante el tensor de marea que actúa sobre las partículas estelares en relación con el centro de masa del sistema. Las simulaciones corren durante 13 Gyr.

Contribuciones Clave y Resultados
El estudio revela que la evolución de los sistemas estelares compactos en halos de ULDM está gobernada por una competencia entre la relajación estelar interna y el calentamiento inducido por la ULDM, lo que conduce a tres regímenes evolutivos distintos:

1. Régimen Dominado por la Relajación (Colapso del Núcleo):

   • En sistemas suficientemente compactos (por ejemplo, S1, S3, S4), la relajación interna de dos cuerpos impulsa el colapso del núcleo.
   • Hallazgo Contraintuitivo: El calentamiento de la ULDM no evita el colapso; en su lugar, lo acelera. Al desprender las estrellas exteriores, el calentamiento de la ULDM reduce la masa estelar ligada y el radio de semimasa. Esta reducción acorta la escala de tiempo de relajación ($t_{rh}$) y reduce la dispersión de velocidades, potenciando así la eficiencia de la relajación de dos cuerpos. En consecuencia, estos sistemas permanecen ligados y se vuelven aún más concentrados de lo que serían en aislamiento.

2. Régimen Dominado por el Calentamiento (Disrupción):

   • En sistemas más extendidos (por ejemplo, S2, S5), la densidad estelar inicial es lo suficientemente baja como para que la escala de tiempo de relajación sea larga. Aquí, el calentamiento de la ULDM domina, inyectando energía que impulsa al sistema hacia un estado de dispersión.
   • Más allá de un umbral crítico, el sistema es completamente disrumpido.

3. El Análogo de "Segue 1":

   • Los autores identifican un régimen cercano a la disrupción cuyos remanentes se asemejan a las UFDs observadas como Segue 1.
   • La simulación S5, que se aproxima a la disrupción a $t \approx 10$ Gyr, produce perfiles de densidad superficial estelar y de dispersión de velocidad que coinciden ampliamente con los datos observacionales de Segue 1.
   • Un sistema con masas estelares más realistas (S6) también reproduce las propiedades de tipo Segue 1 a una edad más temprana ($\sim 6.5$ Gyr). Esto sugiere que para $m_a = 10^{-22}$ eV, Segue 1 podría interpretarse como un sistema actualmente cerca de la disrupción.

Diagrama de Fase y Cuantificación
Para generalizar estos hallazgos, los autores introducen un parámetro adimensional $\alpha$, definido como la relación entre el tiempo de relajación de semimasa ($t_{rh}$) y el tiempo de escala de disrupción característico de la ULDM ($t_{dis}$):
$$\alpha \equiv \frac{t_{rh}}{t_{dis}}$$

• $\alpha < 0.05$: La relajación domina, conduciendo al colapso del núcleo.
• $\alpha > 0.05$: El calentamiento domina, conduciendo a la disrupción.
• $t_{rh} > 4.7$ Gyr (y $t_{dis} > 85$ Gyr): Ambos efectos son insignificantes durante la edad del Universo; los sistemas permanecen estables.

Utilizando simulaciones extensas, los autores construyen un diagrama de fase bidimensional en el plano $(N, M_\star R_h^{-3})$. Este diagrama separa claramente el espacio de parámetros en regiones de colapso del núcleo, disrupción y evolución insignificante, demostrando que el destino evolutivo está determinado por los parámetros iniciales del sistema ($M_\star, N, R_h$) y la masa de la partícula de ULDM.

Significado y Reivindicaciones
El artículo sostiene que la dinámica de los sistemas estelares compactos en halos de ULDM es "rica y no trivial", desafiando la visión simplista de que el calentamiento de la ULDM siempre conduce a la expansión o que los sistemas pequeños son inmunes a los efectos de calentamiento.

• Restricciones Robustas: Los autores argumentan que el modelado preciso de la interacción entre la relajación y el calentamiento es esencial para derivar restricciones robustas sobre la ULDM. Ignorar la aceleración del colapso del núcleo o los umbrales específicos de disrupción podría conducir a límites incorrectos sobre la masa de la partícula de ULDM.
• Interpretación Observacional: Los resultados proporcionan un escenario evolutivo plausible para las UFDs como Segue 1, sugiriendo que pueden ser sistemas en un estado transitorio cerca de la disrupción en lugar de objetos en equilibrio estático.
• Limitaciones: Los autores declaran explícitamente que su análisis se realiza con una masa de partícula de ULDM y una configuración de halo fijas. Señalan que incorporar física adicional, como una función de masa inicial (IMF) estelar realista y la evolución estelar, complicaría aún más el panorama y es necesario para establecer conexiones más robustas con las observaciones.