Multi-species Dark Matter with Warmth and Randomness

Resumen Técnico: Materia Oscura Multi-especie con Calidez y Aleatoriedad

Enunciado del Problema
Las restricciones observacionales actuales sobre el espectro de potencia de la materia (PS) son consistentes con condiciones iniciales adiabáticas casi invariantes de escala y un crecimiento independiente de la escala en escalas subhorizonte, lo que respalda los modelos estándar de Materia Oscura Fría (CDM). Sin embargo, desviaciones de este comportamiento —específicamente la supresión dependiente de la escala (debida a la "calidez" o dispersión de velocidades) o la mejora (debida a fluctuaciones de Poisson/ruido de disparo provenientes de bajas densidades numéricas)— podrían revelar las propiedades microscópicas del sector oscuro. Aunque existen marcos para materia oscura de una sola especie con calidez o ruido de Poisson, y para perturbaciones adiabáticas multi-especie, falta un marco analítico general que maneje simultáneamente:

1. Sistemas multi-especie: Números arbitrarios de componentes de materia oscura con fracciones de masa distintas, distribuciones de velocidad y densidades numéricas.
2. Efectos de discreción: Fluctuaciones de Poisson que surgen de poblaciones escasas (por ejemplo, Agujeros Negros Primordiales, solitones, minicúmulos).
3. Calidez: Dispersión de velocidad finita (flujo libre) que suprime la agrupación a pequeña escala.

Los métodos existentes a menudo tratan estos efectos de forma aislada o dependen en gran medida de simulaciones de $N$-cuerpos, las cuales son computacionalmente costosas y difíciles de generalizar a través de espacios de parámetros arbitrarios.

Metodología
Los autores desarrollan un marco analítico general para la evolución lineal de las perturbaciones de densidad en un sistema de materia oscura multi-especie. El enfoque se basa en los siguientes pilares:

• Jerarquía BBGKY truncada: Partiendo de la ecuación de Liouville para un sistema gravitacionalmente agrupado multi-especie, los autores derivan la evolución del espectro de potencia truncando la jerarquía de Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) a nivel de dos partículas. Esto descarta correlaciones de tres partículas y de orden superior, restringiendo la validez al orden lineal en perturbaciones.
• Ecuaciones integrales de Volterra: La solución de la jerarquía truncada se expresa en términos de tres familias de funciones de transferencia ($T^{(a)}_k, T^{(b)}_k, T^{(c)}_k$) gobernadas por ecuaciones integrales de Volterra acopladas. Estas ecuaciones incorporan los núcleos de flujo libre derivados de las distribuciones de velocidad iniciales de todas las especies.
• Condiciones iniciales generalizadas: El marco acomoda distribuciones iniciales arbitrarias en el espacio de fases. Distingue entre:
  • Modos adiabáticos: Condiciones iniciales correlacionadas (por ejemplo, provenientes de la inflación).
  • Modos isocurvatura: Fluctuaciones de Poisson no correlacionadas que surgen de la naturaleza discreta de los componentes escasos.
• Implementación numérica: Se proporciona un algoritmo eficiente para resolver las ecuaciones de Volterra numéricamente. Los autores también suministran código público para evaluar los espectros de potencia totales, inter-especie e intra-especie.

Resultados Clave
El artículo presenta expresiones analíticas y soluciones numéricas para la evolución de los espectros de potencia en escenarios de materia oscura de dos componentes, validados contra simulaciones de $N$-cuerpos:

1. Espectros de potencia totales y cruzados: El espectro de potencia total del contraste de densidad $P_\delta$ se descompone en contribuciones adiabáticas e isocurvatura, cada una evolucionando mediante funciones de transferencia específicas. El marco también proporciona expresiones para espectros de potencia cruzada ($P_{SS'}$) y espectros intra-especie ($P_{SS}$), permitiendo el análisis de cómo los componentes subdominantes siembran perturbaciones en los componentes dominantes.
2. Estudios de caso:
   • Dominante Fría + Subdominante Caliente/Poisson: Cuando una pequeña fracción de materia oscura es caliente y escasa, genera ruido de Poisson. El componente frío dominante suprime el crecimiento a escala de Jeans de los modos isocurvatura del componente caliente, lo que conduce a una escala $k^{-1}$ en tiempos tardíos en lugar de la supresión más pronunciada observada en modelos de una sola especie caliente. El espectro adiabático muestra una supresión superficial proporcional a la fracción de masa $f_S$.
   • Dominante Caliente + Subdominante Fría/Poisson: Un componente frío subdominante con ruido de Poisson crece de manera invariante de escala dentro del fondo caliente, dominando eventualmente el espectro de potencia a pequeñas escalas donde el componente caliente está suprimido.
3. Validación: Los resultados analíticos fueron validados utilizando simulaciones de $N$-cuerpos (usando Gadget 4) de un modelo de dos componentes (CDM dominante y partículas calientes subdominantes con ruido de Poisson). Las predicciones analíticas coinciden bien con los resultados de la simulación dentro del régimen lineal. Las desviaciones aparecen cuando el componente subdominante entra en el régimen no lineal ($k^3 P_{22}/2\pi^2 > 1$), incluso si el componente dominante permanece lineal.
4. Contexto observacional: El artículo mapea los parámetros teóricos (fracción de masa $f_S$, dispersión de velocidad $\sigma_{eq,S}$ y densidad numérica $\bar{n}_S$) a escalas observables. Demuestra que diferentes combinaciones de calidez y ruido de Poisson pueden compensarse parcialmente, potencialmente enmascarando desviaciones del $\Lambda$CDM, o crear firmas distintivas detectables por sondeos de galaxias y observaciones del bosque Lyman-$\alpha$.

Significado y Afirmaciones
Los autores afirman que este trabajo proporciona el primer marco analítico general capaz de tratar discreción (fluctuaciones de Poisson) y calidez (dispersión de velocidad) simultáneamente en un contexto de materia oscura multi-especie.

• Generalidad: El marco no está restringido a modelos de partículas específicos; se aplica a mezclas de componentes fríos/calientes, agujeros negros primordiales, solitones y estructuras de interferencia de ondas, siempre que se traten como partículas puntuales clásicas en escalas mayores que su estructura interna.
• Eficiencia: El método ofrece una alternativa computacionalmente eficiente a las simulaciones de $N$-cuerpos para explorar el espacio de parámetros de la materia oscura multi-especie, capaz de calcular espectros de potencia en segundos.
• Limitaciones: Los autores señalan explícitamente que el marco está restringido a:
  • Teoría de perturbaciones lineal (escalas subhorizonte).
  • Tratamiento de partículas puntuales clásicas (ignorando efectos de tamaño finito como las escalas de De Broglie, que se abordan en un artículo complementario para materia oscura de onda).
  • Solo interacciones gravitacionales (se ignoran las interacciones no gravitacionales).

El artículo concluye que este marco permite el estudio sistemático de cómo los componentes subdominantes, escasos o calientes influyen en el crecimiento de la estructura cósmica, ofreciendo una herramienta para interpretar las restricciones observacionales actuales y futuras sobre la complejidad del sector oscuro.