Optimal Transport Reconstruction of Biased Tracers in Primordial Non-Gaussian Fields

Resumen Técnico: Reconstrucción mediante Transporte Óptimo de Trazadores Sesgados en Campos No Gaussianos Primordiales

Planteamiento del Problema
El Transporte Óptimo (OT, por sus siglas en inglés) se ha consolidado como un marco eficiente para reconstruir el campo de desplazamiento cósmico, mapeando las posiciones iniciales a las ubicaciones actuales (o viceversa). Trabajos previos demostraron que el OT recupera con precisión los desplazamientos para trazadores sesgados bajo la suposición de condiciones iniciales gaussianas. Sin embargo, una cuestión crítica sigue abierta: ¿sigue siendo válido este marco cuando las condiciones iniciales son no gaussianas? Este es un desafío oportuno dado el creciente esfuerzo por restringir la no gaussianidad primordial (PNG) utilizando los próximos estudios de la era de las grandes encuestas cosmológicas. Una dificultad específica surge porque las implementaciones de OT suelen asumir una distribución espacial uniforme inicial. Si bien esto se cumple para el campo de materia oscura completo, no se cumple para subconjuntos sesgados (por ejemplo, halos). En escenarios gaussianos, el "polvo" (dust) no observado puede modelarse eficazmente; sin embargo, en escenarios no gaussianos, el propio polvo adquiere sesgos dependientes de la escala específicos que deben tenerse en cuenta para asegurar una reconstrucción precisa.

Metodología
Los autores investigan la aplicabilidad de la reconstrucción por OT a campos no gaussianos primordiales "locales", donde el potencial gravitatorio $\phi_{NG}$ se modela como una serie de potencias de un campo gaussiano $\phi_G$ (truncada al segundo orden, $f_{NL} \neq 0$). En este modelo, el contraste de densidad $\delta_{NG}$ contiene un término proporcional al producto de campos gaussianos, lo que conduce a un sesgo dependiente de la escala distintivo en el agrupamiento de trazadores sesgados en grandes escalas (pequeños $k$).

El estudio utiliza el conjunto de simulaciones Quijote-PNG, que consta de 20 cajas periódicas ($L = 1 h^{-1}$ Gpc) con $512^3$ partículas, evolucionando desde $z=127$ hasta $z=0$ bajo una cosmología $\Lambda$CDM plana. Los trazadores sesgados se identifican como halos con masa $> 20 m_p$.

La innovación metodológica central consiste en modelar el "polvo" no observado (la fracción de masa $1-f_b$ que no está en el catálogo de halos) para satisfacer la restricción de que el campo de materia completo debe exhibir un sesgo independiente de la escala. Los autores derivan que si los trazadores sesgados exhiben un sesgo dependiente de la escala caracterizado por los parámetros $b_\delta$ y $b_\phi$ (donde el sesgo escala como $b(k) \approx b_\delta + 2f_{NL}b_\phi \alpha(k)$), el polvo debe poseer un sesgo dependiente de la escala compensatorio con signo opuesto.

Se prueban dos enfoques de reconstrucción:

1. OT-K: Utiliza un modelo de filtro de Wiener para estimar la distribución del polvo basándose en los trazadores observados. Esto requiere suponer la combinación $2f_{NL}b_\phi$ (que puede ser restringida por datos de encuestas) y la fracción de masa $f_b$.
2. OT-W: Utiliza las partículas de polvo reales de la simulación (una muestra aleatoria) combinadas con los halos ponderados por masa para realizar la reconstrucción, sirviendo como un límite superior de rendimiento.

La calidad de la reconstrucción se evalúa comparando las posiciones reconstruidas con las posiciones reales de los protohalos y analizando la dependencia de la escala de la divergencia del campo de desplazamiento ($\theta_{OT}$), que sirve como un proxy para el campo de densidad inicial.

Contribuciones Clave y Resultados

• Caracterización del Polvo: El artículo demuestra que en modelos de PNG locales, el "polvo" no es un fondo pasivo, sino que exhibe su propio sesgo característico dependiente de la escala. Específicamente, la amplitud del sesgo del polvo está determinada por el requisito de que la suma de los campos de trazadores sesgados y de polvo resulte en un campo total independiente de la escala.
• Extensión de OT a la No Gaussianidad: Los autores muestran que la metodología de OT se extiende con éxito a condiciones iniciales no gaussianas. Al incorporar un modelo realista para el polvo que tiene en cuenta el sesgo dependiente de $f_{NL}$, la reconstrucción de los trazadores sesgados sigue siendo precisa.
• Precisión de la Reconstrucción:
  • Posiciones de los Trazadores: La reconstrucción OT-K (usando el modelo de polvo de filtro de Wiener) recupera correctamente el sesgo dependiente de la escala de los protohalos. La reconstrucción OT-W (usando el polvo real) muestra un acuerdo aún mejor, lo que indica que el método es robusto pero sensible a la fidelidad del modelo de polvo.
  • Proxy del Campo Completo: La divergencia del campo de desplazamiento de OT ($\theta_{OT}$) muestra casi ninguna dependencia de la escala en valores pequeños de $k$, reflejando el comportamiento del verdadero campo inicial. El modelo OT-W extiende esta independencia de escala a valores de $k$ más grandes en comparación con OT-K, logrando una mejora de un factor de 2 en el rango donde $b(k) \approx 1$.
• Robustez ante $f_{NL}$: El valor de la media cuadrática del campo reconstruido $\langle \theta_{OT}^2 \rangle$ permanece estable a través de diferentes valores de $f_{NL}$ ($+100, 0, -100$), con variaciones dentro de la desviación estándar, lo que sugiere que el método no está sesgado por la presencia de no gaussianidad en sí misma.

Significancia y Afirmaciones
El artículo afirma que el Transporte Óptimo proporciona un marco útil para reconstruir estructuras cósmicas incluso cuando las condiciones iniciales son no gaussianas, siempre que el modelo para el "polvo" no observado sea suficientemente realista. Los autores enfatizan que el éxito del método radica en el hecho de que la firma distintiva dependiente de la escala de la PNG está presente tanto en los trazadores sesgados como en el polvo, pero con signos opuestos, permitiendo que el campo completo permanezca sin sesgo.

Los autores afirman modestamente que no demuestran una mejora directa en las restricciones sobre $f_{NL}$ a partir de este trabajo por sí solo. En cambio, proporcionan un marco fundamental que permite que los análisis basados en reconstrucción se apliquen de manera consistente en cosmologías de PNG. Sugieren que este trabajo es un paso necesario para incorporar técnicas de reconstrucción en futuros análisis de alta precisión de PNG, ayudando potencialmente a conectar los trazadores sesgados observados más estrechamente con sus contrapartes lagrangianas primordiales y mitigando contribuciones no lineales o estocásticas sutiles relevantes para las encuestas de próxima generación. El artículo concluye que un mejor modelado del polvo (por ejemplo, mediante métodos iterativos o aprendizaje automático) podría producir mayores beneficios.