Cosmic variance or galaxy bias? Disentangling finite-volume and galaxy formation effects in cosmological analysis

Resumen Técnico: Desentrañando los Efectos de la Variación Cósmica y la Formación de Galaxias en el Análisis Cosmológico

Planteamiento del Problema
Los actuales y futuros sondeos de redshift de galaxias a gran escala (p. ej., DESI, Euclid) pretenden determinar los parámetros cosmológicos con una precisión sin precedentes. Sin embargo, lograr resultados no sesgados requiere un control riguroso sobre los sistemas de error teóricos. Un desafío primordial es distinguir entre dos fuentes de incertidumbre fuertemente acopladas: la variación cósmica, que surge de la varianza intrínseca del campo aleatorio gaussiano inicial en un volumen finito, y el sesgo de las galaxias (galaxy bias), que describe los complejos procesos no lineales de la formación de galaxias.

En los marcos perturbativos estándar, los parámetros de sesgo se estiman ajustando las observaciones al espectro de potencia de materia promediado en el conjunto. Este enfoque asume implícitamente que el volumen del sondeo representa el promedio del conjunto, ignorando la variación cósmica específica de la realización observada. En consecuencia, los parámetros de sesgo estimados se vuelven degenerados, confundiendo la verdadera física de la formación de galaxias con las fluctuaciones estadísticas del volumen cósmico específico sondeado. Esta degeneración complica aplicaciones tales como la reconstrucción de las Oscilaciones Acústicas de Baryones (BAO) y la inferencia de parámetros cosmológicos, donde las suposiciones sobre los valores de sesgo son críticas.

Metodología
Los autores desarrollan un nuevo formalismo basado en la Teoría de Perturbaciones Euleriana (EPT) para desentrañar matemáticamente la variación cósmica del sesgo de las galaxias.

1. Reencuadrar la Variación Cósmica como Sesgo: Los autores proponen tratar el campo de densidad de una realización arbitraria de volumen finito ($\delta'$) como un trazador "sesgado" del campo de densidad promedio del conjunto ($\delta$). Al introducir un nuevo conjunto de parámetros de sesgo ($b_{1,dm}, b_{2,dm}$), expresan el campo de volumen finito como:
   $$\delta'(x) = b_{0,dm} + b_{1,dm}\delta(x) + \frac{1}{2}b_{2,dm}\delta^2(x) + \dots$$
   En este marco, a medida que el volumen tiende al infinito, $b_{1,dm} \to 1$ y $b_{2,dm} \to 0$. Esta formulación permite la separación del "sesgo de la variación cósmica" del "sesgo físico de las galaxias" ($b_{1,g}, b_{2,g}$).

2. Aplicación al Desplazamiento de BAO: Los autores aplican este formalismo al desplazamiento no lineal de la escala de BAO. Utilizando las expresiones analíticas derivadas de [22], descomponen el desplazamiento total de BAO ($\Delta \alpha$) en dos componentes distintos:

   • Un componente impulsado por la variación cósmica intrínseca del campo de materia oscura (controlado por $b_{1,dm}, b_{2,dm}$).
   • Un componente impulsado por el sesgo de las galaxias (controlado por $b_{1,g}, b_{2,g}$).
     Esta descomposición revela que la variación cósmica por sí sola puede inducir un desplazamiento de BAO no despreciable, particularmente en volúmenes pequeños, el cual es degenerado con el desplazamiento causado por el crecimiento gravitacional no lineal.

3. Manejo de la Restricción de Integral: El artículo aborda la restricción de integral (el sesgo introducido al forzar que la densidad media de un volumen de sondeo sea cero). Los autores demuestran que este efecto puede ser reabsorbido en el término constante de la expansión del sesgo ($b'_0$), asegurando que no introduzca parámetros de sesgo adicionales ni complique la separación de la variación cósmica.

4. Validación Numérica: Para validar los argumentos analíticos, los autores utilizan la suite de simulaciones UNIT. Analizan cuatro conjuntos distintos de simulaciones de N-cuerpos (volúmenes de $2 \times 10^3$ partículas en $(1 h^{-1}\text{Gpc})^3$):

   • 100 simulaciones estándar.
   • 100 simulaciones de amplitud fija (suprimiendo la varianza de amplitud).
   • 2 $\times$ 50 simulaciones emparejadas (suprimiendo la varianza de fase).
   • 2 $\times$ 50 simulaciones de amplitud fija y emparejadas (suprimiendo ambos).
     Miden la posición del pico de BAO de una manera independiente del modelo (ajustando una Gaussiana a la función de correlación residual) para cuantificar la reducción de la varianza.

Resultados Clave

• Desentrañamiento Analítico: El nuevo formalismo separa con éxito las contribuciones de la variación cósmica y el sesgo de las galaxias al desplazamiento no lineal de BAO. El análisis confirma que la variación cósmica puede inducir un desplazamiento comparable al del sesgo de las galaxias, lo que requiere un cuidado adicional al utilizar enfoques perturbativos en volúmenes finitos.
• Supresión de la Varianza: Los resultados numéricos de las simulaciones UNIT demuestran que controlar las condiciones iniciales reduce significativamente la incertidumbre en la posición del pico de BAO:
  • Fijar (controlar la amplitud) reduce la incertidencia por un factor de $\sim 1.5 - 2.1$.
  • Emparejar (controlar la fase) reduce la incertidumbre por un factor de $\sim 1.5 - 1.7$.
  • Combinado (Fijo y Emparejado): La combinación de ambas técnicas produce un factor de supresión de la varianza de $\sim 3.4 - 5.4$ comparado con los conjuntos estándar.
• Convergencia: Se encuentra que la posición del pico de BAO está bien convergida al nivel de sub-percentual en todos los tipos de simulaciones, sin que las técnicas de fijación o emparejamiento introduzcan sesgos aparentes respecto al conjunto estándar.

Significancia y Estrategias Propuestas
El artículo concluye que los métodos actuales para estimar los parámetros de sesgo incluyen inadvertidamente los efectos degenerados de la variación cósmica. Para abordar esto, los autores proponen dos estrategias prácticas para futuros análisis cosmológicos:

1. Combinación Óptima de Multi-trazadores: En lugar de buscar únicamente maximizar la precisión estadística, los análisis futuros podrían combinar trazadores complementarios (p. ej., galaxias y vacíos cósmicos) específicamente para cancelar la variación cósmica. Esto podría implicar la redefinición de muestras de pre-reconstrucción para trazar regímenes de densidad complementarios o la optimización de los pesos en la función de correlación combinada para imitar el efecto de las realizaciones emparejadas.
2. Reconstrucción con Supresión de Varianza: Los autores sugieren aplicar técnicas de supresión de varianza (fijación y emparejamiento) tanto a los datos observados (mediante métodos de reconstrucción como la inversión hacia atrás o la reconstrucción bayesiana hacia adelante) como a los catálogos de simulaciones (mocks) utilizados para la estimación de la covarianza. Esto permitiría el análisis de una versión "supresora de varianza" del Universo observado y de su "anti-Universo" emparejado. Señalan que esto requiere una implementación cuidadosa para evitar sesgar el desplazamiento de BAO a través de la cosmología adoptada en el proceso de reconstrucción.

Los autores enfatizan que, si bien estas estrategias son prometedoras, requieren más pruebas en catálogos de simulaciones detallados (p. ej., tipo DESI) y su integración con técnicas de análisis de vanguardia como la reconstrucción de BAO y la Teoría de Campos Efectivos.