Core-halo scaling relations in self-interacting scalar field dark matter

Enunciado del Problema
La Materia Oscura de Campo Escalar (SFDM), a menudo denominada Materia Oscura Difusa (FDM) cuando se desprecian las auto-interacciones, postula que la materia oscura consiste en bosones ultraligeros ($m \sim 10^{-22}$ eV/$c^2$). Si bien los modelos estándar de FDM explican con éxito la formación de núcleos solitónicos en halos galácticos, el impacto de las auto-interacciones (SI) en la estructura y evolución de estos halos sigue siendo una pregunta abierta. La literatura previa ha establecido relaciones de escala para FDM no interactuante (por ejemplo, relaciones masa-núcleo), pero no está claro cómo las auto-interacciones repulsivas o atractivas modifican estas relaciones, particularmente en lo que respecta a la estabilidad de los núcleos, la relación masa-radio y el potencial de colapso gravitacional en agujeros negros supermasivos. Además, si bien el régimen de Thomas-Fermi (TF) de SFDM fuertemente repulsivo ha sido estudiado utilizando aproximaciones de fluido, las simulaciones numéricas directas de fusiones de solitones en este régimen sin aproximaciones de fluido han sido escasas.

Metodología
Los autores realizan simulaciones numéricas tridimensionales de múltiples fusiones de solitones para investigar halos de SFDM en una amplia gama de intensidades y signos de auto-interacción.

• Ecuaciones Gobernantes: El estudio resuelve las ecuaciones no relativistas de Gross-Pitaevskii-Poisson (GPP), que describen la evolución de la función de onda del campo escalar $\psi$ acoplada a la gravedad. Las ecuaciones incluyen un término de auto-interacción cúbico (proporcional a la longitud de dispersión $a_s$) y un término quíntico de orden superior.
• Esquema Numérico: Los autores utilizan un método de división temporal pseudo-espectral optimizado para GPUs. Este enfoque trata los términos de auto-interacción como un potencial efectivo, permitiendo el uso de un esquema de golpe-desplazamiento-golpe (KDK). Para interacciones atractivas, se implementa un método KDK de orden superior (O(4)) para garantizar la estabilidad.
• Configuración de la Simulación: Las simulaciones involucran la fusión de $N_{sol}$ solitones (que van de 10 a 50) dentro de una caja de simulación estática y periódica de $L=100$ kpc. La masa de la partícula se fija en $m = 10^{-22}$ eV/$c^2$.
• Espacio de Parámetros: El estudio explora tres regímenes distintos:
  1. Libre (FDM): $a_s = 0$.
  2. SI Repulsiva: $a_s > 0$, incluido el régimen de Thomas-Fermi fuertemente repulsivo ($\Lambda \gg 1$).
  3. SI Atractiva: $a_s < 0$, hasta valores donde el sistema se acerca al colapso gravitacional.
• Análisis: Los autores analizan perfiles de densidad, masa del núcleo ($M_c$), radio del núcleo ($r_c$) y energía total ($E$) para derivar relaciones de escala. Comparan los resultados de la simulación con predicciones analíticas derivadas de la minimización de energía y soluciones auto-similares.

Contribuciones y Resultados Clave

1. Estructura del Núcleo y Densidad:

   • SI Repulsiva: Conduce a la formación de núcleos más masivos y extendidos con densidades centrales más bajas en comparación con el caso de FDM libre. La interacción repulsiva proporciona soporte de presión adicional, evitando que el núcleo se vuelva tan compacto como en el caso no interactuante.
   • SI Atractiva: Aumenta las densidades centrales. Para interacciones atractivas suficientemente fuertes, la densidad del núcleo diverge, indicando un colapso gravitacional. Los autores identifican una escala de masa crítica ($M_{max}$) para modelos atractivos; si la masa del núcleo supera este umbral, se pierde el equilibrio, lo que potencialmente conduce a la formación de agujeros negros supermasivos.
   • Evolución: Los núcleos repulsivos permanecen estables y extendidos durante toda la simulación ($t = 10\tau_{dyn}$), mientras que los núcleos atractivos muestran un aumento monótono en la densidad central hasta el colapso.

2. Relaciones de Escala Masa-Radio y Núcleo-Halo:

   • Relación Masa-Radio: Los resultados de la simulación para la relación masa-radio del núcleo ($M_c$ vs. $r_c$) se alinean bien con las predicciones analíticas (Chavanis 2011; Padilla et al. 2021) incluso en presencia de fuertes auto-interacciones. Las interacciones repulsivas desplazan la relación hacia arriba (permitiendo más masa para un radio dado), mientras que las interacciones atractivas la desplazan hacia abajo.
   • Relación Masa-Núcleo-Halo: Los autores investigan la relación entre la relación de masa núcleo-halo ($M_c/M$) y la cantidad invariante $\Xi = |E| / [M^3(Gm/\hbar^2)]$. Descubren que, si bien la pendiente de esta relación permanece relativamente universal (similar al caso FDM $M_c/M \propto \Xi^{1/3}$), la normalización depende de la intensidad de la auto-interacción. Esto sugiere que el mecanismo que alimenta masa al núcleo es principalmente gravitacional y universal una vez que los halos se han virializado, mientras que la SI determina localmente el soporte de presión y el radio del núcleo.
   • Relación Tamaño-Energía: La relación entre el tamaño del núcleo ($r_c$) y la energía total ($E$) es sensible a la SI. Las interacciones repulsivas resultan en pendientes más suaves (núcleos más grandes para una energía dada), mientras que las interacciones atractivas resultan en pendientes más pronunciadas.

3. Descomposición de Energía:

   • En casos repulsivos, la energía de auto-interacción domina la región del núcleo, explicando la menor densidad y el mayor tamaño.
   • En casos atractivos y libres, la energía potencial gravitacional domina el núcleo.
   • Los autores confirman que las relaciones de escala no son estrictamente universales, sino que dependen de la etapa evolutiva del halo y de los parámetros específicos de la SI.

Significado y Afirmaciones
El artículo afirma que las auto-interacciones proporcionan un mecanismo natural para regular las propiedades de los núcleos de halos de SFDM. El estudio demuestra que:

• Las relaciones de escala núcleo-halo derivadas para FDM libre son casos límite dentro de una familia más amplia de modelos de SFDM que incluyen auto-interacciones.
• Las interacciones repulsivas pueden estabilizar los núcleos contra el colapso y producir estructuras más grandes y difusas, lo cual podría tener implicaciones para el problema del "núcleo-cúspide" en la formación de galaxias.
• Las interacciones atractivas pueden impulsar los núcleos hacia la inestabilidad y el colapso, ofreciendo un canal de formación potencial para semillas de agujeros negros supermasivos en centros galácticos.
• El método pseudo-espectral se valida como una herramienta robusta para simular SFDM en el régimen de Thomas-Fermi sin recurrir a aproximaciones de fluido, capturando con éxito la transición del comportamiento solitónico al comportamiento tipo fluido.

Los autores concluyen que, si bien la dinámica gravitacional fundamental del crecimiento del núcleo es universal, las propiedades estructurales específicas (densidad, tamaño, estabilidad) son altamente sensibles al signo y magnitud de las auto-interacciones y a la historia evolutiva del halo. Estos hallazgos sugieren que las restricciones observacionales sobre los núcleos galácticos podrían utilizarse para restringir la longitud de dispersión y la masa de la partícula de la SFDM.