Large-scale halo velocity correlations and the impact of finite simulation volumes

Resumen Técnico: Correlaciones de velocidad de halos a gran escala y el impacto de los volúmenes de simulación finitos

Planteamiento del problema
Las funciones de correlación de velocidad medidas a partir de los campos de velocidad peculiar de la materia oscura y de los halos en simulaciones de N-cuerpos cosmológicas exhiben una supresión sistemática en su amplitud en escalas grandes en comparación con las predicciones lineales teóricas. Esta discrepancia es particularmente problemática para la cosmología de precisión, donde las estadísticas de velocidad se utilizan para restringir el parámetro de tasa de crecimiento $f(z)\sigma_8(z)$ y para probar la gravedad. La supresión se atribuye al volumen finito de las cajas de simulación, lo que impone un corte duro en los modos de Fourier, eliminando las perturbaciones de longitud de onda larga que están naturalmente presentes en el Universo real. Si bien se sabe que este efecto altera las correlaciones a gran escala, su impacto específico en las estadísticas de velocidad, su dependencia de la masa del halo y sus consecuencias para la inferencia de parámetros en el modelado basado en simulaciones (como los emuladores) no han sido cuantificados completamente.

Metodología
Los autores investigan este problema utilizando series de simulaciones de N-cuerpos del proyecto DARK QUEST, comparando específicamente dos tamaños de caja:

• L1G: Longitud de lado de $1 \, h^{-1}$ Gpc.
• L2G: Longitud de lado de $2 \, h^{-1}$ Gpc.

Para desentrañar los efectos del tamaño de la caja de los efectos de la impronta de los modos de longitud de onda larga, los autores realizan un análisis de subvolúmenes controlado:

1. Subdivisión: Las cajas L2G más grandes se subdividen en volúmenes más pequeños con longitudes de lado de $1 \, h^{-1}$ Gpc ($L\times0.5_{2G}$) y $0.3 \, h^{-1}$ Gpc ($L\times0.15_{2G}$).
2. Comparación: Comparan las correlaciones de velocidad medidas en estos subvolúmenes frente a las de simulaciones L1G independientes. Crucialmente, los subvolúmenes retienen los modos de longitud de onda larga de la caja L2G de la cual proceden, mientras que las simulaciones L1G independientes no lo hacen.
3. Modelado: La función de autocorrelación de velocidad teórica se modela mediante una transformada de Hankel del espectro de potencia de la materia $P_m(k)$. Los autores introducen un parámetro libre, $k_{min}$, que representa el límite inferior de la integral para imitar los modos de longitud de onda larga ausentes.
4. Análisis: Ajustan el modelo a los datos de la simulación para restringir $f\sigma_8$, probando dos escenarios: (a) marginalizando sobre $k_{min}$ como un parámetro libre, y (b) fijando $k_{min} \approx 0$ (ignorando el efecto del volumen finito). También analizan la dependencia de estos efectos con la masa del halo dividiendo los halos en rangos de masa de grupo y de cúster.

Resultados clave

• Origen de la supresión: La supresión de las amplitudes de la correlación de velocidad es impulsada por la ausencia de modos de longitud de onda larga en el volumen de la simulación, no meramente por el tamaño físico de la región analizada. Los subvolúmenes extraídos de la caja L2G más grande ($L\times0.5_{2G}$) exhiben el mismo nivel de supresión que la caja L2G original, a pesar de tener el mismo volumen que la caja L1G independiente. Esto confirma que la "memoria" de los modos de la caja más grande persiste en los subvolúmenes.
• Dependencia de la masa: La supresión depende de la masa. Los halos más masivos (de tamaño de cúster, $M_h > 10^{14} h^{-1} M_\odot$) exhiben una supresión de la amplitud de la correlación de velocidad sistemáticamente más fuerte en comparación con los halos menos masivos (de tamaño de grupo, $M_h \sim 10^{13} h^{-1} M_\odot$).
• Inferencia de parámetros ($f\sigma_8$):
  • Marginalización de $k_{min}$: Cuando $k_{min}$ se trata como un parámetro libre y se marginaliza, el $f\sigma_8$ recuperado es consistente con el valor real de entrada de las simulaciones ($f\sigma_8 = 0.4852$), independientemente del tamaño de la caja.
  • Ignorar el efecto: Si se fija $k_{min}$ a cero (ignorando el efecto del volumen finito), la restricción sobre $f\sigma_8$ se vuelve significativamente sesgada, particularmente cuando no se utiliza la matriz de covarianza completa (es decir, cuando se ignoran las correlaciones entre los bins de separación). Incluso con la matriz de covarianza completa, persiste un desplazamiento en el valor de mejor ajuste, aunque este permanece dentro de $1\sigma$.
• Importancia de la covarianza: El sesgo en $f\sigma_8$ se exacerba cuando se desestiman los elementos fuera de la diagonal de la matriz de covarianza.

Significancia y afirmaciones
El artículo argumenta que los efectos de volumen finito en las simulaciones de N-cuerpos constituyen un sistema crítico, a menudo subestimado, para el modelado de estadísticas de velocidad basado en simulaciones.

• Implicaciones de modelado: Para enfoques basados en simulaciones como los emuladores, el simple uso de la salida bruta de cajas finitas sin corregir los modos de longitud de onda larga ausentes conduce a restricciones cosmológicas sesgadas. Los autores proponen que introducir un corte de escala grande efectivo ($k_{min}$) como un parámetro libre es un paso necesario para recuperar restricciones no sesgadas de $f\sigma_8$.
• Distinción de las observaciones: Los autores destacan una diferencia fundamental entre los artefactos de simulación y las funciones de ventana observacionales. Mientras que las funciones de ventana observacionales pueden modelarse y deconvolverse explícitamente, el corte duro en las simulaciones está impreso en las condiciones iniciales y persiste a través de la evolución no lineal, lo que lo convierte en una limitación más rígida.
• Interpretación física: La supresión dependiente de la masa implica que los efectos de volumen finito y el sesgo de velocidad no pueden separarse limpiamente cuando se analizan trazadores sesgados. Esto complica las pruebas del principio de equivalencia de Einstein utilizando modelos calibrados por simulaciones, ya que los efectos de volumen finito no contabilizados podrían imitar u ocultar señales físicas genuinas.

El trabajo concluye que los análisis futuros de datos de velocidad peculiar y de la dispersión de la temperatura de Sunyaev–Zel'dovich cinético (kSZ) deben contabilizar consistentemente los modos de longitud de onda larga ausentes cuando utilicen simulaciones como plantillas teóricas para asegurar interpretaciones cosmológicas robustas.