Constraining Dark Acoustic Oscillations with the High-Redshift UV Luminosity Function

Resumen Técnico: Restricciones de las Oscilaciones Acústicas Oscuras con la Función de Luminosidad UV de alto desplazamiento al rojo

Planteamiento del Problema
Aunque el modelo estándar $\Lambda$CDM postula que la materia oscura (DM) es fría y sin colisiones, diversas extensiones bien motivadas del sector oscuro predicen "Oscilaciones Acústicas Oscuras" (DAOs, por sus siglas en inglés). Estas surgen cuando una fracción de la DM interactúa con una especie de radiación oscura (DR) en el universo temprano, lo que provoca una supresión y oscilaciones amortiguadas en el espectro de potencia de la materia $P(k)$ en escalas pequeñas. Las restricciones existentes de la Radiación de Fondo de Microondas (CMB) son estrictas en escalas grandes, pero limitan la fracción de DM interactuante a $\lesssim 5\%$ para DAOs fuertes. Por el contrario, las restricciones del bosque Lyman-$\alpha$ son efectivas para DAOs débiles o modelos totalmente interactuantes, pero enfrentan desafíos para distinguir las firmas de las DAOs fuertes de los efectos de la Materia Oscura Cálida (WDM). Existe la necesidad de obtener restricciones agnósticas al modelo sobre la estructura de pequeña escala del universo, particularmente en desplazamientos al rojo altos ($z \sim 3-9$), donde la evolución no lineal de la estructura y la impronta específica de las DAOs en la función de masa de halos (HMF) pueden ser sondeadas.

Metodología
Los autores desarrollan un flujo de trabajo para restringir los parámetros de las DAOs utilizando la Función de Luminosidad Ultravioleta (UVLF) de galaxias de alto desplazamiento al rojo. La metodología consta de cuatro componentes:

1. Modelo de Función de Transferencia Fenomenológica:
   Para evitar el gasto computacional de resolver ecuaciones de Boltzmann para cada punto de un escaneo de Monte Carlo por Cadenas de Markov (MCMC), los autores introducen un nuevo modelo de función de transferencia fenomenológica de cuatro parámetros $T(k)$ relativo a $\Lambda$CDM. Este modelo captura la supresión y las oscilaciones amortiguadas de las DAOs y permite una fracción subunitaria de DM interactuante ($f$). Los parámetros son:

   • $f$: la fracción de DM que experimenta DAOs.
   • $k_{\text{peak}}$: el número de onda del primer pico de la DAO.
   • $k_{\text{damp}}$: la escala de amortiguación de las oscilaciones.
   • $A$: la amplitud de las oscilaciones.
     El modelo es calibrado contra espectros de potencia generados por el código de Boltzmann modificado CLASS-aDM (para Materia Oscura Atómica) y CLASS (para Materia Oscura Interactuante), demostrando que los parámetros se mapean directamente a cantidades físicas como el horizonte de sonido oscuro y la escala de amortiguación de difusión.

2. Formalismo de Press-Schechter Extendido (EPS) y Calibración mediante N-cuerpos:
   Para predecir la HMF a partir del espectro de potencia lineal, los autores emplean el formalismo de Press-Schechter Extendido. Reconociendo que las funciones de ventana estándar (por ejemplo, top-hat en el espacio real) fallan para espectros de potencia truncados, utilizan una función de ventana de "k suave" (smooth-k). Crucialmente, realizan una suite de 13 simulaciones de N-cuerpos de solo DM utilizando el código GÍZMO con condiciones iniciales MUSIC. Estas simulaciones cubren varias fracciones de DM interactuante ($f=0.1, 0.5, 1$) y escalas de DAO ($k_{\text{DAO}} = 5, 20, 66, 190 \, h/\text{Mpc}$). Las simulaciones se utilizan para calibrar los parámetros de EPS ($q, \beta, c$) para asegurar que la HMF semianalítica reproduce con precisión los resultados de las simulaciones en desplazamientos al rojo $z=4$ a $12$.

3. Conexión Halo-Galaxia:
   La conexión entre la HMF y la observable UVLF se modela utilizando el código GALLUMI. Este flujo de trabajo asume que la luminosidad UV es impulsada por tasas de formación estelar (SFR), que en promedio siguen una relación de ley de potencia doble con la tasa de acreción del halo. El modelo incluye ocho parámetros de molestia astrofísicos (que describen la normalización, la escala de pivote y la evolución de la relación SFR-masa del halo, así como la estocasticidad en la magnitud), los cuales son marginalizados para contabilizar las incertidumbres en la física de formación de galaxias de alto desplazamiento al rojo.

