Beyond the Power Spectrum: A New Framework for Non-Stationary Fields with Applications to Light-Cone Effects in Line Intensity Mapping

1. El Problema: La Limitación de la Invariancia Translacional

Las encuestas cosmológicas modernas (como DESI, Euclid, HERA, SKA) cubren volúmenes extremadamente grandes, abarcando escalas de gigaparsecs. En estos volúmenes, no se puede ignorar la evolución cósmica a lo largo de la línea de visión (del extremo cercano al lejano de la encuesta).

• El Efecto del Cono de Luz: A grandes distancias, el volumen observado corresponde a una porción significativa de la historia cósmica. Por ejemplo, en la Época de Reionización (EoR), las regiones más lejanas (alto redshift) muestran propiedades estadísticas diferentes a las regiones cercanas (bajo redshift), ya que la ionización del medio interestelar evoluciona con el tiempo.
• Fallo del Espectro de Potencia Tradicional: El espectro de potencia ($P(k)$) es la estadística resumen estándar, pero se basa en la suposición de invariancia translacional (las propiedades estadísticas son las mismas en todo el campo). Cuando esta suposición se rompe debido al efecto del cono de luz, las diferentes modos de Fourier a lo largo de la línea de visión se correlacionan.
• Pérdida de Información: El espectro de potencia tradicional solo captura la varianza de cada modo individual y asume que los modos son independientes. Al ignorar las covarianzas entre modos (los términos no diagonales en la matriz de covarianza), se pierde una cantidad significativa de información cosmológica y astrofísica contenida en la evolución a lo largo de la línea de visión.

2. Metodología: Una Nueva Base No-Fourier

Los autores proponen un nuevo marco matemático para extraer información de campos no estacionarios, específicamente diseñado para la línea de visión.

• Descomposición Espectral (Eigen Decomposition):
  • En lugar de utilizar la base de Fourier estándar (que diagonaliza la matriz de covarianza solo si el campo es estacionario), los autores calculan la matriz de covarianza a lo largo de la línea de visión ($Cov_{LS}$) utilizando simulaciones cosmológicas (ej. 21cmFAST).
  • Realizan una descomposición en valores propios (eigen decomposition) de esta matriz de covarianza.
  • Los autovectores resultantes forman una nueva base óptima que diagonaliza la matriz de covarianza, capturando así las correlaciones inducidas por el efecto del cono de luz.
• Características de la Nueva Base:
  • Los autovectores se asemejan a modos sinusoidales (como Fourier), pero su amplitud está modulada por una envolvente que depende del redshift (evolución cósmica).
  • El primer autovector corresponde a la señal de temperatura de brillo promedio globalmente.
• Base Híbrida (Solución Práctica):
  • Dado que la convergencia de los autovectores de orden superior (con valores propios pequeños) es lenta y susceptible a fluctuaciones numéricas, proponen una base híbrida.
  • Esta base combina los primeros ~20 autovectores (que capturan la mayor parte de la varianza y la física del cono de luz) con modos de Fourier analíticos para el resto del espectro.
  • Se utiliza el proceso de Gram-Schmidt para asegurar la ortogonalidad estricta entre los autovectores dominantes y los modos de Fourier restantes.
• Nueva Estadística Resumen ($M(k_\perp, i)$):
  • Definen una nueva estadística análoga al espectro de potencia, pero proyectada sobre esta nueva base.
  • $M(k_\perp, i)$ representa la potencia de los coeficientes de expansión en la dirección perpendicular ($k_\perp$) y el índice del autovector ($i$) en la dirección radial.
  • Esta estadística comprime la información de la línea de visión de $N$ frecuencias a $N$ coeficientes, preservando las correlaciones, a diferencia de las funciones de correlación cruzada que requieren $N(N+1)/2$ parámetros.

3. Contribuciones Clave

1. Formalismo Matemático: Derivación rigurosa de una base óptima para campos cosmológicos que no son invariantes bajo traslación a lo largo de la línea de visión.
2. Herramienta de Compresión de Datos: Desarrollo de una estadística resumen eficiente que preserva la información de las correlaciones entre modos de Fourier, superando la compresión subóptima del espectro de potencia tradicional en grandes volúmenes.
3. Guía Práctica: Presentación de un método híbrido viable para observaciones reales, mitigando problemas de convergencia numérica y sobreajuste.
4. Validación Realista: Aplicación del método a un escenario realista de mapeo de intensidad de línea (21 cm) utilizando simulaciones que incluyen ruido térmico instrumental, estrategias de evitación de foregrounds (ventana de la EoR) y la geometría de un interferómetro tipo HERA.

4. Resultados

Los autores realizaron una predicción de Fisher para comparar el rendimiento de su nuevo método frente a dos enfoques de Fourier tradicionales:

• Enfoque 1: Espectro de potencia 2D completo del cono de luz (Full Lightcone).
• Enfoque 2: Espectro de potencia 2D por rebanadas (Slicing), donde el cono de luz se divide en cajas pequeñas estacionarias.
• Enfoque Propuesto: La nueva estadística $M(k_\perp, i)$ basada en la base híbrida.

Hallazgos principales:

• Mejora en Restricciones: El método propuesto logra restricciones en los parámetros astrofísicos aproximadamente un 30% más ajustadas en promedio en comparación con los métodos basados en Fourier.
• Reducción de Degeneraciones: La mejora no proviene principalmente de un aumento en la información diagonal (varianza), sino de una reducción significativa en las degeneraciones entre parámetros (elementos fuera de la diagonal de la matriz de Fisher). Esto permite una inferencia más precisa de parámetros como la eficiencia de formación estelar ($f_*$), la fracción de escape de fotones ionizantes ($f_{esc}$) y la luminosidad de rayos X.
• Robustez: Se demostró que la base es cualitativamente robusta ante variaciones en los parámetros astrofísicos, lo que sugiere que se puede construir la base utilizando estimaciones preliminares de parámetros obtenidas mediante espectros de potencia tradicionales.

5. Significado e Impacto

• Optimización de Futuras Observaciones: A medida que los telescopios como HERA, SKA y telescopios de mapeo de intensidad (LIM) recopilen datos de volúmenes cada vez más grandes, la capacidad de extraer información de la evolución a lo largo de la línea de visión será crucial. Este método permite explotar al máximo la alta resolución en la línea de visión que ofrecen estos instrumentos.
• Superación de Limitaciones Actuales: Proporciona una alternativa superior al enfoque de "rebanado" (que pierde escalas grandes) y al enfoque de cono de luz completo (que mezcla información de diferentes épocas de manera subóptima).
• Generalidad: Aunque se aplica aquí a la línea de 21 cm durante la EoR, el formalismo es general y puede aplicarse a cualquier campo cosmológico o observable donde la invariancia translacional se rompa a lo largo de la línea de visión (ej. Época del Amanecer Cósmico, Edad Oscura).
• Eficiencia Computacional: Ofrece una compresión de datos mucho más eficiente que calcular funciones de correlación cruzada completas, facilitando el manejo de conjuntos de datos masivos en la cosmología de precisión.

En conclusión, este trabajo establece un nuevo estándar para el análisis de datos en grandes volúmenes cosmológicos, demostrando que considerar explícitamente la evolución temporal (efecto del cono de luz) mediante una base de datos adaptada puede revelar información astrofísica oculta que los métodos tradicionales pasan por alto.