Assembly bias and local Primordial non-Gaussianity from DESI DR1 Quasars

Título: Sesgo de ensamblaje y no-gaussianidad primordial local a partir de los Cuásares de DESI DR1

1. El Problema

La medición de la no-gaussianidad primordial local (PNG), parametrizada por el parámetro $f_{\text{NL}}$, es una de las pruebas más potentes para distinguir entre diferentes modelos de inflación cósmica. La señal clave para restringir $f_{\text{NL}}$ en estudios de estructura a gran escala (como los de quasares) es el sesgo dependiente de la escala inducido en la agrupación de halos de materia oscura.

Sin embargo, existe una degeneración completa entre $f_{\text{NL}}$ y el parámetro de respuesta del sesgo, $b_\phi$. La relación teórica estándar (relación de universalidad) se expresa como:
$$b_\phi = 2\delta_c(b_1 - p)$$
donde $b_1$ es el sesgo lineal, $\delta_c \approx 1.686$ es el umbral de colapso esférico y $p$ es un parámetro que depende de la historia de formación del halo.

• El problema central es que el valor de $p$ no es universal para todas las muestras de trazadores. Para una selección basada puramente en masa, $p=1$. Sin embargo, los cuásares se forman preferentemente en halos que han sufrido fusiones recientes, lo que implica un valor de $p > 1$ (típicamente $p \approx 1.6$ en modelos simplificados).
• En el análisis previo de DESI (DR1), esta incertidumbre obligó a reportar dos restricciones separadas para $f_{\text{NL}}$ (una asumiendo $p=1$ y otra $p=1.6$), dejando la medición final degenerada y menos precisa.

2. Metodología

El objetivo del trabajo es romper esta degeneración construyendo un prior físico robusto para el parámetro $p$ basado en simulaciones hidrodinámicas de alta resolución, en lugar de asumir valores arbitrarios.

• Simulaciones Utilizadas:

  • IllustrisTNG (TNG300-1): Se utilizó para crear una muestra de cuásares que imita los criterios de selección de DESI. Se seleccionaron halos con masa similar a la inferida para DESI ($M \sim 4.6 \times 10^{12} h^{-1}M_\odot$) y se identificaron los cuásares basándose en su tasa de acreción (relación de Eddington, $\chi$).
  • CAMELS: Se utilizó la suite de simulaciones CAMELS (conjuntos CV y 1P) para validar la robustez del resultado frente a la varianza cósmica y a las incertidumbres en la física subgrid (modelado de retroalimentación de AGN, crecimiento de agujeros negros, etc.).

• Procedimiento de Calibración:

  1. Se seleccionaron cuásares en las simulaciones siguiendo la densidad numérica y los cortes de masa de DESI.
  2. Se rastrearon los árboles de fusión de los halos anfitriones para determinar su redshift de formación ($z_f$).
  3. Se aplicó el formalismo de Press-Schechter Extendido (ePS) para calcular el parámetro $p$ para cada halo, utilizando la variable $\omega_f$ que codifica la historia de acreción.
  4. Se construyó una distribución empírica $\Pi(p)$ combinando los resultados de IllustrisTNG y CAMELS.
  5. Se probó la robustez de esta distribución introduciendo "ruido" en la selección (variando los umbrales de luminosidad) para simular sesgos observacionales.

• Aplicación a Datos Reales:

  • Se aplicó el prior $\Pi(p)$ a las cadenas de Markov (MCMC) existentes del análisis de DESI DR1.
  • Se realizó una marginalización sobre la distribución de $p$ para obtener una única restricción posterior para $f_{\text{NL}}$, en lugar de fijar $p$ a un valor único.

3. Contribuciones Clave

• Determinación Física de $p$: Por primera vez, se proporciona una distribución de probabilidad para el parámetro de sesgo de ensamblaje $p$ específica para la muestra de cuásares de DESI, derivada de simulaciones hidrodinámicas realistas en lugar de suposiciones teóricas simplificadas.
• Validación de Robustez: Se demuestra que el valor medio de $p$ es estable frente a:
  • Variaciones en los criterios de selección (ruido en la selección de objetivos).
  • Cambios en la física subgrid (parámetros de retroalimentación de AGN y crecimiento de agujeros negros) mediante el uso de CAMELS.
  • Varianza cósmica.
• Unificación de Restricciones: Se logra transformar la doble restricción degenerada de DESI DR1 en una única medición física y coherente.

4. Resultados Principales

• Distribución de $p$: La distribución de $p$ para los cuásares de DESI no es un valor fijo, sino una distribución sesgada hacia valores altos.
  • Valor medio: $\bar{p} \simeq 1.41$ (en IllustrisTNG) y $\bar{p} \simeq 1.38$ (en CAMELS tras variar parámetros astrofísicos).
  • Esto confirma que los cuásares residen en halos con historias de fusión reciente, pero el valor no alcanza el límite superior teórico de $1.6$ para fusiones extremas, sino que se sitúa en un valor intermedio debido a la mezcla de historias de formación.
  • La dependencia con el redshift es débil en el rango relevante para DESI ($z \sim 1.5 - 2.1$).
• Nueva Restricción sobre $f_{\text{NL}}$: Al aplicar este prior al conjunto de datos DESI DR1, se obtiene una nueva restricción unificada:
  $$f_{\text{NL}} = -3.3 \pm 9.2$$
• Comparación con DESI DR1 original:
  • La nueva restricción es ligeramente más ajustada (menor incertidumbre) que las restricciones individuales reportadas anteriormente cuando se asumen $p=1$ o $p=1.6$, debido a que el valor de $b_\phi$ se optimiza correctamente.
  • Se demuestra que fijar $p$ a su valor medio ($1.4$) durante la inferencia produce resultados estadísticamente indistinguibles de marginalizar sobre toda la distribución $\Pi(p)$, dada la precisión actual de los datos.

5. Significado e Impacto

• Avance en Cosmología Observacional: Este trabajo elimina una de las mayores fuentes de incertidumbre sistemática en la búsqueda de no-gaussianidad primordial usando trazadores de alta redshift. Al proporcionar un prior basado en física real, se mejora la credibilidad de las restricciones sobre los modelos de inflación.
• Preparación para Futuros Datos: Aunque el impacto en la precisión de DR1 es moderado (debido a que los errores estadísticos aún dominan), la metodología es crucial para las futuras liberaciones de datos de DESI y otros sondeos (como Euclid o el LSST/Rubin), donde la precisión estadística aumentará y los errores sistemáticos en el modelado de trazadores serán el factor limitante.
• Validación de Modelos Astrofísicos: Los resultados confirman que los modelos de formación de cuásares en simulaciones como IllustrisTNG (donde los cuásares son impulsados por fusiones) son consistentes con las observaciones de agrupación a gran escala, validando el uso de estas simulaciones para calibrar sesgos cosmológicos.

En resumen, el artículo establece un nuevo estándar para el tratamiento del sesgo de ensamblaje en cuásares, transformando una fuente de degeneración en un parámetro calibrado físicamente, lo que permite una medición más precisa y robusta de la física del universo temprano.