Control variates from Eulerian and Lagrangian perturbation theory: Application to the bispectrum

Resumen Técnico: Variables de Control de la Teoría de Perturbaciones Euleriana y Lagrangiana: Aplicación al Bispectro

Planteamiento del Problema
Las simulaciones cosmológicas N-cuerpo son esenciales para modelar la estructura no gaussiana de la red cósmica, particularmente para estadísticas de orden superior como el bispectro. Sin embargo, estas simulaciones sufren de "varianza cósmica" debido a la aleatoriedad de sus condiciones iniciales. Para lograr la alta precisión requerida por los sondeos de galaxias de Fase IV, típicamente se necesita promediar sobre muchas realizaciones de simulación, lo cual es computacionalmente prohibitivo. Aunque existen técnicas de reducción de varianza como el "emparejamiento y fijación", son menos efectivas para los bispectros que para los espectros de potencia. El método de variables de control ofrece una alternativa prometedora al explotar simulaciones más económicas y correlacionadas (sustitutas) para reducir la varianza muestral de los observables. Un desafío crítico es identificar modelos sustitutos que sean tanto analíticamente tratables (para determinar su media sin ensambles costosos) como altamente correlacionados con el bispectro completo N-cuerpo en un amplio rango de escalas.

Metodología
Los autores investigan variables de control perturbativas derivadas tanto de la Teoría de Perturbaciones Euleriana (EPT) como de la Teoría de Perturbaciones Lagrangiana (LPT), utilizando el bispectro de la materia en $z=0.5$ como estudio de caso. El estudio utiliza la suite de simulaciones N-cuerpo Quijote ($N=1000$ realizaciones) como verdad fundamental.

1. Análisis de Modelo de Juguete: Los autores analizan primero un modelo simplificado donde el campo "verdadero" es generado por la Aproximación de Zel'dovich (ZA). Demuestran que las variables de control "lineales" (Gaussianas), aunque tienen media cero, exhiben una descorrelación exponencial con el campo no lineal proporcional a $\exp[-\frac{1}{2}\Sigma^2(k_1^2 + k_2^2 + k_3^2)]$, donde $\Sigma$ es la dispersión de desplazamiento.
2. Variables de Control Eulerianas: Los autores construyen variables de control utilizando campos EPT hasta quinto orden en la densidad ($O(\delta^5)$). Emplean un parámetro de conteo $\epsilon$ para aislar las contribuciones a nivel de árbol ($O(\epsilon)$) y los diagramas de orden superior de media cero ($O(\epsilon^2)$) de la realización completa de la red. Abordan el desafío de igualar las predicciones analíticas con las mediciones de red desarrollando esquemas de "triángulo efectivo" y promediado discreto para corregir los efectos de agrupamiento (binning).
3. Variables de Control Lagrangianas: Los autores evalúan los bispectros de la ZA y la Teoría de Perturbaciones Lagrangiana de Segundo Orden (2LPT). Muestran que, aunque la 2LPT captura la estructura correcta a nivel de árbol y los flujos a gran escala, su bispectro completo no es analíticamente calculable, requiriendo estimación empírica de la media. Por el contrario, el bispectro de la ZA es analíticamente calculable pero posee una estructura a nivel de árbol incorrecta en comparación con la EPT, lo que conduce a una pobre correlación en bajo-$k$.
4. Variables de Control Desplazadas: Para combinar las ventajas de ambos enfoques, los autores introducen una "variable de control desplazada". Esta se construye utilizando la base de "operador desplazado", donde el campo se define como un traza sesgada lagrangiana con solo el desplazamiento lineal (Zel'dovich) exponenciado. Esta construcción asegura que el campo retenga el bispectro a nivel de árbol correcto (igualando a la EPT) mientras reordena los desplazamientos lineales (manteniendo las propiedades de correlación de la LPT).
5. Cálculo Analítico: Para el bispectro de la ZA, los autores desarrollan un método para calcular predicciones analíticas promediadas por bin utilizando pesos derivados de cálculos a nivel de árbol y muestreo Monte-Carlo para bins anchos, logrando un acuerdo sub-1% con las mediciones de red.

Contribuciones y Resultados Clave

• Límites de Descorrelación: El estudio confirma que las variables de control basadas en EPT sufren de descorrelación exponencial en escalas menores que la dispersión de desplazamiento ($k \gtrsim \Sigma^{-1}$), limitando su utilidad a grandes escalas.
• ZA vs. 2LPT: La aproximación de Zel'dovich, aunque analíticamente tratable, es una variable de control subóptima para el bispectro en comparación con el espectro de potencia porque su bispectro a nivel de árbol difiere de la predicción de la EPT. La 2LPT ofrece mejor correlación pero carece de tratabilidad analítica para la media.
• Rendimiento de la Variable de Control Desplazada: La "variable de control desplazada" propuesta supera a todos los demás métodos perturbativos. Combina la estructura correcta a nivel de árbol de la EPT con la reordenación de desplazamientos lineales de la LPT.
  • Reducción de Varianza: Aplicar la variable de control desplazada a una sola simulación N-cuerpo de $V = 1 (\text{Gpc}/h)^3$ reduce la varianza del bispectro de la materia en un factor equivalente a promediar sobre $10^4$ simulaciones para los triángulos de menor-$k$.
  • Precisión: Esto permite mediciones de precisión sub-2% del bispectro de la materia para todas las configuraciones de triángulos en el rango $k_{\text{min}} = 0.04$ a $k_{\text{max}} = 0.47 \, h\text{Mpc}^{-1}$ utilizando una sola simulación.
• Tratabilidad Analítica: Los autores demuestran que la variable de control desplazada es, en principio, analíticamente calculable a todos los órdenes, aunque la implementación numérica eficiente para el bispectro completo sigue siendo una tarea para trabajos futuros. Actualmente, estiman su rendimiento numéricamente.

Significado y Afirmaciones
El artículo afirma que el desarrollo de emuladores precisos del bispectro es ahora factible sin el costo prohibitivo de promediar miles de simulaciones N-cuerpo. Al utilizar la variable de control desplazada, la varianza inherente a la modelización basada en simulaciones puede suprimirse a un nivel donde la precisión está limitada únicamente por la medición de configuraciones de triángulos a pequeña escala, en lugar de la varianza muestral.

Los autores posicionan este trabajo como sentar las bases teóricas para extender las variables de control a toda la jerarquía de $n$-puntos. Afirman explícitamente que la metodología es suficiente para extenderse a trazas sesgadas (galaxias) y espacio de redshift, aunque estas aplicaciones se dejan para trabajos futuros. El artículo no afirma haber resuelto el cálculo analítico del bispectro desplazado completo, sino establece su viabilidad teórica y su rendimiento superior en pruebas numéricas. El significado principal radica en permitir la modelización cosmológica basada en simulaciones de alta precisión, adecuada para los próximos sondeos de estructura a gran escala.