DESI DR2 Reference Mocks: Clustering results from UCHUU ELGs and QSOs

Resumen Técnico: Mocks de Referencia para DESI DR2 a partir de Uchuu ELGs y QSOs

Problema y Motivación
El Instrumento de Espectroscopía de Energía Oscura (DESI) tiene como objetivo mapear la distribución de materia del universo utilizando más de 50 millones de espectros de galaxias y cuásares para restringir la historia de la expansión y el crecimiento de la estructura. Para utilizar plenamente estos datos, particularmente para los próximos análisis cosmológicos de la Data Release 2 (DR2), se requieren catálogos simulados (mocks) realistas. Estos mocks deben emular las funciones de selección y las características estadísticas de DESI para probar los procesos (pipelines), evaluar las incertidumbres sistemáticas y validar las predicciones teóricas.

Un desafío específico abordado en este trabajo es el modelado de las Galaxias con Líneas de Emisión (ELGs) y los Cuásares (QSOs). A diferencia de las Galaxias Elípticas Luminosas (LRGs), que se cree que son trazadores completos en ocupación de masa de halo, las ELGs y los QSOs son trazadores incompletos en espacios de parámetros conocidos (como luminosidad o masa estelar). En consecuencia, las técnicas tradicionales de Ajuste de Abundancia de Subhalos (SHAM), que dependen de un mapeo uno a uno entre las propiedades del halo y los observables de las galaxias, son insuficientes. Los autores pretenden construir mocks de referencia que reproduzcan las estadísticas de agrupamiento (clustering) de las ELGs y QSOs de DESI DR1 y DR2 en los redshifts $0.8 < z < 2.1$ (extendiéndose hasta $z=3.5$ para QSOs) para mejorar la construcción de conos de luz simulados y potenciar la comprensión de la conexión galaxia-halo.

Metodología
El estudio utiliza la simulación N-body Uchuu, una simulación de alta resolución y gran volumen ($2 h^{-1}$Gpc de caja) que contiene $12,800^3$ partículas de materia oscura, ejecutada dentro de una cosmología $\Lambda$CDM plana basada en los parámetros de Planck 15. La simulación proporciona instantáneas (snapshots) desde $z=14$ hasta $z=0$, con halos y subhalos de materia oscura identificados mediante el buscador de halos Rockstar.

Para poblar estos halos con trazadores, los autores emplean un marco de SHAM modificado:

1. Modelado de QSOs:

   • Los QSOs se asignan a los halos basándose en la velocidad circular máxima ($V_{peak}$), definida como la $V_{max}$ máxima medida a través de todas las salidas de redshift.
   • El modelo asume una distribución gaussiana de $V_{peak}$ tanto para los QSOs centrales como para los satélites, caracterizada por una media ($V_{mean}$) y una desviación estándar ($\sigma_V$).
   • Un parámetro libre, la fracción de satélites ($f_{sat}$), controla la normalización relativa de las distribuciones central y de satélite.
   • Para dar cuenta de las grandes incertidumbres de redshift inherentes a la espectroscopía de QSOs (debido a las líneas de emisión anchas), se aplican errores de redshift gaussianos ($\sigma \approx 500$ km/s) a los catálogos simulados.
   • Los parámetros se ajustan minimizando el estadístico $\chi^2$ entre el monopolo del mock de la función de correlación de dos puntos (2PCF) y los datos de DESI en escalas de $3$a$70 h^{-1}$Mpc.

2. Modelado de ELGs:

   • Reconociendo que la formación estelar en galaxias satélite puede ser inhibida (quenched) por el barrido de presión por impacto (ram pressure stripping) si se mueven demasiado rápido en relación con su halo central, los autores introducen criterios de selección específicos para las satélites de las ELGs.
   • Las satélites se seleccionan solo si cumplen con: (i) $V_{max}/V_{peak} > 0.9$ (para evitar fusiones mayores) y (ii) bajas velocidades relativas con respecto al halo central ($V_{xyz,pec} - V_{central} < 300$ km/s) y un $V_{peak}$ central bajo ($< 300$ km/s).
   • De manera similar a los QSOs, el modelo utiliza una distribución de $V_{peak}$ gaussiana con parámetros $V_{mean}$, $\sigma_V$ y $f_{sat}$, pero los criterios de selección se aplican antes del ajuste de abundancia.

3. Construcción de Conos de Luz:

   • Los conos de luz de los mocks se construyen replicando las cajas de la simulación Uchuu y cortándolas en capas esféricas para cubrir el volumen del estudio DESI.
   • Los mocks se enmascaran para coincidir con la huella de DESI DR2 y se ponderan para reproducir la densidad numérica comóvil observada $n(z)$ en bins de redshift de $\Delta z = 0.01$.
   • Los mejores parámetros se determinan independientemente para diferentes bins de redshift y luego se parametrizan como funciones de redshift (por ejemplo, evolución lineal para el $f_{sat}$ de QSOs y evolución cuadrática para el $V_{mean}$ de ELGs).

