The impact of source and survey modelling on the connection between [O III] emitters and Ly $\alpha$ forest transmission at z ~ 6

Resumen Técnico: El impacto del modelado de fuentes y de sondeos en la conexión entre los emisores de [O iii] y la transmisión del bosque de Ly $\alpha$ a $z \sim 6$

Planteamiento del problema
Observaciones recientes del Telescopio Espacial James Webb (JWST) de galaxias emisoras de [O iii] han revelado una correlación cruzada con la transmisión del bosque de Ly $\alpha$ que presenta un desafío para los modelos numéricos existentes. Específicamente, las mediciones indican un exceso en la transmisión de Ly $\alpha$ intergaláctica en escalas de $\sim 10\text{--}40$ cMpc de las galaxias. Aunque los modelos teóricos predicen tal exceso impulsado por una ionización mejorada alrededor de las galaxias, reproducir la escala y la amplitud observadas ha resultado difícil. El trabajo teórico previo a menudo carecía de un modelado detallado de la población de emisores de [O iii] o de las restricciones geométricas específicas y la profundidad de los sondeos de JWST (como ASPIRE y EIGER), lo que potencialmente derivó en discrepancias entre la simulación y la observación.

Metodología
Los autores utilizan la suite de simulaciones hidrodinámicas cosmológicas de alta resolución Sherwood-Relics, que emplea un esquema de transferencia de radiación híbrido para modelar la reionización fragmentada y la reionización completa a $z=5.3$. Para cerrar la brecha entre los halos de la simulación y las galaxias observadas, los autores construyen un nuevo modelo empírico que incorpora los siguientes pasos:

1. Emparejamiento de abundancia (Abundance Matching): Mapean la función de masa de halos de materia oscura de la simulación a la función luminosidad UV observada de las galaxias a $z \sim 6$ (basándose en Bouwens et al. 2021) utilizando la técnica de emparejamiento de abundancia.
2. Conversión de luminosidad: Las luminosidades UV se convierten en luminosidades de [O iii] utilizando una relación de escala ajustada a los datos de JWST (Matthee et al. 2023). Esta relación tiene en cuenta la dispersión en la razón $L_{[\text{O III}]}/L_{\text{UV}}$ y aplica un límite de sensibilidad de flujo ($1.0 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$) para imitar los límites de detección de los sondeos.
3. Modelado de sondeo hacia adelante (Survey Forward Modeling): Los autores generan catálogos de sondeos simulados que incorporan explícitamente la geometría y la profundidad de los sondeos ASPIRE y EIGER. Esto incluye la simulación de líneas de visión de cuásares específicos, restricciones de campo de visión y el número de emisores detectados por línea de visión.
4. Análisis estadístico: Generan un conjunto de 1,024 realizaciones para sondeos tipo ASPIRE y 2,500 realizaciones para sondeos tipo Zhu et al. (2026) para cuantificar la varianza en las mediciones de correlación cruzada. Comparan estos modelos contra datos observacionales para:
   • Agrupamiento de galaxias (autocorrelación proyectada promediada por volumen).
   • Correlación cruzada unidimensional de la transmisión de la galaxia–Ly $\alpha$.
   • Densidad superficial de galaxias en dos dimensiones frente a la profundidad óptica efectiva de Ly $\alpha$ ($\tau_{\text{eff}}^{50}$).

Resultados clave

• Agrupamiento de galaxias: El modelo reproduce el agrupamiento de los emisores de [O iii] (Huang et al. 2026; Eilers et al. 2024) con un acuerdo razonable ($< 1\text{--}2\sigma$) al aplicar un corte de luminosidad de $\log_{10}(L_{[\text{O III}]}/\text{erg s}^{-1}) \ge 42.4$. Las discrepancias en escalas pequeñas ($r_\perp \lesssim 0.4$ cMpc) se atribuyen al "término de un solo halo" (one-halo term) ausente en el buscador de halos Friends-of-Friends de la simulación, el cual no puede albergar múltiples galaxias dentro de un solo halo.
• Correlación cruzada unidimensional: El modelado empírico de la fuente y del sondeo conduce a un ligero aumento en la masa media del halo anfitrión ($\langle \log_{10}(M_h/M_\odot) \rangle \approx 10.7$) en comparación con los enfoques previos de corte de masa más simples. Sin embargo, esto no resuelve la discrepancia de escala: el pico de la correlación cruzada en los modelos ocurre a $\sim 20$ cMpc, mientras que las observaciones muestran un pico a $\sim 30$ cMpc.
• Gran dispersión: Un hallazgo principal es la muy grande dispersión en las mediciones de los sondeos simulados de la correlación cruzada, particularmente en los redshifts más altos ($z=6.4$). Esta dispersión surge de la rareza de los picos de transmisión de Ly $\alpha$ y el número limitado de líneas de visión en los sondeos reales. En consecuencia, la gran varianza permite que el modelo logre un ajuste estadísticamente bueno con los datos observacionales (incluyendo la alta amplitud vista en Kakiichi et al. 2025) a pesar del desajuste en la ubicación del pico.
• Correlación cruzada bidimensional: El modelo muestra un buen acuerdo estadístico con la relación entre la densidad superficial de galaxias y la profundidad óptica efectiva de Ly $\alpha$ (Zhu et al. 2026). Los resultados indican que tanto las líneas de visión más transparentes ($\tau_{\text{eff}}^{50} \lesssim 3$) como las más opacas ($\tau_{\text{eff}}^{50} \gtrsim 5$) pueden coincidir con subdensidades de galaxias.
  • Las líneas de visión transparentes corresponden a regiones subdensas y calientes que han sido reionizadas recientemente.
  • Las líneas de visión opacas corresponden a regiones más frías con una mayor probabilidad de contener islas neutras ($x_{\text{HI}} > 0.5$ sobre $\ge 2$ cMpc).

Significancia y conclusiones
El artículo argumenta que la gran dispersión inherente a las observaciones actuales de galaxias e IGM, impulsada por tamaños de muestra limitados y la geometría de los sondeos, implica que los datos presentes pueden tener dificultades para descartar una amplia gama de modelos de fuentes ionizantes. Si bien el modelado detallado de la selección de fuentes y la geometría del sondeo es crítico para evaluar con precisión la varianza de la medición, esto no altera fundamentalmente la forma de la correlación cruzada ni resuelve la discrepancia de escala entre modelos y observaciones.

Los autores concluyen que el progreso futuro requiere:

1. Aumento de los tamaños de muestra observacional de los programas en curso como EIGER y ASPIRE.
2. Simulaciones en volúmenes mayores ($> 250$ cMpc) para capturar mejor el extremo de alta masa de la función de masa de halos y reducir la varianza cósmica.
3. Modelos de fuentes diversos, incluyendo núcleos galácticos activos tenues, para probar la robustez de las restricciones de la historia de la reionización.

El estudio refuerza que la relación entre las galaxias y el medio intergaláctico al final de la reionización es compleja, donde la densidad de galaxias no escala simplemente con el estado de ionización del gas.