Evolution of Compact Stellar Systems in Ultralight Dark Matter Halos: Dependence on Stellar and Dark Matter Parameters

Resumen Técnico: Evolución de Sistemas Estelares Compactos en Halos de Materia Oscura Ultraligera

Planteamiento del Problema
Los sistemas estelares compactos, como las galaxias enanas ultra débiles (UFD), se utilizan frecuentemente para restringir la masa de partícula de la materia oscura ultraligera (ULDM). La base teórica para estas restricciones se basa en el "efecto de calentamiento", donde la interferencia de ondas en el halo de ULDM transfiere energía a los sistemas estelares incrustados, provocando su expansión y potencialmente contradiciendo los datos observacionales. Sin embargo, la literatura reciente presenta conclusiones contradictorias respecto a la viabilidad de los modelos de ULDM. Algunos estudios sugieren que el despojo de marea (tidal stripping) por parte de la Vía Láctea (MW) suprime este calentamiento, permitiendo masas de ULDM más bajas, mientras que otros argumentan que el calentamiento sigue siendo significativo para masas más altas. Además, el papel de la dinámica estelar interna —específicamente la relajación de dos cuerpos y la evolución estelar (SE)— en la mediación de estas interacciones ha sido ampliamente ignorado. Este trabajo investiga el destino a largo plazo de las poblaciones estelares compactas en halos de ULDM, centrándose en la interacción entre la metalicidad estelar, el campo de marea de la MW y la masa de la partícula de ULDM, particularmente en el régimen donde el tamaño del sistema ($R_h$) es mucho menor que la longitud de onda de de Broglie de la ULDM ($\lambda_{dB}$).

Metodología
Los autores emplean una estrategia de simulación de dos etapas utilizando simulaciones directas de $N$-cuerpos acopladas con la evolución de la función de onda de la ULDM:

1. Evolución del Halo de ULDM: El halo de ULDM se modela mediante el método de descomposición de autoestados y se evoluciona a través del método pseudo-espectral (paquete PyUltraLight) en una caja cúbica de 40 kpc. El sistema estelar se trata inicialmente como una partícula puntual para registrar el tensor de marea de la ULDM, $\nabla\nabla V_{\text{ULDM}}$, en su posición.
2. Evolución del Sistema Estelar: La evolución interna del sistema estelar se simula utilizando el código de $N$-cuerpos directos PeTar. El sistema sigue un perfil Plummer con $N=5 \times 10^4$ partículas y un radio de semiluz inicial $R_h = 5$ pc.
3. Mecanismo de Acoplamiento: La influencia de la ULDM se incorpora mediante una aproximación de tensor de marea de primer orden. La aceleración de marea $\Delta a(r, r_c) \simeq -(r - r_c) \cdot \nabla\nabla V_{\text{ULDM}}(r_c)$ se añade a la aceleración de cada partícula estelar relativa al centro de masa del sistema.
4. Conjuntos de Simulación: El estudio realiza 11 ejecuciones de simulación variando:
   • Metalicidad Estelar ($Z$): $10^{-4}$ vs. $10^{-2}$.
   • Evolución Estelar (SE) e IMF: SE completa con IMF, SE sin IMF, sin SE, y partículas de masa igual.
   • Contexto Ambiental: Evolución aislada vs. evolución dentro del campo de marea de la MW (incluyendo el despojo de marea del halo externo de ULDM).
   • Masa de Partícula de ULDM ($m_{22} \equiv m_a/10^{-22}$ eV): Valores que oscilan entre 1 y 3.
   • Dinámica Orbital: Centro de masa fijo vs. movimiento orbital libre.

Contribuciones Clave y Resultados

• Impacto de la Metalicidad y Evolución Estelar:
  El estudio encuentra que la metalicidad estelar altera significativamente la resistencia del sistema a la disrupción. Los sistemas con mayor metalicidad ($Z=10^{-2}$) producen agujeros negros (BH) de menor masa debido a vientos estelares más fuertes. Esto reduce la eficiencia de la segregación de masas, retrasando la transferencia de estrellas de baja masa hacia las periferias del sistema y demorando la disrupción. Por el contrario, los sistemas donde la SE se desactiva exhiben una segregación de masas intensificada y una disrupción rápida. Los sistemas con partículas de masa igual (sin SE, sin IMF) no se disrupcion, sino que experimentan colapso del núcleo.

• Papel del Campo de Marea de la Vía Láctea y el Movimiento Orbital:
  Contrario a los hallazgos en el régimen $\lambda_{dB} \lesssim R_h$ (donde el despojo de marea de la MW suprime el calentamiento al eliminar estados de alta energía de la ULDM), los autores encuentran que en el régimen $\lambda_{dB} \gg R_h$, el campo de marea de la MW potencia la disrupción. El campo de marea de la MW despoja el halo externo de ULDM, alterando el movimiento orbital del sistema estelar dentro del halo restante. Esta modificación aumenta la frecuencia de los pasos por el pericentro (el punto de mayor cercanía al solitón), lo que conduce a eventos de calentamiento por marea más frecuentes e intensos. Los sistemas que evolucionan en aislamiento (S4) permanecen estables más tiempo que aquellos sujetos al despojo de marea de la MW (S0).

• Dominancia del Calentamiento Inducido por el Solitón:
  Al comparar las simulaciones con fluctuaciones completas de ULDM frente a aquellas que utilizan solo un potencial de solitón estático (S0-soliton, S4-soliton), los autores demuestran que el efecto de calentamiento es impulsado primordialmente por la influencia de marea del solitón, más que por las estructuras granulares fluctuantes del halo de ULDM. Se produce una pérdida significativa de partículas estelares cuando el sistema estelar es desplazado más allá del radio del solitón ($d > r_c$), donde experimenta el potencial gravitatorio fuertemente inhomogéneo del solitón.

• Dependencia de la Masa de la Partícula de ULDM:
  En el régimen $\lambda_{dB} \gg R_h$, el efecto de calentamiento se vuelve más fuerte con el aumento de la masa de la partícula de ULDM. Esta tendencia es opuesta a la observada en el régimen $\lambda_{dB} \lesssim R_h$. Los autores atribuyen esto a la relación núcleo-halo del solitón: a medida que la masa de la partícula aumenta, la masa del solitón aumenta y su radio disminuye, lo que resulta en un campo de marea más fuerte ($GM_c/r_c^3 \propto m_{22}^2$). En consecuencia, los sistemas estelares se disuelven más rápidamente para masas de partículas de ULDM mayores en este régimen.

Significancia
El artículo concluye que la evolución de los sistemas estelares compactos en halos de ULDM es altamente compleja y cualitativamente diferente dependiendo de la relación entre el tamaño del sistema y la longitud de onda de de Broglie. Los hallazgos sugieren que las restricciones existentes sobre la masa de la partícula de ULDM derivadas de sistemas compactos (como Segue 1) pueden requerir una revisión cuidadosa. Específicamente, la interacción entre la metalicidad estelar, los procesos de relajación interna y la historia orbital específica dentro del campo de marea de la MW puede alterar significativamente los tiempos de supervivencia de estos sistemas. El estudio enfatiza que los regímenes $\lambda_{dB} \gg R_h$ y $\lambda_{dB} \lesssim R_h$ deben tratarse por separado al derivar restricciones astrofísicas de la ULDM, ya que los mecanismos físicos que impulsan el calentamiento y la disrupción difieren fundamentalmente entre ellos.