4. Análisis de Datos:
   Los autores combinan mediciones de la UVLF del Telescopio Espacial Hubble (HST, campos tanto en blanco como con lentes), el Telescopio Espacial James Webb (JWST), el Telescopio Subaru y el Canadá-Francia-Hawái (CFHT) cubriendo desplazamientos al rojo $z=3$ a $9$. Realizan escaneos Bayesianos MCMC y análisis de verosimilitud de perfil frecuentista para restringir los parámetros de las DAOs, fijando los parámetros de $\Lambda$CDM a los valores óptimos de Planck 2018.

Contribuciones Clave

• Nuevo Modelo Fenomenológico: Introducción de un modelo de función de transferencia flexible que acomoda fracciones subdominantes de DM interactuante (como la Materia Oscura Atómica), cerrando la brecha entre modelos específicos de física de partículas y restricciones generales.
• Calibración de Simulaciones: Primera calibración del formalismo EPS para modelos de DAO utilizando un conjunto dedicado de simulaciones de N-cuerpos, validando el enfoque semianalítico para predecir la HMF en presencia de oscilaciones amortiguadas.
• Restricciones Agnósticas al Modelo: Derivación de restricciones en el espacio de parámetros de las DAOs ($f$ vs. $k_{\text{peak}}$) que son independientes de la física de partículas específica del sector oscuro, permitiendo la reinterpretación fácil en varios modelos (por ejemplo, IDM, Materia Oscura Atómica).

Resultados

• Restricciones sobre $f$ y $k_{\text{peak}}$: El análisis arroja una fuerte correlación entre la fracción de DM interactuante y la escala de las DAOs.
  • Si la fracción de DM interactuante es significativa ($f \gtrsim 0.07$), el número de onda del primer pico de la DAO debe ser $k_{\text{peak}} \gtrsim 50 \, h/\text{Mpc}$ para permanecer consistente con los datos de la UVLF.
  • Por el contrario, si $k_{\text{peak}}$ es pequeño (indicando un desacoplamiento más temprano o interacciones más fuertes), la fracción de DM interactuante está restringida a $f \lesssim 0.07$.
• Comparación con Otros Sondas:
  • Para fracciones altas de DM interactuante ($f \sim 1$), las restricciones de la UVLF mejoran las cotas derivadas de la CMB por aproximadamente un orden de magnitud en $k_{\text{peak}}$.
  • Para escalas muy pequeñas ($k_{\text{peak}} \lesssim 1 \, h/\text{Mpc}$), las restricciones de la CMB siguen siendo más estrictas sobre la abundancia de DM interactuante ($f \lesssim 5\%$).
  • Las restricciones sobre la sección eficaz de dispersión DM-barión en modelos de Materia Oscura Interactuante (IDM) resultan ser más de cinco órdenes de magnitud más fuertes que las restricciones de la CMB y competitivas con los límites derivados de las abundancias de satélites de la Vía Láctea.
• Parámetros No Restringidos: El análisis no encuentra restricciones significativas sobre la amplitud de la DAO ($A$) o la relación de la escala de amortiguación ($k_{\text{damp}}/k_{\text{peak}}$) al marginalizar sobre otros parámetros, ya que la sensibilidad de la UVLF está impulsada principalmente por la supresión inicial del espectro de potencia más que por las características oscilatorias.

Significancia
Este artículo demuestra que las mediciones de la función de luminosidad ultravioleta de alto desplazamiento al rojo proporcionan una sonda poderosa y complementaria a los datos de la CMB y del bosque Lyman-$\alpha$ para restringir la física de la materia oscura a pequeña escala. Al calibrar con éxito el formalismo EPS contra simulaciones de N-cuerpos y emplear un modelo flexible de conexión halo-galaxia, los autores establecen límites robustos sobre la fracción de la materia oscura que puede experimentar oscilaciones acústicas. Estos resultados empujan la época de desacoplamiento de la materia oscura interactuante a tiempos más tempranos (mayores $k_{\text{peak}}$) o requieren que la fracción interactuante sea subdominante ($f < 7\%$) para escalas de interacción específicas. El trabajo destaca la utilidad de los observables de alto desplazamiento al rojo para distinguir entre diferentes paradigmas del sector oscuro y establece una base para futuras restricciones utilizando los datos provenientes del Observatorio Vera Rubin, el Telescopio Espacial Nancy Grace Roman y el Square Kilometer Array.