Contribuciones Clave

• Mocks de Referencia: El artículo presenta el primer conjunto de mocks de referencia para DESI DR2 específicamente ajustados para reproducir el agrupamiento de ELGs y QSOs utilizando la simulación Uchuu.
• Implementación de SHAM Modificado: Los autores adaptan con éxito la técnica SHAM para trazadores incompletos. Para las ELGs, introducen un criterio de selección basado en la velocidad para dar cuenta de la inhibición (quenching) ambiental de la formación estelar en las satélites, una modificación que no está presente en el SHAM estándar.
• Evolución de Redshift: El estudio identifica y parametriza la evolución de redshift de los parámetros de SHAM. Específicamente, la fracción de satélites ($f_{sat}$) para los QSOs aumenta linealmente con el redshift, mientras que la velocidad media ($V_{mean}$) para las ELGs sigue una tendencia cuadrática. Esto permite una parametrización más eficiente que el ajuste de parámetros independientes en cada bin.
• Validación: Los mocks son rigurosamente validados contra los datos de DESI DR1 y DR2, reproduciendo las densidades numéricas y las estadísticas de agrupamiento en escalas de $1$a$100 h^{-1}$Mpc.

Resultados

• Estadísticas de Agrupamiento: Los mocks de SHAM modificado muestran un excelente acuerdo con los datos de DESI DR2 tanto para el monopolo como para el cuadrupolo de la función de correlación de dos puntos (2PCF) en el espacio de redshift y la función de correlación proyectada. El acuerdo se mantiene en el rango de ajuste ($3$a$70 h^{-1}$Mpc) y se extiende a la escala BAO ($\sim 100 h^{-1}$Mpc).
• Espectro de Potencia: Los mocks reproducen razonablemente bien el monopolo y el cuadrupolo del espectro de potencia, con desviaciones menores en números de onda ($k$) altos para las ELGs, lo cual es esperado ya que el espectro de potencia no fue utilizado en el proceso de ajuste.
• Distribución de Ocupación de Halo (HOD):
  • ELGs: La HOD consiste en un componente gaussiano para las centrales (con un pico en $M_{halo} \approx 10^{12} h^{-1}M_{\odot}$ con una ocupación de 0.095) y un componente de ley de potencia para las satélites (que domina en $M_{halo} > 10^{12.4} h^{-1}M_{\odot}$).
  • QSOs: La HOD muestra de manera similar un componente central gaussiano (con un pico en $M_{halo} \approx 10^{12.4} h^{-1}M_{\odot}$ con una ocupación de 0.319) y un componente de satélite de ley de potencia (que domina en $M_{halo} > 10^{12.8} h^{-1}M_{\odot}$).
• Sesgo (Bias) a Gran Escala: El sesgo a gran escala ($b$) medido de los mocks coincide con el sesgo derivado de los datos de DESI DR2 en el rango de $10$a$80 h^{-1}$Mpc. Ambos trazadores muestran un aumento monotónico del sesgo con el redshift, consistente con la ocupación de halos más masivos en tiempos cósmicos tempranos.

Significancia y Reivindicaciones
El artículo afirma que estos resultados proporcionan un puente crítico entre la teoría y la observación para los análisis de DESI. Al reproducir las densidades numéricas observadas y los patrones de agrupamiento, estos mocks sirven como herramientas esenciales para:

1. Estimación de la Matriz de Covarianza: Pueden utilizarse para generar matrices de covarianza para futuros análisis cosmológicos de DESI.
2. Pruebas de Sistemas: Permiten la prueba de los procesos (pipelines) cosmológicos y la evaluación de las incertidumbres sistemáticas inherentes al estudio.
3. Perspectiva Física: La aplicación exitosa de las técnicas de SHAM modificado ofrece una mejor comprensión de la conexión galaxia-halo, particularmente respecto a la dependencia ambiental de la formación estelar en las galaxias satélites (para las ELGs) y la evolución de las poblaciones de QSOs.

Los autores enfatizan que, aunque los modelos están ajustados a los datos de DESI, están construidos sobre la cosmología de Planck y la simulación Uchuu, asegurando que sean físicamente consistentes con los estándares cosmológicos actuales. El trabajo no pretende descubrir nueva física, sino proporcionar la infraestructura necesaria (el mock realista) para permitir una inferencia cosmológica precisa a partir del estudio DESI.